收敛域在圆外,对应原始时间序列为右边序列。
上代码:
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个元素相同。