24小时热门版块排行榜    

Znn3bq.jpeg
查看: 772  |  回复: 4

micronano

金虫 (正式写手)


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

求解如下微分方程时,出现警告
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];
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
micronano(金币+2): 2011-03-14 16:35:06
引用回帖:
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.
...

好像是微分方程奇异,求解不下去了。
2楼2011-03-11 15:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

free--fire

新虫 (小有名气)


micronano(金币+2): 2011-03-14 16:35:19
引用回帖:
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.
...

方程好像写得有问题
另外就是最好将数据处理一下,不要一个方程里面系数差别太大
3楼2011-03-14 09:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jswoo

铁杆木虫 (著名写手)


★ ★
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提示,也不表示一定是错误
4楼2011-03-14 10:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jswoo

铁杆木虫 (著名写手)



xiegangmai(金币+1): 鼓励讨论交流! 2011-03-14 12:37:21
micronano(金币+1): 2011-03-14 16:36:14
引用回帖:
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.
...

这个方程有点怪
dy(1)=Y(2);      
dy(2)=(1+Y(2)^2)/Y(1)+9.4597e+08*(1+Y(2)^2)^1.5;  
方程组里居然连自变量都没
再核实核实方程组可靠否?
5楼2011-03-14 10:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 micronano 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0859,337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考博] 材料方向考博,求推荐 +3 言语aaa 2026-04-05 4/200 2026-04-08 22:22 by nxgogo
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 男生,一志愿沪9生物学071000,初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
信息提示
请填处理意见