bycj.net
当前位置:首页 >> mAtlAB oDE45 >>

mAtlAB oDE45

fun_u=@(t,u) u-2*t/u; [t,u]=ode45(fun_u, 0:0.01:1, 1); plot(t,u);

ode45可以用来解微分方程,基本用法如下: 一、常用格式:[t,y]=ode45(odefun,tspan,y0) 参数说明: odefun:用以表示f(t,y)的函数句柄或inline函数,t是标量,y是标量或向量。 tspan:如果是二维向量[t0,tf],表示自变量初值t0和终值tf;如果是高...

为ode45求解器设置event触发事件,当事件触发的时候就可以停止求解器继续求解 具体参考我的回答例子 http://zhidao.baidu.com/question/1046647733744158819.html?oldq=1

ode45是用4阶方法提供候选解,5阶方法控制误差,是一种自适应步长的方法。而我们平时用的4阶和5阶龙格库塔法的公式中步长是给定的。具体算法和原理你可以看Dormand, J. R. and P. J. Prince, "A family of embedded Runge-Kutta formulae," J. C...

引用一下matlab论坛里大神的回答: 总得来说:二者算法相似,只不过ode45比ode23精度要高一点,其它没什么差别。 具体ode是Matlab专门用于解微分方程的功能函数;solver有变步长(variable-step)和定步长(fixed-step)两种类型,不同类型有着...

箭头的图叫向量场或叫矢量场图,使用quiver,但一般需要和求梯度的gradient函数一起使用。 你可以参考http://zhidao.baidu.com/question/262953572739405925.html?oldq=1

这张图来自于MATLAB 技术论坛,这是一个非常有价值的网站,里面高手云集,希望你可以进去看看,里面有非常详细的关于ODE的介绍 http://www.matlabsky.com/thread-528-1-1.html

把 “[t,C]=ode45(@func1,tspan,C0,[],k)” 改成 [t,C]=ode45(@(t,C)func1(t,C,k),tspan,C0) 你原来那句是老用法,新版不支持了,改用@.

相关常数都没给,这里随便设置一些数值: m=1;c=2;k=3;y0=4;yp0=5;tf=10;dy=@(t,y)[y(2);-(c*y(2)+k*y(1))/m];[t,y]=ode45(dy,[0 tf],[y0 yp0]);plot(t,y(:,1))

ode45是最常用的求解非刚性微分方程(组)的指令,它采用变步长四、五阶Runge-Kutta法,求解精度比ode23、ode15要高。但ode45的计算量比较大,而ode23计算量小,且误差大 调用格式:[t,y]=ode45(odefun,tspan,y0) odefun 用以表示f(t,y)...

网站首页 | 网站地图
All rights reserved Powered by www.bycj.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com