| 4 | 1/1 | 返回列表 |
| 查看: 1370 | 回復(fù): 3 | ||||
ghost118木蟲 (小有名氣)
|
[求助]
matlab龍貝格求積 已有1人參與
|
|
小弟按照書上的算法程序輸入了matlab,被積函數(shù)是ln(1+x)/x,積分區(qū)間是0到1,下面是算法程序: a=0;b=1;eps=0.5*10^-6; M=1;h=b-a;err=1;J=0;R=zeros(4,4); R(1,1)=h*(1+f(b))/2; while(err>eps) J=J+1; h=h/2; x=a+h:2*h:b-h; R(J+1,1)=R(J,1)/2+h*sum(f(x)); for K=1:min(3,J) R(J+1,K+1)=R(J+1,K)+(R(J+1,K)-R(J,K))/(4^K-1); end if(J>3) err=abs(R(J+1,4)-R(J,4)); end end quad=R(J,4); 運(yùn)行出來的結(jié)果和實(shí)際相差很大,而且用這個(gè)程序做書上的例題,可以做出來,是不是被積函數(shù)是反常積分的緣故呢,希望能有大神從中指點(diǎn),不勝感激 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
自適應(yīng)龍貝格積分: quadl(@(x) log(1+x)./x,0,1) ans = 0.822467034106751 自適應(yīng)Gauss-Kronrod 積分: quadgk(@(x) log(1+x)./x,0,1) ans = 0.822467033424113 |

木蟲 (小有名氣)
木蟲 (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|---|---|---|---|---|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 307求調(diào)劑 +9 | 超級(jí)伊昂大王 2026-03-24 | 10/500 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 6/300 |
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +6 | 101次希望 2026-03-28 | 6/300 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +7 | WW.' 2026-03-23 | 9/450 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|