| 查看: 973 | 回复: 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]; |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
【长期招聘】复旦大学 AI 制药实验室 博士后招聘
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有5人回复
-
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解二阶非线性常微分方程！！！解决之后另有重谢！！！
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
香港岭南大学郭瑛课题组博士招生 (2027年秋季入学）
+4/120
-
东北大学应用经济学博士招生-2027级
+1/97
-
香港理工大学化学系李光超课题组2027Fall博士招聘
+1/85
-
博后出站想去海外工作，植物分子生物学领域
+1/46
-
海南大学徐月山老师招收2027年第一批博士研究生1名，欢迎联系咨询。
+1/41
-
涂层低温高湿（雨淞覆冰）测试
+1/37
-
慢蹲一个与子偕老的爱人
+1/32
-
东北大学管理科学与工程专业博士招生（2027年度）
+1/25
-
浙江大学-郭强兵课题组诚招博士后/博士生/科研助理
+1/17
-
您感兴趣的遗传病是哪些？您认为产前诊断在阻断遗传病传递作用有多大？
+1/16
-
香港科技大学（广州）收水文学方向科研助理
+1/9
-
中国石油大学（北京）石油工程专业-开发方向课题组 2027博士招生
+1/7
-
西安交通大学前沿院/机械学院招收2027级硕博研究生（微纳制造与光学方向）
+1/6
-
蔚蓝海岸大学招2027届CSC有机方法学博士生
+1/5
-
海南大学海洋技术与装备学院导师招收材料专业博士
+1/4
-
【全奖博士招生】澳洲莫纳什大学（DECRA资助 | 顶尖团队联培 | 尽快入学优先）
+1/3
-
Circular Energy & Resources期刊征稿
+1/3
-
PLA检测蛋白质翻译后修饰：单分子级灵敏度与原位空间解析
+1/2
-
上海应用物理研究所-2027年博士-申请考核 -还有1个名额（散尽金币再发一次）
+1/1
-
请问手中有可以提取3微米，6微米的球形颗粒，提取设备，这个技术有市场吗？
+1/1
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
5楼2011-03-14 10:13:04