《DSP using MATLAB》示例Example4.9

时间:2023-01-19 17:03:12

《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)

  运行结果:

《DSP using MATLAB》示例Example4.9

《DSP using MATLAB》示例Example4.9

查表找逆z变换,

《DSP using MATLAB》示例Example4.9

《DSP using MATLAB》示例Example4.9

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