ms6r2.m

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

n=0:25;
N=length(n);
x=[n>=3];          %产生输入信号序列
%x=[n==0]-[n==5];
%x=(exp(0.1*n)).*[n>=3];
%x=(0.5).^n.*[n>=0];
subplot(2,2,1);stem(n,x,'filled','k');            
title('x(n)');
a=[1,-0.5,0,0,0,0,1,-0.5];
b=[1,-1,1,0,0,0,0,0];
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)');
x01=0;y01=0;
xi=filtic(b,a,x01,y01);
y2=filter(b,a,x);         %filter求系统的响应
subplot(2,2,4);stem(y2,'filled','k');
title('y2(n)');