一种基于非线性变换的调相信号频偏估计及跟踪方法技术

本发明专利技术涉及一种基于非线性变换的调相信号频偏估计及跟踪方法，通过二相BPSK、四相平衡QPSK和非平衡四相UQPSK调制信号的理论分析，首次提出了利用一次或二次非线性变换，得到了去除符号影响的射频单载波信号，为傅立叶FFT频谱分析，实现快速、准确地载波频偏估计创造了先决条件。与现有卫星技术中对调相信号直接进行FFT频谱分析等措施相比，本发明专利技术解决了包络谱峰值点分散，频偏估计误差大等技术缺点，大大提高了调相信号频偏估计及跟踪的准确性和实时性，在实现卫星不同链路调相信号的载波频偏估计及跟踪方面有着广阔的应用前景。

【技术实现步骤摘要】
一种基于非线性变换的调相信号频偏估计及跟踪方法
本专利技术属于卫星测控与通信
，涉及一种针对调相信号的基于非线性变换的载波频偏估计及跟踪方法。
技术介绍
在卫星测控和通信过程中，调相信号一般分为两相BPSK调制、功率平衡四相QPSK调制、功率非平衡四相UQPSK调制等等。针对两相BPSK和平衡四相QPSK，目前大多采取直接FFT频谱分析方法，实现射频载波的捕获过程。因载波调相后包络频谱呈开口向下的“抛物线”特性，包络谱峰值点分散(较平坦频谱顶部)使得很难通过查找谱线最大值方法完成准确地载波频偏估计，频偏估计误差大。同时，直接FFT频谱分析针对的是去调制符号前宽带包络谱，奈奎斯特带通采样定理决定了采样频率fs必须不小于2倍的符号速率fb。较高的采样频率决定了FFT频谱分析的信噪比较低，捕获精度差。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提出了一种基于一次或多次非线性变换的调相信号载波频偏最佳估计方案，通过非线性变换达到去符号的目的，使得包络谱变为单载波信号，然后通过直接FFT谱分析即可满足快速、准确地载波频偏估计及跟踪的型号需求。本专利技术的技术解决方案是：一种基于非线性变换的调相信号频偏估计及跟踪方法，步骤如下：(1)进行中频调相信号的带通数据采样，当采样频率fs满足阈值范围时，记录此时采样频率fs对应的采样信号Xa；(2)利用数控振荡器NCO将采样频率fs划分为2路正交的本地载波信号，分别和采样信号Xa乘积混频后，形成数字下变频后的零中频信号，即Xi和Xq；(3)滤除零中频信号Xi和Xq所含带外系统时钟、谐波高阶的频率分量，形成低通滤波后的零中频信号，即Yi和Yq；(4)判断调制信号的类型，若为两相BPSK调制信号和非平衡四相UQPSK调制信号，则对低通滤波后的零中频信号Yi和Yq进行一次非线性变换，得到线性变换后的信号Zi和Zq；若为平衡四相QPSK调制信号，则对低通滤波后的零中频信号Yi和Yq进行二次非线性变换，得到线性变换后的信号Zi和Zq；所述一次非线性变换的具体公式为：所述二次非线性变换的具体公式为：(5)对非线性变换后信号Zi和Zq进行Nb次累加求和，得到FFT频偏估计用零中频信号，即Vi＝∑Zi/Nb和Vq＝∑Zq/Nb，其采样频率降至fg＝fs/Nb；Nb为正整数且满足使得降速后的采样频率fg不小于2倍的载波多普勒频偏；(6)对多次累加后得到的信号Vi和Vq进行傅立叶FFT频谱分析，得到频率峰值最大点Mmax，则载波频偏估计值△f1＝Mmax*△fm,其中频率分辨率△fm＝fg/Nf，且△fm小于跟踪环路带宽B，采样频率fs与FFT点数Nf乘积小于系统捕获时间T0。利用载波频偏估计值△f1对数控振荡器NCO进行补偿和校正；(7)利用△f1补偿和校正后的数控振荡器NCO重复步骤(1)至(3)得到低通滤波后的零中频信号Yi和Yq，对其进行逆调制环完成信号的载波鉴相，得到载波鉴相结果W，具体公式为：二相BPSK：W＝Sign(Yi)*Yq/Sqrt(Yi2+Yq2)；平衡四相QPSK：W＝Sign(Yi)*Yq-Sign(Yq)*Yi/Sqrt(Yi2+Yq2)；非平衡四相UQPSK：W＝Sign(Yi)*Yq-ξ*Sign(Yq)*Yi/Sqrt(Yi2+Yq2)；其中ξ为测距/数传功率比；(8)将载波鉴相结果W通过低通滤波器后得到基于逆判决的相位差估计值△f2，直接补偿至数控振荡器NCO；利用相位差估计值△f2补偿和校正后的数控振荡器NCO重复步骤(7)至(8)直至数控振荡器NCO输出的两路正交的本地载波信号和输入信号Xa达到同频同相。所述步骤(1)中采样频率fs满足阈值范围在实信号采样的情况下确定方法为：其中Nc为满足上式且满足fs>2*(fH-fL)的最大正整数，fc为中频载波频率，fb为基带调制符号速率，频率带限信号X频率限制在(fL，fH)范围内，fL为左边带频率点，fL＝fc-fb，fH为右边带频率点，fH＝fc+fb。本专利技术与现有技术相比的有益效果是：1)、针对调制信号进行一次或两次非线性变换使得宽带调制包络谱变为单载波，峰值明显，频偏估计及跟踪误差小。2)、经非线性变换后的单载波采样信号不受基带调制符号速率的限制，可通过多次累加方式进一步提高输入信号的信噪比，谱线清晰，频偏估计及跟踪精度高。3)、基于平衡四相调制的逆判决反馈环的载波鉴相公式，提出了采取功率因子进行补偿的逆判决反馈环修正算法，首次解决了非平衡四相UQPSK调制的载波跟踪问题。附图说明图1为本专利技术涉及原理框图；图2为本专利技术的流程图。具体实施方式下面结合附图对本专利技术的具体实施方式进行进一步的详细描述。本专利技术的原理框图如图1所示，频偏估计方法由数字下变频、低通滤波器、非线性变换、FFT频谱分析、NCO补偿与闭环、载波鉴相等模块组成。其中A为输入信号幅度，I和Q为输入信号的调制基带数据，A0为本地产生的两路正交信号的幅度，ωi和ω0分别是输入信号和本地载波的角频率。本专利技术提出的一种基于非线性变换的调相信号频偏估计方法，计算流程如图2所示，实施步骤如下：(1)根据中频载波频率fc(42Mhz)、基带调制符号速率fb(I路长码3.069Mcps，Q路高速数传5Mbps)，利用奈奎斯特带通采样定理，结合系统时钟f0(36Mhz)，确定AD采样频率fs(36Mhz)，并完成中频调相信号的带通数据采样，得到采样信号Xa；(2)采样频率fs利用数控振荡器NCO产生2路正交的本地载波信号，和采样信号Xa乘积混频后，形成数字下变频后的零中频信号，即Xi和Xq；所述数控振荡器NCO是指根据各个相位角度θ计算相应的正弦值和余弦值，并以相位角度θ作为地址存储各相位的正弦和余弦样本，其中本振信号频率fn(6Mhz)用相位步进控制字Ka实现，满足以下公式：fn＝Ka*fs/2Na，数控振荡器的相位累加宽度为Na(32位)。(3)利用低通滤波器，滤除零中频信号Xi和Xq所含带外系统时钟、谐波高阶的频率分量，形成低通滤波后的零中频信号，即Yi和Yq；经过步骤(1)～(3)，星上设备实现带通采样、数字下变频和低通滤波等，为载波频偏估计和相位跟踪完成了前端信号调理。(4)进行载波频偏估计。若背景型号为两相BPSK和非平衡四相UQPSK调制信号，对低通滤波后的零中频信号进行一次非线性变换，得到一次非线性变换信号Zi和Zq；若平衡四相QPSK调制信号，则对低通滤波后的零中频信号进行两次非线性变换，得到二次非线性变换信号Zi和Zq。经非线性变换后，宽带调制包络谱变为射频“单载波”信号，或“单载波+IQ调制谱”信号，信号峰值明显；所述非线性变换是指对低通滤波后的零中频信号Yi和Yq，进行乘积或平方等运算，即得到Zi和Zq，区别于加减线性运算，称之为非线性变换：此公式运算几次，即为几次非线性变换，也就是说上述公式应用两次即为二次非线性变换。本背景型号由于是非平衡四相UQPSK调制，进行一次非线性变换，变换后信号特征呈现“单载波+IQ调制谱”特性。(5)为了提高信噪比，降低采样速率，对非线性变换后信号Zi和Zq进行Nb次(Nb＝16)累加求和，其中Nb为正整数且须满足使得降速后的采样频率fg不小于2倍的载波多普勒频偏，本文档来自技高网...

【技术保护点】
一种基于非线性变换的调相信号频偏估计及跟踪方法，其特征在于步骤如下：(1)进行中频调相信号的带通数据采样，当采样频率fs满足阈值范围时，记录此时采样频率fs对应的采样信号Xa；(2)利用数控振荡器NCO将采样频率fs划分为2路正交的本地载波信号，分别和采样信号Xa乘积混频后，形成数字下变频后的零中频信号，即Xi和Xq；(3)滤除零中频信号Xi和Xq所含带外系统时钟、谐波高阶的频率分量，形成低通滤波后的零中频信号，即Yi和Yq；(4)根据调制信号的类型，对低通滤波后的零中频信号Yi和Yq进行非线性变换，得到线性变换后的信号Zi和Zq；(5)对非线性变换后信号Zi和Zq进行Nb次累加求和，得到FFT频偏估计用零中频信号，即Vi＝∑Zi/Nb和Vq＝∑Zq/Nb；(6)对多次累加后得到的信号Vi和Vq进行傅立叶FFT频谱分析，得到频率峰值最大点Mmax，则载波频偏估计值△f1＝Mmax*△fm,其中△fm为频率分辨率，利用载波频偏估计值△f1对数控振荡器NCO进行补偿和校正；(7)利用△f1补偿和校正后的数控振荡器NCO重复步骤(1)至(3)得到低通滤波后的零中频信号Yi和Yq，对其进行逆调制环完成信号的载波鉴相，得到载波鉴相结果W；(8)将载波鉴相结果W通过低通滤波器后得到基于逆判决的相位差估计值△f2，直接补偿至数控振荡器NCO；利用相位差估计值△f2补偿和校正后的数控振荡器NCO重复步骤(7)至(8)直至数控振荡器NCO输出的两路正交的本地载波信号和输入信号Xa达到同频同相。...

【技术特征摘要】
1.一种基于非线性变换的调相信号频偏估计及跟踪方法，其特征在于步骤如下：(1)进行中频调相信号的带通数据采样，当采样频率fs满足阈值范围时，记录此时采样频率fs对应的采样信号Xa；(2)利用数控振荡器NCO将采样频率fs划分为2路正交的本地载波信号，分别和采样信号Xa乘积混频后，形成数字下变频后的零中频信号，即Xi和Xq；(3)滤除零中频信号Xi和Xq所含带外系统时钟、谐波高阶的频率分量，形成低通滤波后的零中频信号，即Yi和Yq；(4)根据调制信号的类型，对低通滤波后的零中频信号Yi和Yq进行非线性变换，得到线性变换后的信号Zi和Zq；(5)对非线性变换后信号Zi和Zq进行Nb次累加求和，得到FFT频偏估计用零中频信号，即Vi＝∑Zi/Nb和Vq＝∑Zq/Nb；(6)对多次累加后得到的信号Vi和Vq进行傅立叶FFT频谱分析，得到频率峰值最大点Mmax，则载波频偏估计值△f1＝Mmax*△fm,其中△fm为频率分辨率，利用载波频偏估计值△f1对数控振荡器NCO进行补偿和校正；(7)利用△f1补偿和校正后的数控振荡器NCO重复步骤(1)至(3)得到低通滤波后的零中频信号Yi和Yq，对其进行逆调制环完成信号的载波鉴相，得到载波鉴相结果W；(8)将载波鉴相结果W通过低通滤波器后得到基于逆判决的相位差估计值△f2，直接补偿至数控振荡器NCO；利用相位差估计值△f2补偿和校正后的数控振荡器NCO重复步骤(7)至(8)直至数控振荡器NCO输出的两路正交的本地载波信号和输入信号Xa达到同频同相。2.根据权利要求1所述的一种基于非线性变换的调相信号频偏估计方法，其特征在于：所述步骤(1)中采样频率fs满足阈值范围使用奈奎斯特带通采样定理确定。3.根据权利要求2所述的一种基于非线性变换的调相信号频偏估计方法，其特征在于使用奈奎斯特带通采样定理的具体方法为：其中Nc为满足上式且满足fs>2*(fH-fL)的最大正整数，fc为中频载波频率，fb为基带调制符号速率，频率带限信号X频率限制在(fL，f...

【专利技术属性】
技术研发人员：赵鸿，李静澜，余晓川，孙重磊，
申请(专利权)人：西安空间无线电技术研究所，
类型：发明
国别省市：陕西,61