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

micronano

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 524.7
帖子: 304
在線: 182.7小時(shí)
蟲號: 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ū)間，請幫忙分析一下。

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è)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請回稿延期一個(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à)值貼推薦，對您同樣有幫助:

求解二階非線性常微分方程�。�！解決之后另有重謝�。�！已經(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的回帖

jswoo

鐵桿木蟲 (著名寫手)

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

★ ★
xiegangmai(金幣+2): 謝謝應(yīng)助！ 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:
求解如下微分方程時(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.
...

我也不知道你那常微分方程組寫得對不對？再核實(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的回帖

查看全部 5 個(gè)回答

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小時(shí)
蟲號: 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í)
蟲號: 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ā)表

[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 296求調(diào)劑 +3	www_q 2026-03-18	6/300	2026-03-19 22:28 by zhq0425
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 本人考085602 化學(xué)工程專碩 +17	不知道叫什么！ 2026-03-15	19/950	2026-03-19 15:06 by 盡舜堯1
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 268求調(diào)劑 +6	簡單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