版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時(shí)熱門版塊排行榜

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

返回列表

文刀星

木蟲 (小有名氣)

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

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

對一個(gè)函數(shù)積分，給定初始值，積分到某一時(shí)刻，出現(xiàn)了NaN的情況，最后一項(xiàng)是一個(gè)常數(shù)，不用管；然后以出現(xiàn)NaN之前的一個(gè)數(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））的下一個(gè)節(jié)點(diǎn)0.0357時(shí)，出現(xiàn)NaN情況，見下面部分?jǐn)?shù)據(jù)；

x1=x(17,

;
[t,x,te,xe,ie]=ode45(@myGAzrfun,[t(17),2],x1,options)
%然后將t=0.0324的時(shí)刻所對應(yīng)的的初值及時(shí)間重新積分，卻又能夠繼續(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ǐn)?shù)據(jù)；
%下面是積分到最后的部分?jǐn)?shù)據(jù)，而且也是我想要的數(shù)據(jù)，說明應(yīng)當(dāng)是正確的，可是為啥會(huì)出現(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

回復(fù)此樓

» 猜你喜歡

一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)
327求調(diào)劑已經(jīng)有3人回復(fù)

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

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

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

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

wuhanfeng

新蟲 (小有名氣)

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

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

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

贊一下

回復(fù)此樓

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

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

文刀星

木蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學(xué)生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時(shí)
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(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就不行了。這是怎么回事��？？？請各位大神指點(diǎn)指點(diǎn)！給點(diǎn)意見！

贊一下

回復(fù)此樓

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

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

wuhanfeng

新蟲 (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

文刀星

木蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

allengjn

新蟲 (正式寫手)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

文刀星

木蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(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就消失了。因?yàn)槲业谋环e方程中有一個(gè)函數(shù)不穩(wěn)定，所以提高精度就好點(diǎn)。我猜測你的應(yīng)該也是在128步時(shí)出現(xiàn)了對負(fù)數(shù)開方的情況 ...

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

贊一下

回復(fù)此樓

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

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

文刀星

木蟲 (小有名氣)

應(yīng)助: 13 (小學(xué)生)
金幣: 3651.2
散金: 500
紅花: 5
帖子: 60
在線: 100.9小時(shí)
蟲號: 3693767
注冊: 2015-02-22
性別: GG
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(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,...);

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

贊一下

回復(fù)此樓

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

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

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

返回列表

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

[考研] 材料與化工272求調(diào)劑 +16	阿斯蒂芬2004 2026-03-28	16/800	2026-03-29 22:57 by tantai88888
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9	景隨風(fēng) 2026-03-29	14/700	2026-03-29 21:22 by chixmc
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +10	辰沐5211314 2026-03-26	10/500	2026-03-29 21:01 by 唐沐兒
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +10	leaflight 2026-03-24	10/500	2026-03-29 20:58 by goldfish51
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +9	張zz111 2026-03-27	10/500	2026-03-29 14:57 by 唐沐兒
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 330分求調(diào)劑 +5	qzenlc 2026-03-29	5/250	2026-03-29 07:37 by 無際的草原
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 343求調(diào)劑 +5	愛羈絆 2026-03-28	5/250	2026-03-28 20:53 by 唐沐兒
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 化學(xué)調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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