版塊導(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小時熱門版塊排行榜

返回列表

tianaylk

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線: 9小時
蟲號: 865837
注冊: 2009-10-08
專業(yè): 電力系統(tǒng)

[求助] matlab 積分計算算不出結(jié)果，求幫助�。�！

q =1.0e+008 *[2.3916; -3.6135; 1.8065; -0.6109];
syms y ;
F2=int(exp(q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3)),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)
可以計算出F4 = 149272173.68701522388239715123227，但是當(dāng)q =1.0e+009 *[2.3916; -3.6135; 1.8065; -0.6109] 時，運(yùn)算不出結(jié)果，其實(shí)結(jié)果大概為149272173*10，懷疑是不是因?yàn)镕2數(shù)值太大，matlab處理不了？q的取值為1.0e+008數(shù)量級時F2=2*10^64828081，q的取值為1.0e+009數(shù)量級時F2就計算不出來了，提示 numeric::int(1/exp(152725000*y^4 - (5051340117333333*y^3)/8388608 + 890400000*y^2 - (4908174540800001*y)/8388608 - 1348575000), y = 9/10..11/10) 但其實(shí)取對數(shù)以后得到的F3并不大，請各位大俠幫忙解決一下~~~非常感激~~~

回復(fù)此樓

» 猜你喜歡

一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武理材料305分求調(diào)劑已經(jīng)有6人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有7人回復(fù)

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

求matlab 解方程的問題已經(jīng)有5人回復(fù)
關(guān)于Matlab算符號積分已經(jīng)有6人回復(fù)
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
求助使用MATLAB計算綠色葉片的覆蓋度的問題已經(jīng)有24人回復(fù)
matlab 計算過程中工作空間變量保存及讀取求解決內(nèi)存不足的方法已經(jīng)有5人回復(fù)
【求助】matlab 計算一個二重積分，高手們幫幫忙吧已經(jīng)有6人回復(fù)
求一定積分的解析解或數(shù)值解(數(shù)值解則要求給出matlab程序) 已經(jīng)有5人回復(fù)
怎樣用MATLAB計算出圖中白色圓形的面積已經(jīng)有6人回復(fù)
求助matlab積分的問題已經(jīng)有4人回復(fù)
【求助】催化反應(yīng)動力學(xué)matlab計算各基元反應(yīng)的速率常數(shù)時，該如何避免較小量被忽略？已經(jīng)有3人回復(fù)
【求助】matlab中保存計算數(shù)據(jù)的問題已經(jīng)有8人回復(fù)

1樓 2012-03-22 15:57:01

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝應(yīng)助 2012-03-23 02:57:02

q =1.0e+009 *[2.3916; -3.6135; 1.8065; -0.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))/1e10),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e10
Warning: Explicit integral could not be found.

F4 =

9999999998.5398333654497435804148

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2012-03-22 16:42:16

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

ldoop

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 3485.8
散金: 3151
紅花: 16
沙發(fā): 1
帖子: 2225
在線: 285.5小時
蟲號: 533053
注冊: 2008-03-25
專業(yè): 群落生態(tài)學(xué)

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝應(yīng)助 2012-03-23 02:57:07

應(yīng)該是值溢出了，建議稍微處理下表達(dá)式，比如分子都是 10^9乘以一個數(shù)，然后相加后除以10^10，那么直接10就可以了

贊一下

回復(fù)此樓

3樓2012-03-22 17:16:39

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

tianaylk

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線: 9小時
蟲號: 865837
注冊: 2009-10-08
專業(yè): 電力系統(tǒng)

引用回帖:

2樓: Originally posted by dbb627 at 2012-03-22 16:42:16:
q =1.0e+009 *;
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))/1e10),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e10
Warning: Explicit integral could not b ...

你好，非常感謝你的回復(fù)。但是你的變形是不等價的，因此得出的結(jié)果是有問題的，請再幫我考慮一下，謝謝！

贊一下

回復(fù)此樓

4樓2012-03-22 18:47:28

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

tianaylk

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線: 9小時
蟲號: 865837
注冊: 2009-10-08
專業(yè): 電力系統(tǒng)

引用回帖:

3樓: Originally posted by ldoop at 2012-03-22 17:16:39:
應(yīng)該是值溢出了，建議稍微處理下表達(dá)式，比如分子都是 10^9乘以一個數(shù)，然后相加后除以10^10，那么直接10就可以了

你好，非常感謝你的回復(fù)。你看到的可能是一樓作過變形后的表達(dá)式，但是那個變形和我的表達(dá)式是不等價的，原函數(shù)的形式請看主貼，貌似不太好處理...請再幫我考慮一下，謝謝！

贊一下

回復(fù)此樓

5樓2012-03-22 18:50:14

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★
jjdg: 金幣+1, 感謝應(yīng)助 2012-03-23 02:57:21

引用回帖:

4樓: Originally posted by tianaylk at 2012-03-22 18:47:28:
你好，非常感謝你的回復(fù)。但是你的變形是不等價的，因此得出的結(jié)果是有問題的，請再幫我考慮一下，謝謝！

