自控-二阶系统Matlab仿真

自动制原二阶系统能分析Matlab仿真大作业附题目+完整报告容

1tnt1tnt设阶制统图1所，中环递数(s)

10(s(sR()

-

n(s)n

C()123要：1分用图2和图3示测反控和例分制种式善统性，果求善系的尼=则分取多

G(

ss

n

10,

10,

1210

G

stn

)

)

s2

nK)stn

K时K2

1dn21dn2

G)

Tsn(s

dns2Tsdnnn

)

T=(0.707n

2请用分别画出第1题的系对位跃入响应；

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('图一单位阶跃响应曲线');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('图二单位阶跃响应曲线');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('图三单位阶跃响应曲线');xlabel('t/s');ylabel('c(t)');3分求在位坡入，3系的态差

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

是系统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

是系统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个系统动性和态能并较析速馈控和例分制改系性的同处求峰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调量5试绘图3对应的系中数化的轨图分变对系性的响用画出别0，1的统位跃应，较动性。

G()

Ts10Td(s(s

()

dd

2

Ts101d2

Matlab代clcclearnum=[100];den=[1110];G=tf(num,den);rlocus(G);title('Td变化的参数根轨迹');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],