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