版塊導航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

[求助] ode45積分出現(xiàn)NaN 已有1人參與

對一個函數(shù)積分，給定初始值，積分到某一時刻，出現(xiàn)了NaN的情況，最后一項是一個常數(shù)，不用管；然后以出現(xiàn)NaN之前的一個數(shù)據(jù)重新積分，卻可以積分到最后，什么情況？

options=odeset('events',@eventRe);
x0=[0.9922,-0.1396,0,1.0157,0,0,0.1667];
[t,x,te,xe,ie]=ode45(@myGAzrfun,[0,2],x0,options)
%上面的積分積分到t=0.0324（t（17））的下一個節(jié)點0.0357時，出現(xiàn)NaN情況，見下面部分數(shù)據(jù)；

x1=x(17,

;
[t,x,te,xe,ie]=ode45(@myGAzrfun,[t(17),2],x1,options)
%然后將t=0.0324的時刻所對應(yīng)的的初值及時間重新積分，卻又能夠繼續(xù)積分下去；

0.9939 -0.1396       0 1.0139 0.0125       0 0.1667
0.9948 -0.1396       0 1.0130 0.0188       0 0.1667
0.9957 -0.1396       0 1.0121 0.0252       0 0.1667
0.9965 -0.1395       0 1.0112 0.0315       0 0.1667
   NaN    NaN    NaN    NaN    NaN    NaN 0.1667
   NaN    NaN    NaN    NaN    NaN    NaN 0.1667
   NaN    NaN    NaN    NaN    NaN    NaN 0.1667
   NaN    NaN    NaN    NaN    NaN    NaN 0.1667
   NaN    NaN    NaN    NaN    NaN    NaN 0.1667
%上面是積分遇到問題時的部分數(shù)據(jù)；
%下面是積分到最后的部分數(shù)據(jù)，而且也是我想要的數(shù)據(jù)，說明應(yīng)當是正確的，可是為啥會出現(xiàn)上面的情況？
0.4418 -0.1594       0 1.0001 0.1092       0 0.1667
0.4221 -0.1612       0 1.0001 0.1096       0 0.1667
0.4030 -0.1631       0 1.0000 0.1100       0 0.1667
te =
0.1213
xe =
0.4030 -0.1631       0 1.0000 0.1100       0 0.1667
ie =
   1

360截圖20150516124651022.jpg

回復此樓

» 猜你喜歡

310求調(diào)劑已經(jīng)有5人回復
材料與化工272求調(diào)劑已經(jīng)有17人回復
材料工程專碩求調(diào)劑已經(jīng)有3人回復
337求調(diào)劑已經(jīng)有5人回復
材料與化工328分調(diào)劑已經(jīng)有9人回復
318一志愿吉林大學生物與醫(yī)藥求調(diào)劑已經(jīng)有5人回復
一志愿南開大學0710生物學359求調(diào)劑已經(jīng)有4人回復
求調(diào)劑已經(jīng)有6人回復
085600材料與化工調(diào)劑已經(jīng)有4人回復
327求調(diào)劑已經(jīng)有4人回復

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

微分方程ode45求解，最小二乘法優(yōu)化微分方程參數(shù)，程序運行求助已經(jīng)有5人回復
matlab 運行ode45出錯，不知道什么原因已經(jīng)有4人回復
matlab編程 ODE45 已經(jīng)有4人回復
Matlab中ode45出現(xiàn)大誤差，關(guān)于用ODE45解偏微分方程的疑問，有待解決，大家共同討論！已經(jīng)有4人回復
matlab 計算后顯示NAN 已經(jīng)有5人回復
求助利用matlab的ODE45求解微分方程已經(jīng)有16人回復
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復
【求助】急求各位大俠，關(guān)于用ode45解微分方程已經(jīng)有8人回復

1樓 2015-05-16 12:56:14

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

wuhanfeng

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 113.6
紅花: 2
帖子: 64
在線: 94.9小時
蟲號: 1612956
注冊: 2012-02-12
專業(yè): 傳熱傳質(zhì)學

可能在函數(shù)在這個點產(chǎn)生奇性，或者值比較大，matlab算不出來

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復此樓

2樓2015-05-16 14:37:27

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

引用回帖:

2樓: Originally posted by wuhanfeng at 2015-05-16 14:37:27
可能在函數(shù)在這個點產(chǎn)生奇性，或者值比較大，matlab算不出來

先謝謝你的回復！
那為什么以上一個點代入繼續(xù)積分卻可以呢？而且感覺我的函數(shù)中沒有產(chǎn)生奇異點的可能，因為只有x（4）做分母，但第四列沒有出現(xiàn)零啊。

贊一下

回復此樓

3樓2015-05-16 14:45:41

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

現(xiàn)在發(fā)現(xiàn)把該語句[t,x,te,xe,ie]=ode45(@myGAzrfun,[0,2],x0,options)
改為[t,x,te,xe,ie]=ode45(@myGAzrfun,[0,1],x0,options)就可以正常積分出來了；也就是把[0,2]改為[0,1]就可以了；后來試了幾組數(shù)據(jù)，發(fā)現(xiàn)1.2還行，1.3就不行了。這是怎么回事啊？？？請各位大神指點指點！給點意見！

贊一下

回復此樓

4樓2015-05-16 17:03:21

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

wuhanfeng

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 113.6
紅花: 2
帖子: 64
在線: 94.9小時
蟲號: 1612956
注冊: 2012-02-12
專業(yè): 傳熱傳質(zhì)學

不知道你的函數(shù)具體形式，總之ode45是數(shù)值計算方法，只要在區(qū)間內(nèi)的值都能算出來，積分就能積的出來

贊一下

回復此樓

5樓2015-05-16 18:56:59

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

引用回帖:

5樓: Originally posted by wuhanfeng at 2015-05-16 18:56:59
不知道你的函數(shù)具體形式，總之ode45是數(shù)值計算方法，只要在區(qū)間內(nèi)的值都能算出來，積分就能積的出來

搞了一天，確實是在積分過程中出現(xiàn)了一個Inf，但是現(xiàn)在知道問題了。但是還是不明白為啥改變積分區(qū)間就會有不同的結(jié)果。謝謝你啦！

贊一下

回復此樓

6樓2015-05-16 22:10:22

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機械動力學

引用回帖:

6樓: Originally posted by 文刀星 at 2015-05-16 22:10:22
搞了一天，確實是在積分過程中出現(xiàn)了一個Inf，但是現(xiàn)在知道問題了。但是還是不明白為啥改變積分區(qū)間就會有不同的結(jié)果。謝謝你啦！...

我也遇到這個問題了，我的程序如果積分終點時間是0.000151就能得到151步的結(jié)果，如果積分時間是大于等于0.000152，就在128步的位置出現(xiàn)復數(shù)，但是129步的實部與128步還是連續(xù)的，請問樓主知道為什么了么？求解惑啊！

贊一下

回復此樓

7樓2015-07-18 11:46:55

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

引用回帖:

7樓: Originally posted by allengjn at 2015-07-18 11:46:55
我也遇到這個問題了，我的程序如果積分終點時間是0.000151就能得到151步的結(jié)果，如果積分時間是大于等于0.000152，就在128步的位置出現(xiàn)復數(shù)，但是129步的實部與128步還是連續(xù)的，請問樓主知道為什么了么？求解惑啊 ...

具體問題的原因我還是不太清楚。我后來把ode45的abtol和retol設(shè)到了e-10甚至更小，然后inf就消失了。因為我的被積方程中有一個函數(shù)不穩(wěn)定，所以提高精度就好點。我猜測你的應(yīng)該也是在128步時出現(xiàn)了對負數(shù)開方的情況，可以試一下將精度提高。ode45是變步長積分，它會根據(jù)積分區(qū)間取步長。但問啥會出現(xiàn)我們遇到的問題，我還是不太明白。如果你找到了，希望你也能告訴我一聲，謝謝！

贊一下

回復此樓

8樓2015-07-19 13:09:04

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

binnya

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 1
在線: 48分鐘
蟲號: 4412656
注冊: 2016-02-18

引用回帖:

8樓: Originally posted by 文刀星 at 2015-07-19 13:09:04
具體問題的原因我還是不太清楚。我后來把ode45的abtol和retol設(shè)到了e-10甚至更小，然后inf就消失了。因為我的被積方程中有一個函數(shù)不穩(wěn)定，所以提高精度就好點。我猜測你的應(yīng)該也是在128步時出現(xiàn)了對負數(shù)開方的情況 ...

樓主樓主我也是這樣的情況，你能具體說下怎么改的abtol和retol嗎？不是很明白，拜托了�。�！

贊一下

回復此樓

9樓2016-02-18 13:06:43

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導航、制導與傳感技術(shù)

引用回帖:

9樓: Originally posted by binnya at 2016-02-18 13:06:43
樓主樓主我也是這樣的情況，你能具體說下怎么改的abtol和retol嗎？不是很明白，拜托了�。�！...

用odeset函數(shù)。具體參見help文件，大概是這樣的：options=odeset('Abtol',1e-10,'Retol',1e-8);[....]=ode45(@fun,tspan,X0,options,...);

提高有時候提高精度也不一定能積出來，實在不行，個人覺得，如果精度要求不高，可以自己編寫4階龍格庫塔法求解。

贊一下

回復此樓

10樓2016-02-19 08:38:57

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主文刀星的主題更新

返回列表

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

[考研] 085600材料與化工調(diào)劑 +4	kikiki7 2026-03-30	4/200	2026-03-30 06:46 by zhyzzh
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 343求調(diào)劑 +6	愛羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無際的草原
[考研] 348求調(diào)劑 +5	小懶蟲不懶了 2026-03-28	5/250	2026-03-29 10:34 by 唐沐兒
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 086502化學工程342求調(diào)劑 +6	阿姨復古不過 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 331環(huán)境科學與工程求調(diào)劑 +3	熠然好運氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當當
[考研] 0703化學338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷 +10	材料學257求調(diào)劑 2026-03-23	11/550	2026-03-27 10:13 by YCIT- LHL
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 材料考研調(diào)劑生 +3	黃粱一夢千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪�。�！

24小時熱門版塊排行榜

[求助] ode45積分出現(xiàn)NaN 已有1人參與

» 猜你喜歡

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

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