龍格—庫(kù)塔法分析lorenz方程課程設(shè)計(jì)-在線瀏覽
2024-11-01 19:54本頁(yè)面
【正文】 [0， 5， 0] [0， 15， 0] [0， 20， 0] 從 上面可以看出 ，隨著初始 y 值的增大，奇怪吸引子中曲線在其附近來(lái)回跳躍的兩個(gè)位置中的 一個(gè)吸引力變?nèi)酰硪粋€(gè)吸引力變強(qiáng) 。 這說(shuō) 明在空間存在一些區(qū)域，當(dāng)初始位置位于這些區(qū)域外時(shí)解將出現(xiàn)奇怪吸引子的性質(zhì)，而在這些區(qū)域以內(nèi)解將呈現(xiàn)普通吸引子的性質(zhì)。當(dāng) x=0,y=0時(shí)在方程中 dx/dt=0,dy/dt=0， x， y 方向的值不發(fā)生變化，僅 z 方向的值變化，因此解 為一直線。 當(dāng)只調(diào)整 b 且將初始變量取為 [0， eps， 0]。當(dāng) b達(dá)到某一值 時(shí)，個(gè)位置喪失吸引力，另一位置則將曲線完全吸引過(guò)來(lái)變成普通吸引子。 五、 實(shí)驗(yàn)結(jié)論 本實(shí)驗(yàn)利用龍格 — 庫(kù)塔法對(duì) Lorenz 方程進(jìn)行了分析，從實(shí)驗(yàn)中我們得出， Lorenz 方程的解對(duì)初始變量和參數(shù) r、 s、 b 具有很強(qiáng)的敏感性。 U=zeros(n)。ok39。 end for r=2:n L(r,1)=A(r,1)/U(1,1)。 for r=1:i1 z=z+L(i,r)*U(r,j)。 end if abs(U(i,i))eps flag=39。 return。 for q=1:i1 m=m+L(k,q)*U(q,i)。 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)，做最小二乘擬合，得出均方誤差 方誤差。 end x=[1 0 1]。 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]。 end sum(E) 輸出結(jié)果為 a = f1 = 4878849915647781/1298074214633706907132624082305024*x^31600/2639*x^25348847520430703/649037107316853453566312041152512*x+1514/2639 （擬合的多項(xiàng)式）
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1