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

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

2、只要運行matlab，利用matlab軟件中的simulink仿真功能，以及matlab編程功能，可以完成常見的控制系統(tǒng)典型環(huán)節(jié)動態(tài)響應。1 ）比例環(huán)節(jié)： 仿真過程及其結果如下：2 ）慣性環(huán)節(jié)： 仿真過程及其仿真結果如下：3 ）積分環(huán)節(jié)： 仿真過程及其結果：4 微分環(huán)節(jié)： 仿真過程及其結果：2. 典型控制系統(tǒng)matlab仿真分析（1）驗證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)脈沖響應：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;單位階躍響應：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