matlab 數值計算課 二階微分方程-龍格庫塔方法 & ODE45

詳見mathworkshtml 龍格庫塔方法 寫成矩陣（狀態方程）的形式更簡潔一點（其實這兩種方法結果是同樣的，若是C是[1,0,0]的話，就很明顯了） 例如：求系統在0-5s內的單位階躍相應，已知傳遞函數： 狀態方程： web clear %狀態方程 A=[0,1 -100,-10]; B=[0; 1]; C=[100,0]; %初值 x=[0;0];y=0;t0=0;tf=5;

>>阅读原文<<