Data-Aided and Non-Data-Aided SNR Estimators for CPM Signals in Ka-Band Satellite Communications

Adaptive coded and modulation (ACM) is an effective measure to resist rain attenuation in Ka-band satellite communications. The accuracy of the estimator for signal-to-noise ratio (SNR) is one of the main factors that affect ACM performance. This paper establishes the channel model of satellite communication in Ka-band and derives and analyzes data-aided (DA) maximum likelihood (ML) and method-of-moments estimators for SNR of continuous phase modulation (CPM) in Ka-band fading channel. Simulations and analysis indicate that the normalized mean square error (NMSE) of the DA ML estimator is closer to the Cramer-Rao bounds (CRB) at low SNR, but the performance worsens with the SNR increases from medium to high values. The M2M4 estimator performs poorly at low SNR and best at medium SNR. When the SNR is low, the performance of the DA ML estimator is better than the M2M4 estimator; however, the two estimators have similar performance when SNR is high. The performances of the two estimators will become increasingly better the greater the lengths of the observation signals become. However, the influence of the signal length will become increasingly smaller as the SNR becomes larger.


Introduction
Many bands can be used for satellite communication; these bands include L, S, C, Ku, and Ka bands.Adopting a high-frequency band in the construction of a broadband satellite communication network is necessary to address the issue of the limited bandwidth of satellite communication.Current high-throughput satellite (HTS) systems for broadband-distributed user access are designed following two main concepts: The use of Ka-band radio frequency (RF) links both for the forward and for the return link; The use of multispot coverage: this technique is largely applied to increase the system throughput through frequency reuse and system reconfigurability [1].While with the exploitation of extremely high-frequency (EHF) bands (30-300 GHz) for broadband transmission over satellite links, the Q-V band (30-50 GHz) and W-band (75-110 GHz) seem to offer very promising perspectives [2].This paper use Ka-band as the operating frequency band for satellite communication because of the advantages of the wideband, little disturbance, and small-sized terminals [3].
At frequencies beyond 10 GHz, rainfall is one of the important atmospheric parameters (hydrometeors) which causes signal degradation along both terrestrial and satellite tracks [4][5][6].Therefor the rain attenuation is considered to be the main factor in Ka-band satellite communications [7][8][9].Rain rate data was measured through OTT Parsivel, which shows the rain rate of about 50 mm/h and attenuation of 10.7, 11.6, and 11.3 dB for 12. 25, 19.8, and 20.73 GHz, respectively, for 0.01% of the time for the combined values of rain rate and rain attenuation statistics [10,11].The design of millimeter wave communication links and the study of propagation impairments Information 2017, 8, 75; doi:10.3390/info8030075www.mdpi.com/journal/informationat higher frequencies due to a hydrometeor, particularly rain, require the knowledge of 1-min.rainfall rate data [12,13].ITU-R series play an important role in prediction of rainfall attenuation.The effectiveness of ITU-RP.618-12method is seen in which the better prediction for rain-induced attenuation and the propagation data and prediction methods required for the design of Earth-space telecommunication systems.ITU-RP.311-16shows the data acquisition, representation and analysis of radio wave propagation studies [14].The adaptive-coded and modulation (ACM) system is an effective measure to resist rain attenuation because it can guarantee communication quality and overcome the problem of system resource waste caused by rainfall reserve.Hence, ACM could be one of the technologies with the most potential and value.
The basic principle of ACM technology in Ka-band satellite communication is to estimate the signal-to-noise ratio (SNR) at the receiver and feed the value back to the transmitter.The transmitter changes the coding and modulation mode according to the signal-to-noise ratio (SNR).The common modulation used for satellite communications mainly includes multiple phase shift keying (M-PSK) and Multiple Quadrature Amplitude Modulation (M-QAM) [15].The limitations of these modulation techniques in the frequency band utilization restrict the development of satellite communications because of the continuous improvement of the channel bandwidth requirements of satellite communications.Continuous phase modulation (CPM) is a constant envelope continuous phase modulation signal with a wide bandwidth, high power utilization, and a certain coding gain compared with phase shift keying (PSK) modulation [16,17].CPM also has a larger signal state under some specific parameters, which makes a larger matched filters, leading to a high implementation complexity of optimal receiver.However, this problem can be solved with the effective choice of parameters and the use of large capacity devices.Hence, CPM has become one of the most important modulation methods in recent years.This paper uses the CPM signal to replace the PSK modulation as the main modulation mode of ACM system in Ka-band satellite communication.
The accuracy of channel state estimation is one of the main factors that affect the performance of ACM, and the SNR is used to characterize channel state information.The SNR estimation method consists of time domain and frequency domain methods.The time domain method can be divided into data-aided (DA) and non-data aided (NDA) methods [18].The DA method has higher estimation accuracy than NDA method but involves the insertion of a periodic pilot sequence, which is inefficient.A significant number of studies have been conducted on SNR estimation technology in recent years.The DA and NDA ML SNR estimators for CPM signals were studied by Rice, and the effects of the frequency and phase shifts were also considered.However, the simulation is based on the additive Gaussian white noise (AWGN) channel [19].Chang studied the SNR estimation method for MPSK signals in Nakagami fading channels and proposed a joint estimation algorithm for SNR and fading factor; however, the algorithm is suitable only for MPSK [20].Based on the second generation digital satellite TV broadcasting standard (DVB-S2), Wu proposed an SNR estimation algorithm based on eigenvalue decomposition in the AWGN channel; however, it is not used in the Ka band [21].Li proposed the ML SNR estimation algorithm for M-APSK signals in slow time-varying fading channels, which improve the estimation accuracy under low SNR, but the CPM signal is not introduced [22].
The SNR estimation algorithm in Ka-band satellite communication is based mostly on PSK modulation, and no relevant literature applies the CPM signal to the ACM technology in Ka-band satellite communication.This study analyzes the channel model in Ka-band satellite communication and focuses on studies on SNR estimation algorithm for CPM signals in the Ka-band fading channel.DA ML and M 2 M 4 estimators are derived and analyzed.The SNR fitting curve, normalized mean square error (NMSE), and CRLB are used to analyze the performance of the estimation algorithm, and the effects of several parameters on SNR estimation are discussed.
The paper is organized as follows.Section 2 establishes the channel model of Ka-band satellite communications.Section 3 presents the schematic of the system and the expression of the received signal.Section 4 develops and analyzes two SNR estimators and provides the expression of CRLB.Section 5 presents the simulation results.Section 6 gives the conclusions.

