Znn3bq.jpeg
²é¿´: 1382  |  »Ø¸´: 3

ghost118

ľ³æ (СÓÐÃûÆø)

[ÇóÖú] matlabÁú±´¸ñÇó»ý ÒÑÓÐ1È˲ÎÓë

СµÜ°´ÕÕÊéÉϵÄËã·¨³ÌÐòÊäÈëÁËmatlab£¬±»»ýº¯ÊýÊÇln(1+x)/x,»ý·ÖÇø¼äÊÇ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);
ÔËÐгöÀ´µÄ½á¹ûºÍʵ¼ÊÏà²îºÜ´ó£¬¶øÇÒÓÃÕâ¸ö³ÌÐò×öÊéÉϵÄÀýÌ⣬¿ÉÒÔ×ö³öÀ´£¬ÊDz»ÊDZ»»ýº¯ÊýÊÇ·´³£»ý·ÖµÄÔµ¹ÊÄØ£¬Ï£ÍûÄÜÓдóÉñ´ÓÖÐÖ¸µã£¬²»Ê¤¸Ð¼¤

[ ·¢×ÔÊÖ»ú°æ http://muchong.com/3g ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ghost118: ½ð±Ò+30, ¡ïÓаïÖú 2014-11-15 21:00:15
×ÔÊÊÓ¦Áú±´¸ñ»ý·Ö£º
quadl(@(x)  log(1+x)./x,0,1)

ans =

   0.822467034106751

×ÔÊÊÓ¦Gauss-Kronrod »ý·Ö£º
quadgk(@(x)  log(1+x)./x,0,1)

ans =

   0.822467033424113
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2014-11-12 18:49:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ghost118

ľ³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by ÔÂÖ»À¶ at 2014-11-12 18:49:39
×ÔÊÊÓ¦Áú±´¸ñ»ý·Ö£º
quadl(@(x)  log(1+x)./x,0,1)

ans =

   0.822467034106751

×ÔÊÊÓ¦Gauss-Kronrod »ý·Ö£º
quadgk(@(x)  log(1+x)./x,0,1)

ans =

   0.822467033424113

¸ç£¬´ð°¸ÊǶԵ쬵«ÊÇÓдúÂëô£¬»òÕß°ïÎÒÕâ´úÂë¸Ä¸Ä£¬¿ì½»ÁË£¬¹òÇó°¡

[ ·¢×ÔСľ³æ¿Í»§¶Ë ]
3Â¥2014-11-13 00:16:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ghost118

ľ³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
3Â¥: Originally posted by ghost118 at 2014-11-13 00:16:51
¸ç£¬´ð°¸ÊǶԵ쬵«ÊÇÓдúÂëô£¬»òÕß°ïÎÒÕâ´úÂë¸Ä¸Ä£¬¿ì½»ÁË£¬¹òÇó°¡
...

ÎÒÃÇÀÏʦҪ¿´´úÂë

[ ·¢×ÔСľ³æ¿Í»§¶Ë ]
4Â¥2014-11-13 00:17:33
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ghost118 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 211±¾¿Æ²ÄÁÏ»¯¹¤Çóµ÷¼Á +15 YHLAH 2026-04-11 16/800 2026-04-12 12:44 by BruceLiu320
[¿¼ÑÐ] Ò»Ö¾Ô¸°²´óÉúÎïѧ07³õÊÔ322¡¢±¾¿Æ¶þ±¾¡¢µ÷¼ÁÇóÖú +5 Àî¶àÃ×lee. 2026-04-12 6/300 2026-04-12 09:59 by ѩɽ·Éºü7233
[¿¼ÑÐ] 310Çóµ÷¼Á +12 666ÕæºÃ 2026-04-11 13/650 2026-04-12 08:40 by Sealedwind
[¿¼ÑÐ] 271Çóµ÷¼Á +20 2261744733 2026-04-11 22/1100 2026-04-11 23:14 by labixiaoqiao
[¿¼ÑÐ] 295Çóµ÷¼Á +6 ?ÒªÉϰ¶? 2026-04-05 7/350 2026-04-11 19:02 by laoshidan
[¿¼ÑÐ] µ÷¼ÁÇóÖú +6 ¹ûÈ»ÓÐÎÒ 2026-04-11 7/350 2026-04-11 16:22 by Ã÷Ô´ËʱÓÐ
[¿¼ÑÐ] »¹Óл¯¹¤¶þÂÖµ÷¼ÁµÄѧУÂð 5+14 »¯¹¤ÈË999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[¿¼ÑÐ] ±¾¿Æ211 ¹¤¿Æ085400 280·ÖÇóµ÷¼Á ¿É¿çרҵ +11 LZH£¨µÈ´ýµ÷¼ÁÖÐ 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[¿¼ÑÐ] 085506-Çóµ÷¼Á-285·Ö +3 À×Å··ÉÌß 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[¿¼ÑÐ] 071000ÉúÎïѧµ÷¼ÁÇóÖú +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by Çóµ÷¼Ázz
[¿¼ÑÐ] »·¾³×¨Ë¶µ÷¼Á +16 »á˵»°µÄÖâ×Ó 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[¿¼ÑÐ] 347²ÄÁÏר˶Çóµ÷¼Á +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[¿¼ÑÐ] ¿¼Ñе÷¼Á-²ÄÁÏÀà-284 +28 Ïë»»ÊÖ»ú²»Ïë½âÊ 2026-04-08 28/1400 2026-04-09 20:08 by µ¹Êý321?
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖпÆÔº105500רҵ×Ü·Ö315Çóµ÷¼Á +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[¿¼ÑÐ] µ÷¼Á +12 ÔÂ@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[¿¼ÑÐ] 283µç×ÓÐÅÏ¢Çóµ÷¼Á +4 ÈýʯWL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[¿¼ÑÐ] ¶þ´Îµ÷¼ÁÇóÀÏʦÊÕÁô +3 ЦЦԬ 2026-04-08 3/150 2026-04-08 23:50 by ×íÔÚ·çÀï
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +4 ö©??? 2026-04-08 4/200 2026-04-08 21:44 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 285Çóµ÷¼Á +12 AZMK 2026-04-05 18/900 2026-04-08 20:43 by ÄæË®³Ë·ç
[¿¼ÑÐ] 0854Çóµ÷¼Á +9 ºàÊÏ·¬ÇÑɳ˾ 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û