Estimate frequency offset for QAM signal
The QAMCoarseFrequencyEstimator System object™ estimates frequency offset for a QAM signal.
To estimate frequency offset for a QAM signal:
H = comm.QAMCoarseFrequencyEstimator creates a rectangular QAM coarse frequency offset estimator object, H. This object uses an open-loop, FFT-based technique to estimate the carrier frequency offset in a received rectangular QAM signal.
H = comm.QAMCoarseFrequencyEstimator(Name,Value) creates a rectangular QAM coarse frequency offset estimator object, H, with the specified property Name set to the specified Value. You can specify additional name-value pair arguments in any order as (Name1,Value1,...,NameN,ValueN).
Desired frequency resolution (Hz)
Specify the desired frequency resolution for offset frequency estimation as a positive, real scalar of data type double. This property establishes the FFT length that the object uses to perform spectral analysis. The value for this property must be less than or equal to half the SampleRate property. The default is 0.001.
Sample rate (Hz)
Specify the sample rate in samples per second as a positive, real scalar of data type double. The default is 1.
|clone||Create QAMCoarseFrequencyEstimator object with same property values|
|isLocked||Locked status for input attributes and nontunable properties|
|release||Allow property value and input characteristics changes|
|reset||Reset states of the QAMCoarseFrequencyEstimator object|
|step||Estimate frequency offset for QAM signal|
Estimate and correct for a -250 Hz frequency offset in a 16-QAM signal using the QAM Coarse Frequency Estimator System object™.
Create a rectangular QAM modulator System object using name-value pairs to set the modulation order to 16 and the constellation to have an average power of 1 W.
hMod = comm.RectangularQAMModulator('ModulationOrder',16, ... 'NormalizationMethod','Average power', ... 'AveragePower',1);
Create a square root raised cosine transmit filter System object.
hTxFilter = comm.RaisedCosineTransmitFilter;
Create a phase frequency offset object, where the FrequencyOffset property is set to -250 Hz and SampleRate is set to 4000 Hz using name-value pairs.
hPFOError = comm.PhaseFrequencyOffset(... 'FrequencyOffset',-250, ... 'SampleRate',4000);
Create an AWGN channel object with a 25 dB signal-to-noise ratio.
hAWGN = comm.AWGNChannel(... 'NoiseMethod','Signal to noise ratio (SNR)', ... 'SNR',25);
Create a QAM coarse frequency estimator System object with a sample rate of 4 kHz and a frequency resolution of 1 Hz.
hFreqEst = comm.QAMCoarseFrequencyEstimator(... 'SampleRate',4000, ... 'FrequencyResolution',1);
Create a second phase frequency offset object to correct the offset. Set the FrequencyOffsetSource property to Input port so that the frequency correction estimate is an input argument to the step function.
hPFOCorrect = comm.PhaseFrequencyOffset(... 'FrequencyOffsetSource','Input port', ... 'SampleRate',4000);
Create a spectrum analzyer object to view the frequency response of the signals.
hSA = dsp.SpectrumAnalyzer('SampleRate',4000);
Generate a 16-QAM signal, filter the signal, apply the frequency offset, and pass the signal through the AWGN channel.
modData = step(hMod,randi([0 15],4096,1)); % Generate QAM signal txFiltData = step(hTxFilter,modData); % Apply Tx filter offsetData = step(hPFOError,txFiltData); % Apply frequency offset noisyData = step(hAWGN,offsetData); % Pass through AWGN channel
Plot the frequency response of the noisy, frequency-offset signal using the spectrum analyzer. The signal is shifted 250 Hz to the left.
hSA.Title = 'Received Signal'; step(hSA,noisyData);
Estimate the frequency offset using the step function associated with hFreqEst. Observe that the estimate is close to the -250 Hz target.
estFreqOffset = step(hFreqEst,noisyData)
estFreqOffset = -249.7559
Correct for the frequency offset using hPFOCorrect and the inverse of the estimated frequency offset.
compensatedData = step(hPFOCorrect,noisyData,-estFreqOffset);
Plot the frequency response of the compensated signal using the spectrum analyzer. The signal is now properly centered.
hSA.Title = 'Frequency-Compensated Signal'; step(hSA,compensatedData);
 Nakagawa, T., Matsui, M., Kobayashi, T., Ishihara, K., Kudo, R., Mizoguchi, M., and Y. Miyamoto. "Non-data-aided wide-range frequency offset estimator for QAM optical coherent receivers", Optical Fiber Communication Conference and Exposition (OFC/NFOEC), 2011 and the National Fiber Optic Engineers Conference , March, 2011, pp. 1–3.
 Wang, Y., Shi. K., and E. Serpedin. "Non-Data-Aided Feedforward Carrier Frequency Offset Estimators for QAM Constellations: A Nonlinear Least-Squares Approach", EURASIP Journal on Advances in Signal Processing, Vol. 13, 2004, pp. 1993–2001.