基于离散傅里叶变换的新型数字解调器

基于离散傅里叶变换的新型数字解调器

0数字视频监控系统对于时间连续的信号，时间间隔转换后的信号可以转换为频率范围，从而在频域中进行分析。对于时间离散的信号,信号是用序列表示的,相应的有离散的傅里叶变换算法,经过变换后的结果也是一个序列,该序列的各元素均为虚数。如果对该序列的每个元素取模,得到一组实数序列,这组实数序列与信号含有的各次谐波的频率是一一对应的,即为某一谐波的幅值,当前这种技术多用在对信号的分析上。在通信技术中,由于数字通信有抗干扰性强,易于加密等优点而得到广泛的应用。数字通信中信号的调制方式有2ASK、2FSK和2PSK(二进相移键控)三种,其中的2ASK和2FSK可以认为与两种频率的交替变化有关,如果用离散傅里叶变换先对调制信号进行分析,然后再判断不同频率对应的幅值是否满足一定条件,进而判断出该位的数字逻辑,就可以很轻松的实现对信号的解调。1打造参数频率和旋转因子一个周期为2l的函数f(x)可用傅里叶基数展开为:将连续函数的傅里叶基数展开式(1)离散化。为了离散化式(1),在周期区间(0,2l)上等间隔的取N个点,取样间隔为∆t,那么,这里要注意dx→∆t。则f(x)的离散化序列为{x0,x1,x2,,xN-1},且,由此式(1)的离散化形式为:在这里对做一个变换,对其分子分母同乘以∆t后变为:,由此可得出第kx项为一个正弦和一个余弦周期函数之和,其频率同为:,其中T为所取序列总的时间长度。随着k的增大,三角函数的频率逐渐增加,周期逐渐减小,其周期为:时,谐波的频率最大为:,该频率称为Nyquist频率,当k从N/2取到N时,其结果与k从0取到N/2是镜像对称的。现在将式(2)的各次谐波写成如下形式:其中:,为k次谐波的振幅;,为k次谐波的初相。在这里如果将ak和bk分别表示为一个虚数的实部和虚部,结合欧拉公式,则可直接得到离散傅里叶变换的复数形式:其中,称为旋转因子。X(k)是一个虚数,与kc的关系为:。用式(4)求离散傅里叶变换时有很多快速傅里叶变换(FFT)算法,如时域抽取法基2FFT,这里不再详细分析快速算法。从以上的分析可以看出,离散傅里叶变换的计算量与N的大小有关,N取值增大,计算量会成倍的增加。快速傅里叶变换算法与一般离散傅里叶变换算法的运算量相比较,它极大的减少了运算量,使离散傅里叶变换的计算时间大大缩短,为数据的实时处理提供了可能性,两种算法运算量的比较如图1。从图1中可以看出,采用FFT算法可以极大地减少运算量,所以在后面的仿真中也是采用快速傅里叶变换函数来做离散傅里叶变换,图1中纵坐标数值要乘以104。22函数相乘相乘对于2ASK来说,其调制的数学原理为:0e(t)=s(t)cosωt,其中0e(t)就是调制后的信号,s(t)是随时间变化的二进制序列,就是调制信号;cosωt是时间连续的余弦函数,其频率由ω决定,是载波信号;因此2ASK原理可理解为:在不同时刻用二进制数字信号与一个余弦函数相乘,也就是用数字信号去调制模拟的载波信号。2FSK的数学原理为:e0(t)=s(t)cos(ω1t+θn)+cos(ω2t+ψn),可以看出有两个频率为1ω和ω2的余弦函数,如果某时刻数字信号为二进制“1”,则e0(t)=cos(ω1t+θn),则输出为角频率为1ω的余弦信号,如果某时刻数字信号为二进制“0”,则e0(t)=cos(ω2t+ψn),则输出为角频率为ω2的余弦信号,可见移频键控就是用两个幅度相同,而频率不同的余弦信号来表示二进制位信号,也就是用二进制位信号调制两个余弦信号。从以上的分析可以看出,2ASK和2FSK都是由二进制数字信号的位信号来控制两个不同频率的余弦信号输出的过程。对于2ASK来说,可以理解为另一个余弦信号频率为零。这样一来就可以利用离散傅里叶变换具有对信号频谱分析的特点来分析信号中所含有的频率及其振幅,确定该段信号是被“0”或“1”所调制,从而就可以把数字信号解调出来。3fsk信号相乘2ASK和2FSK信号的解调,总的来说有两种基本方法,相干解调法和非相干解调法(又叫包络检波法)。2ASK的相干解调法要对信号进行带通滤波,全波整流,低通滤波,抽样判决,这样才能得到二进制数字信号。非相干解调法不用全波整流,而是用本地产生一个频率、相位都与载波频率相同的余弦信号与接收到的信号相乘,其他部分相同。2FSK因为有两个频率,所以要把2ASK中解调的方法分两路进行。可见要对2ASK和2FSK信号进行解调,硬件电路还是比较复杂的,如果用离散傅里叶变换来解调,就会简化很多。它的基本原理很简单,原理如图2所示。首先要对接收到的信号放大,然后再进行采样(A/D转换),采样的点数就是第二部分离散傅里叶变换分析中的N,对N个点做离散傅里叶变换,最终将得到一组与频率对应的振幅序列ck。只要知道载波信号的频率fk和采样周期∆t,利用公式就可以计算出k(k取整数),在振幅序列中找到第k个数ck,判断ck是否大于给定的判决值,如果ck大于该值,说明采样的这段信号中含有载波频率fk,然后根据在调制时约定的是用“0”还是“1”调制的fk,就可以判断这段采样信号所承载的数字逻辑。从上面的设计可以看出,用离散傅里叶变换进行数字解调不需要对信号滤波,因为傅里叶变换本身就具有对信号频谱分析的功能,而且整个过程主要是用程序来实现,易于维护和改动,适用性强。4随机干扰及仿真通过对2ASK和2FSK信号解调仿真的实例来说明用离散傅里叶变换实现数字解调的效果。2ASK信号的载波频率为fk=1kHz,2FSK的载波频率为1f=1kHz,f2=1kHz数字调制信号都是随机产生的,采样周期∆t取0.00002s,采样点数对于2FSK信号来说N=128,对于2ASK信号来说N=256(取128个点完全可以解调出来,这里为了使仿真图形更直观,而选256个点),加入随机干扰信号,干扰信号最大振幅为0.8。这里只解调出8位数字二进制位以说明问题。按照以上的设计通过MATLAB做仿真,仿真结果如图3和图4。图3(a)、图4(a)为调制后信号,不含干扰;图3(b)、图4(b)为二进制数字调制信号,是随即产生的;图3(c)、图4(c)为加入随机干扰后的调制信号;图3(d)、图4(d)就是用傅里叶变换解调出来的数字信号,与原信号相同,而且经过多次仿真,并加长仿真数据,仿真结果都没有出现误码。5仿真的可行性分析通过以上的原理分析可以看出,应用离散傅里叶变换来解调2ASK和2FSK信号,原理上是完全可行的,通过仿