實(shí)驗(yàn)一典型環(huán)節(jié)的模擬研究

1、實(shí)驗(yàn)一 典型環(huán)節(jié)的模擬研究一實(shí)驗(yàn)?zāi)康?通過(guò)搭建控制系統(tǒng)典型環(huán)節(jié)模型，熟悉并掌握自動(dòng)控制仿真的方法。2通過(guò)對(duì)典型環(huán)節(jié)的軟件仿真研究，熟悉并掌握Matlab軟件的使用方法。3了解并掌握各典型環(huán)節(jié)的傳遞函數(shù)及其特性，觀察和分析各典型環(huán)節(jié)的響應(yīng)曲線，掌握電路模擬和軟件仿真研究方法。二實(shí)驗(yàn)內(nèi)容1搭建各種典型環(huán)節(jié)的模擬電路，觀測(cè)并記錄各種典型環(huán)節(jié)的階躍響應(yīng)曲線。2調(diào)節(jié)模擬電路參數(shù)，研究參數(shù)變化對(duì)典型環(huán)節(jié)階躍響應(yīng)的影響。3運(yùn)行Matlab軟件中的simulink仿真功能，完成各典型環(huán)節(jié)階躍特性的軟件仿真研究，并與理論計(jì)算的結(jié)果作比較。三實(shí)驗(yàn)步驟1. 典型環(huán)節(jié)的simulink仿真分析在實(shí)驗(yàn)中觀測(cè)實(shí)驗(yàn)結(jié)果時(shí)，

2、只要運(yùn)行Matlab，利用Matlab軟件中的simulink仿真功能，以及Matlab編程功能，可以完成常見(jiàn)的控制系統(tǒng)典型環(huán)節(jié)動(dòng)態(tài)響應(yīng)。1 ）比例環(huán)節(jié)： 仿真過(guò)程及其結(jié)果如下：2 ）慣性環(huán)節(jié)： 仿真過(guò)程及其仿真結(jié)果如下：3 ）積分環(huán)節(jié)： 仿真過(guò)程及其結(jié)果：4 微分環(huán)節(jié)： 仿真過(guò)程及其結(jié)果：2. 典型控制系統(tǒng)Matlab仿真分析（1）驗(yàn)證Laplace變換表：1)syms s t w a bF1=1/s;f1=ilaplace(F1);f1=simple(f1) f1 =12）syms s t w a bF1=1/s2;f1=ilaplace(F1);f1=simple(f1) f1 =t3）

