求助!怎样用MATLAB编程实现四阶龙格库塔法求解一阶常微分方程组_百度知 ...

http://wenku.baidu.com/view/fa7ee0ebaeaad1f346933fb7.html 这里有个三阶龙格库塔的程序，可以作为参考

function theta=danbai(t,X)x=X(1);dx=X(2);ddx=-sin(x);theta=[dx;ddx];--- 命令窗口输入 >> [t,Y]=ode45(@danbai,[0 6],[pi/3 -1/2]);>> plot(t,Y(:,1),'ro-',t,Y(:,2),'bv-');>> legend('\theta-t','d\theta-t')自编龙格库塔 --- function [y,z]...

用Matlab四阶龙格库塔法求常微分方程可以按照以下方法去实现。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...

dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')y=(3*sin(5*x))/exp(2*x)一楼有错误

function R=rk4(f,a,b,ya,N)y'=f(x,y)a,b左右端点 N为迭代步长 h为步长 ya为初值 h=(b-a)/N;T=zeros(1,N+1);Y=zeros(1,N+1);T=a:h:b;Y(1)=ya;for j=1:N k1=h*feval(f,T(j),Y(j));k2=h*feval(f,T(j)+h/2,Y(j)+k1/2);k3=h*feval(f,T(j)+h/...

四阶龙格-库塔是求解常微分方程（常微分方程组）精度最高的一种数值方法。四阶龙格-库塔迭代公式为 根据四阶龙格-库塔迭代公式，可用matlab软件编写得到如下计算代码：调用方法：clc;close all;y0=1;h=0.1;a=0;b=1;[x,y] = runge_kutta(@func,y0,h,a,b);disp(' x y')[x',y']...

function [t,x]=rk4(funname,t0,t1,x0,dt)t=[];x=[];while t0<t1 t=[t;t0];x=[x;x0];k1=funname(t0,x0);t0=t0+dt/2;k2=funname(t0,x0+dt*k1/2);k3=funname(t0,x0+dt*k2/2);t0=t0+dt/2;k4=funname(t0,x0+dt*k3);x0=x0+dt/6*(k1+2*k2+2*k3+k4);t...

用matlab的四、五阶龙格库塔方法（ode45）可以求得其数值解。设定其微分方程组的初始条件X(0)=0，S(0)=0。解决方法：1、常微分方程组的自定义函数odefun（t，y），其主要内容 dy(1)=3*y(2)/(60+y(2))*y(1)-0.16*y(1);dy(2)=4*250-4*y(2)-5*y(2)/(60+y(2))*y(1)-...

以上是4阶龙格库塔法的代码:自己写函数，存为f.m function dxdt = f (t,x)dxdt(1)=exp(x(1)*sin(t))+x(2);dxdt(2)=exp(x(2)*cos(t))+x(1); % x(1)是你的f，x(2)是你的g dxdt=dxdt(:);end 自己给出t0,x0,h的值（初始时间，初值，步长）如果求t0到t1的轨迹的话...

功能：用四阶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:...