Attenuation Factor
Many factors in the atmosphere cause the additional decline of the signal in the Ka-band satellite communication, precipitation, water vapor, clouds, oxygen, and flickers.The most important factor is rain attenuation, which is related to frequency, weather conditions, and system availability.These factors can be considered in the satellite communication link by modifying the free space propagation formula.
Signal fading caused by rainfall is the most serious propagation loss in Ka-band satellite communication.Measured data show that rain attenuation is a function of carrier frequency and system feasibility.In the case of large elevation [22], L r can be expressed as where C 1 , C 2 , δ 1 , and δ 2 are the functions of the system availability, and f is the carrier frequency.When the elevation is small, the formula can be simplified as where θ is the elevation of the earth station, and θ 0 is the reference elevation.Atmospheric absorption includes mainly oxygen and water vapor absorption.The International Telecommunication Union Radio communications sector (ITU-R) provides the expression of oxygen absorption loss under a slant path as follows where h o is the effective height of dry air, γ o is the oxygen loss coefficient, and θ is the elevation.ITU-R also gives the expression of water vapor absorption loss where h w is the effective height of water vapor, γ w is the water vapor loss coefficient, and θ is the elevation.The size of cloud attenuation is related to the content and temperature of liquid water.ITU-R gives the expression of cloud attenuation where h refers to the cloud layer thickness, r and i are the real and imaginary parts of the dielectric constant of water, respectively, and θ is the elevation.The refractive index of the propagation path rises randomly with position and time.A flicker exhibits rapid fluctuation when the radio waves pass through the atmosphere, causing fast fluctuation of the amplitude and phase of the radio wave.A flicker can attenuate and enhance the electric wave signal.
ITU-R suggests that the total attenuation value L t can be expressed as The free space propagation formula in the design of the Ka-band satellite communication system should be amended with the additional attenuation caused by rainfall, atmosphere, cloud, and flicker.It can be expressed as where P R and P T are the average power of the receiver and transmitter, respectively, (w) is the unit, f is the carrier frequency, l is the distance from satellite to earth terminal, c is the speed of light in vacuum, A R and A T are the effective contact areas of receiving antenna and transmitting antenna, respectively, L t is the total attenuation for link, and p refers to system unavailability.Formula (7) shows that the multiplicative factor m( f , p) = f 2 /L t increases with the increase of f .L t increases rapidly with the increase of f and the feasibility of the system.m( f , p) decreases with the increase of f .Numerical calculations show that the multiplicative factor m( f , p) remains unchanged in a wide frequency range when system feasibility is low (e.g., the system feasibility is less than 99.7%) [23].The total fading of the channel can be considered as non-frequency selective fading given that the feasibility of the Ka-band satellite communication system is usually less than 99.7%.The entire channel can be viewed as flat fading and AWGN; flat fading is slow over time and is related only to the weather [24].

