课题介绍2.PPT

二、DTMF信号的解码 1.FFT算法 (1)信号接收； (2)用FFT对信号进行频谱分析； (3)还原为数字键； (2)中： 几个问题： 1）取样频率为何取8KHz？ 答：语音信号的最高频率fc =4KHz,因为取样频率fs≥2fc才能保证取样后的信号不失真， 所以fs=2×4k=8kHz(工业标准) 2）为何取200点为一帧做FFT? 答：为在频谱图中分辨出不同的频率分量，于是对信号取200点为一帧，则频谱分辨率F=fs/N=8000/200=40Hz73Hz(表1中任意两频率的最小间隔)。 每个信号(+间隔)占1600/200=8帧 3）为了对每帧(200点）信号作FFT， 问:对每帧信号作多少点FFT? 即N=? 4)每帧信号幅频谱仅画64点(N/4)。(存入r) r:8*n行，64列 (n---信号的个数,一帧占一行) 答：因为信号x为实数序列，所其幅频谱|y|具有偶对称性，于是，幅频谱可以仅画N/2点，其中第N/2点对应实际频率为fs/2=4KHz,(书P59 图3-19) 又因为DTMF 信号中最高频率为1633Hz，小于2KHz（fs/4），因此，这里只画N/4=64点。 5)频谱横坐标为频率点k k=fk/F=(fk/fs)*N; f----实际频率。 DTMF信号是两个正弦波的迭加，它的幅频谱就是两根谱线，谱线的横坐标就是该信号的两个频率点KL和KH 。 例：按键 “1”对应DTMF信号的频率 fL=697Hz,fH=1209Hz。(fs=8kHz,N=256) 问:(1)KL=? KH=? (2)运行该程序，幅频谱图中产生什么 现象？ 6)用阀值法消除频谱泄漏现象。(存入c) 答：由于信号x是有限长的，这就相当于对无限长的信号加矩形窗，所以在频谱图中必然会出现频谱泄漏现象,使信号能量散布到其他谱线位置。为此，在程序中应选择一适当阀值，将出现在这两条谱线周围的幅度较小的谱线消除（置0），最后，将处理后的幅频谱数据存入数组c中。 (3)中： 1)用查表(sm)法将频率点转换为对应数字键。 在幅频谱图中，频率轴的定标方式为频率点K而不是实际频率f,转换关系为：K=f/F,因此，数字键 0-9对应频率点如下表所示： 频率点矩阵sm= [31,44;23,40;23,44; 23,48;26,40;26,44; 26,48;28,40;28,44; 28,48] (对应数字键0~9) 数组c中不等于0的下标就是各信号的频率点，KL,KH,查表sm，即可将各DTMF信号还原为相应的数字键。 *用到的函数： i)find(c)----找出c中≠0的数据的下标。 ii)nnz(c)----找出c中≠0的数据的个数。 2)查找过程。 从sm的第一行开始查，查到, 则数字键AN=这时sm的下标 - 1; (以数字键“1”为例)跳出本级循环。 FFT算法解码程序: A=wavread(‘D2.wav’); %接收到的DTMF信号 A=A’; %A转置 subplot(212); plot(A);%绘图 N=256; for s=1:8*n %对每帧信号作N=256点 的FFT R=A(200*(s-1)+1:200*s); y=fft(R,N); c(s,:)=abs(y(1:64)); %幅频谱取64点，存入c r(s,:)=c(s,:); %r=c z=find(c(s,:)40); %消除频谱泄漏现象（阀 值=40），结果再存入c c(s,z)=zeros(size(z)); end sm=[31 44;23 40;23 44;23 48;26 40;26 44;26 48;28 40;28 44;28 48]; %0-9对应的频率点表矩阵 for i3=1:8*n b=nnz(c(i3,:)); %b:c中≠0数据的个数； if b==2 %若b=2,则c为信号幅频谱，其 ≠0的下标q1即为频率点。 q1=find(c(i3,:)); for i4=1:10