| ²é¿´: 2047 | »Ø¸´: 7 | ||
longwen8ľ³æ (ÎÄ̳¾«Ó¢)
|
[ÇóÖú]
΢·Ö·½³Ì×éµÄËĽ×Áú¸ñ¿âËþ¹«Ê½Çó½âmatlab°æ ÒÑÓÐ1È˲ÎÓë
|
|
΢·Ö·½³Ì×éµÄËĽ×Áú¸ñ¿âËþ¹«Ê½£¬Çó½âmatlab°æµÄ´úÂë¡£ÏÂÃæÊÇÎҵĴúÂ룬ÔËÐв»³öÀ´¡£ function varargout=rungekutta(varargin) clc,clear x0=0;xn=1.0;y0=1;h=0.05;%hΪ²½³¤ [y,x]=rungekutta4(x0,xn,y0,h); Function z=f(x,y); z=y-2*x/y; function [y,x]=rungekutta4(x0,xn,y0,h) x=x0:h:xn; n=(xn-x0)/h; y1=x; y1(1)=y0; for i=1:n %Áú¸ñ¿âËþ·¨µÄËã·¨ K1=f(x(i),y1(i)); K2=f(x(i)+h/2,y1(i)+h/2*K1); K3= f(x(i)+h/2,y1(i)+h/2*K2); K4= f(x(i)+h,y1(i)+h*K3); y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4); yy1(i)=(1+2*x(i+1)) ^0.5; error(i)=y1(i+1)-yy1(i); yy(i+1)=yy1(i); e(i+1)=error(i); end y=y1; n=(xn-x0)/h; fprintf(¡®i x(i) yiÊýÖµ½â y(i)ÕæÖµ Îó²î\n¡¯); fprintf(¡®-------------------------------------------------------------\n¡¯); for i=1:n fprintf(¡®%2d %12.4f %14.8f %14.8f %12.8f\n¡¯,i, x(i+1), y(i+1), yy(i+1), error(i)); end plot(x,e,¡®r.¡¯);%»³öÎó²îÇúÏß xlable(¡®xÖᡯ),ylable(¡®Îó²î¡¯); title(¡®²½³¤Îª0.05ʱµÄÎó²îÇúÏß¡¯); |
» ²ÂÄãϲ»¶
°ïÎÒµÄÓ¢Óï¿ÚÓïÀÏʦÕÒѧÉú
ÒѾÓÐ4È˻ظ´
¼ÆËãÈý½×Á¦³£ÊýÓÃthirdorderÀ©°ûʱ²ÎÊýÈçºÎÑ¡Ôñ£¿ÁÙ½üÔ×ÓÊýÔõôȡֵ
ÒѾÓÐ7È˻ظ´
ÎïÀíѧIÂÛÎÄÈóÉ«/·ÒëÔõôÊÕ·Ñ?
ÒѾÓÐ227È˻ظ´
2026ÄêÑ»·¾¼Ã¹¦ÄܲÄÁϹú¼Ê»áÒ飨ICFMCE 2026£©
ÒѾÓÐ0È˻ظ´
2026ÄêµÚÎå½ìµçÆø¡¢µç×ÓÓëÐÅÏ¢¹¤³Ì¹ú¼Ê»áÒ飨ISEEIE 2026£©
ÒѾÓÐ0È˻ظ´
¹þ¶û±õÀí¹¤´óѧÎïÀíϵÕÐÊÕÎïÀíѧ¿¼Ñе÷¼Á
ÒѾÓÐ4È˻ظ´
0702Ò»Ö¾Ô¸¼ª´óBÇøÇóµ÷¼Á
ÒѾÓÐ5È˻ظ´
Çóµ÷¼Á
ÒѾÓÐ0È˻ظ´
0702Ò»Ö¾Ô¸¼ª´óBÇøÇóµ÷¼ÁÓÐÂÛÎÄ
ÒѾÓÐ0È˻ظ´
ÇëÎÊ»¹ÓÐûÓÐÓÃLatexдÎÄÕµÄС»ï°éÃÇ£¿
ÒѾÓÐ0È˻ظ´
hzlhm
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- Ó¦Öú: 387 (˶ʿ)
- ½ð±Ò: 17808.5
- ºì»¨: 53
- Ìû×Ó: 2879
- ÔÚÏß: 606.3Сʱ
- ³æºÅ: 1517335
- ×¢²á: 2011-11-30
- ÐÔ±ð: GG
- רҵ: ³£Î¢·Ö·½³ÌÓ붯Á¦ÏµÍ³
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¶À¹ÂÉñÓî: ½ð±Ò+5, ¹ÄÀø½»Á÷ 2020-12-18 19:25:10
longwen8: ½ð±Ò+20, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸ 2020-12-18 19:42:37
¶À¹ÂÉñÓî: ½ð±Ò+5, ¹ÄÀø½»Á÷ 2020-12-18 19:25:10
longwen8: ½ð±Ò+20, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸ 2020-12-18 19:42:37
|
Ò»¡¢¼¸¸ö×Ô¶¨Ò庯ÊýλÖôíÎó¡£Ó¦ÎªÕâÑù·ÅÖà 1¡¢Ö÷³ÌÐò £¬ function rungekutta() clc,clear x0=0;xn=1.0;y0=1;h=0.05;%hΪ²½³¤ [y,x]=rungekutta4(x0,xn,y0,h); end 2¡¢ ËĽ×Áú¸ñ¿âËþº¯Êý³ÌÐò function [y,x]=rungekutta4(x0,xn,y0,h) ¡£¡£¡£ end 3¡¢zº¯Êý³ÌÐò function z=f(x,y); z=y-2*x/y; end ¶þ¡¢xlable(¡®xÖᡯ),ylable(¡®Îó²î¡¯);Õâ¾ä´úÂëÖк¯ÊýÊéд´íÎó¡£Ó¦Îª xlabel('xÖá'),ylabel('Îó²î'); Èý¡¢Ð޸ĺóÔËÐеõ½ÈçϽá¹û |