Statistical Channel Model of Ka-Band Satellite Communication
When transmitting a digital modulation signal s(t) = Re[s 1 (t) exp(j2π f c t)] through the channel, s 1 (t) is the time domain complex baseband, and s 1 ( f ) is the frequency domain complex baseband.The equivalent lowpass signal (excluding AWGN) can be expressed as where C( f , t) is the time-frequency transfer function of the channel that can be viewed as a complex value independent frequency in the signal bandwidth because the Ka-band fading channel is a non-frequency selective channel.Considering the spectrum of s 1 ( f ) concentrated in the vicinity of f = 0 Hz, C( f , t) = C(0, t) can be substituted into Formula ( 8) The fading characteristics can be regarded as a constant in a symbol interval because of the slow time-varying characteristics of multiplicative interference.Formula (9) can be simplified as where T is the cycle of modulation symbol and α and θ are the envelope and phase of the equivalent lowpass channel, respectively.The researches show that the characteristics of a Ka-band satellite communication channel are determined mainly by the weather.The envelope and phase obey the Gauss distribution [24], and the probability density function can be expressed as where, σ 1 , σ 2 , m 1 , and m 2 are available in Table 1.
Table 1 shows the probability distribution parameters of the signal envelope and phase under various weather conditions in a Ka-band satellite communication channel [24].Where, σ 1 and σ 2 are the variances of the signal envelope and phase, respectively, m 1 and m 2 are the mean of the signal envelope and phase, respectively.Previous studies have accumulated a considerable amount of data from many propagation characteristic measurement experiments in Olympus-star, Italsat-star and ACTS-star [24,25], the values in Table 1 is thus obtained.Through the above analysis, we can establish the Ka-band satellite communication channel statistical model as shown in Figure 1.The digital modulation signal can be expressed as a result of the frequency independent complex multiplicative interference factor and AWGN after the non-frequency selective fading channel.C(t) = α exp(jθ) is the channel multiplication vector, and z(t) is the AWGN.
Information 2017, 8, 75 5 of 15 Table 1 shows the probability distribution parameters of the signal envelope and phase under various weather conditions in a Ka-band satellite communication channel [24].Where, 1 σ and 2 σ are the variances of the signal envelope and phase, respectively, 1 m and 2 m are the mean of the signal envelope and phase, respectively.Previous studies have accumulated a considerable amount of data from many propagation characteristic measurement experiments in Olympus-star, Italsat-star and ACTS-star [24,25], the values in Table 1 is thus obtained.Through the above analysis, we can establish the Ka-band satellite communication channel statistical model as shown in Figure 1.The digital modulation signal can be expressed as a result of the frequency independent complex multiplicative interference factor and AWGN after the nonfrequency selective fading channel.( ) exp( ) C t j α θ = is the channel multiplication vector, and ( ) z t is the AWGN.

