²é¿´: 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ʱµÄÎó²îÇúÏß¡¯);
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

hzlhm

ÖÁ×ðľ³æ (ÖøÃûдÊÖ)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¶À¹ÂÉñÓî: ½ð±Ò+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('Îó²î');
Èý¡¢Ð޸ĺóÔËÐеõ½ÈçϽá¹û
QQ£º2120156492
2Â¥2020-12-18 18:12:48
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

longwen8

ľ³æ (ÎÄ̳¾«Ó¢)

ÒýÓûØÌû:
2Â¥: Originally posted by hzlhm at 2020-12-18 18:12:48
Ò»¡¢¼¸¸ö×Ô¶¨Ò庯ÊýλÖôíÎó¡£Ó¦ÎªÕâÑù·ÅÖÃ
1¡¢Ö÷³ÌÐò £¬
function rungekutta()
clc,clear
x0=0;xn=1.0;y0=1;h=0.05;%hΪ²½³¤
=rungekutta4(x0,xn,y0,h);
end
2¡¢ ËĽ×Áú¸ñ¿âËþº¯Êý³ÌÐò
function =rungeku ...

ʲô½á¹û£¿´óÉñ¡£
3Â¥2020-12-18 19:43:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

hzlhm

ÖÁ×ðľ³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
3Â¥: Originally posted by longwen8 at 2020-12-18 19:43:08
ʲô½á¹û£¿´óÉñ¡£...

ÔËÐнá¹û£º
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
ͼÏñÉÏ´«²»ÉÏÀ´¡£
QQ£º2120156492
4Â¥2020-12-18 22:02:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

longwen8

ľ³æ (ÎÄ̳¾«Ó¢)

ÒýÓûØÌû:
4Â¥: Originally posted by hzlhm at 2020-12-18 22:02:15
ÔËÐнá¹û£º
i       x(i)       yiÊýÖµ½â        y(i)ÕæÖµ      Îó²î
-------------------------------------------------------------
1       0.0500     1.04880886     1.04880885   0.00000001
2      ...

Õâô¶à½á¹û£¬Ñ¡ÄÄÒ»¸ö¡£»¹ÊÇÕâÖÖÔËË㣬¾ÍûÓо¡Í·°¡

·¢×ÔСľ³æAndroid¿Í»§¶Ë
5Â¥2020-12-19 00:00:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

longwen8

ľ³æ (ÎÄ̳¾«Ó¢)

ÒýÓûØÌû:
2Â¥: Originally posted by hzlhm at 2020-12-18 18:12:48
Ò»¡¢¼¸¸ö×Ô¶¨Ò庯ÊýλÖôíÎó¡£Ó¦ÎªÕâÑù·ÅÖÃ
1¡¢Ö÷³ÌÐò £¬
function rungekutta()
clc,clear
x0=0;xn=1.0;y0=1;h=0.05;%hΪ²½³¤
=rungekutta4(x0,xn,y0,h);
end
2¡¢ ËĽ×Áú¸ñ¿âËþº¯Êý³ÌÐò
function =rungeku ...

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
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

hzlhm

ÖÁ×ðľ³æ (ÖøÃûдÊÖ)

QQ£º2120156492
7Â¥2020-12-21 12:41:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

longwen8

ľ³æ (ÎÄ̳¾«Ó¢)

ÒýÓûØÌû:
7Â¥: Originally posted by hzlhm at 2020-12-21 12:41:45
Ôõô·¢¸øÄ㣿

°Ù¶ÈÍøÅÌ»òÕßÓÊÏ䣺709923654@qq.com

·¢×ÔСľ³æAndroid¿Í»§¶Ë
8Â¥2020-12-21 13:58:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ longwen8 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 3/150 2026-03-15 17:32 by СÎïÀí»¯Ñ§
[¿¼ÑÐ] Öпƴó²ÄÁÏר˶319Çóµ÷¼Á +3 ÃÏöβÄÁÏ 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[»ù½ðÉêÇë] ÃæÉϺÍÇà»ùÒ»ÑùÏÞ30Ò³²»ºÏÀí +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤£¨0856£©304ÇóBÇøµ÷¼Á +7 Çñgl 2026-03-10 11/550 2026-03-14 12:18 by Çñgl
[¿¼ÑÐ] 308Çóµ÷¼Á +4 ÊÇLupa°¡ 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[¿¼ÑÐ] Çóµ÷¼Á +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[¿¼ÑÐ] bÇø»·¾³¹¤³ÌÇóµ÷¼Á +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[¿¼ÑÐ] 304Çóµ÷¼Á +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[¿¼ÑÐ] 0703»¯Ñ§Ò»Ö¾Ô¸211 ×Ü·Ö320Çóµ÷¼Á +5 Â꿨°Í¿¨°¡¹þ 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[¿¼ÑÐ] 329Çóµ÷¼Á +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[¿¼ÑÐ] 285»¯¹¤Ñ§Ë¶Çóµ÷¼Á£¨081700£© +6 ²ñ¿¤Ã¨_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[¿¼ÑÐ] 26µ÷¼Á/²ÄÁÏ¿ÆÑ§Ó빤³Ì/×Ü·Ö295/ÇóÊÕÁô +9 2026µ÷¼ÁÏÀ 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[¿¼ÑÐ] ¡¾¿¼Ñе÷¼ÁÇóÊÕÁô¡¿ +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[¿¼ÑÐ] ÇóbÇøÑ§Ð£µ÷¼Á +3 ÖÜ56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[¿¼ÑÐ] 310Çóµ÷¼Á +3 ¡¾ÉÏÉÏÇ©¡¿ 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[¿¼ÑÐ] 308Çóµ÷¼Á +3 ÊÇLupa°¡ 2026-03-12 3/150 2026-03-13 14:30 by Çóµ÷¼Ázz
[¿¼ÑÐ] 277Çóµ÷¼Á +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by °×Ò¹ÓÆ³¤
[¿¼ÑÐ] 321Çóµ÷¼Á£¨Ê³Æ·/ר˶£© +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[¿¼ÑÐ] 341Çóµ÷¼Á +4 µ·µ°ÖíÖí 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[¿¼ÑÐ] 279Çóµ÷¼Á +3 Īxiao 2026-03-10 4/200 2026-03-11 08:06 by Õ¶»êµÎÍÃ×Ó£¡
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û