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

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

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

返回列表

下一頁(yè)

tianaylk

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 865837
注冊(cè): 2009-10-08
專(zhuān)業(yè): 電力系統(tǒng)

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

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專(zhuān)碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)

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

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

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

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

dbb627

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

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

【答案】應(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

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 7 (幼兒園)
金幣: 3485.8
散金: 3151
紅花: 16
沙發(fā): 1
帖子: 2225
在線(xiàn): 285.5小時(shí)
蟲(chóng)號(hào): 533053
注冊(cè): 2008-03-25
專(zhuān)業(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乘以一個(gè)數(shù)，然后相加后除以10^10，那么直接10就可以了

贊一下

回復(fù)此樓

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

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

tianaylk

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 865837
注冊(cè): 2009-10-08
專(zhuān)業(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à)的，因此得出的結(jié)果是有問(wèn)題的，請(qǐng)?jiān)賻臀铱紤]一下，謝謝！

贊一下

回復(fù)此樓

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

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

tianaylk

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 865837
注冊(cè): 2009-10-08
專(zhuān)業(yè): 電力系統(tǒng)

引用回帖:

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

你好，非常感謝你的回復(fù)。你看到的可能是一樓作過(guò)變形后的表達(dá)式，但是那個(gè)變形和我的表達(dá)式是不等價(jià)的，原函數(shù)的形式請(qǐng)看主貼，貌似不太好處理...請(qǐng)?jiān)賻臀铱紤]一下，謝謝！

贊一下

回復(fù)此樓

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

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

dbb627

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

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

【答案】應(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à)的，因此得出的結(jié)果是有問(wèn)題的，請(qǐng)?jiān)賻臀铱紤]一下，謝謝！

確實(shí)是有問(wèn)題
下面變形應(yīng)該等價(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))-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

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 865837
注冊(cè): 2009-10-08
專(zhuān)業(yè): 電力系統(tǒng)

引用回帖:

6樓: Originally posted by dbb627 at 2012-03-22 19:34:19:
確實(shí)是有問(wèn)題
下面變形應(yīng)該等價(jià)
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)問(wè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);

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

贊一下

回復(fù)此樓

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

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

dbb627

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

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

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

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

引用回帖:

7樓: Originally posted by tianaylk at 2012-03-22 20:28:48:
感謝版主解答！可是又發(fā)現(xiàn)問(wè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
這個(gè)是你計(jì)算的數(shù)字實(shí)在很大，matlab沒(méi)法表示，而你想通過(guò)減去一個(gè)數(shù)把它控制100左右，比較困難，我暫時(shí)還想不出較好的方法，也許你考慮將你的參數(shù)歸一化，然后積分，或者考慮換個(gè)軟件

贊一下

回復(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

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 141.5
帖子: 16
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 865837
注冊(cè): 2009-10-08
專(zhuān)業(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
這個(gè)是你計(jì)算的數(shù)字實(shí)在很大，matlab沒(méi)法表示，而你想通過(guò)減去一個(gè)數(shù)把它控制100左右，比較困難，我暫時(shí)還想不出較好的方法 ...

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

贊一下

回復(fù)此樓

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

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

397341529

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1
帖子: 3
在線(xiàn): 56分鐘
蟲(chóng)號(hào): 1712227
注冊(cè): 2012-03-23
性別: GG
專(zhuān)業(yè): 計(jì)算機(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 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專(zhuān)家開(kāi)始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥(niǎo)于飛Chao
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 求調(diào)劑（材料與化工327） +4	愛(ài)吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky

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

[求助] matlab 積分計(jì)算 算不出結(jié)果，求幫助！！！

» 猜你喜歡

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

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

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

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

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

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

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