北邮通信工程MATLAB双音多频.docx

N=205; a1=0.5204; b1=1,-1.7077,1; a2=0.5686; b2=1,-1.6453,1; a3=0.6203; b3=1,-1.5687,1; a4=0.6736; b4=1,-1.4782,1; a5=0.8131; b5=1,-1.1641,1; a6=0.8671; b6=1,-0.9964,1; a7=0.9168; b7=1,-0.7986,1; a8=0.9652; b8=1,-0.5685,1; x=zeros(1,N); x(1)=1; n=0:N-1; y1=filter(a1,b1,x); y2=filter(a2,b2,x); y3=filter(a3,b3,x); y4=filter(a4,b4,x); y5=filter(a5,b5,x); y6=filter(a6,b6,x); y7=filter(a7,b7,x); y8=filter(a8,b8,x); z1=y1+y5; %数字1 z2=y1+y6; %数字2 z3=y1+y7; %数字3 z4=y2+y5; %数字4 z5=y2+y6; %数字5 z6=y2+y7; %数字6 z7=y3+y5; %数字7 z8=y3+y6; %数字8 z9=y3+y7; %数字9 z10=y4+y2; %数字0 z11=y4+y5; %符号* z12=y4+y7; %符号# z13=y1+y8; %A z14=y2+y8; %B z15=y3+y8; %C z16=y4+y8; %D stem(n,z1) stem(n,z2) stem(n,z3) stem(n,z4) stem(n,z5) stem(n,z6) stem(n,z7) stem(n,z8) stem(n,z9) k=18,20,22,24,31,34,38,42; p=goertzel(z1,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z2,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z3,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z4,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z5,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z6,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z7,k+1); val=abs(p); stem(k,val) k=18,20,22,24,31,34,38,42; p=goertzel(z8,k+1); val=abs(p); stem(k,val) k