自動制御原理と信号およびシステムにおいて、MATLAB を使用して単位ステップ応答を見つけ、過渡特性に対する減衰比と固有振動角周波数の影響を解析します。
Xia Deqian の『Automatic Control Theory』第 4 版のタイトル例を出発点として取り上げます。
コードは次のとおりです。
w = 1; %Wn的取值
%w = 5;%Wn的另一个取值
num = [w^2];
p = [0 0.2 0.5 0.707 1.0 1.25]; %ζ的取值
den1 = [1 2*p(1)*w w^2 ];
den2 = [1 2*p(2)*w w^2 ];
den3 = [1 2*p(3)*w w^2 ];
den4 = [1 2*p(4)*w w^2 ];
den5 = [1 2*p(5)*w w^2 ];
den6 = [1 2*p(6)*w w^2 ];
sys1=tf(num, den1);
sys2=tf(num, den2);
sys3=tf(num, den3);
sys4=tf(num, den4);
sys5=tf(num, den5);
sys6=tf(num, den6);
t = 20/w;%限定t时间范围，以获得合适的图像。
step(sys1,sys2,sys3,sys4,sys5,sys6,t);
title('ω=1');
%title('ω=5');
legend('ζ=0','ζ=0.2','ζ=0.5','ζ=0.707','ζ=1.0','ζ=1.25');
効果は次の図に示すとおりです。
w を切り替える必要がある場合は、最初の行の 1 を 5 または他の数字に変更するだけです。