24小时热门版块排行榜    

Znn3bq.jpeg
查看: 776  |  回复: 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的回帖

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的回帖
查看全部 5 个回答
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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工求调剂! +35 RichLi_ 2026-04-06 35/1750 2026-04-11 11:02 by zhq0425
[考研] 284求调剂 +8 archer.. 2026-04-09 9/450 2026-04-11 10:35 by zhq0425
[考研] 296求调剂 +13 汪!?! 2026-04-10 15/750 2026-04-11 10:31 by 逆水乘风
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 362求调剂 +10 我要考大 2026-04-06 14/700 2026-04-10 17:00 by luoyongfeng
[考研] 0856专硕求调剂 希望是a区院校 +21 好好休息好不好 2026-04-09 24/1200 2026-04-10 16:58 by luoyongfeng
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +27 努力奋斗112 2026-04-07 30/1500 2026-04-10 15:06 by Kilig0317
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 0854电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] (调剂)一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 8/400 2026-04-06 16:52 by momo皓
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
信息提示
请填处理意见