System Model
The Ka-band satellite channel can be simulated by a multiplicative interference vector ( ) C t and ( ) z t from Section 2. The schematic of the system can be established as shown in Figure 2.

System Model
The Ka-band satellite channel can be simulated by a multiplicative interference vector C(t) and z(t) from Section 2. The schematic of the system can be established as shown in Figure 2. Table 1 shows the probability distribution parameters of the signal envelope and phase under various weather conditions in a Ka-band satellite communication channel [24].Where, 1 σ and 2 σ are the variances of the signal envelope and phase, respectively, 1 m and 2 m are the mean of the signal envelope and phase, respectively.Previous studies have accumulated a considerable amount of data from many propagation characteristic measurement experiments in Olympus-star, Italsat-star and ACTS-star [24,25], the values in Table 1 is thus obtained.Through the above analysis, we can establish the Ka-band satellite communication channel statistical model as shown in Figure 1.The digital modulation signal can be expressed as a result of the frequency independent complex multiplicative interference factor and AWGN after the nonfrequency selective fading channel.( ) exp( ) C t j α θ = is the channel multiplication vector, and ( ) z t is the AWGN.

System Model
The Ka-band satellite channel can be simulated by a multiplicative interference vector ( ) C t and ( ) z t from Section 2. The schematic of the system can be established as shown in Figure 2.  The data source produces a binary sequence.The data symbol generated at any time is independent of the previous symbol.The probability of the two symbols is the same.The information was sequenced through the CPM modulator and transmitted to the Ka-band fading channel.Gaussian Processes 1 and 2 in Figure 2 are considered as two independent stochastic processes with a mean of zero, and the variances are σ 1 and σ 2 .Gaussian process 1 and the real constant generator 1 are added to generate the Gaussian stochastic process α.Gaussian stochastic process θ can be generated in the same way.θ is sent to the exponential generator and multiplied with α to obtain the multiplicative interference vector C(t) = α exp(jθ).Different weather conditions only need to change the variance and the mean of α and θ.The received signal can be represented as r(t) (the expression is shown in follows) after the Ka-band fading channel.The SNR is estimated according to the received signal and the SNR estimators are shown in Section 4. A suitable coded modulation scheme is selected in the feedback link according to the value of SNR and the ACM schemes switching standards (Since the emphasis of this paper is on SNR estimators, we do not introduce too many coding process).The data is transmitted to the Ka-band fading channel according to the selected coded modulation scheme and the received signal is processed by Viterbi algorithm.
Transmitting terminal uses a complex-valued low-pass CPM signal, which is given by where, for nT where T is the symbol time, a = (a 0 , a 1 • • • ) is the transmitted information sequence, h = K/P is the modulation index, where K and P are the coprime integer, and q(t) is the phase pulse and is a continuous monotone function The CPM signal can be obtained by integrating the rectangular pulse (LREC) with a duration of LT, the raised cosine pulse (LRC), or the Gaussian minimum shift keying (GMSK) pulse and expressed as where g(t) is the pulse shaping function and is nonzero in 0 ≤ t ≤ LT. g(t) is a full response when L = 1 and a partial response when L > 1.The received signal is given by where A is the true amplitude of the received signal, s(t) is the CPM modulation signal, and z(t) is additive Gaussian white noise with the variance of δ 2 .This paper aims to estimate the SNR according to the received signal.The signal fading caused by the envelope plays a decisive role, and signal fading caused by the phase is almost negligible.Therefore, the mean value of the envelope can be approximated as the multiplicative fading in the channel.In practical applications, the multiplication factor m 1 of the channel when the weather condition is different can be determined according to Table 1.The power of the transmission signal can be expressed as The average SNR for Ka-band fading channels can be expressed as

SNR Estimators
After sampling the received signal in Formula ( 16), the output sequence is expressed as where K is the observation point after sampling, r k obeys the Gauss distribution, As k is the mean, and δ 2 is the variance.

