龍格—庫(kù)塔法分析lorenz方程課程設(shè)計(jì)(存儲(chǔ)版)
【正文】 1。 for i=0:n1 x=a+h*i。 end y(1,:)=y0(:)39。 f3=h*feval(fun,x(k)+h/2,y(k,:)+f2/2,varargin{:})。*h。9*y/(1+2*x)39。 [x1,yr]=runge_kutta(fun,x0,xt,y0,pointnum) 一． 運(yùn)行結(jié)果： 1） Runge_kutta 方法： x1 = 0 yr = 小行星軌道問題：一天文學(xué)家要確定一顆小行星繞太陽(yáng)運(yùn)行的軌道，他在軌道平面內(nèi)建立 以太陽(yáng)為原點(diǎn)的直角坐標(biāo)系，在五個(gè)不同的對(duì)小行星作了五次觀察，測(cè)得軌道上五個(gè)點(diǎn)的坐標(biāo)數(shù)據(jù)（單位：萬公里）如下表所示： P1 P2 P3 P4 P5 X 坐標(biāo) 53605 58460 62859 66662 68894 Y 坐標(biāo) 6026 11179 16954 23492 68894 由開普勒第一定律知，小行星軌道為一橢圓，橢圓的一般方程可表示為： a1x2+2a2xy+a3y2+2a4x+2a5y+1=0 現(xiàn)需要建立橢圓的方程以供研究。 b=[1,1,1,1,1]。 y0=1。 xt=1。 end h=(xtx0)/pointnum。 f2=h*feval(fun,x(k)+h/2,y(k,:)+f1/2,varargin{:})。 , ( 0 ) 2 。 i1=i2。 for 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。 五個(gè)節(jié)點(diǎn)時(shí)， matlab編碼為： 首先建立 M文件，并保存 function y=f(x) y=1/(1+25*x^2)。 end U(i,j)=A(i,j)z。 《數(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)。 z 分量： [0， 0， 20] 從上圖 可以看出解的曲線為一直線，這可以從方程的角度來解釋。10。 b=。x0+h/2])。 ％ 返回初始輸入矩陣的行列數(shù)，并將值賦給 s(1) r=r+s。《數(shù)值分析》課程設(shè)計(jì)實(shí)驗(yàn)報(bào)告 龍格 — 庫(kù)塔法分析 Lorenz 方程 一、 問題敘述 考慮著名的 Lorenz 方程 ()dx s y xdtdy rx y xzdtdz xy bzdt? ????? ? ? ???? ????
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1