ms6_2.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 18 行

a=[1,-0.9];                        %输入差分方程的系数
b=[1,0.9];
x01=0;y01=0;                       %输入初始条件
xi=filtic(b,a,x01,y01);            %计算初始状态
N=40;n=0:N-1;
x=(exp((-0.05+j*0.4)*n)).*[n>=2];  %建立输入信号x(n)
y=filter(b,a,x,xi);                %求系统的完全响应
subplot(2,2,1),stem(n,real(x),'filled');
title('输入信号x(n)的实部');
subplot(2,2,2),stem(n,imag(x),'filled');
title('输入信号x(n)的虚部');
subplot(2,2,3),stem(n,real(y),'filled');
title('系统响应y(n)的实部');
subplot(2,2,4),stem(n,imag(y),'filled');
title('系统响应y(n)的虚部');