DA ML Estimator
The probability density function of the received sequence in Formula ( 20) is Given the statistical independence in sampling, Formula ( 21) can be expressed as and the corresponding log-likelihood function is The partial derivative of A and δ 2 is respectively set as and where, ln p(r) has a maximum value if the result is zero.A and δ 2 can be expressed, respectively as and The SNR estimate is computed using

M 2 M 4 Estimator
The DA ML estimator requires the receiver and prior information sequences of the transmitter.However, obtaining prior information is impossible in some cases of practical application.This paper presents the M 2 M 4 estimation method without prior information sequence.The two-order moments and four-order moments of the received sequence in Formula (20) can be expressed as and: where, the conjugate of r k , s k and z k are denoted r * k , s * k and z * k .Given that the signal and noise are independent of each other, and the mean is zero, Equations ( 29) and (30) can be simplified as follows and The estimated value of A and δ 2 can be expressed as The corresponding SNR estimator is M 2 and M 4 can be obtained in practical applications through the average value of the receiving sequence and expressed respectively as and (37)

Performance Index
Two estimators are simulated in the MATLAB simulation platform and the SNR fitting curve is drawn.The abscissa is the SNR value of the theory.The ordinate is the SNR value of the estimate, which can determine whether the SNR estimators are unbiased.The main index to measure the performance of SNR estimators in the simulation is NMSE, which can be expressed as where ρ is the SNR value of the theory and ρ is the SNR value of the estimate.The estimator performs better when the NMSE is smaller.The CRLB provides absolute lower limit of the estimation.The NMSE can be compared with the CRLB to illustrate the absolute performance of the estimator.A parameter vector β = [S N] T is defined as the relationship between SNR, and the parameter vector is where S = ( A m 1 ) 2 is the signal power and N = δ 2 is the noise power.The formula given by Kay in Reference [26] shows that where is the Fisher Information Matrix, which can be expressed as When the SNR is low (0-6 dB), the NDA CRLB is higher than DA CRLB; however, the two lower bounds tend to coincide with the increase of SNR (greater than 6 dB) [27].The CRLB of the DA SNR estimation is simple and unaffected by the shaping filter and sampling.Thus, this study uses DA CRLB as the reference of absolute performance.The CRLB of the complex signal can be expressed as [28] N MATLAB tools are used to simulate and analyze the above theoretical models.The data source produces a binary sequence and sequenced through the CPM modulator.Transmitted the modulation signal to the Ka-band fading channel in a number of frames in the range of SNR.Seven weather conditions are generated according to the numerical information in Table 1.After the fading factor is obtained, the fading channel simulation model is established by using the theoretical channel model in Section 2. The SNR is estimated according to Formulas ( 28) and (35) in Section 4. Numerical results and simulation performance are shown in Section 5.