2Â¥2020-12-18 18:12:48
longwen8
ľ³æ (ÎÄ̳¾«Ó¢)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 5385.8
- É¢½ð: 13079
- ºì»¨: 48
- ɳ·¢: 83
- Ìû×Ó: 10804
- ÔÚÏß: 1496.5Сʱ
- ³æºÅ: 4254310
- ×¢²á: 2015-11-29
- ÐÔ±ð: GG
- רҵ: »·¾³¹¤³Ì
3Â¥2020-12-18 19:43:08
hzlhm
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- Ó¦Öú: 387 (˶ʿ)
- ½ð±Ò: 17808.5
- ºì»¨: 53
- Ìû×Ó: 2879
- ÔÚÏß: 606.3Сʱ
- ³æºÅ: 1517335
- ×¢²á: 2011-11-30
- ÐÔ±ð: GG
- רҵ: ³£Î¢·Ö·½³ÌÓ붯Á¦ÏµÍ³
|
ÔËÐнá¹û£º i x(i) yiÊýÖµ½â y(i)ÕæÖµ Îó²î ------------------------------------------------------------- 1 0.0500 1.04880886 1.04880885 0.00000001 2 0.1000 1.09544514 1.09544512 0.00000002 3 0.1500 1.14017546 1.14017543 0.00000004 4 0.2000 1.18321600 1.18321596 0.00000005 5 0.2500 1.22474493 1.22474487 0.00000006 6 0.3000 1.26491113 1.26491106 0.00000007 7 0.3500 1.30384056 1.30384048 0.00000008 8 0.4000 1.34164088 1.34164079 0.00000009 9 0.4500 1.37840498 1.37840488 0.00000011 10 0.5000 1.41421368 1.41421356 0.00000012 11 0.5500 1.44913781 1.44913767 0.00000014 12 0.6000 1.48323985 1.48323970 0.00000015 13 0.6500 1.51657526 1.51657509 0.00000017 14 0.7000 1.54919353 1.54919334 0.00000019 15 0.7500 1.58113904 1.58113883 0.00000021 16 0.8000 1.61245178 1.61245155 0.00000023 17 0.8500 1.64316793 1.64316767 0.00000026 18 0.9000 1.67332034 1.67332005 0.00000028 19 0.9500 1.70293895 1.70293864 0.00000031 20 1.0000 1.73205115 1.73205081 0.00000034 ͼÏñÉÏ´«²»ÉÏÀ´¡£ |

