second order system analysis 自动控制原理二阶系统的matlab仿真分析

本文主要是介绍second order system analysis 自动控制原理二阶系统的matlab仿真分析，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

二阶系统的matlab仿真分析

二阶系统的matlab仿真分析如上图。

根据二阶函数对阶跃函数的响应函数，我们对参数epsilon进行分析讨论

由于临界阻尼和无阻尼的情况在现实生活中比较难出现，二阶方程的根几乎不可能恰好，实部为0，或者两个实部相同且虚部为0. 于是，并为对以上两种较特殊的情况进行讨论。

选择欠阻尼和过阻尼两种情况进行分析讨论。

可以看出，当epsilon比较小的时候，响应时间短，且伴随有明显的超调。

随着epsilon的增大，超调明显降低，epsilon在0.7(恰巧工程上的最佳阻尼系数是0.7！)之后就没有超调了。随着epsilon的增大响应时间变得越来越长。

Wn = (1/T);

二阶方程的根

实部：X = -epsilon_0*Wn

虚部：Y = j*(Wn ).*sqrt(1-(epsilon_0).^2)

下图是上面十条曲线对应的根分布

可以看出图中所有的点均位于Y轴左侧，说明最终响应都将收敛

虚部为0的点很好的对应了过阻尼状态(同样颜色的是一对实根)，图中关于X轴对称的点，系统处于欠阻尼状态

下图是上升时间和阻尼系数epsilon之间的关系(欠阻尼状态下)

可以看出，随着阻尼系数的增大，上升时间变长！

matlab和本文相关代码：

%%*************************************************************
% code writer: EOF
% code date:2014.03.18
% e-mail: jasonleaster@gmail.com
% code purpose : 
%           I just want to share with someone who is interesting
% in adaptive control. This code is to help people to understand
% second order system.
%%**************************************************************
clear all
clc
syms s f t m;
K01 = 1;
K02 = 1;
K0 = (K01*K02)./(1+K01*K02);
hold on;
figure(1);
T0 = 1;
for epsilon_0 = 0.1:0.2:2T = T0./(1+K0);
%     epsilon_0 = 0.5*(1/(K01*K02*T0));epsilon = epsilon_0./(1+K0);K = K0/(1+K0);f = (K./((T.^2).*(s.^2)+2*epsilon_0.*T.*s+K0)).*(1./s);m = ilaplace(f);ezplot(m,[0,120]);axis([0 60 0 1.2]);
end
legend('0.1','0.3','0.5','0.7','0.9','1.1','1.3','1.5','1.7','1.9');
hold off;
figure(2);
hold on;
for epsilon_0 = 0.1:0.2:2if epsilon_0 <1plot(-epsilon_0.*(1./T),(1./T).*sqrt(1-(epsilon_0).^2),'*');plot(-epsilon_0.*(1./T),-(1./T).*sqrt(1-(epsilon_0).^2),'*');elseif abs(epsilon_0-1.1) < 0.1  plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','r');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','r');elseif abs(epsilon_0-1.3) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','g');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','g');elseif abs(epsilon_0-1.5) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','b');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','b');elseif abs(epsilon_0-1.7) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','y');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','y');elseif abs(epsilon_0-1.9) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','k');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','k');end
end
axis([-6 1 -3 3]);
hold off;
figure(3);
hold on;
Wn = (1./T);
for epsilon_0 = 0.1:0.2:2 Wd = Wn.*(1-(epsilon_0).^2);if epsilon_0 < 1belta = acos(epsilon_0);tr = (pi-belta)./Wd;plot(epsilon_0*10,tr,'*');end
end
legend('epsilon = 0.1','epsilon = 0.3','epsilon = 0.5','epsilon = 0.7','epsilon = 0.9');

这篇关于second order system analysis 自动控制原理二阶系统的matlab仿真分析的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！