资源描述：
《利用Matlab画频谱图.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、5.7n=0:30;xn=sin(0.25*pi*n)+sin(0.075*pi*n);B=[1,-2*cos(0.075*pi),1];A=1;K=1/3.9447;yn=K*filter(B,A,xn);n=0:length(xn)-1;subplot(2,2,1);stem(n,xn,'.');line([0,30],[0,0]);title('aÊä³öx(n)');xlabel('n');ylabel('x(n)');n=0:length(yn)-1;subplot(2,2,2);stem(n,yn,'.');line([0,30],[
2、0,0]);title('bÊä³öy(n)');xlabel('n');ylabel('y(n)');5.10a1=[1,-0.8*exp(j*2*pi/9)];a2=[1,-0.8*exp(-j*2*pi/9)];a3=[1,-0.875*exp(j*4*pi/9)];a4=[1,-0.875*exp(-j*4*pi/9)];A=conv(conv(a1,a2),conv(a3,a4));b1=[1,-1];b2=[1,1];b3=[1,-exp(j*3*pi/4)];b4=[1,-exp(-j*3*pi/4)];B=conv(conv(b1
3、,b2),conv(b3,b4));[H,w]=freqz(B,A,'whole');subplot(2,2,1);zplane(B,A);title('a');xlabel('RE[z]');ylabel('IM[z]');subplot(2,2,2);plot(w/pi,abs(H)/max(abs(H)));title('b');xlabel('w(rad/s)');ylabel('H(ejw)');
