bin.rar_dtmf_信号检测_读取八位DTMF

function DSP(k) N=8;alpha=(N-1)/2; wc=2*697*pi/8000; w0=2*1477*pi/8000; n=(0:7); i=n-alpha; i=i+(i==0)*eps; h1=2*(sin((i).*wc)).*(cos((i).*w0))./((i).*pi); y=0; biaozhi=0; [M,N]=size(k); for m=0:N:205 for j=1:200 y=y+k(1,j+m)^2; end y=0.5*y; if(y>=175)&(biaozhi==0) biaozhi=1; t=[18 20 22 24 31 34 38 42]; g1=goertzel(x(m:m+205),t+1); g1=abs(g1); if((g1(18)>=10)&(g1(31)>=10)) disp(1); elseif((g1(18)>=10)&(g1(34)>=10)) disp(2); elseif((g1(18)>=10)&(g1(38)>=10)) disp(3); elseif((g1(20)>=10)&(g1(31)>=10)) disp(4); elseif((g1(20)>=10)&(g1(34)>=10)) disp(5); elseif((g1(20)>=10)&(g1(38)>=10)) disp(6); elseif((g1(22)>=10)&(g1(31)>=10)) disp(7) ; elseif((g1(22)>=10)&(g1(34)>=10)) disp(8) ; elseif((g1(22)>=10)&(g1(38)>=10)) disp(9) ; elseif((g1(24)>=10)&(g1(34)>=10)) disp(0) end if(y<175) biaozhi=0; end end end end