實驗一經(jīng)典的連續(xù)系統(tǒng)仿真建模方法的值積分法仿真程序。1.編寫四階Runge_Kutta公式的計算程序，對非線性模型(3)式進行仿真。(1)將閥位u增大10％和減小10％，觀察響應(yīng)曲線的形狀;0100200300400500600700800900100001002003004005006007008009001000開大或關(guān)小閥位之后，穩(wěn)態(tài)值會相應(yīng)的從原液位上升或下降，這是符合實際的。(2)研究仿真步長對穩(wěn)定性的影響，仿真步長取多大時RK4算法變得不穩(wěn)定？由(1)可知，當(dāng)步長為40時，仿真結(jié)果是穩(wěn)定的0時的圖像0123456700200400600800100012001400160018002000h(1,1)的數(shù)值穩(wěn)定，但是并不是實際求得的穩(wěn)態(tài)值。h(1,2)的值顯然發(fā)散。進一步取小步長，取hstep=42時，圖像出現(xiàn)偏差，但是穩(wěn)態(tài)值不變020040060080010001200Hstep=65時，圖像偏差明顯32.5210.5020040060080010001200140016001800501015020040060080010001200140016001800(3)利用MATLAB中的ode45()函數(shù)進行求解，比較與(1)中的仿真結(jié)果有何區(qū)別。Ode45調(diào)用的函數(shù)：A=2;將ode45與編寫的龍格庫塔算法畫到同一個坐標(biāo)系中(點表示ode45)：01002003004005006007008009001000可以發(fā)現(xiàn)，ode45更快達到穩(wěn)定值，實際上，在縮小了編寫算法中的步長后，兩種算法的曲010020030040050060070080090010002.編寫四階Runge_Kutta公式的計算程序，對線性狀態(tài)方程(18)式進行(1)將閥位增大10％和減小10％，觀察響應(yīng)曲線的形狀;050100150200250050100150200250(2)研究仿真步長對穩(wěn)定性的影響，仿真步長取多大時RK4算法變得不穩(wěn)定？hstep就不穩(wěn)定了(3)閥位增大10％和減小10％，利用MATLAB中的ode45()函數(shù)進行求解階躍響050100150200250增大10%：050100150200250%DX被ode45調(diào)用的微分方程函數(shù)A=2;Adh=AMTRIX*x+BMATRIX*u;實驗完成后，要寫出實驗報告，內(nèi)容包括：2．實驗所用的仿真程序清單，以及程序結(jié)構(gòu)的簡單說明；3．實驗結(jié)果曲線及分析，穩(wěn)態(tài)值是多少；一般來說，仿真步長越長，系統(tǒng)越傾向于不穩(wěn)定，仿真步長越短