《DSP using MATLAB》示例Example4.9

收敛域在圆外，对应原始时间序列为右边序列。

上代码：

b = 1; a = poly([0.9, 0.9, -0.9]);               %  compute the polynomials coefficients given its roots

[R,p,C] = residuez(b,a)

[delta, n] = impseq(0, 0, 7); x_check = filter(b, a, delta);       % check sequence

x_answer = (0.75)*(0.9).^n + (0.5)*n.*(0.9).^n + (0.25)*(-0.9).^n   % answer swquence

%[b,a] = residuez(R,p,c)

　　运行结果：

查表找逆z变换，

 

    由脉冲序列激励，filter函数得到的x_check和查表得到的x_answer前8个元素相同。