matlab 写入代码:求H(z1)=(z*z-2z-1)/(2*z*z*z-1),H(z2)=(z+1)/(z*z*z-1),H(z3)=(z*z+2)/(z*z*z+2z*z-4z+1),H(z4)=(z*z*z)/(z*z*z+0.2z*z+0.3z+0.4)这些离散系统函数分别对应的系统的零极点、零极点图和系统的稳定性。
时间: 2023-07-08 20:55:09 浏览: 105
好的,下面是 MATLAB 代码,可以求解离散系统函数的零极点、零极点图和稳定性:
1. H(z1) = (z*z-2z-1)/(2*z*z*z-1)
```matlab
% 系统函数
num = [1 -sqrt(2)+1];
den = [2 -1 0 1];
% 求解系统的零极点
z = roots(num);
p = roots(den);
% 绘制零极点图
figure;
zplane(z,p);
title('H(z1)的零极点图');
% 判断系统的稳定性
if all(abs(p) < 1)
disp('H(z1)是稳定系统');
else
disp('H(z1)是不稳定系统');
end
```
2. H(z2) = (z+1)/(z*z*z-1)
```matlab
% 系统函数
num = [1 1];
den = [1 0 -1 0];
% 求解系统的零极点
z = roots(num);
p = roots(den);
% 绘制零极点图
figure;
zplane(z,p);
title('H(z2)的零极点图');
% 判断系统的稳定性
if all(abs(p) < 1)
disp('H(z2)是稳定系统');
else
disp('H(z2)是不稳定系统');
end
```
3. H(z3) = (z*z+2)/(z*z*z+2z*z-4z+1)
```matlab
% 系统函数
num = [1 0 2];
den = [1 2 -4 1];
% 求解系统的零极点
z = roots(num);
p = roots(den);
% 绘制零极点图
figure;
zplane(z,p);
title('H(z3)的零极点图');
% 判断系统的稳定性
if all(abs(p) < 1)
disp('H(z3)是稳定系统');
else
disp('H(z3)是不稳定系统');
end
```
4. H(z4) = (z*z*z)/(z*z*z+0.2z*z+0.3z+0.4)
```matlab
% 系统函数
num = [1 0 0];
den = [1 0.3 0.2 0.4];
% 求解系统的零极点
z = roots(num);
p = roots(den);
% 绘制零极点图
figure;
zplane(z,p);
title('H(z4)的零极点图');
% 判断系统的稳定性
if all(abs(p) < 1)
disp('H(z4)是稳定系统');
else
disp('H(z4)是不稳定系统');
end
```
以上代码可以求解离散系统函数的零极点、零极点图和稳定性。
阅读全文