| ²é¿´: 2852 | »Ø¸´: 13 | ||
value_1987гæ (³õÈëÎÄ̳)
|
[ÇóÖú]
MatlabÓÃÁú¸ñ¿âËþ·¨Çó¶þ½×·ÇÏßÐÔ·½³Ì
|
|
|
ÇóÖú½âÉÏʽ·½³Ì£¬ÓÃÁú¸ñ¿âËþ·½·¨£¬ÓÃ[t£¬y]=ode45(¡®fname¡¯,tspan,y0)£¬Èç¹ûu£¨t£©Êdz£Êý»òÕýÏұ仯µÈ£¬ÎÒ»á±à³ÌÐò£¬µ«ÊÇu£¨t£©ÊÇÒ»¸öËæ»úµÄÊý£¬ÓëtÓйصģ¬ÎÒ±àÁ˸öС³ÌÐòÀ´±íÊ¾Ëæ»úÊý£¬µ«ÊDz»ÖªµÀÔõô·Åµ½MÎļþÀïÃæ¡£ |
» ²ÂÄãϲ»¶
2026Äê8ÔÂ25ÈÕ¹ú×ÔÈ»·Å°ñǰͻȻÊÕµ½ÁÐÈëÆÀÉóר¼ÒÓʼþ£¬ÓйØÏµÂð£¿
ÒѾÓÐ9È˻ظ´
½ñÌìÎñί»á¿ªÍêÁË£¬Ã÷Ìì³ö½á¹ûÂð
ÒѾÓÐ18È˻ظ´
ÔÚ¼á±ù»¹¸Ç×ű±º£µÄʱºò£¬ÎÒ¿´µ½ÁËŷŵÄ÷»¨¡£
ÒѾÓÐ9È˻ظ´
Ã÷ÌìÓ¦¸Ã¿É²éÁË£¡£¿
ÒѾÓÐ6È˻ظ´
ÓÐûÓдóÉñ°ïÎÒ¿´¿´»ù½ð´úÂë
ÒѾÓÐ10È˻ظ´
·Å°ñǰµÄ²»µ¶¨
ÒѾÓÐ19È˻ظ´
ijЩ»ú¹¹£¬ÒÔЧÂʵÍΪÈÙ£¬ÒÔЧÂʵÍ×÷Ϊ´æÔÚ¸Ð
ÒѾÓÐ10È˻ظ´
2026¹ú×ÔÈ»º¯ÆÀ·Ñµ½ÕË
ÒѾÓÐ23È˻ظ´
½ñÈÕ²»·Å°ñ£¿Íø´«¹ú×ÔȻԤ¼Æ 8 Ô 27 Èտɲé½á¹û
ÒѾÓÐ22È˻ظ´
ʲôʱºò¿ª½±£¿
ÒѾÓÐ13È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
Çómatlab ÓÃËĽ×Áú¸ñ-¿âËþ·¨Çó½â΢·Ö·½³Ì
ÒѾÓÐ4È˻ظ´
Áú¸ñ¿âËþ·¨ matlab
ÒѾÓÐ6È˻ظ´
ÒþʽŷÀ·¨Çó½âÒ»½×³£Î¢·Ö·½³Ì
ÒѾÓÐ7È˻ظ´
¡¾ÇóÖú¡¿MatlabÖÐÀûÓÃËĽ×Áú¸ñ-¿âËþ·¨Çó½â΢·Ö·½³Ì£¡£¡£¡£¡
ÒѾÓÐ9È˻ظ´
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
5Â¥2013-03-03 16:37:39
ϦÑôÎ÷ÏÂ
½ð³æ (ÖøÃûдÊÖ)
- Ó¦Öú: 138 (¸ßÖÐÉú)
- ¹ó±ö: 0.041
- ½ð±Ò: 9247.1
- É¢½ð: 73
- ºì»¨: 28
- ɳ·¢: 1
- Ìû×Ó: 1268
- ÔÚÏß: 1164.4Сʱ
- ³æºÅ: 220750
- ×¢²á: 2006-03-20
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
- ¹ÜϽ: ÎïÀí
6Â¥2013-03-04 10:14:56
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
change0618
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
·½ÕÉ´óʦ
- Ó¦Öú: 44 (СѧÉú)
- ½ð±Ò: 17724.5
- ºì»¨: 17
- Ìû×Ó: 2413
- ÔÚÏß: 546.7Сʱ
- ³æºÅ: 496517
- ×¢²á: 2008-01-19
- רҵ: »¯Ñ§·´Ó¦¹¤³Ì
7Â¥2013-03-04 13:03:10
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
9Â¥2013-03-04 15:45:22
value_1987
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17
- Ìû×Ó: 25
- ÔÚÏß: 5.1Сʱ
- ³æºÅ: 1417825
- ×¢²á: 2011-09-26
- רҵ: ¹¤³ÌÈÈÎïÀíÏà¹Ø½»²æÁìÓò
10Â¥2013-03-04 15:56:50









»Ø¸´´ËÂ¥
20