Numerical Results and Analysis
The length of the source information sequence is set to 1000 bits, and the sampling points per symbol are set to IPOINT = 20.For the CPM modulator, the modulation index is set to h = 1/4, the pulse shaping function is set to GMSK, L = 2, and the modulation order is set to M = 8.The average value of the SNR estimators is calculated with 1000 frames.
The SNR fitting curve of the DA ML and M 2 M 4 estimators are shown in Figure 3.The horizontal axis in the figure is shown as the actual SNR, which is the assumed value in the MATLAB simulation and is the channel environment in the actual communication.The ordinate is the SNR estimated by the Formulas ( 28) and (35) in Section 4. The horizontal and ordinate of the truth values are all the actual values of SNR and the true values are the standard used to measure the offset of estimate SNR and actual SNR.The DA ML estimator can be approximated as unbiased when SNR is less than 5 dB.The estimation error is large in the case of medium and high SNR.The estimation bias increases with the increase of SNR.The M 2 M 4 estimator can be approximated as unbiased when SNR is medium (near 5 dB), but the estimation error is significantly large at low and high SNRs.The NMSE of the DA ML and M2M4 estimators are shown in Figure 4.The DA ML estimator has the best performance when SNR is approximately 0 dB, but the performance becomes worse with the increase of SNR.The M2M4 estimator has the best estimation performance when SNR is near 7 dB.The performance curve of the M2M4 estimator is similar to a parabola, and the performance is poor at low SNR.The performance becomes worse with the increase of SNR after 7 dB.The two estimators cannot achieve the Cramer-Rao bounds and the existence of the SNR node.The performance worsens with the increase of SNR after the node because the non-CM characteristics introduced by the fading factor.When the SNR is low, the noise power is the main factor affecting the estimation performance due to the large noise power, while the signal power is the main factor affecting the estimation performance in the medium to high SNR.In the process of signal power estimation, it is necessary to do the division of the mean of the fading factor as in Formulas (18).The difference between the estimated value and actual value of the signal power is greater with the increase of SNR, which will lead to the worse performance.The DA ML estimator displays better performance than the M2M4 estimator, especially in the low SNR.The performance gap between the two estimators is reduced gradually when the SNR is medium or high.The two estimators are close to each other in practical applications when the channel environment is good.Thus, the M2M4 estimator can be used to reduce the complexity.The DA ML estimator can be used to improve the estimation performance when the channel environment is poor, and the transmitter information sequence is known.The NMSE of the DA ML and M 2 M 4 estimators are shown in Figure 4.The DA ML estimator has the best performance when SNR is approximately 0 dB, but the performance becomes worse with the increase of SNR.The M 2 M 4 estimator has the best estimation performance when SNR is near 7 dB.The performance curve of the M 2 M 4 estimator is similar to a parabola, and the performance is poor at low SNR.The performance becomes worse with the increase of SNR after 7 dB.The two estimators cannot achieve the Cramer-Rao bounds and the existence of the SNR node.The performance worsens with the increase of SNR after the node because the non-CM characteristics introduced by the fading factor.When the SNR is low, the noise power is the main factor affecting the estimation performance due to the large noise power, while the signal power is the main factor affecting the estimation performance in the medium to high SNR.In the process of signal power estimation, it is necessary to do the division of the mean of the fading factor as in Formulas (18).The difference between the estimated value and actual value of the signal power is greater with the increase of SNR, which will lead to the worse performance.The DA ML estimator displays better performance than the M 2 M 4 estimator, especially in the low SNR.The performance gap between the two estimators is reduced gradually when the SNR is medium or high.The two estimators are close to each other in practical applications when the channel environment is good.Thus, the M 2 M 4 estimator can be used to reduce the complexity.The DA ML estimator can be used to improve the estimation performance when the channel environment is poor, and the transmitter information sequence is known.The NMSE of the DA ML and M2M4 estimators are shown in Figure 4.The DA ML estimator has the best performance when SNR is approximately 0 dB, but the performance becomes worse with the increase of SNR.The M2M4 estimator has the best estimation performance when SNR is near 7 dB.The performance curve of the M2M4 estimator is similar to a parabola, and the performance is poor at low SNR.The performance becomes worse with the increase of SNR after 7 dB.The two estimators cannot achieve the Cramer-Rao bounds and the existence of the SNR node.The performance worsens with the increase of SNR after the node because the non-CM characteristics introduced by the fading factor.When the SNR is low, the noise power is the main factor affecting the estimation performance due to the large noise power, while the signal power is the main factor affecting the estimation performance in the medium to high SNR.In the process of signal power estimation, it is necessary to do the division of the mean of the fading factor as in Formulas (18).The difference between the estimated value and actual value of the signal power is greater with the increase of SNR, which will lead to the worse performance.The DA ML estimator displays better performance than the M2M4 estimator, especially in the low SNR.The performance gap between the two estimators is reduced gradually when the SNR is medium or high.The two estimators are close to each other in practical applications when the channel environment is good.Thus, the M2M4 estimator can be used to reduce the complexity.The DA ML estimator can be used to improve the estimation performance when the channel environment is poor, and the transmitter information sequence is known.Three parameters affect the estimation performance by setting the different values of the parameters.These parameters include the number of sampling points per symbol, the length of the transmitted information sequence, and the modulation order.The common parameters are the same as the previous settings.
The NMSEs of DA ML and M 2 M 4 for IPOINT = 10, 20, 40 are plotted in Figure 5.The estimation performance improves with the increase in the number of sampling points per symbol.
Three parameters affect the estimation performance by setting the different values of the parameters.These parameters include the number of sampling points per symbol, the length of the transmitted information sequence, and the modulation order.The common parameters are the same as the previous settings.
The NMSEs of DA ML and M2M4 for IPOINT = 10, 20, 40 are plotted in Figure 5.The estimation performance improves with the increase in the number of sampling points per symbol.The NMSEs of DA ML and M2M4 for the length of transmission information for 1000, 2000, and 4000 bits are plotted in Figure 6.The estimation performance improves the longer the transmitted information.Three parameters affect the estimation performance by setting the different values of the parameters.These parameters include the number of sampling points per symbol, the length of the transmitted information sequence, and the modulation order.The common parameters are the same as the previous settings.
The NMSEs of DA ML and M2M4 for IPOINT = 10, 20, 40 are plotted in Figure 5.The estimation performance improves with the increase in the number of sampling points per symbol.The NMSEs of DA ML and M2M4 for the length of transmission information for 1000, 2000, and 4000 bits are plotted in Figure 6.The estimation performance improves the longer the transmitted information.The three figures show that the three parameters have no influence on the estimation performance when SNR is greater than 8 dB.The observation points of the receiver can be expressed as 2 ( )* log ( ) We can obtain the conclusion by considering the three parameters.The more observation points the receiver has, the better the two estimators performs.The number of observation points has no influence on the estimation performance when the SNR is greater than 8 dB.When the channel environment is poor in practical applications, better estimation performance can be obtained by increasing the number of sampling points per symbol, increasing the length of the transmitted information sequence, and reducing the CPM modulation order.
SNR estimators are the main research content of this paper, and we only select the modulation orders as the example to illustrate the principle of ACM control loop.The length of the source information sequence is set to 1000 bit, and the sampling points per symbol are set to IPOINT = The BER for different modulation orders are shown in Figure 8.For the same SNR, the BER of the system decreases as the modulation order becomes larger.The target BER is assumed to be 10 −4 .When the SNR is less than 10.6 dB, all the three modulation orders can't achieve the target BER, we can choose a low BER method as the modulation order of the system, that is M = 8.When the SNR is within the range of 10.6 dB to 11.4 dB, only M = 8 can achieve the target BER, so we still choose M = 8 as the modulation order of the system.When SNR is in the range of 11.4 dB to 14 dB, both M = 4 and M = 8 can achieve the target BER.Because of the higher frequency band utilization rate of M = 4, we select M = 4 as the modulation order of the system.When the SNR is greater than 14 dB, all the three modulation orders can achieve the target BER.Because of the highest frequency band utilization in M = 2, we select M = 2 as the modulation order of the system.When SNR grows the modulation order is decreased by the ACM control loop, when SNR drops down higher modulations have to be used.The three figures show that the three parameters have no influence on the estimation performance when SNR is greater than 8 dB.The observation points of the receiver can be expressed as We can obtain the conclusion by considering the three parameters.The more observation points the receiver has, the better the two estimators performs.The number of observation points has no influence on the estimation performance when the SNR is greater than 8 dB.When the channel environment is poor in practical applications, better estimation performance can be obtained by increasing the number of sampling points per symbol, increasing the length of the transmitted information sequence, and reducing the CPM modulation order.
SNR estimators are the main research content of this paper, and we only select the modulation orders as the example to illustrate the principle of ACM control loop.The length of the source information sequence is set to 1000 bit, and the sampling points per symbol are set to IPOINT = 10 * log 2 (M).For the CPM modulator, the modulation index is set to h = 1/2, the pulse shaping function is set to RC, L = 2, and the modulation orders are set to M = 2, 4, 8, respectively.The average value of the SNR estimators is calculated with 1000 frames.The BER for different modulation orders are shown in Figure 8.For the same SNR, the BER of the system decreases as the modulation order becomes larger.The target BER is assumed to be 10 −4 .When the SNR is less than 10.6 dB, all the three modulation orders can't achieve the target BER, we can choose a low BER method as the modulation order of the system, that is M = 8.When the SNR is within the range of 10.6 dB to 11.4 dB, only M = 8 can achieve the target BER, so we still choose M = 8 as the modulation order of the system.When SNR is in the range of 11.4 dB to 14 dB, both M = 4 and M = 8 can achieve the target BER.Because of the higher frequency band utilization rate of M = 4, we select M = 4 as the modulation order of the system.When the SNR is greater than 14 dB, all the three modulation orders can achieve the target BER.Because of the highest frequency band utilization in M = 2, we select M = 2 as the modulation order of the system.When SNR grows the modulation order is decreased by the ACM control loop, when SNR drops down higher modulations have to be used.

