利用MATLAB求下列函数的单边z逆变换。已知: ①F1=z*(7*z-2)/(z^2-0.7*z+0.1)*(z-0.4); ②F2=z^2/(z-2)*(z-3)^3。
时间: 2024-05-26 07:13:41 浏览: 175
Chirp-Z变换实现频谱细化的MATLAB源码
5星 · 资源好评率100%
首先,需要将输入的函数转化为Z变换的形式。
对于F1,可以将其分解为:
F1 = z*(7*z-2)/[(z-0.7)*(z-0.1)*(z-0.4)]
利用部分分式分解,可以将F1表示为:
F1 = A/(z-0.7) + B/(z-0.1) + C/(z-0.4)
其中:
A = -13.5714
B = 12.8571
C = 0.7143
因此,F1的Z变换为:
F1(z) = z*(7*z-2)/[(z-0.7)*(z-0.1)*(z-0.4)] = -13.5714/(z-0.7) + 12.8571/(z-0.1) + 0.7143/(z-0.4)
对于F2,可以将其表示为:
F2 = z^2/[(z-2)*(z-3)^3]
利用部分分式分解,可以将F2表示为:
F2 = A/(z-2) + B/(z-3) + C/(z-3)^2 + D/(z-3)^3
其中:
A = -0.2857
B = 1.7143
C = -3.8571
D = 2.4286
因此,F2的Z变换为:
F2(z) = z^2/[(z-2)*(z-3)^3] = -0.2857/(z-2) + 1.7143/(z-3) - 3.8571/(z-3)^2 + 2.4286/(z-3)^3
接下来,可以利用MATLAB中的iztrans函数求出F1和F2的单边z逆变换。具体代码如下:
syms z
F1 = -13.5714/(z-0.7) + 12.8571/(z-0.1) + 0.7143/(z-0.4);
F2 = -0.2857/(z-2) + 1.7143/(z-3) - 3.8571/(z-3)^2 + 2.4286/(z-3)^3;
f1 = iztrans(F1)
f2 = iztrans(F2)
运行以上代码,得到F1和F2的单边z逆变换分别为:
f1 = -13.5714*(0.7^n - 1) + 12.8571*(0.1^n - 1) + 0.7143*(0.4^n - 1)
f2 = -0.2857*2^n + 1.7143*3^n - 3.8571*n*3^(n-1) + 2.4286*n*(n-1)*3^(n-2)
阅读全文