7-2 第七章FIR设计-2要点.pptVIP

　　这就是Goertzel算法。每输入一个样值，分别用迭代法计算出各个vk(n)，每次计算只需进行两次实数乘和两次实数加。 当n=N时，用式（8-3）计算出各yk(n)的值，这些值就是所要求的各点的DFT值。 　　另外，考虑到最后对信号进行判断时需要幅度值|X(k)|或幅度的平方值|X(k)|2，因此对式（8-3）进行变换： 这样， 算法中就避免了复数运算，最后对每个值点再进行四次乘法和二次加法运算即可。 　　3. DTMF信号产生与检测算法的实现 　　在检测信号时，为避免信号前后的0值带来的对DFT计算的影响，同时考虑到 DTMF 信号的高频和低频分量的最小间隔为73 Hz， 故选取每200点为一段， 分别用Goertzel算法计算每一段的8个频率点的DFT的平方值，然后设定一个阈值对得到的值进行判断（若某一频率点的幅度平方值大于阈值， 判为有效，否则为无效），再根据有效频率点的组合给出最后检测到的数字。 　　为了模拟接收到DTMF信号在数据流中的开始位置的检测效果，可以在产生的采样数据前加上不同个数的0。 　　具体的程序如下（首先输入要产生和检测的信号个数， 然后输入加在抽样数据前的0的个数，最后输入单个数字）： function DTMF %用来模拟产生和检测DTMF信号 s=input(′请输入检测的数字个数： ＼n′); l=input(′请输入加在数据前面的0的个数： ＼n′); out=zeros(1600*s+1,1); %设置输入样值数组 tab=［941 1336;697 1209;697 1336;697 1447;770 1209;770 　　　　　1336;770 1477;852 1209；852 1336;852 1477; for I=1:s k=input(′请输入数字（0~9）： ＼n′); 　for j=1:800 x=tab(k+1,1,)/8000; y=tab(k+1,2)/8000; z=sin(x*j*2*pi)+sin(y*j*2*pi); 　out(1600*(i-1)+j+1,1)=14000*z; %合成信号采样值 　　out(1600*(i-1)+j+800+1,1)=0 end end out=out./28000; wavwrite(out,8000,16,DTMF1); %写入DTMF1.WAV文件 w=zeros(8); w(1)=2*cos(2*pi*17/200); w(2)=2*cos(2*pi*19/200); w(3)=2*cos(2*pi*21/200); w(4)=2*cos(2*pi*23/200); w(5)=2*cos(2*pi*130/200); w(6)=2*cos(2*pi*33/200); w(7)=2*cos(2*pi*36/200); w(8)=2*cos(2*pi*40/200); result=zeros(8*s,8); a=wavread(′DTMF′); 读出DTMF1.WAV文件中的数据 for i=1:8*s; b=zeros(8,3); for j=1:200; for k=1:8 b(k,3)=w(k)*b(k,2)-b(k,1)+a((i-1)*200+j); b(k,1)=b(k,2); b(k,2)=b(k,3); end end for k=1:8 result(i,k)=b(k,2)^2+b(k,1)^2-w(k)*b(k,1)*b(k,2)-600; if(result(i,k)0) result(i,k)=0; end end end for i=1:8*s 判断数字值 a=nnz(result(i,:)) if(a= =2) q=find(result(i,:)); if(q(1)= =1) if(q(2)= =5) res(i)=1; elseif(q(2)= =6) res(i)=2; elseif(q(2)= =7) res(i)=3; end elseif(q(1)= =2) if(q(2)= =5) res(i)=4; elseif(q(2)= =6)