龍格—庫塔法分析lorenz方程課程設(shè)計(jì)-文庫吧資料
2024-09-06 19:54本頁面
【正文】 r i=1:5 E1(i)=(f(x(i))(a(1)*x(i)^3+a(2)*x(i)^2+a(3)*x(i)+a(4)))^2。 for i=1:9 y(i)=f(x(i))。 for i=1:5 E(i)=(f(x(i))(a(1)*x(i)^3+a(2)*x(i)^2+a(3)*x(i)+a(4)))^2。 for i=1:5 y(i)=f(x(i))。 五個(gè)節(jié)點(diǎn)時(shí)， matlab編碼為： 首先建立 M文件，并保存 function y=f(x) y=1/(1+25*x^2)。 end L(k,i)=(A(k,i)m)/U(i,i)。 end for k=i+1:n m=0。failure39。 end U(i,j)=A(i,j)z。 end for i=2:n for j=i:n z=0。 for i=1:n U(1,i)=A(1,i)。 flag=39。 《數(shù)值分析》課程設(shè)計(jì)實(shí)驗(yàn)報(bào)告 數(shù)值分析試驗(yàn)報(bào)告 矩陣的 LU 分解 1. 題目：求 4 階矩陣??????????????401815618962156946242的 LU 分解 2. 方法：杜里特爾分解法 3. 程序： function f=LU_de(A) [m,n]=size(A) L=eye(n)。 改變s 和 r 的值也有類似的現(xiàn)象。具體情況如下： s=， r=， b=， r=， b= s=， r=， b=， r=， b= 增大 b 值時(shí)， Lorenz 曲線在其附近來回跳躍的兩個(gè)位置會(huì)一個(gè)加強(qiáng)，一個(gè)減弱。 （ 3） 調(diào)整參數(shù) r、 s、 b 對圖像的影響 為便于分析 ，我們只調(diào)整 r、 s、 b 三個(gè)參數(shù)中的任意一個(gè)。 z 分量： [0， 0， 20] 從上圖 可以看出解的曲線為一直線，這可以從方程的角度來解釋。初始 y 繼續(xù)增大到某一特定值，情況又會(huì)變回來。初始變量值取為 [10， 10， 10] ， [20， 20， 20]時(shí)，依然如此。 ％ 調(diào)用前面的算法程序 plot3(T(1,:),T(2,:),T(3,:))。10。fun39。 dy=dy39。 ％ Lorenz方程表達(dá)式 dy(2)=(r*x(1)x(3)*x(1)x(2))。 b=。 end （ 2） 主程序 function dy=fun(x) ％ 定義函數(shù) s=。 T(:,t)=[T(1:r1,t1)+(k1+k2*2+k3*2+k4)*(h/6)。x0+h])。x0+h/2])。x0+h/2])。 for t=2:n+1 ％ 以下是 具體的求解過程 k1=feval(f,T(1:r1,t1))。 T(:,1)=[y0。
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1