| 5 | 1/1 | 返回列表 |
| 查看: 1654 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
求助被積函數(shù)有絕對(duì)值號(hào)的問題。
|
|||
|
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)該怎么改寫這個(gè)函數(shù)? ================================================== 還有一個(gè)問題,對(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),感激涕零! |
木蟲 (知名作家)
迷糊領(lǐng)域榮譽(yù)領(lǐng)主
|
-5+exp(2)+exp(3)-exp(1) ans = 19.7563 已經(jīng)手算出來(lái),樓主第一種方法是對(duì)的。 第二種方法分塊錯(cuò)了,要分成3塊。 |

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

至尊木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑推薦 +6 | 清酒714 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 08工科275分求調(diào)劑 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +6 | GENJIOW 2026-03-31 | 6/300 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 307分求調(diào)劑 +6 | (o~o) 2026-03-31 | 6/300 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考博] 材料專業(yè)申博 +5 | 杜雨婷dyt 2026-03-29 | 5/250 |
|
|
[考研] 一志愿大連理工大學(xué),機(jī)械工程學(xué)碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|