| 5 | 1/1 | 返回列表 |
| 查看: 1652 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
求助被積函數(shù)有絕對值號的問題。
|
||
|
f=@(z1,z2)exp(abs(z1-z2));%對該函數(shù)進(jìn)行積分,積分范圍0<z1<3,1<z2<3 quad2d(@(z1,z2)f(z1,z2),0,3,1,3) ans = 19.7563 >> %去掉絕對值號計算,對f函數(shù)的積分轉(zhuǎn)化成對f1,f2兩個函數(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)確定第二種方法是錯誤的,因為對f1的積分quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)來說,并不表示1<z2<c,c有可能小于1 但不知道第一種方法是否正確? 如果要去掉絕對值號,應(yīng)該怎么改寫這個函數(shù)? ================================================== 還有一個問題,對函數(shù)f4積分 f4=@(z1,z2)exp(z1-z2); quad2d(@(z1,z2)f4(z1,z2),2,5,3,8) ans = 6.9739 但我用分段求和的方法來積分這個函數(shù),結(jié)果如下 sum1=0; N=1000; %分段數(shù) d1=(5-2)/N; %變量z1的步長 d2=(8-3)/N; %變量z2的步長 for i=1:N sum2=0; z1=(2+(i-1)*d1+2+i*d1)/2; %先固定z1,對z2積分,sum2表示對z2的積分結(jié)果,sum1表示對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 這是什么原因,哪里有錯誤嗎,還是精度產(chǎn)生的誤差?末學(xué)對matlab非常外行,請大賢指點,感激涕零! |
木蟲 (知名作家)
迷糊領(lǐng)域榮譽領(lǐng)主

木蟲 (知名作家)
迷糊領(lǐng)域榮譽領(lǐng)主

至尊木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +9 | 小黃鴨寶 2026-03-30 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +17 | 荷包蛋hyj 2026-03-26 | 17/850 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 看看我!301求調(diào)劑 081700 +3 | 一個小明同學(xué) 2026-03-24 | 3/150 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
[考研] 英一數(shù)一408,總分284,二戰(zhàn)真誠求調(diào)劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|