4.差分方程的求解

题目：设一系统的常系数差分方程为y(n)-ay(n-1)=x(n),若系统的输入序列为x(n)=δ(n),试求系统的输出序列y(n)【代码不是很懂】

a=0.7;%差分方程系数
ys=0;%初始条件y(-1)=0
xn=zeros(1,20);%生成1行20列的零矩阵
xn(1,1)=1;%在一行一列的值为1
B=1;%差分方程x的系数，若有多个，[1 2 ...]表示
A=[1,-a];
xi=filtic(B,A,ys);%xi是等效初始条件的输入序列
%{
filtic(B,A,Y,X)
输入X，输出Y；B是x的系数，A是y的系数
filtic(B,A,Y)表示在过去X=0
%}
yn=filter(B,A,xn,xi);
%{
filter(B,A,xn,xi);一维滤波函数，用来解差分方程
B是x的系数，A是y的系数
xn是输入信号，xi是等效初始条件的输入序列
%}
N=length(yn);%求序列y(n)的序列长度
n=0:N-1;%输出有y(n)的时间序列
stem(n,yn,'.');
grid;
xlabel('n');
ylabel('h(n)');