【正文】 0 是初始函數(shù)值， h 是步長， n 為步數(shù) if nargin5 n=100。 x0=x0+h。 ％顯示三分量的關系圖 axis([20 20 50 50 0 50]) ％ 定義坐標軸長度 view(3) ％設定觀察角度 四、 實驗數(shù)據(jù)結果及分析 （ 1）各初始變量相同時 的圖像分析 各初始變量取相同的值 [10， 10， 10]， 運行上述程序 后 ， 得到如下圖像： 從圖中可以看出， 各初始變量相同時， 曲線總是被吸引回奇怪吸引子附近作來回跳躍。ok39。 end a=polyfit(x,y,3) syms x f1=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。 x1=x+h。 f3=f3(:)39。,39。 68894^2 2*68894*68894 68894^2 2*68894 2*68894]。 end Runge_kutta 方法 x0=0。 f1=f1(:)39。 h=(ba)/n。 end end L U end 4. 結果 LU_de(A) m = 4 n = 4 L = 1 0 0 0 2 1 0 0 1 2 1 0 3 3 2 1 U = 2 4 2 6 0 1 2 3 0 0 3 6 0 0 0 1 在編寫 程序過程中由于角標較多因此在運行過程中出現(xiàn)了不少角標不對的錯誤 題目；給出函數(shù) f(x)=1/(1+25x^2)，求 f（ x 在 [1,1]上取 5 個和 9個等距節(jié)點，做最小二乘擬合，得出均方誤差 方誤差。 五、 實驗結論 本實驗利用龍格 — 庫塔法對 Lorenz 方程進行了分析，從實驗中我們得出， Lorenz 方程的解對初始變量和參數(shù) r、 s、 b 具有很強的敏感性。,0,[10。 k3=feval(f,[k2*(h/2)+T(1:r1,t1)。《數(shù)值分析》課程設計實驗報告 龍格 — 庫塔法分析 Lorenz 方程 一、 問題敘述 考慮著名的 Lorenz 方程 ()dx s y xdtdy rx y xzdtdz xy bzdt? ????? ? ? ???? ???? 其中 s， r， b 為變化區(qū)域 內 有一定限制的實參數(shù)，該方程形式簡單，表面上看并無驚人之處，但由該方程揭示出的許多現(xiàn)象，促使“混沌”成為數(shù)學研究的嶄新領域，在實際應用中也產生了巨大的影響。x0+h/2])。10。 《數(shù)值分析》課程設