借助MATLAB分析系统稳定性

%% 

% 借助MATLAB分析系统稳定性

 

% 1.控制系统闭环传递函数为(3*s^4+2*s^3+s^2+4*s+2)/(3*s^5+5*s^4+s^3+2*s^2+2*s+1),分析其稳定性。

num=[3,2,1,4,2];

den=[3,5,1,2,2,1];

[z,p]=tf2zp(num,den) %[z,p,k]=tf2zp(b,a)从传递函数参数b和a中找到零点矩阵z,极点向量p,

% 以及相关的增益向量 k。该函数将多项式传递函数表示转换为零极点增益形式。

pzmap(num,den) %零极点图

ii=find(real(p)>0) %k=find(X)返回一个包含数组X中每个非零元素的线性索引的向量。极点向量实部大于0

n1=length(ii)

if(n1>0),disp("System is unstable,with'int2str(n1)'unstable poles");

else disp("System is stable");

end

disp("The unstable poles are:"),disp(p(ii))

% 2.系统开环传递函数为G(s)=40/(s*(0.1*s+1)*(0.05*s+1)),绘制其乃式图及伯德图，并求出系统的相角裕量和幅值裕量。

sys=tf([40],[0.005,0.15,1,0]) %使用 tf 创建实数值或复数值传递函数模型，或将动态系统模型转换为传递函数形式。

nyquist(sys)

bode(sys)

[Gm,Pm,Wcg,Wcp]=margin(sys) %求稳定裕量，Gm为增益裕量，Pm为相位裕量，Wcg为相位-180°处的频率，Wcp为剪切频率

%% 
% 借助MATLAB分析系统稳定性
% 1.控制系统闭环传递函数为(3*s^4+2*s^3+s^2+4*s+2)/(3*s^5+5*s^4+s^3+2*s^2+2*s+1),
% 分析其稳定性。
num=[3,2,1,4,2];
den=[3,5,1,2,2,1];
[z,p]=tf2zp(num,den)   %[z,p,k]=tf2zp(b,a)从传递函数参数b和a中找到零点矩阵z,极点向量p,
% 以及相关的增益向量 k。该函数将多项式传递函数表示转换为零极点增益形式。
pzmap(num,den)               %零极点图
ii=find(real(p)>0)   %k=find(X)返回一个包含数组X中每个非零元素的线性索引的向量。极点向量实部大于0
n1=length(ii)
if(n1>0),disp("System is unstable,with'int2str(n1)'unstable poles");
else disp("System is stable");
end
disp("The unstable poles are:"),disp(p(ii))
% 2.系统开环传递函数为G(s)=40/(s*(0.1*s+1)*(0.05*s+1)),绘制其乃式图及伯德图，并求出系统的相角裕量和幅值裕量。
sys=tf([40],[0.005,0.15,1,0])  %使用 tf 创建实数值或复数值传递函数模型，或将动态系统模型转换为传递函数形式。
nyquist(sys)
bode(sys)
[Gm,Pm,Wcg,Wcp]=margin(sys)   %求稳定裕量，Gm为增益裕量，Pm为相位裕量，Wcg为相位-180°处的频率，Wcp为剪切频率