3、syms s t w a bF1=1/(s+a);f1=ilaplace(F1);f1=simple(f1) f1 =exp(-a*t)4）syms s t w a bF1=1/(s2+w2);f1=ilaplace(F1);f1=simple(f1) f1 =sin(w2)(1/2)*t)/w5）syms s t w a bF1=s/(s2+w2);f1=ilaplace(F1);f1=simple(f1) f1 =cos(csgn(w)*w*t)6）syms s t w a bF1=s/(s+a)/(s+b);f1=ilaplace(F1);f1=simple(f1) f1 =1/(b-a

4、)*(b*exp(-b*t)-a*exp(-a*t)7）syms s t w a bF1=1/(s+a)2/s;f1=ilaplace(F1);f1=simple(f1)f1 =1/a2*(1-(1+a*t)*exp(-a*t)8）syms s t w a b nF1=1/(s+a)n;f1=ilaplace(F1);f1=simple(f1) f1 =ilaplace(s+a)(-n),s,t) syms s t w a b nF1=1/(s+a)5;f1=ilaplace(F1);f1=simple(f1) f1 =1/24*t4/exp(a*t)9）syms s t w a b nF1=

5、(s+b)/(s+a)/(s2+w2);f1=ilaplace(F1);f1=simple(f1) f1 =(exp(-a*t)*b*w-exp(-a*t)*a*w-b*cos(w*t)*w+a*cos(w*t)*w+w2*sin(w*t)+a*b*sin(w*t)/(w3+w*a2)10）syms s t w a b nF1=(s+a)/(s+b)2+w2);f1=ilaplace(F1);f1=simple(f1) f1 =exp(-b*t)*(cos(w*t)*w-b*sin(w*t)+a*sin(w*t)/w1）p34 例25 syms silaplace(s+2)/(s2+4*s+3

6、)ezplot(ilaplace(s+2)/(s2+4*s+3),0 5);grid;set(gca,ytick,0:.05:1.2); 2.） 例26syms sf1=ilaplace(s-3)/(s2+2*s+2)ezplot(ilaplace(s-3)/(s2+2*s+2),0 10);grid;set(gca,ytick,0:.05:1.2); 3.）.p34 例27syms sezplot(ilaplace(s+2)/(s2+2*s+1)/(s2+3*s),0 20);grid;set(gca,ytick,0:.05:1.2);4）P76 例226G1=tf(1,1 10);G2=t

7、f(1,1 1);G3=tf(1 0 1,1 4 4);numg4=1 1;deng4=1 6;G4=tf(numg4,deng4);H1=zpk(-1,-2,1);numh2=2;denh2=1;H3=1;nh2=conv(numh2,deng4);dh2=conv(denh2,numg4);H2=tf(nh2,dh2);sys1=series(G3,G4);sys2=feedback(sys1,H1,+1);sys3=series(G2,sys2);sys4=feedback(sys3,H2);sys5=series(G1,sys4);sys=feedback(sys5,H3)Zero/p

8、ole/gain: 0.083333 (s+1)2 (s+2) (s2 + 1)-(s+10.12) (s+2.44) (s+2.349) (s+1) (s2 + 1.176s + 1.023)脈沖響應(yīng)：syms sezplot(ilaplace( 0.083333*(s+1)2*(s+2)*(s2+1)/(s+10.12)/(s+2.44)/(s+2.349)/(s+1)/(s2 + 1.176*s + 1.023),0 20);grid;單位階躍響應(yīng)：syms sezplot(ilaplace( 0.083333*(s+1)2*(s+2)*(s2+1)/s/(s+10.12)/(s+2.4

9、4)/(s+2.349)/(s+1)/(s2 + 1.176*s + 1.023),0 20); grid;set(gca,ytick,0:.05:1.2);5）P78，26 syms t s;r=1*sym(Heaviside(t);c=1-exp(-2*t)+exp(-t);R=laplace(r);C=laplace(c);G=simple(C/R);GN,GD=numden(simple(C/R)g=ilaplace(G)GN =4*s+s2+2GD =(s+2)*(1+s)g =Dirac(t)+2*exp(-2*t)-exp(-t)所以系統(tǒng)的傳遞函數(shù)是： G（s）= GN/ GD=

10、(4*s+s2+2)/(s+2)*(1+s)6 P78，27（1）微分方程拉普拉斯變換：syms t s A B C D y x;z0=A*diff(sym(y(t),2)+B*diff(sym(y(t)+C*sym(y(t);z=laplace(z0);y=laplace(D*sym(Heaviside(t);S=z-yS =A*(s*(s*laplace(y(t),t,s)-y(0)-D(y)(0)+B*(s*laplace(y(t),t,s)-y(0)+C*laplace(y(t),t,s)-D/s（2）微分方程的解：syms t s Y X A B C D;A=1;B=3;C=2;D=

11、2;S=A*s2*Y+s+B*s*Y+B+C*Y-D/s;F=solve(S,Y);f0=ilaplace(F);y=simple(factor(f0)y =1+2/exp(t)2-4/exp(t)syms ty =1+2/exp(t)2-4/exp(t)ezplot(y),0,15)7 P41，213syms t s A B C D y x Y;z0=A*diff(sym(y(t),2)+B*diff(sym(y(t)+C*sym(y(t);z=laplace(z0);y=laplace(D*sym(Heaviside(t);S=z-yA=1;B=1;C=1;D=1;S=A*s2*Y+s+B*s*Y+B+C*Y-D/s;F=solve(S,Y);f0=ilaplace(F);y=simple(factor(f0)ezplot(y),0,15) S =A*(s*(s*laplace(y(t),t,s)-y(0)-D(y)(0)+B*(s*laplace(y(t),t,s)-y(0)+C*laplace(y(t),t,s)-D/sy =1-2*exp(-1/2*t)*cos(1/2*3(1/2)*t)-2/3*exp(-1/2*t)*3(1/2)*sin(1/2*3(1/2)*t