4Â¥2020-12-18 22:02:15
longwen8
ľ³æ (ÎÄ̳¾«Ó¢)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 5385.8
- É¢½ð: 13079
- ºì»¨: 48
- ɳ·¢: 83
- Ìû×Ó: 10804
- ÔÚÏß: 1496.5Сʱ
- ³æºÅ: 4254310
- ×¢²á: 2015-11-29
- ÐÔ±ð: GG
- רҵ: »·¾³¹¤³Ì
5Â¥2020-12-19 00:00:51
longwen8
ľ³æ (ÎÄ̳¾«Ó¢)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 5385.8
- É¢½ð: 13079
- ºì»¨: 48
- ɳ·¢: 83
- Ìû×Ó: 10804
- ÔÚÏß: 1496.5Сʱ
- ³æºÅ: 4254310
- ×¢²á: 2015-11-29
- ÐÔ±ð: GG
- רҵ: »·¾³¹¤³Ì
|
function rungekutta() clc,clear x0=0;xn=1.0;y0=1;h=0.05;%hΪ²½³¤ [y,x]=rungekutta4(x0,xn,y0,h); End function [y,x]=rungekutta4(x0,xn,y0,h) x=x0:h:xn; n=(xn-x0)/h; y1=x; y1(1)=y0; for i=1:n %Áú¸ñ¿âËþ·¨µÄËã·¨ K1=f(x(i),y1(i)); K2=f(x(i)+h/2,y1(i)+h/2*K1); K3= f(x(i)+h/2,y1(i)+h/2*K2); K4= f(x(i)+h,y1(i)+h*K3); y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4); yy1(i)=(1+2*x(i+1))^0.5; error(i)=y1(i+1)-yy1(i); yy(i+1)=yy1(i); e(i+1)=error(i); end function z=f(x,y); z=y-2*x/y; end y=y1; n=(xn-x0)/h; fprintf('i x(i) yiÊýÖµ½â y(i)ÕæÖµ Îó²î\n'); fprintf('-------------------------------------------------------------\n'); for i=1:n fprintf(¡®%2d %12.4f %14.8f %14.8f %12.8f\n¡¯,i, x(i+1), y(i+1), yy(i+1), error(i)); end plot(x,e, 'r. ');%»³öÎó²îÇúÏß xlabel('xÖá'),ylabel('Îó²î'); title('²½³¤Îª0.05ʱµÄÎó²îÇúÏß'); ´óÉñ£¬ÎÒÐÞ¸ÄÍ꣬»¹ÊÇÔËÐв»³öÀ´°¡¡£ÄãÄܲ»ÄÜÖ±½Ó°ÑÄãµÄÍêÕû´úÂë·¢¸øÎÒ£¬ÎÒ×Ô¼ºÕÒÕÒ´íÎó¡£ |
6Â¥2020-12-21 12:10:48
hzlhm
ÖÁ×ðľ³æ (ÖøÃûдÊÖ)
- Ó¦Öú: 387 (˶ʿ)
- ½ð±Ò: 17808.5
- ºì»¨: 53
- Ìû×Ó: 2879
- ÔÚÏß: 606.3Сʱ
- ³æºÅ: 1517335
- ×¢²á: 2011-11-30
- ÐÔ±ð: GG
- רҵ: ³£Î¢·Ö·½³ÌÓ붯Á¦ÏµÍ³

7Â¥2020-12-21 12:41:45
longwen8
ľ³æ (ÎÄ̳¾«Ó¢)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 5385.8
- É¢½ð: 13079
- ºì»¨: 48
- ɳ·¢: 83
- Ìû×Ó: 10804
- ÔÚÏß: 1496.5Сʱ
- ³æºÅ: 4254310
- ×¢²á: 2015-11-29
- ÐÔ±ð: GG
- רҵ: »·¾³¹¤³Ì
8Â¥2020-12-21 13:58:39













»Ø¸´´ËÂ¥