micronano

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 524.7
帖子: 304
在線: 182.7小時(shí)
蟲號(hào): 698144

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

求解如下微分方程時(shí)，出現(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ū)間，請(qǐng)幫忙分析一下。

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

回復(fù)此樓

» 猜你喜歡

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

查看全部散金貼

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

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

xiegangmai

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

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

free--fire

新蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 543.5
帖子: 199
在線: 51.7小時(shí)
蟲號(hào): 1051983

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

引用回帖:

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

贊一下

回復(fù)此樓

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

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

jswoo

鐵桿木蟲 (著名寫手)

仿真EPI: 2
應(yīng)助: 37 (小學(xué)生)
金幣: 7511.9
帖子: 2221
在線: 1210.1小時(shí)
蟲號(hào): 770961

★ ★
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-03-14 12:37:04
micronano(金幣+5): 2011-03-14 16:35:29

引用回帖:

我也不知道你那常微分方程組寫得對(duì)不對(duì)？再核實(shí)下。

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

贊一下(1人)

回復(fù)此樓

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

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

jswoo

鐵桿木蟲 (著名寫手)

仿真EPI: 2
應(yīng)助: 37 (小學(xué)生)
金幣: 7511.9
帖子: 2221
在線: 1210.1小時(shí)
蟲號(hào): 770961

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

引用回帖:

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

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +18	rare12345 2026-03-18	18/900	2026-03-19 22:36 by 學(xué)員8dgXkO
[考研] 生物學(xué)調(diào)劑招人�。。� +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +5	吃吃吃才有意義 2026-03-19	5/250	2026-03-19 16:18 by 30660438
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[基金申請(qǐng)] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu