一种适用于数字调制信号的载波频率估计方法.docxVIP

一种适用于数字调制信号的载波频率估计方法 1 单承载频率估计 在瞬态分析信号时，为了计算信号的瞬态相位和瞬时频率，有必要精确测量波形频率的精确值。在调制识别中,许多特征是从瞬时幅度、瞬时相位和瞬时频率中提取的,当然也需要载波频率的精确值。因此,在进行瞬态分析和调制识别之前必须知道载波频率的精确值,但是作为通信的第三方是不可能预先知道载波频率的精确值的,只能通过载波频率估计来获取载波频率。目前估计载波频率的方法有4种。其中周期图法把信号周期图最高峰的位置作为载波频率,只适用于载波分量大的信号,如低调制深度的AM和MASK信号,不适用其它调制方式的信号;频率居中法把信号频谱的中心频率作为载波频率,对于对称信号效果尚好,而对于非对称信号效果不好。其缺点是精度较低,特别是在有噪声和信道中有信道衰落时;零交叉法基于射频信号的零交叉,其缺点是对噪声特别敏感,但可以通过定位微弱信号区对该方法进行改进;最小均方算法对于相位连续的信号效果不错,而对于有相位跳变点的信号则误差很大。 为对数字调制信号进行实时瞬态分析,本文从实际应用出发,从含有线性相位分量的瞬时相位求出的瞬时频率中存在载波频率分量,从中估计出载波频率。在信噪比较高的情况下,对于各种数字调制方式的估计精度均可接受。 2 单通道总频率和同化作用 已知信号的瞬时频率fis(t)由其瞬时相位?is(t)的微分求得 fis(t)=12π?d?is(t)dt(1) 该微分可由中心有限差分求得 fis(i)=fs4π[?is(i+1)-?is(i-1)](2) 式中:fs——采样率。若用?ic(i)表示载波的相位分量,?in(i)表示由窄带高斯白噪声引入的干扰相位分量(下面所提噪声均指窄带高斯白噪声),?uw(i)表示混合无卷叠瞬时相位, 如果载波频率fc是精确已知的,则 ?uw(i)=?is(i)+2πfcfs?i+?in(i)(3) 对式(3)求中心有限差分,并把式(2)代入得 fi(i)=fs4π[?uw(i+1)-?uw(i-1)]=fis(i)+fc+fs4π[?in(i+1)-?in(i-1)](4) 令fin(i)=fs4π[?in(i+1)-?in(i-1)](5) 把式(5)代入式(4),得 fi(i)=fis(i)+fc+fin(i) (6) 式中:fis(i)——信号的瞬时频率,fin(i)——由噪声引入的干扰瞬时频率。由于?in(i)服从(-π,π)区间上的均匀分布,因此fin(i)的均值为零,即序列{fin(i)}中不存在直流分量,而且对于数字调制信号,序列{fis(i)}中也不存在直流分量,所以序列{fi(i)}中的直流分量即为所要求的载波频率。 设序列{fi(i)}的长度为N,利用最小二乘法,令平行于x轴的直线y=fc代表未知的载波频率,使平方和 Q=Ν∑i=1[fi(i)-fc]2(7) 最小,则Ω达到最小值必满足条件?Q/?fc=0,而 ?Q?fc=2Ν∑i=1[fc-fi(i)]=2Νfc-2Ν∑i=1fi(i)(8) 从而得 fc=1ΝΝ∑i=1fi(i)=E[{fi(i)}](9) 由上式可得序列{fi(i)}的数学期望即所求的载波频率。 3 数字调制方法 在实际信号中,存在瞬时相位和频率估计对噪声更加敏感的微弱信号区,这相当于接收到一个信噪比很低的信号,而且对微弱信号区进行载波频率估计,估计的精度很低。因此,载波频率估计要尽量在非微弱信号区中进行。这就需要对微弱信号样品定位,方法很简单, 对每个样品检验不等式 ai(i)≤at(10) 式中:ai(i)——归一化瞬时幅度,at——信号段的微弱区所选定的门限。由此定位了非微弱信号区和微弱信号区。 另外,数字调制会在码元转换处引起瞬时频率序列{fi(i)}的跳变,如图1所示,从而降低了载波频率的估计精度。对于这一问题,通过对瞬时频率序列{fi(i)}进行中值滤波,可以消除ASK和PSK信号的瞬时频率序列{fi(i)}中的跳变,使滤波后的瞬时频率序列{fi(i)}中只含载波频率从而显著提高了ASK和PSK信号的载波频率估计精度。但是对于FSK信号,瞬时频率序列{fi(i)}的变化反应了调制信号的变化,所以中值滤波不能消除且不应消除FSK信号的瞬时频率序列{fi(i)}中的跳变,因此中值滤波对于FSK信号的载波频率估计精度的提高效果一般。 选择非微弱信号区进行载波频率估计,并对序列{fi(i)}进行中值滤波,可以显著地提高载波频率的估计精度。 4 i相位卷叠并进行中值滤波器设计 步骤1求实信号s(i)的解析信号z(i) z(i)=s(i)+j?s(i)(11) 式中:s?(i)——实信号的希尔伯特变换。 步骤2计算z(i)的归一化瞬时幅度ai(i) ai(i)=s(i)2+s?(i)2/amax(12