龍格—庫塔法分析lorenz方程課程設計-預覽頁
【正文】 [0， 20， 0] 從 上面可以看出 ，隨著初始 y 值的增大，奇怪吸引子中曲線在其附近來回跳躍的兩個位置中的 一個吸引力變弱，另一個吸引力變強 。當 x=0,y=0時在方程中 dx/dt=0,dy/dt=0， x， y 方向的值不發(fā)生變化，僅 z 方向的值變化，因此解 為一直線。當 b達到某一值 時，個位置喪失吸引力，另一位置則將曲線完全吸引過來變成普通吸引子。 U=zeros(n)。 end for r=2:n L(r,1)=A(r,1)/U(1,1)。 end if abs(U(i,i))eps flag=39。 for q=1:i1 m=m+L(k,q)*U(q,i)。 end x=[1 0 1]。 end sum(E) 輸出結果為 a = f1 = 4878849915647781/1298074214633706907132624082305024*x^31600/2639*x^25348847520430703/649037107316853453566312041152512*x+1514/2639 （擬合的多項式） ans = （均方誤差） 九個點的時候， matlab 編碼為： x=[1 0 1]。 end sum(E1) 輸出結果為： a = f2 = 728732707776039/2535301202056458802993406410752*x^31263030908712921/2251799813685248*x^2+4915246442354361/20282409603651670423947251286016*x+1093229300764671/2251799813685248 （最小二乘擬合多項式） ans = （均方誤差） 用復合梯形公式求積分 120 41 dxx??的值。 i2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h。 i2=0。 end end 四階龍格 庫塔法分別求解下列初值問題； 20 . 9(1 ) 39。 [ 0 , 1 ]1yy y xxxyy y xx? ? ??? ? ? ?? function[x,y]=runge_kutta(fun,x0,xt,y0,pointnum,varargin) if nargin3 y0=0。*h。 f2=f2(:)39。 f4=f4(:)39。 x=x0+[0:pointnum]39。 f=f(:)39。 fun=inline(39。,39。 pointnum=10。 62859^2 2*62859*16954 16954^2 2*62859 2*16954。 x=A\b 二．運行結果：
黨政相關相關推薦
鄂ICP備17016276號-1