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

求教用4阶龙格库塔解微分方程

求教用4阶龙格库塔解微分方程求助各位高手，帮忙编程，用4阶龙格库塔解微分方程，万分感谢！！

dX/dT=cX(1-X)-k1*k2*m*X/(k1*x+k2)

wuliping1 发表于 2008-11-20 11:36

这是四阶龙4阶龙格库塔解微分方程

这是四阶龙4阶龙格库塔解微分方程的MATLAB子程序，另外在百度上面搜索下，请了解详细过程。
function x=sr(a,b,h,x1)
x=zeros(1,length(x1));
for i=1:length(x1)-1
    k1=h*(a*x(i)-b*x(i).^3+x1(i));
    k2=h*(a*(x(i)+k1/2)-b*(x(i)+k1/2).^3+x1(i));
    k3=h*(a*(x(i)+k2/2)-b*(x(i)+k2/2).^3+x1(i+1));
    k4=h*(a*(x(i)+k3)-b*(x(i)+k3).^3+x1(i+1));
    x(i+1)=x(i)+(1/6)*(k1+2*k2+2*k3+k4);
end

gwcrow 发表于 2008-11-20 18:52

十分感谢！！
但是能不能帮我写一下dX/dT=cX(1-X)-k1*k2*m*X/(k1*x+k2)这个方程的matlab程序，本人初学matlab，不胜感激

sogooda 发表于 2008-11-20 21:37

回复 板凳 gwcrow 的帖子

建议楼主自己也稍微想一想，matlab初学是可以的，但是以前是不是学过c或者fortran？楼上的提供的程序应该比较容易读懂了，不要太依靠别人，这样不好。

ch_j1985 发表于 2008-11-20 21:51

回复 板凳 gwcrow 的帖子

搜索一下版面，然后自己尝试着编写吧

dj-100 发表于 2008-11-21 11:20

建议参考清华大学出版社的《C常用算法》，有详细的求解方法！
页: [1]
查看完整版本: 求教用4阶龙格库塔解微分方程