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

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

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

分段函数的一个简单方法就是分了几段就写几个函数: 定义一个主调函数、两个子函数 最后一个主调程序调用 function dx=fun(t,x) dx=zeros(2,1); f=; dc(t)=diff('f','t');%f就是你的c(t) dx(1)=x(2)*dc(t); dx(2)=x(1)+6*t^2;%0

相关常数都没给,这里随便设置一些数值: 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))

MATLAB提供了7个常微分方程求解器(solver),分别是ode45,ode23,ode113,ode15s,ode23s,ode23t,ode23tb,其中前3个适用于求解非刚性(Nonstiff)问题,后4个适用于刚性问题。所谓刚性问题,简单点说,就是系统包含多个相互作用但变化速度相差十...

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

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

t>=0即可用ode45,它的物理意义类似时间,但也不局限于时间。你的问题比较模糊。 满意请采纳。

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