matlab數值分析 常微分方程初邊值問題數值解?標準龍格庫塔四階四段公式?歐拉法

1、標準龍格庫塔四階四段公式

function y=rk4(fun,a,b,y0,n) h=(b-a)/n; y(1)=y0; for k=1:nx=a+(k-1)*h;k1=h*feval(fun,x,y(k));k2=h*feval(fun,x+h/2,y(k)+k1/2);k3=h*feval(fun,x+h/2,y(k)+k2/2);k4=h*feval(fun,x+h,y(k)+k3);y(k+1)=y(k)+(k1+2*k2+2*k3+k4)/6; end

函數文件

?

function u=frk4(x,y)u=y-2*x/y;

調用程序

y=rk4(@frk4,0,1,1,10)

2 歐拉法

function y=euler(f,a,b,y0,h) n=(b-a)/h; y(1)=y0; for i=1:nx(i)=a+(i-1)*h;y(i+1)=h*feval(f,x(i),y(i)); end

函數文件

?

function u=feuler(x,y)u=x^3-y/x;

調用程序

y=euler(@feuler,1,2,0.4,0.2)

?

?

總結

以上是生活随笔為你收集整理的2021-01-07 matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。