micronano

金蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 524.7
帖子: 304
在線: 182.7小時
蟲號: 698144

[交流] 【求助】求解常微分方程出現(xiàn)警告，求助

求解如下微分方程時，出現(xiàn)警告
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.
無法得到要求的解區(qū)間，請幫忙分析一下。

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];

回復此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復
265求調(diào)劑已經(jīng)有4人回復
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有12人回復
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復
329求調(diào)劑已經(jīng)有5人回復
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復已經(jīng)有4人回復
材料學碩318求調(diào)劑已經(jīng)有5人回復
一志愿中國海洋大學，生物學，301分，求調(diào)劑已經(jīng)有6人回復
081700化工學碩調(diào)劑已經(jīng)有3人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

求解二階非線性常微分方程�。�！解決之后另有重謝�。�！已經(jīng)有5人回復
哪位高手可以幫我解答，常微分方程，解的延拓課后習題求解已經(jīng)有11人回復
求解二階變系數(shù)常微分方程已經(jīng)有7人回復
求助非線性常微分方程和常系數(shù)雙曲偏微分方程求解已經(jīng)有7人回復
常微分方程求解公式已經(jīng)有7人回復
求解方程，在線等……求大神指點！已經(jīng)有7人回復
Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復
Matlab數(shù)值求解非線性常微分方程已經(jīng)有5人回復
求助，matlab求解二元二階的常微分方程組已經(jīng)有9人回復
求助ANSYS求解總是警告已經(jīng)有6人回復
一個二階常微分方程的求解方法？已經(jīng)有5人回復
matlab求解常微分方程已經(jīng)有6人回復
二階常微分方程求解已經(jīng)有5人回復
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復
高手進來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復
隱式歐拉法求解一階常微分方程已經(jīng)有7人回復
1stOpt求解常微分方程邊值問題已經(jīng)有5人回復
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復
【求助】ansys求解時出現(xiàn)警告：有4個小支點調(diào)解方程求解，原因是？已經(jīng)有3人回復
【求助】出現(xiàn)以下警告一般是哪些設置出現(xiàn)了問題？已經(jīng)有4人回復
【求助】常微分方程組求解中系數(shù)與某變量值關聯(lián)的問題已經(jīng)有12人回復
【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】已經(jīng)有9人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

free--fire

新蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 543.5
帖子: 199
在線: 51.7小時
蟲號: 1051983

micronano(金幣+2): 2011-03-14 16:35:19

引用回帖:

Originally posted by micronano at 2011-03-11 15:44:43:
求解如下微分方程時，出現(xiàn)警告
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.
...

方程好像寫得有問題
另外就是最好將數(shù)據(jù)處理一下，不要一個方程里面系數(shù)差別太大

贊一下

回復此樓

3樓2011-03-14 09:07:47

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個回答

xiegangmai

版主 (職業(yè)作家)

micronano(金幣+2): 2011-03-14 16:35:06

引用回帖:

好像是微分方程奇異，求解不下去了。

贊一下

回復此樓

2樓2011-03-11 15:55:53

已閱回復此樓關注TA 給TA發(fā)消息送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

引用回帖:

我也不知道你那常微分方程組寫得對不對？再核實下。

按照你的方程組形式和初始條件幫你簡化了下，兩三行就可以解決問題了。
---------------------------------------------------------------------------------------------------
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))
---------------------------------------------------------------------------------------------------
但結(jié)果同樣遇到問題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提示，也不表示一定是錯誤

贊一下(1人)

回復此樓

4樓2011-03-14 10:01:32

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jswoo

鐵桿木蟲 (著名寫手)

仿真EPI: 2
應助: 37 (小學生)
金幣: 7511.9
帖子: 2221
在線: 1210.1小時
蟲號: 770961

★
xiegangmai(金幣+1): 鼓勵討論交流！ 2011-03-14 12:37:21
micronano(金幣+1): 2011-03-14 16:36:14

引用回帖:

這個方程有點怪
dy(1)=Y(2);
dy(2)=(1+Y(2)^2)/Y(1)+9.4597e+08*(1+Y(2)^2)^1.5;
方程組里居然連自變量都沒
再核實核實方程組可靠否？

贊一下(1人)

回復此樓

5樓2011-03-14 10:13:04

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085700資源與環(huán)境308求調(diào)劑 +6	墨墨漠 2026-03-18	6/300	2026-03-20 00:02 by 23Postgrad
[考研] 294求調(diào)劑材料與化工專碩 +14	陌の森林 2026-03-18	14/700	2026-03-19 22:38 by 學員8dgXkO
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5	Doleres 2026-03-19	5/250	2026-03-19 20:14 by 制度的
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿北京化工大學0703化學318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 0817調(diào)劑 +3	沒有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 0703化學調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper