控制系統(tǒng)數(shù)字仿真13794

1、.現(xiàn)代工程控制理論實驗報告實驗名稱：控制系統(tǒng)數(shù)字仿真技術實驗時間： 2015/5/3精品.目錄一、實驗目的3二、實驗內容3三、實驗原理3四、實驗方案61、分別離散法；62、整體離散法；73、歐拉法94、梯形法105、龍格庫塔法11五、實驗結論12小結：14一、 實驗目的精品.1、 探究多階系統(tǒng)狀態(tài)空間方程的求解；2、 探究多種控制系統(tǒng)數(shù)字仿真方法并對之進行精度比較；二、 實驗內容1、 對上面的系統(tǒng)進行仿真，運用分別離散法進行分析；2、 對上面的系統(tǒng)進行仿真，運用整體離散法進行分析；3、 對上面的系統(tǒng)進行仿真，運用歐拉法進行分析；4、 對上面的系統(tǒng)進行仿真，運用梯形法進行分析；5、 對上面的系統(tǒng)

2、進行仿真，運用龍澤庫塔法進行分析；6、 對上面的幾種方法進行總計比較，對他們的控制精度分別進行分析比較；三、 實驗原理1、 控制系統(tǒng)狀態(tài)空間方程整體離散法的求解；控制系統(tǒng)的傳遞函數(shù)一般為有兩種控制框圖簡化形式如下：精品.ki控制器可以用框圖表示如下：慣性環(huán)節(jié)表示如下：高階系統(tǒng)的框圖如下對于上面的框圖可以簡寫傳遞函數(shù)根據各環(huán)節(jié)間的關系可以列寫出式子中出現(xiàn)的系數(shù)a、b、c和d，下面進行整體離散法求傳遞函數(shù)的推導精品.這樣，如果知道系數(shù)，就可以知道高階系統(tǒng)的傳遞函數(shù)和狀態(tài)空間方程。2、 在控制系統(tǒng)的每一個環(huán)節(jié)都加一個采樣開關，構成分別離散法求解系統(tǒng)的狀態(tài)空間方程；采樣開關其實是一個零階保持器精品.比

3、例環(huán)節(jié)：積分環(huán)節(jié)：慣性環(huán)節(jié)：四、 實驗方案1、 分別離散法；系統(tǒng)框圖根據上面提到的分別離散法得到仿真的公式已知系數(shù)：k1=0.93;k2=2.086;t1=73.3;t2=96.1;n1=2;n2=4;kp1=0.32;ki1=0.0018;kp2=2;ki2=0.00008;精品.慣性環(huán)節(jié)的系數(shù)：fai1=exp(-dt/t1);faim1=1-fai1;fai2=exp(-dt/t2);faim2=1-fai2;pid 控制環(huán)節(jié)：up1=e*kp1;x(1)=x(1)+ki1*dt*e；up2=e1*kp2;x(2)=x(2)+ki2*dt*e1;慣性環(huán)節(jié)：x(3)=fai1*x(3)+k

4、1*faim1*u1;x(4)=fai1*x(4)+faim1*x(3);x(5)=fai2*x(5)+k2*faim2*x(4);x(6)=fai2*x(6)+faim2*x(5);x(7)=fai2*x(7)+faim2*x(6);x(8)=fai2*x(8)+faim2*x(7);2、 整體離散法；將系統(tǒng)框圖拆開系統(tǒng)的狀態(tài)空間方程為：精品.可以得到此時狀態(tài)方程的系數(shù)精品.由上面的推導可知求出就可以得到系統(tǒng)的狀態(tài)空間方程在matlab中仿真時為for i=1:n1*n2 faim=faim+(dti)*(a(i-1)/factorial(i);endfai=faim*a+eye(n1*n2

5、);faim=faim*b;for j=1:lp x=fai*x+faim*r; y=c*x+d*r; y1=y1 y; t=t j*dt;end3、 歐拉法由上面已經求出系統(tǒng)的狀態(tài)空間方程，所以這里直接引用，歐拉法的求解過程如下：在matlab中的仿真程序如下：for i=1:lp精品. xk=a*x+b*r; x=x+xk*dt; y=c*x+d*r; y1=y1 y; t=t dt*i;end4、 梯形法類似于歐拉法，梯形法的推導如下在matlab中仿真的程序如下：for i=1:lp xk=a*x+b*r; xk1=x+dt*xk; xk2=a*xk1+b*r; e=(xk+xk2)/

6、2; x=x+dt*e; y=c*x+d*r; y1=y1 y; t=t dt*i;end5、 龍格庫塔法推導如下：精品.在matlab中的仿真程序如下：for i=1:lp e1=a*x+b*r; xk1=x+dt*e1/2; e2=a*xk1+b*r; xk2=x+dt*e2/2; e3=a*xk2+b*r; xk3=x+dt*e3/2; e4=a*xk3+b*r; e=(e1+e2+e3+e4)/6; x=x+dt*e; y=c*x+d*r; y1=y1 y; t=t dt*i;end五、 實驗結論5種方法仿真圖形精品.放大后的圖像此時，可以看出，分別離散已經開始遠離其他的線繼續(xù)放大精品.此時分別離散已經明顯遠離其他，并且歐拉法也開始遠離其他的線精品.最終可以看出，龍格庫塔法與整體離散法得到的仿真曲線最接近。小結：利用不同的方法對多階系統(tǒng)的狀態(tài)方程進行求解，分別離散法，因為零階保持器的緣故，所以誤差比較大；歐拉法通過簡單的取切線的端點作為下一步的起點，提升了精確性，但是本身也存在缺點，當步數(shù)增加時，誤差在逐漸累積；詳細實例見附件；梯形法是歐拉法的升級版，首先可以由歐拉法求得下一時刻的值，再代入校正得到一個更精準的值，這樣，可以較歐拉法得到更精準的值；龍格庫塔法是至