Conclusions
ACM technology is one of the effective measures to resist rain attenuation in Ka-band satellite communication.The accuracy of SNR estimation is one of the main factors that affect the performance of the ACM system.This paper establishes the channel model of Ka-band satellite communication, and derives and analyzes the DA ML and M 2 M 4 estimators for the SNR of CPM in Ka-band fading channel.The SNR fitting curve, NMSE, and CRLB were used as indices for analysis of the performance of the two estimators.The two estimators cannot achieve the CRLB, and the estimation performance in the medium and high SNR worsens with the increase of SNR.When the SNR is low, the performance of the DA ML estimator is better than the M 2 M 4 estimator.The higher the number of observation points, the better the estimation performance of the two.When the SNR is high, the two estimators have similar estimation performance, and the observation points have no influence on the estimation performance.The two estimators in the commonly used SNR range have a better estimation performance and can be approximated as unbiased.Therefore, CPM can be a competitive modulation mode of the ACM technology in Ka-band satellite communication.

Figure 2 .
Figure 2. Schematic of the system.

Figure 2 .
Figure 2. Schematic of the system.

Figure 2 .
Figure 2. Schematic of the system.

Figure 3 .
Figure 3.Estimated SNR comparison of DA ML and M2M4 in Ka-band fading channel.

Figure 4 .
Figure 4. NMSE comparison of DA ML and M2M4 in Ka-band fading channel.