確實(shí)是有問題
下面變形應(yīng)該等價
q =1.0e+008  *[2.3916; -3.6135; 1.8065; -0.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e8),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e8
Warning: Explicit integral could not be found.

F4 =

149272173.68701522388239715123227

>> q =1.0e+009  *[2.3916; -3.6135; 1.8065; -0.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e8),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e8
Warning: Explicit integral could not be found.

F4 =

log(numeric::int(1/exp(152725000*y^4 - (5051340117333333*y^3)/8388608 + 890400000*y^2 - (4908174540800001*y)/8388608 - 1248575000), y = 9/10..11/10)) + 100000000

>> q =1.0e+009  *[2.3916; -3.6135; 1.8065; -0.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e9),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e9
Warning: Explicit integral could not be found.

F4 =

1492721786.6977684879699111309597

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2012-03-22 19:34:19

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

tianaylk

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線: 9小時
蟲號: 865837
注冊: 2009-10-08
專業(yè): 電力系統(tǒng)

引用回帖:

6樓: Originally posted by dbb627 at 2012-03-22 19:34:19:
確實(shí)是有問題
下面變形應(yīng)該等價
q =1.0e+008 *;
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e8),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e8
Warnin ...

感謝版主解答！可是又發(fā)現(xiàn)問題了，如果把q的數(shù)值改變，例如把q(4)從-0.6109改為-5.6109，就又算不出結(jié)果了
q =1.0e+009  *[2.3916; -3.6135; 1.8065; -5.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e9),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e9
Warning: Explicit integral could not be found.
> In sym.int at 64

F4 =

log(numeric::int(1/exp(1402725000*y^4 - (5051340117333333*y^3)/8388608 - 6609600000*y^2 - (4908174540800001*y)/8388608 + 3401425000), y = 9/10..11/10)) + 1000000000

q =1.0e+009  *[2.3916; -3.6135; 1.8065; -5.6109];
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))-1e10),y,0.9,1.1);
F3=log(F2);
F4=vpa(F3)+1e10
Warning: Explicit integral could not be found.
> In sym.int at 64
??? Error using ==> mupadfeval at 28
Error: Singularity [ln];
during evaluation of 'mlfloat'

Error in ==> sym.vpa at 42
   r = mupadfeval('mlfloat',s);

也就是說隨著q的變化，指數(shù)中應(yīng)該減去10的多少次方不是固定的，應(yīng)該怎么辦呢？還請版主再幫忙解答一下，謝謝！

贊一下

回復(fù)此樓

7樓2012-03-22 20:28:48

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★
jjdg: 金幣+1, 感謝應(yīng)助 2012-03-23 02:57:31
余澤成: 金幣+3, 謝謝應(yīng)助，歡迎常來程序語言版！ 2012-03-23 16:58:09

引用回帖:

7樓: Originally posted by tianaylk at 2012-03-22 20:28:48:
感謝版主解答！可是又發(fā)現(xiàn)問題了，如果把q的數(shù)值改變，例如把q(4)從-0.6109改為-5.6109，就又算不出結(jié)果了
q =1.0e+009 *;
syms y ;
F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y ...

log(exp(100))

ans =

100

>> log(exp(1000))

ans =

Inf
這個是你計算的數(shù)字實(shí)在很大，matlab沒法表示，而你想通過減去一個數(shù)把它控制100左右，比較困難，我暫時還想不出較好的方法，也許你考慮將你的參數(shù)歸一化，然后積分，或者考慮換個軟件

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

8樓2012-03-22 21:33:10

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

tianaylk

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線: 9小時
蟲號: 865837
注冊: 2009-10-08
專業(yè): 電力系統(tǒng)

引用回帖:

8樓: Originally posted by dbb627 at 2012-03-22 21:33:10:
log(exp(100))

ans =

100

>> log(exp(1000))

ans =

Inf
這個是你計算的數(shù)字實(shí)在很大，matlab沒法表示，而你想通過減去一個數(shù)把它控制100左右，比較困難，我暫時還想不出較好的方法 ...

是啊，計算的確有點(diǎn)大，但其實(shí)取對數(shù)以后并不大，不知道該如何解決這個問題，已經(jīng)采用符號積分，但還是只能計算到q為1e8數(shù)量級，希望版主有空再幫忙考慮一下，非常感謝~

贊一下

回復(fù)此樓

9樓2012-03-22 22:11:09

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

397341529

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1
帖子: 3
在線: 56分鐘
蟲號: 1712227
注冊: 2012-03-23
性別: GG
專業(yè): 計算機(jī)應(yīng)用技術(shù)

看不懂啊

回復(fù)此樓

10樓2012-03-23 19:30:00

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

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

返回列表

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

[考研] 求材料調(diào)劑 +7	隔壁陳先生 2026-03-12	7/350	2026-03-18 16:14 by 楓橋ZL
[考研] 344求調(diào)劑 +5	knight344 2026-03-16	6/300	2026-03-18 15:17 by knight344
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian

24小時熱門版塊排行榜

[求助] matlab 積分計算 算不出結(jié)果，求幫助�。�！

» 猜你喜歡

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

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] matlab 積分計算算不出結(jié)果，求幫助�。�！

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