二程序框图clf; n=0:299; x=cos(20*pi*n/256)+cos(200*pi*n/256); num=0.5 0.27 0.77; den=1 0 0; y=filter(num,den,x); subplot(2,1,1); plot(n,x); xlabel(时间信号n); ylabel(信号幅度 ); title( 输入信号 ); subplot(2,1,2); plot(y); xlabel(时间信号n); ylabel(信号幅度 ); title( 输出信号 ); 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 7 页clf; n=0:299; x=cos(20*pi*n/256)+cos(200*pi*n/256); num=0.45 0.5 0.45; den=1 -0.53 0.46; y=filter(num,den,x); subplot(2,1,1); plot(n,x); xlabel(时间信号n); ylabel(信号幅度 ); title( 输入信号 ); subplot(2,1,2); plot(y); xlabel(时间信号n); ylabel(信号幅度 ); title( 输出信号 ); 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 7 页学习文档仅供参考n=40; num=0.5 0.27 0.77; den=1 0 0; y=impz(num,den,n); stem(y); xlabel(时间信号n); ylabel(信号幅度 ); title( 冲激响应 ); grid; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 7 页n=40; num=0.45 0.5 0.45; den=1 -0.53 0.46; y=impz(num,den,n); stem(y); xlabel(时间信号n); ylabel(信号幅度 ); title( 冲激响应 ); grid; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 7 页学习文档仅供参考n=0:60; a=2; b=3; x1=cos(20*pi*n/256); x2=cos(200*pi*n/256); x=a*x1+b*x2; num=0.45 0.5 0.45; den=1 -0.53 0.46; y1=filter(num,den,x1); y2=filter(num,den,x2); y=filter(num,den,x); yt=a*y1+b*y2; subplot(2,1,1); stem(n,y); ylabel(信号幅度 ); subplot(2,1,2); stem(n,yt); ylabel(信号幅度 ); 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 7 页n=0:60; D=10; a=1; b=1; x=a*cos(20*pi*n/256)+b*cos(200*pi*n/256); xd=zeros(1,D) x; num=0.45 0.5 0.45; den=1 -0.53 0.46; ic=0 0; y=filter(num,den,x,ic); yd=filter(num,den,xd,ic); N=length(y); d=y-yd(1+D:N+D); subplot(3,1,1); stem(n,y); ylabel(信号幅度 ); title( 输出 yn); grid; subplot(3,1,2); stem(n,yd(1:length(y); ylabel(信号幅度 ); subplot(3,1,3); stem(n,d); xlable(时间序号n); ylable(信号幅度 ); title( 差值信号 ); grid; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 7 页学习文档仅供参考精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 7 页