基于df的一种多模式调解调算法.docxVIP

基于df的一种多模式调解调算法 0 信号解调方法 近年来，作为软件无线电技术的重要组成部分，数字解调理论的研究引起了人们的关注。为了满足实时数字化解调技术对算法要达到计算量小,误码率低的要求,信号的解调方法上也有了很大的改进。对于实际中应用较为广泛的FSK信号,利用短时离散傅里叶变换(STFT)来实现其数字解调成为一个研究热点。文献介绍了使用DSTFT解调FSK信号的方法,同原有的方法相比,具有算法简单,抗干扰能力强的优点。但它们都是以理想的FSK信号为解调对象,没有考虑到实际中接收到的FSK信号码元波形存在着过渡区,并且计算量也很大(由于使用移动窗函数等原因),所以上述文献中介绍的解调方法在实际应用中都会受到限制。 1 基于dft的数字解调算法 理想的FSK信号在码元切换时可能发生载频跳变,使得旁瓣功率较大,需要经过限带滤波后才能发射。而滤波作用使得码元跳变处的高频成分被滤除,所以实际接收到FSK信号的码元波形分为两种区域,码元中间部分是稳定区,前、后部分为过渡区。稳定区的波形比较接近于载波频率的正弦波,过渡区波形的频率(相位)不稳定,且幅度明显降低。显然,稳定区非常适合对码元进行解调判决。数字化解调算法不应忽略码元波形过渡区的存在。 DFT解调原理如下:对接收到的FSK信号进行采样后,通过两路数字带通滤波器,分别滤出两个载波频率成分;然后对每个码元分别用波形稳定区内若干个载波周期的采样值进行DFT,就可以得到两个幅值,通过比较它们的大小来恢复数字码元。 为方便讲解,设FSK信号的载频为1f和f2,采样频率为fs,采样后的样值序列为x(i),n(i)是高斯白噪声序列,则加噪声信号序列为s(i)=x(i)+n(i)。如果FSK信号的第i个样值属于1码元,则x(i)=cos(2πf1*i/fs),否则x(i)=cos(2πf2*/fs),基于DFT的数字化解调算法如图1所示。 图1中,数字带通滤波器BPF1、BPF2的中心频率分别为1f和f2。对滤波输出的两路信号分别进行DFT,求出同一个码元的两个幅值1A和2A。求1A的计算公式如下: 其中s(i)是位于码元稳定区中的采样值。求2A时只需要把式(1)中的1f改写为f2。 码元判决:若1A2A,码元判为1,否则判为0。由于是比较1A和2A的相对大小,可以省去式(2)中的开平方计算,从而减少计算量。 2 同频同相插值a 位同步又称为时钟同步。为了确定码元的切换时刻,可以分别对两个数字带通滤波器输出信号的每一个(或几个)载波周期的采样值进行DFT,得到两个幅值序列1B(k)和2B(k)。为了克服噪声的影响,可以对1B(k)和2B(k)进行适当的滤波。为了细化位同步的时间刻度,对1B(k)和2B(k)序列进行插值,求出每一个采样时刻的1B(k)和2B(k)。在1码元,1B(i)应大于2B(i),在0码元,1B(i)应小于2B(i)。当B1(i)-B2(i)的值由正数转变为负数或者由负数转变为正数时,就可以认为此处为码元的同步点,即B1(i)-B2(i)的极性改变时刻对应于码元切换时刻。由求得的码元切换时刻就可以调整再生时钟的相位,从而解决再生时钟与接收码元之间的同频同相问题,使收到的码元可以获得正确的判决。在调整再生时钟的相位时,可以采取适当的保护措施(与基带传输系统的位同步调整问题相同,此处不再讨论)。用整数个载波周期的采样值来进行DFT是为了减少“功率谱泄漏”。如果一个载波周期内的采样点数不是整数,则可以用连续几个载波周期内的采样值来进行DFT。比如当f1/fs=8.5时,可以取N=17。实际上在位同步维持状态下,只需要在预定的位同步时刻前后一段时间内进行1B和2B的计算并判断B1(i)-B2(i)的极性改变时刻即可。 3 带通滤波器同位dft 为了验证该解调方案的正确性,采用Matlab对整个解调及位同步过程进行仿真。设FSK信号的码速率为1k Baud,载波频率1f=10k Hz,f2=15k Hz,采样频率fs=120k Hz,数字带通滤波器1和2的带宽均取为2 k Hz,码元个数取为10 000个。图2中的(a)和(b)分别是经过带通滤波器1和2滤波后的波形。由图2可知,两路滤波后的波形幅度都存在过渡区和稳定区,根据该特点,对两路信号采样序列分别进行DFT运算后所提取的幅值如图3所示,这说明按照上文(2位同步方法)中所提到的方法即可确定出码元的位同步点。 改变噪声的功率值并进行码元判决,最终可得到这10000个码元在不同信噪比情况下的误码率,如图4所示。图4中给出了对2FSK信号采用DFT进行解调的误码率与传统的相干解调及非相干解调的误码率比较,由图4可知该解调方案提高了通信系统的抗噪声性能,这是因为解调时只使用了码元波形稳定区内的采样数据,更有利判决,且采用DFT进行解