数字信号处理：matlab解差分方程

1. 验证全响应

%验证全响应=零状态响应+零输入响应
%y(n)=4y(n-1)+x(n),其中x(n)=δ(n),y(-1)=2.

clc;%清屏
clear all;%清除所有变量的值
b=[1];
a=[1,-4];
ys=[2];
xs=[0];%没有初始值，就是0
xn=[1, zeros(1,4)];%输入序列，假设长度是5，则输出长度也是5
n=0:length(xn)-1;%序列横坐标

y0=filtic(b,a,ys,xs);%由初始状态得到的y0值,
yn1=filter(b,a,zeros(1, length(xn)), y0);%零输入响应,由y0递推可得。
yn2=filter(b,a,xn, 0);%零状态响应
yn=filter(b,a,xn, y0);%全响应

结果：

2. 验证系统的时不变

%验证时不变
%y(n)=4y(n-1)+x(n),其中x(n)=δ(n),y(-1)=0.

clc;%清屏
clear all;%清除所有变量的值
b=[1];
a=[1,-4];
ys=[0];
xs=[0];
xn=[1, zer