资源描述:
《北邮通信工程MATLAB双音多频.docx》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、>>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];
2、>>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>>
3、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,z
4、4)>>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=goe
5、rtzel(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=ab
6、s(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)