Figure 3 .
Figure 3.Estimated SNR comparison of DA ML and M 2 M 4 in Ka-band fading channel.

Figure 4 .
Figure 4. NMSE comparison of DA ML and M2M4 in Ka-band fading channel.

Figure 4 .
Figure 4. NMSE comparison of DA ML and M 2 M 4 in Ka-band fading channel.

Figure 5 .
Figure 5. NMSE comparison of DA ML and M2M4 in different IPOINT.

Figure 6 .
Figure 6.NMSE comparison of DA ML and M2M4 in different lengths of information.

Figure 5 .
Figure 5. NMSE comparison of DA ML and M 2 M 4 in different IPOINT.

Figure 5 .
Figure 5. NMSE comparison of DA ML and M2M4 in different IPOINT.

Figure 6 .
Figure 6.NMSE comparison of DA ML and M2M4 in different lengths of information.

Figure 6 .
Figure 6.NMSE comparison of DA ML and M 2 M 4 in different lengths of information.

Figure 7 .
Figure 7. NMSE comparison of DA ML and M2M4 in different M.

2 10 *
log ( ) M .For the CPM modulator, the modulation index is set to h = 1/2, the pulse shaping function is set to RC, L = 2, and the modulation orders are set to M = 2, 4, 8, respectively.The average value of the SNR estimators is calculated with 1000 frames.

Figure 7 .
Figure 7. NMSE comparison of DA ML and M 2 M 4 in different M.

Table 1 .
Ka-band satellite channel envelope and phase model.

Table 1 .
Ka-band satellite channel envelope and phase model.

Table 1 .
Ka-band satellite channel envelope and phase model.