24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

micronano

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 524.7
帖子: 304
在线: 182.7小时
虫号: 698144

[交流] 【求助】求解常微分方程出现警告，求助

求解如下微分方程时，出现警告
Warning: Failure at t=4.893789e-005. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.738623e-019) at time t.
无法得到要求的解区间，请帮忙分析一下。

Y0=[0.0001; 15.2571];
tspan=[0.4894e-4 0];
[t,YY]=ode23s('DYdt',tspan,Y0);
plot(YY(:,1),t)

function Yd=DYdt(t,Y)
Yd=[Y(2); (1+Y(2)^2)/Y(1)+ 9.4597e+008*(1+Y(2)^2)^1.5];

回复此楼

» 猜你喜欢

有时候真觉得大城市人没有县城人甚至个体户幸福已经有9人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有4人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有3人回复
退学或坚持读已经有28人回复
面上项目申报已经有3人回复
酰胺脱乙酰基已经有9人回复
博士延得我，科研能力直往上蹿已经有7人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求解二阶非线性常微分方程！！！解决之后另有重谢！！！已经有5人回复
哪位高手可以帮我解答，常微分方程，解的延拓课后习题求解已经有11人回复
求解二阶变系数常微分方程已经有7人回复
求助非线性常微分方程和常系数双曲偏微分方程求解已经有7人回复
常微分方程求解公式已经有7人回复
求解方程，在线等……求大神指点！已经有7人回复
Matlab数值求解二阶常微分方程已经有9人回复
Matlab数值求解非线性常微分方程已经有5人回复
求助，matlab求解二元二阶的常微分方程组已经有9人回复
求助ANSYS求解总是警告已经有6人回复
一个二阶常微分方程的求解方法？已经有5人回复
matlab求解常微分方程已经有6人回复
二阶常微分方程求解已经有5人回复
1stopt4阶非线性常微分方程，帮小弟看一下，跪求。已经有7人回复
高手进来用matlab求解常微分方程，十分感谢！已经有17人回复
隐式欧拉法求解一阶常微分方程已经有7人回复
1stOpt求解常微分方程边值问题已经有5人回复
matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助已经有3人回复
【求助】ansys求解时出现警告：有4个小支点调解方程求解，原因是？已经有3人回复
【求助】出现以下警告一般是哪些设置出现了问题？已经有4人回复
【求助】常微分方程组求解中系数与某变量值关联的问题已经有12人回复
【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】已经有9人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-11 15:44:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jswoo

铁杆木虫 (著名写手)

仿真EPI: 2
应助: 37 (小学生)
金币: 7511.9
帖子: 2221
在线: 1210.1小时
虫号: 770961

★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-03-14 12:37:04
micronano(金币+5): 2011-03-14 16:35:29

引用回帖:

Originally posted by micronano at 2011-03-11 15:44:43:
求解如下微分方程时，出现警告
Warning: Failure at t=4.893789e-005. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.738623e-019) at time t.
...

我也不知道你那常微分方程组写得对不对？再核实下。

按照你的方程组形式和初始条件帮你简化了下，两三行就可以解决问题了。
---------------------------------------------------------------------------------------------------
function main
fun = @(t,Y)[Y(2);(1+Y(2)^2)/Y(1)+9.4597e+08*(1+Y(2)^2)^1.5];
[t,Y] = ode45(fun,[0.4894e-4 0],[0.0001 15.2571]);
plot(Y(:,1),Y(:,2))
---------------------------------------------------------------------------------------------------
但结果同样遇到问题Warning: Failure at t=4.893789e-005. Unable to meet integration tolerances
without reducing the step size below the smallest value allowed
(1.738623e-019) at time t.
-------------------------------------------
分析：
1. 按提示需要减小step size
2. 检查你的方程和初值条件，如果都没问题，那么即使warning提示，也不表示一定是错误