| 查看: 687 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求解常微分方程出现警告，求助
|
|||
|
求解如下微分方程时，出现警告 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]; |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解二阶非线性常微分方程！！！解决之后另有重谢！！！
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
散金
+5/655
-
坐标济南，来碰碰运气
+1/447
-
骨生物材料与侗药调控类器官再生湖南省普通高等学校重点实验室主任团队招聘了
+5/150
-
上海交通大学任垭萌课题组招聘博士后
+1/80
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
+1/68
-
湖南科技大学资安学院管青军教授2026年招收审核制博士生
+1/66
-
留学导师避雷——望传播
+1/58
-
北京-89175-事业单位-诚征女友
+1/58
-
北京—征老婆
+1/51
-
昆士兰科技大学（QUT）博士招生信息 导师：李志勇教授
+1/38
-
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1-2名
+1/27
-
坐标济南，山东农科院招 有机合成 or 药物化学 联培硕士研究生
+1/24
-
上海海洋大学水生态修复团队2026年博士招生需求
+1/8
-
北京林业大学林学院——昆虫病毒/病毒组——博士招生
+1/8
-
南京邮电大学材料科学与工程学院柔性电子研究所2026年招聘公告
+1/4
-
武汉双一流高校干细胞与肿瘤生物学团队招聘2026级申请考核制博士生
+1/4
-
武汉双一流高校干细胞与肿瘤生物学团队招聘2026级申请考核制博士生
+1/4
-
2026 博士自荐-机器人机构学方向
+1/3
-
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+1/3
-
液晶拓扑光子学博士招生（电子科技大学）
+1/2
5楼2011-03-14 10:13:04
2楼2011-03-11 15:55:53
3楼2011-03-14 09:07:47
★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-03-14 12:37:04
micronano(金币+5): 2011-03-14 16:35:29
xiegangmai(金币+2): 谢谢应助！ 2011-03-14 12:37:04
micronano(金币+5): 2011-03-14 16:35:29
|
我也不知道你那常微分方程组写得对不对？再核实下。 按照你的方程组形式和初始条件帮你简化了下，两三行就可以解决问题了。 --------------------------------------------------------------------------------------------------- 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提示，也不表示一定是错误 |
4楼2011-03-14 10:01:32