北邮通信工程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=[18,20,22,24,31,34,38,42];>> p=goertzel(z9,k+1);>> val=abs(p);>> stem(k,val)。