急!!求matlab 用二、三、四阶龙格-库塔法求解常微分方程 y'=-2xy,y...

看不到了

广州市风扬信息科技有限公司，专注于技术创新和解决方案提供，我们致力于将人工智能的无限潜力转化为实际的解决方案，为各行各业的客户提供智能化的产品和服务。使命与愿景我们的使命是通过前沿的人工智能技术，推动社会进步和商业发展。我们的...

1、首先建立自定义微分方程函数 function f = ode_fun(x,y)f=y+2*x/y^2;end 2、然后用四阶龙格库塔法求其数值解 figure(2)y0=[1]; %初值y（0）=1 h=0.1;a=0;b=5;[x,y] = runge_kutta(@(x,y)ode_fun(x,y),y0,h,a,b);disp(' x y')A=[x',y']plot(x,y...

K(1,1)=f1(x(i-1),y(i-1),z(i-1));K(2,1)=f2(x(i-1),y(i-1),z(i-1));K(1,2)=f1(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2);K(2,2)=f2(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2);K(1,3)=f1(x(i...

初始值向量，步长，时间起点，时间终点n=floor((b-a)/h);%求步数 x(1)=a;%时间起点 y(:,1)=y0;%赋初值%按龙格库塔方法进行求解 for ii=1:n x(ii+1)=x(ii)+h; k1=ufunc(x(ii),y(:,ii));

四阶龙格-库塔法 y_rk4(1)=y0;for i=1:n rk4_K1=myode(x(i),y_rk4(i));rk4_K2=myode(x(i)+h/2,y_rk4(i)+(h/2)*rk4_K1);rk4_K3=myode(x(i)+h/2,y_rk4(i)+(h/2)*rk4_K2);rk4_K4=myode(x(i)+h,y_rk4(i)+h*rk4_K3);y_rk4(i+1)=y_rk4(...

三阶龙格—库塔法的计算公式为：三阶龙格—库塔公式的Matlab程序代码：functiony=DELGKT3_kuta(f,h,a,b,y0,varvec)formatlong;N=(b-a)/h;y=zeros(N+1,1);y(1)=y0;x=a:h:b;var=findsym(f);fori=2:N+1K1=Funval(f,varvec,[x(i-1)y(i-1)]);K2=Funval(f,varvec,[x(i-...

主程序执行以下步骤：定义自定义函数1、2、3、4、5，分别对应欧拉法、二阶、四阶龙格库塔算法实现。通过自定义函数，对微分方程求解，计算t从0到2的数值解。比较不同方法求得的y(t)与理论解y(t)=exp(-t)+t的差异，观察数值解的精度。结果展示：将t从0修改至10，展示不同方法求解结果。作者：...

//例子求y'=y-2*x/y(0<x<1）；y0=1;/ int main(){ double x[6],y[6];printf("用二阶龙格-库塔方法\n");RungeKutta（1,0,1,5,x,y,2,function);for(int i=0;i<6;i++)printf("x[%d]=%f,y[%d]=%f\n",i,x[i],i,y[i]);printf("用三阶龙格-库塔方法\n");Rung...

功能：用四阶Runge-Kutta 法求解常微分方程 --- function R=Rungkuta4(f, a, b, n, ya)f:微分方程右端函数句柄 a,b:自变量取值区间的两个端点 n:区间等分的个数 ya:函数初值y(a)R=[x',y']:自变量X 和解Y 所组成的矩阵 h=(b-a)/n;x=zeros(,n+1);y=zeros(1,n+1);x=a:...

该二阶微分方程用龙格库塔法可以这样来求解。第一步，根据该二阶微分方程，自定义微分方程函数，func(t,y)第二步，根据初始条件，确定y和y'的初值，即y0=[0,0]第三步，使用ode45函数求解【t,y】的数值解，即 [t,y] = ode45(@func,[0 0.0005],y0);第四步，根据t、y、y'值，绘制t—y...