closeall;clearall;A=[1,-];B=[,]; x1n=[11111111zeros(1,50)];x2n=ones(1,128); hn=impz(B,A,58); %subplot(2,2,1);y='h(n)';tstem(hn,y);%subplot(2,2,1);y='x2(n)';tstem(x2n,x);subplot(2,2,1);y='h(n)';stem(hn);title('(a)系统单位脉冲响应h(n)')y1n=filter(B,A,x1n); subplot(2,2,2);y='y1(n)';stem(y1n);title('(b)系统对R8(n)的响应y1(n)')y2n=filter(B,A,x2n); subplot(2,2,4);y='y2(n)';stem(y2n);title('(c)系统对u(n)的响应y2(n)')x1n=[11111111]; h1n=[ones(1,10)zeros(1,10)];h2n=[(1,10)];y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)subplot(2,2,1);y='h1(n)';stem(h1n);title('(d)系统单位脉冲响应h1(n)')subplot(2,2,2);y='y21(n)';stem(y21n);title('(e)h1(n)与R8(n)的卷积y21(n)')subplot(2,2,3);y='h2(n)';stem(h2n);title('(f)系统单位脉冲响应h2(n)')subplot(2,2,4);y='y22(n)';stem(y22n);title('(g)h2(n)与R8(n)的卷积y22(n)')un=ones(1,256); n=0:255;xsin=sin(*n)+sin(*n);A=[1,-,];B=[1/,0,-1/]; y31n=filter(B,A,un); y32n=filter(B,A,xsin);figure(3)subplot(2,1,1);y='y31(n)';stem(y31n)title('(h)谐振器对u(n)的响应y31(n)')subplot(2,1,2);y='y32(n)';stem(y32n);title('(i)谐振器对正弦信号的响应y32(n)')%clearall;closeallx1n=[ones(1,4)];%产生矩阵序列R4M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb];subplot(2,2,3);stem(x2n);x3n=[xb,xa];X1k8=fft(x1n,8);X1k16=fft(x1n,16);X2k8=fft(x2n,8);X2k16=fft(x2n,16);X3k8=fft(x3n,8);X3k16=fft(x3n,16);subplot(2,2,1);mstem(abs(X1k8));title('(1a)8点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');axis([0,2,0,*max(abs(X1k8))])%axis([0,8,0,4])subplot(2,2,2);mstem(X1k16);title('(1b)16点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');axis([0,2,0,*max(abs(X1k16))])figure(2)subplot(2,2,1);mstem(X2k8);title('(2a)8点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');axis([0,2,0,*max(abs(X2k8))])subplot(2,2,2);mstem(X2k16);title('(2b)16点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');axis([0,2,0,*max(abs(X2k16))])subplot(2,2,3);mstem(X3k8);title('(3a)8点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');axis([0,2,0,*max(abs(X3k8))])subplot(2,2,4);mstem(X3k16);title('(3b)16点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度'
低通滤波器幅频特性 来自淘豆网www.taodocs.com转载请标明出处.