當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

micronano

金蟲 (正式寫手)

應(yīng)助: 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];

回復(fù)此樓

» 猜你喜歡

287求調(diào)劑已經(jīng)有5人回復(fù)
22408 344分求調(diào)劑一志愿華電計算機(jī)技術(shù) 已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有3人回復(fù)
317求調(diào)劑已經(jīng)有10人回復(fù)
廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機(jī)專業(yè)調(diào)劑）已經(jīng)有3人回復(fù)
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？已經(jīng)有8人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有25人回復(fù)
085600材料與化工已經(jīng)有6人回復(fù)
265求調(diào)劑已經(jīng)有8人回復(fù)

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

求解二階非線性常微分方程�。�！解決之后另有重謝！��！已經(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ù)
一個二階常微分方程的求解方法？已經(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 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
【求助】ansys求解時出現(xiàn)警告：有4個小支點(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的回帖

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:
求解如下微分方程時，出現(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.
...

這個方程有點(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的回帖

查看全部 5 個回答

xiegangmai

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

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

引用回帖:

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

贊一下

回復(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小時
蟲號: 1051983

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

引用回帖:

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

贊一下

回復(fù)此樓

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

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

jswoo

鐵桿木蟲 (著名寫手)

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

引用回帖:

我也不知道你那常微分方程組寫得對不對？再核實(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提示，也不表示一定是錯誤

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 295復(fù)試調(diào)劑 +7	簡木ChuFront 2026-03-19	7/350	2026-03-20 14:11 by 星空星月
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +4	西南交通專材355 2026-03-19	4/200	2026-03-20 11:39 by 花開富貴幸福人?/a>
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7	Doleres 2026-03-19	7/350	2026-03-20 10:39 by EBSD
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 0817調(diào)劑 +3	沒有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易