自控-二階系統(tǒng)Matlab仿真

自動制原二階系統(tǒng)能分析Matlab仿真大作業(yè)附題目+完整報告容

1tnt1tnt設(shè)階制統(tǒng)圖1所，中環(huán)遞數(shù)(s)

10(s(sR()

-

n(s)n

C()123要：1分用圖2和圖3示測反控和例分制種式善統(tǒng)性，果求善系的尼=則分取多

G(

ss

n

10,

10,

1210

G

stn

)

)

s2

nK)stn

K時K2

1dn21dn2

G)

Tsn(s

dns2Tsdnnn

)

T=(0.707n

2請用分別畫出第1題的系對位躍入響應(yīng)；

n

1210Matlab代clcclearwn=sqrt(10);zeta=1/(2*sqrt(10));t=0::12;Gs=tf(wn^2,[1,2*zeta*wn,wn^2]);step(Gs,t)title('圖一單位階躍響應(yīng)曲線');xlabel('t/s');ylabel('c(t)');

s

2

n2(Kst

10,n

0.707tMatlab代clcclearwn=sqrt(10);zeta=;t=0::12;Gs=tf(wn^2,[1,2*zeta*wn,wn^2]);step(Gs,t)title('圖二單位階躍響應(yīng)曲線');xlabel('t/s');ylabel('c(t)');)

s

2

dndnn

n

dMatlab代clcclearwn=sqrt(10);zeta=;t=0::12;Gs=tf([*wn^2,wn^2],[1,2*zeta*wn,wn^2]);

vvssvvvssvstep(Gs,t)title('圖三單位階躍響應(yīng)曲線');xlabel('t/s');ylabel('c(t)');3分求在位坡入，3系的態(tài)差

r()

G

ss

是Iess

lim

s0

s

11(s)H(2

lim

s0

11,KsG(s)H(s)Kv

lim

s

sG(s)H()

lim

Ks0s中K=10，

ess

110

r()

G

2.237s)s(3.47)(0.224st

是系統(tǒng)ess

lim

s0

s

11(s)H(2

lim

s0

11,KsG(s)H(s)Kv

lim

s

sG(s)H()

lim

s0

Ks中

1e

r()

Ts10(0.374sG(nss

是系統(tǒng)ess

lim

s0

s

11(s)H(2

lim

s0

11,KsG(s)H(s)Kv

lim

s

sG(s)H()

lim

Ks0s中K=10，

ess

110

4列比3個系統(tǒng)動性和態(tài)能并較析速饋控和例分制改系性的同處求峰clcclearwn=sqrt(10);zeta=1/(2*sqrt(10));t=0::12;G=tf(wn^2,[1,2*zeta*wn,wn^2]);C=dcgain(G);[y,t]=step(G);plot(t,y);[Y,k]=max(y);timetopeak=t(k)percentovershoot=100*(Y-C)/Cn=1;whiley(n)<Cn=n+1;endrisetime=t(n)i=length(t)while(y(i)>*C)&(y(i)<*C)i=i-1;end

settingtime=t(i)

000調(diào)量5試繪圖3對應(yīng)的系中數(shù)化的軌圖分變對系性的響用畫出別0，1的統(tǒng)位躍應(yīng)，較動性。

G()

Ts10Td(s(s

()

dd

2

Ts101d2

Matlab代clcclearnum=[100];den=[1110];G=tf(num,den);rlocus(G);title('Td變化的參數(shù)根軌跡');xlabel('實軸');ylabel('虛軸');)

dns2T)sdnn

10,

12

在分別取和1，可語Matlab代clcclearwn=;zeta=;t=0::12;Td=[0,,,,1];holdon;fori=1:length(Td)

Gs=tf([Td(i)*wn^2,wn^2],