参考代码:
a1=[1 -1.845 0.850586]; b=1;
a2=[1 -1.85 0.85]; b=1;
N=100;
x=[1 zeros(1,N)];
y1=filter(b,a1,x);
y2=filter(b,a2,x);
stairs(0:N, [y1; y2].')
也可以用离散传递函数:
a1=[1 -1.845 0.850586];
a2=[1 -1.85 0.85];
G1=tf([1 0 0],a1,1);
G2=tf([1 0 0],a2,1);
impulse(G1,G2)
从响应结果看,系统2不稳定(临界稳定),这也可以通过求出系统特征根进行验证:
>> roots(a2)
ans =
1.0000
0.8500
可以看到,有一个根刚好落在单位圆上,故系统临界稳定。