微分方程數(shù)值解實驗-wenkub
2023-05-01 23:19:06 本頁面
【正文】 end[t,y]=ode45(Rk4,[0,30],[12,2,9])suptitle(39。ktz=x1*y18/3*z1。 ksy=x2*(28z2)y2。%eval_r(aby)。ktz=0。ksy=0。 dy(1)=10*(y(1)+y(2))。我們要在考慮精度的基礎上，選擇合適的。關鍵詞：微分方程數(shù)值解、MATLAB 二、常微分方程數(shù)值解 基本思路 常微分方程數(shù)值解法(numerical methods forordinary differential equations)，實用上許多很有價值的常微分方程的解不能用初等函數(shù)來表示，.常微分方程初值問題的數(shù)值解法是求方程（1）的解在點列上的近似值，這里是到的步長，一般略去下標記為。其中，常微分方程求解是微分方程的重要基礎內(nèi)容。但是，對于許多的微分方程，往往很難得到甚至不存在精確的解析表達式，這時候，數(shù)值解提供了一個很好的解決思路。 （1） 經(jīng)典的方法是一個四階的方法，它的計算公式是： （2）方法的優(yōu)點是：單步法、精度高，計算過程便于改變步長，缺點是計算量較大，每前進一步需要計算四次函數(shù)值。 算法步驟 、四階龍格－庫塔（RK）方法流程圖：輸入待求微分方程、求解的自變量范圍、初值以及求解范圍內(nèi)的取點數(shù)等。 dy(2)=28*y(1)y(2)y(1)*y(3)。kty=0。 kfx=10*(y3x3)。kfz=x3*y38/3*z3。ksz=x2*y28/3*z2。 x=x3+h/12*(23*kfx16*ksx+5*ktx)。RungeKutta4階法39。legend(39。t39。ylabel(39。,14)。legend(39。t39。ylabel(39。,14)。z關于t 的變化關系圖39。,39。z39。subplot(2,2,4)plot3(y(:,1),y(:,2),y(:,3))。,1)。FontSize39。,39。y39。view(40,60)。t=1:1:3000。 x2=。 y3=。b(i)=y。z1=z2。x3=x。x=%f39。z=%f\n39。) %總標題subplot(2,2,1)。x關于t 的變化關系圖39。,39。x39。subplot(2,2,2)。y關于t 的變化關系圖39。,39。y39。subplot(2,2,3)。z關于t 的變化關系圖39。,39。z39。subplot(2,2,4)plot3(a,b,c)。,1)。FontSize39。,39。y39。view(40,60)。此時，存在混沌和一個奇怪吸引子。利用Matlab程序計算：命令窗口輸入： RungeKutta4階法： exe11Adams4階外插法： exe12結(jié)果輸出： 由于迭代次數(shù)的關系，結(jié)果數(shù)據(jù)量很大，故這里以圖片來展示結(jié)果。 用matlab編寫源程序 Matlab程序源代碼：function [u,tt,uu]=parabola(k,q,l) %q=theta l=taoh=。 %初始值endu=u39。 for i=2:9 A(i,i)=1+2*q*t/h/h。 B(1,1)=(1q)*t/h/h*(u(2)2*u(1))+u(1)。 %求解矩陣v w=u。 endtt=k*t。fprintf(39。τ=，θ= 時 ：%f\n39。subplot(2,2,1)。grid on。[u,tt,uu]=parabola(100,)。fprintf(39。,abs(uu(4)u(4)))x=::。)。)。,uu(4))[u,tt,uu]=parabola(100,)。誤差為：%f\n\n\n39。o39。τ=，θ=39。τ=，θ= 時 ：%f\n39。,u(4))fprintf(39。plot(x,u,x,uu,39。ti
教學教案相關推薦
鄂ICP備17016276號-1