初始化

clc
clear

已知開環傳遞函數

Gs=tf([40*40],[1 0.6*40 40*40]);

求其單位負反饋的閉環傳遞函數

G=feedback(Gs,1);
figure(1)

閉環傳遞函數的分子

num=cell2mat(G.num);

閉環傳遞函數的分母

den=cell2mat(G.den);

常數K（學過控制工程的知道其含義）

finalvalue=polyval(num,0)/polyval(den,0); %y(∞)=lim sY(s)=limG(s)
[y,t]=step(G);
[yp,k]=max(y);

峰值時間tp

tp=t(k);

超調量Mp

Mp=100*(yp-finalvalue)/finalvalue;

上升時間tr

for i=1:length(y)
    if y(i)>finalvalue
        tr=(t(i)+t(i-1))/2;
        break
    end
end

調節時間ts，%2的誤差帶

len=length(t);
while (y(len)>0.98*finalvalue)&(y(len)<1.02*finalvalue)< span="">
    len=len-1;
end;
ts=t(len);

階躍響應

step(G)
disp(['PD調節前系統的上升

總結

以上是生活随笔為你收集整理的MATLAB从入门到精通-控制系统动态特性的计算以及阶跃响应波德图等的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。