ms6r3.m

《数字信号处理实验(MATLAB版)》这

n=0:25;
N=length(n);
%x=[n==3];          %产生输入信号序列
%x=[ones(1,5),zeros(1,N-5)];
%x=cos(2*pi/3*n)+sin(3*pi/10*n);
x=0.6.^n.*[n>=3]
subplot(2,2,1);stem(n,x,'filled','k');            
title('x(n)');
b=[0.187632 ,0.241242 ,0.241242 ,0.1876321];
a=[1,-0.602012 ,0.495684 ,-0.0359244];
h=impz(b,a,n);        %求系统的单位冲激响应
subplot(2,2,2);stem(n,h,'filled','k');
title('h(n)');
y1=dlsim(b,a,x);          %dlsim求系统的响应
subplot(2,2,3);stem(y1,'filled','k');
title('y1(n)');
y01=[5,5];
xi=filtic(b,a,y01);
y2=filter(b,a,x);         %filter求系统的响应
subplot(2,2,4);stem(y2,'filled','k');
title('y2(n)');