gwcrow 发表于 2008-10-10 10:48

求助：关于用matlab求解龙格库塔的方法

请问怎么用四阶龙格库塔求解微分方程，用matlab的方法求dV/dt=KV。
谢谢各位高手啦

蜜雪儿 发表于 2008-10-10 12:10

function =rk4(dfun,xspan,y0,h)
x=xspan(1):h:xspan(2);
y(1)=y0;
for k=1:length(x)-1
    k1=feval(dfun,x(k),y(k));
    k2=feval(dfun,x(k)+h/2,y(k)+h/2*k1);
    k3=feval(dfun,x(k)+h/2,y(k)+h/2*k2);
    k4=feval(dfun,x(k+1),y(k)+h*k3);
    y(k+1)=y(k)+h*(k1+2*k2+2*k3+k4)/6;
end
x=x';y=y';

gwcrow 发表于 2008-10-10 16:52

万分感谢！！！！:victory: :victory:

happy 发表于 2008-10-11 06:06

http://forum.vibunion.com/forum/viewthread.php?tid=17615

相较二楼的要完善很多

gwcrow 发表于 2008-10-11 11:20

确实很完善，万分感谢！！！！！！！！
页: [1]
查看完整版本: 求助：关于用matlab求解龙格库塔的方法