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

li_qun

銅蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 61
散金: 759
紅花: 3
帖子: 1240
在線: 419.7小時(shí)
蟲(chóng)號(hào): 692320
注冊(cè): 2009-01-14

[求助] 求助被積函數(shù)有絕對(duì)值號(hào)的問(wèn)題。

f=@(z1,z2)exp(abs(z1-z2));%對(duì)該函數(shù)進(jìn)行積分，積分范圍0<z1<3,1<z2<3

quad2d(@(z1,z2)f(z1,z2),0,3,1,3)

ans =

19.7563

>> %去掉絕對(duì)值號(hào)計(jì)算，對(duì)f函數(shù)的積分轉(zhuǎn)化成對(duì)f1，f2兩個(gè)函數(shù)的積分
f1=@(z1,z2)exp(z1-z2);%積分范圍0<z1<3,1<z2<z1
f2=@(z2,z1)exp(z2-z1);%積分范圍1<z2<3,0<z1<z1

c=@(z1)z1;
d=@(z2)z2;
quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)+quad2d(@(z2,z1)f2(z2,z1),1,3,0,d)

ans =

19.3884
末學(xué)已經(jīng)確定第二種方法是錯(cuò)誤的，因?yàn)閷?duì)f1的積分quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)來(lái)說(shuō)，并不表示1<z2<c，c有可能小于1
但不知道第一種方法是否正確？
如果要去掉絕對(duì)值號(hào)，應(yīng)該怎么改寫(xiě)這個(gè)函數(shù)？

==================================================
還有一個(gè)問(wèn)題，對(duì)函數(shù)f4積分

f4=@(z1,z2)exp(z1-z2);
quad2d(@(z1,z2)f4(z1,z2),2,5,3,8)

ans =

6.9739

但我用分段求和的方法來(lái)積分這個(gè)函數(shù)，結(jié)果如下
sum1=0;
N=1000;    %分段數(shù)
d1=(5-2)/N;  %變量z1的步長(zhǎng)
d2=(8-3)/N;  %變量z2的步長(zhǎng)
for i=1:N
sum2=0;
z1=(2+(i-1)*d1+2+i*d1)/2; %先固定z1，對(duì)z2積分，sum2表示對(duì)z2的積分結(jié)果，sum1表示對(duì)z1的積分結(jié)果
for j=1:1000
z2=(3+(i-1)*d2+3+i*d2)/2;
sum2=sum2+f4(z1,z2)*d2;
end
sum1=sum1+sum2*d1;
end
>> sum1

sum1 =

2.3857
這是什么原因，哪里有錯(cuò)誤嗎，還是精度產(chǎn)生的誤差？末學(xué)對(duì)matlab非常外行，請(qǐng)大賢指點(diǎn)，感激涕零！

回復(fù)此樓

» 猜你喜歡

【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有3人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有3人回復(fù)
一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑已經(jīng)有3人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有8人回復(fù)
291求調(diào)劑已經(jīng)有25人回復(fù)

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

一個(gè)絕對(duì)值求積分問(wèn)題已經(jīng)有8人回復(fù)
求助 matlab隱函數(shù)畫(huà)圖問(wèn)題已經(jīng)有4人回復(fù)
求助一個(gè)inline函數(shù)的問(wèn)題已經(jīng)有5人回復(fù)
高數(shù)：導(dǎo)函數(shù)與原函數(shù)問(wèn)題求助已經(jīng)有12人回復(fù)
請(qǐng)教一個(gè)積分用留數(shù)定理計(jì)算的問(wèn)題已經(jīng)有15人回復(fù)
求助 Origin 自定義函數(shù)輸入問(wèn)題已經(jīng)有4人回復(fù)
關(guān)于勢(shì)函數(shù)的求助已經(jīng)有23人回復(fù)
求助excel2007中normdist函數(shù)的使用問(wèn)題已經(jīng)有1人回復(fù)
【求助】linux信號(hào)處理函數(shù)中，一定不能用malloc這類(lèi)不可重用函數(shù)嗎？已經(jīng)有4人回復(fù)
【求助】delta函數(shù)平方后積分已經(jīng)有5人回復(fù)
【求助】一個(gè)求數(shù)組絕對(duì)值的小問(wèn)題已經(jīng)有7人回復(fù)
【求助】用graythresh函數(shù)遇到的問(wèn)題已經(jīng)有8人回復(fù)
【求助】相關(guān)函數(shù)歸一化問(wèn)題已經(jīng)有6人回復(fù)
【求助】請(qǐng)教mathematica中關(guān)于Plot函數(shù)的問(wèn)題【已完成】已經(jīng)有5人回復(fù)
【求助】求積分程序中被積函數(shù)問(wèn)題已經(jīng)有15人回復(fù)
【求助】Mittag-leffler函數(shù)問(wèn)題已經(jīng)有3人回復(fù)
【求助】一個(gè)求解反函數(shù)的問(wèn)題已經(jīng)有9人回復(fù)
【求助】只知道有限個(gè)函數(shù)值求積分已經(jīng)有10人回復(fù)
【求助】請(qǐng)教非平衡格林函數(shù)的問(wèn)題已經(jīng)有5人回復(fù)

1樓 2013-08-18 23:13:38

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

abacus

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 10 (幼兒園)
金幣: 15994.5
散金: 100
紅花: 11
帖子: 2273
在線: 639.9小時(shí)
蟲(chóng)號(hào): 553166
注冊(cè): 2008-05-02
專(zhuān)業(yè): 建筑環(huán)境與結(jié)構(gòu)工程

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與！ 2013-08-24 21:53:16

如圖所示，將積分區(qū)域分為兩部分即可。
代碼如下：

CODE:

%% Method 1 

f=@(z1,z2)exp(abs(z1-z2));%對(duì)該函數(shù)進(jìn)行積分，積分范圍0<z1<3,1<z2<3 an=quad2d(f,0,3,1,3) 

%% Method 2 

% 積分區(qū)域分為 兩部分 

bn=quad2d(@(z1,z2)exp(z1-z2),1,3,1,@(z1)z1); % z1>z2 cn=quad2d(@(z2,z1)exp(z2-z1),1,3,0,@(z2)z2); % z2<z1 

dn=bn+cn %最終結(jié)果

無(wú)標(biāo)題.png

贊一下

回復(fù)此樓

5樓2013-08-19 16:20:36

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

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

tgvictor99

木蟲(chóng) (知名作家)

迷糊領(lǐng)域榮譽(yù)領(lǐng)主

應(yīng)助: 40 (小學(xué)生)
金幣: 2667.4
散金: 1907
紅花: 49
沙發(fā): 1
帖子: 5238
在線: 923.8小時(shí)
蟲(chóng)號(hào): 2524003
注冊(cè): 2013-06-27
專(zhuān)業(yè): 流體力學(xué)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與！ 2013-08-24 21:53:02

樓主，你分段搞錯(cuò)了。
f2=@(z2,z1)exp(z2-z1);%積分范圍1<z2<3,0<z1<z1
f2的積分范圍應(yīng)該是1<z2<3,0<z1<z2

贊一下

回復(fù)此樓

每個(gè)人的人生都是時(shí)代的產(chǎn)物。每個(gè)人的血統(tǒng)、命運(yùn)和夢(mèng)想組成了現(xiàn)在這個(gè)時(shí)代。

2樓2013-08-19 14:44:48

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

li_qun

銅蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 61
散金: 759
紅花: 3
帖子: 1240
在線: 419.7小時(shí)
蟲(chóng)號(hào): 692320
注冊(cè): 2009-01-14

引用回帖:

2樓: Originally posted by tgvictor99 at 2013-08-19 14:44:48
樓主，你分段搞錯(cuò)了。
f2=@(z2,z1)exp(z2-z1);%積分范圍1<z2<3,0<z1<z1
f2的積分范圍應(yīng)該是1<z2<3,0<z1<z2

