| ²é¿´: 2858 | »Ø¸´: 13 | |||
| µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû | |||
value_1987гæ (³õÈëÎÄ̳)
|
[ÇóÖú]
MatlabÓÃÁú¸ñ¿âËþ·¨Çó¶þ½×·ÇÏßÐÔ·½³Ì
|
||
|
ÇóÖú½âÉÏʽ·½³Ì£¬ÓÃÁú¸ñ¿âËþ·½·¨£¬ÓÃ[t£¬y]=ode45(¡®fname¡¯,tspan,y0)£¬Èç¹ûu£¨t£©Êdz£Êý»òÕýÏұ仯µÈ£¬ÎÒ»á±à³ÌÐò£¬µ«ÊÇu£¨t£©ÊÇÒ»¸öËæ»úµÄÊý£¬ÓëtÓйصģ¬ÎÒ±àÁ˸öС³ÌÐòÀ´±íÊ¾Ëæ»úÊý£¬µ«ÊDz»ÖªµÀÔõô·Åµ½MÎļþÀïÃæ¡£ |
» ²ÂÄãϲ»¶
·¶½øÖоÙÒ»ÎĵÄÖÐÐÄ˼Ïë
ÒѾÓÐ9È˻ظ´
2026¹ú×ÔÈ»º¯ÆÀ·Ñµ½ÕË
ÒѾÓÐ24È˻ظ´
ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.55.1.O.5.4,¿ÆÄ¿È«,¿É+¼±
ÒѾÓÐ5È˻ظ´
2026Äê8ÔÂ25ÈÕ¹ú×ÔÈ»·Å°ñǰͻȻÊÕµ½ÁÐÈëÆÀÉóר¼ÒÓʼþ£¬ÓйØÏµÂð£¿
ÒѾÓÐ25È˻ظ´
ÊÛSCIÒ»ÇøÎÄÕ£¬ÎÒ:8.O.551.O.5.4,¿ÆÄ¿È«,¿ÉÙ¤¼±
ÒѾÓÐ6È˻ظ´
ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.55.1.O.5.4,¿ÆÄ¿È«,¿É+¼±
ÒѾÓÐ5È˻ظ´
ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8O.55.1.O.5.4,¿ÆÄ¿ÆëÈ«,¿É+¼±
ÒѾÓÐ4È˻ظ´
½ñÌìÎñί»á¿ªÍêÁË£¬Ã÷Ìì³ö½á¹ûÂð
ÒѾÓÐ22È˻ظ´
ÔÚ¼á±ù»¹¸Ç×ű±º£µÄʱºò£¬ÎÒ¿´µ½ÁËŷŵÄ÷»¨¡£
ÒѾÓÐ9È˻ظ´
Ã÷ÌìÓ¦¸Ã¿É²éÁË£¡£¿
ÒѾÓÐ6È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
Çómatlab ÓÃËĽ×Áú¸ñ-¿âËþ·¨Çó½â΢·Ö·½³Ì
ÒѾÓÐ4È˻ظ´
Áú¸ñ¿âËþ·¨ matlab
ÒѾÓÐ6È˻ظ´
ÒþʽŷÀ·¨Çó½âÒ»½×³£Î¢·Ö·½³Ì
ÒѾÓÐ7È˻ظ´
¡¾ÇóÖú¡¿MatlabÖÐÀûÓÃËĽ×Áú¸ñ-¿âËþ·¨Çó½â΢·Ö·½³Ì£¡£¡£¡£¡
ÒѾÓÐ9È˻ظ´
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
|
function [y,z]=Runge_kutta(a,b,y0,z0,h) x=a:h:b; y(1)=y0; z(1)=z0; n=(b-a)/h+1; for i=2:n K(1,1)=f1(x(i-1),y(i-1),z(i-1)); K(2,1)=f2(x(i-1),y(i-1),z(i-1)); K(1,2)=f1(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(2,2)=f2(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(1,3)=f1(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(2,3)=f2(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(1,4)=f1(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); K(2,4)=f2(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); y(i)=y(i-1)+h/6*(K(1,1)+2*K(1,2)+2*K(1,3)+K(1,4)); z(i)=z(i-1)+h/6*(K(2,1)+2*K(2,2)+2*K(2,3)+K(2,4)); end y(2); plot(y,'r') ×Ô¼º±àдÁËËæ»úÊýµÄMÎļþ£¬µ«ÊÇÿ¸öK¶¼Êǵ÷ÓÃÒ»´ÎMÎļþ£¬¾Í»á²úÉú²»Í¬Ëæ»úÊý£¬Ó¦¸ÃÒª²úÉúÏàÍ¬Ëæ»úÊý²ÅÐа¡ |
8Â¥2013-03-04 15:43:55
mechaniclmw
ľ³æ (СÓÐÃûÆø)
- Ó¦Öú: 9 (Ó×¶ùÔ°)
- ½ð±Ò: 2548
- É¢½ð: 197
- ºì»¨: 2
- Ìû×Ó: 292
- ÔÚÏß: 205.8Сʱ
- ³æºÅ: 2001977
- ×¢²á: 2012-09-15
- ÐÔ±ð: GG
- רҵ: ¶¯Á¦Ñ§Óë¿ØÖÆ
2Â¥2013-03-03 14:18:14
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
3Â¥2013-03-03 14:45:31
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
4Â¥2013-03-03 14:47:54









»Ø¸´´ËÂ¥
20