| 6 | 1/1 | 返回列表 |
| 查看: 1651 | 回復(fù): 5 | ||
[求助]
求助被積函數(shù)有絕對值號的問題。
|
|
f=@(z1,z2)exp(abs(z1-z2));%對該函數(shù)進行積分,積分范圍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ā)一次吧。結(jié)果說明,你的第一次計算是正確的。 %% Method 1 f=@(z1,z2)exp(abs(z1-z2));%對該函數(shù)進行積分,積分范圍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é)果 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 085404 22408 315分 +4 | zhuangyan123 2026-03-31 | 4/200 |
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +4 | zzll03 2026-03-31 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 研s. 2026-03-25 | 4/200 |
|
|
[考研] 調(diào)劑申請 +5 | 張張張張zy 2026-03-31 | 5/250 |
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑,不限專業(yè),有學(xué)上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|
[考研] 286分調(diào)劑 +10 | Faune 2026-03-30 | 11/550 |
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[碩博家園] 招收生物學(xué)/細胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|