第二章 DFT和其快速算法.pptVIP

hk(n)的Z变换为 写出系统输出 由于其中的系数为复数，使用不便。对Hk(z)做进一步变换 可以推导出递推计算公式 这就是Goertzel算法。 此式展开后就是对应Hk(z)的差分方程。 计算步骤： 1、每输入一个样值，分别用迭代法计算出各个vk(n)。 3、对应的幅度平方值性为： 2、当n=0,1,...N-1时，即可计算出对应k的DFT值。 DTMF数字信号的生成 在DTMF信号合成时，每输入一个数字，合成100ms的DTLF信号，后面再加上100ms的间隔。采样频率规定为8kHz，这样，每输入—个数字，计算800点合成抽样值，在后面再加上800点的0值。 DTMF音调编码器 DTMF音调频率可由数字谐振器产生的正弦波构成。该谐振器由双极点滤波器来实现，该滤波器可以用以下差分方程描述： 其中k是一个常数，其定义为： k = 2 cos(2π×音调频率/ 采样速率) 由于DTMF拨号盘只需要少量音调，k的8个数值可以预先计算。例如，在8kHz采样速率下产生第一列音调(770Hz)所需的常数为： k = 2 cos(2π × 770 / 8000) = 2 cos(0.60) = 1.65 还有一个数值必须计算：振荡器开始运行所需的初始脉冲。很明显，若y(n-1)与y(n-2)都为0，随后的每个y(n)都将为0。为了启动振荡器，将y(n-1)设置为0，并将y(n-2)设置为 y(n-2) = -A ×sin(2π × 音调频率/ 采样速率) 在我们的实例中，假定需要单位正弦波，则上式可简化为： y(n-2) = -1 × sin(2π × 770 / 8000) = -sin(0.60) = -0.57 close all;clear all y=zeros(1,8000);fs=8000;n=8000 tone_f=[697,770,852,941,1209,1336,1477,1633]; k = 2*cos(2*pi*tone_f/fs); y_2=-1*sin(2*pi*tone_f/fs); for m=1:8 y(1)=-1*y_2(m);y(2)=k(m)*y_2(m); for n=3:8000; y(n)=k(m)*y(1)-y(2); a=y(1);y(2)=a;b=y(n);y(1)=b; end h=fftshift(fft(y)/n); x_co=[-pi:2*pi/n:pi-2*pi/n];%最大值π对应频率为最高采样频率(单位为弧度)的一半 subplot(2,4,m);stem(x_co.*fs/2/pi,abs(h));grid on;%频率坐标显示(单位为Hz) sound(y,fs); end 信号的高频和低频分量的最小间隔为73Hz，选取以每200点为一段，分别用Goertzel算法。 计算每一段的8个频率点的DFT的平方值。这样，就可以将相邻频率点分开。 在计算出每—段的8个频率点的幅度平方值后，设定一个阈值t，对得到的值进行判断， 苦某一频率点的幅度平方值大于阀值t，判为有效，否则判为无效。再根据有效频率点的组合，给出最后检测出的数字。 最小频差：770-697=73Hz，因此，最短采样时间应大于1/73=13.7毫秒。为了留有余量，一般选40ms，分辨率为25Hz。 K的选择：ωk＝2πk/N，因为ω=ΩT，所以对应模拟频率点为fk＝采样频率×k/N，选择合适的N得到规定的各个fk。反过来讲，对需分析的8个fk而言对应一个共同的最接近的整数N。 结论：采样频率＝8000Hz；采样点数＝205；采样时间＝0.04秒 8个基频 对应的k值 最接近的整数k值 绝对误差 697 17.861 18 0.139 770 19.531 20 0.269 852 21.833 22 0.167 941 24.133 24 0.133 1209 30.981 31 0.019 1336 34.235 34 0.235 1477 37.848 38 0.152 1633 41.846 42 0.154 当N=205时计算出的k值 close all;clear all; tm=[1,2,3,65; 4,5,6,66; 7,8,9,67; 42,0,35,68]; n=205;k=[18,20,22,24,31,34,38,42]; f1=[697,770,852,941];%行向量频率 f2=[1209,1336,1477,1633];%列向量频率 tn=input(输入6位电话号码=); tnr=0;%接收端电话号码初值为0 for m=1:6