謝謝，注釋錯(cuò)了，程序里沒(méi)錯(cuò)

贊一下

回復(fù)此樓

3樓2013-08-19 14:53:48

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

tgvictor99

木蟲(chóng) (知名作家)

迷糊領(lǐng)域榮譽(yù)領(lǐng)主

應(yīng)助: 40 (小學(xué)生)
金幣: 2667.4
散金: 1907
紅花: 49
沙發(fā): 1
帖子: 5238
在線: 923.8小時(shí)
蟲(chóng)號(hào): 2524003
注冊(cè): 2013-06-27
專(zhuān)業(yè): 流體力學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
li_qun: 金幣+10, ★★★很有幫助, 謝謝，的確錯(cuò)了 2013-08-19 18:00:54
li_qun: 金幣+10, ★★★很有幫助 2013-08-23 08:48:58
xiegangmai: 金幣+2, 鼓勵(lì)討論交流！ 2013-08-24 21:53:09

引用回帖:

3樓: Originally posted by li_qun at 2013-08-19 14:53:48
謝謝，注釋錯(cuò)了，程序里沒(méi)錯(cuò)...

-5+exp(2)+exp(3)-exp(1)

ans =

19.7563

已經(jīng)手算出來(lái)，樓主第一種方法是對(duì)的。
第二種方法分塊錯(cuò)了，要分成3塊。

贊一下

回復(fù)此樓

每個(gè)人的人生都是時(shí)代的產(chǎn)物。每個(gè)人的血統(tǒng)、命運(yùn)和夢(mèng)想組成了現(xiàn)在這個(gè)時(shí)代。

4樓2013-08-19 15:25:03

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

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

[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑 +12	12.27 2026-03-30	14/700	2026-04-02 00:18 by 欣喜777
[考研] 085600 295分求調(diào)劑 +15	W55j 2026-03-30	18/900	2026-04-01 23:00 by 來(lái)看流星雨10
[考研] 304求調(diào)劑 +12	素年祭語(yǔ) 2026-03-31	15/750	2026-04-01 22:41 by peike
[考研] 材料求調(diào)劑 +8	呢呢妮妮 2026-04-01	10/500	2026-04-01 22:24 by salamander`
[考研] 江蘇科技大學(xué)招材料研究生 +4	Su032713. 2026-04-01	5/250	2026-04-01 22:03 by cccchenso
[考研] 286求調(diào)劑 +16	PolarBear11 2026-03-26	16/800	2026-04-01 21:31 by 七度不信任
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +6	pipiver 2026-03-30	6/300	2026-04-01 16:11 by oooqiao
[考研] 材料調(diào)劑 +10	Eujd1 2026-03-31	11/550	2026-04-01 11:23 by ivanqyq
[考研] 環(huán)境工程 085701，267求調(diào)劑 +15	minht 2026-03-29	16/800	2026-04-01 10:13 by li_sujuan99
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +6	篤行致遠(yuǎn). 2026-03-28	6/300	2026-04-01 09:28 by oooqiao
[考研] 085701環(huán)境工程，267求調(diào)劑 +17	minht 2026-03-26	17/850	2026-04-01 09:11 by xiayizhi
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 340求調(diào)劑 +4	希望如此i 2026-03-31	4/200	2026-03-31 16:40 by 690616278
[考研] 266分，求材料相關(guān)專(zhuān)業(yè)調(diào)劑 +10	哇呼哼呼哼 2026-03-30	12/600	2026-03-31 11:00 by 熊一刀
[考研] 085600，專(zhuān)業(yè)課化工原理，320分求調(diào)劑 +6	大饞小子 2026-03-29	6/300	2026-03-31 10:03 by 氯化亞硝酰
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8	披星河 2026-03-30	8/400	2026-03-31 08:53 by lbsjt
[考研] 297求調(diào)劑 +17	田洪有 2026-03-26	18/900	2026-03-30 18:32 by nothing投稿中
[考研] 332求92調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid