| ²é¿´: 689 | »Ø¸´: 1 | |||
| µ±Ç°Ö÷ÌâÒѾ´æµµ¡£ | |||
[½»Á÷]
¡¾ÇóÖú¡¿¹ØÓÚÊýѧ½¨Ä£¶öÀÇ×·ÍÃÎÊÌâµÄ½â´ð·½·¨
|
|||
|
ËÄ¡¢¸ß½×³£Î¢·Ö·½³ÌÄ£ÐÍ¡ª¶öÀÇ×·ÍÃÎÊÌâ ÏÖÓÐÒ»Ö»ÍÃ×Ó¡¢Ò»Æ¥ÀÇ£¬ÍÃ×ÓλÓÚÀǵÄÕýÎ÷100Ã×´¦£¬¼ÙÉèÍÃ×ÓÓëÀÇͬʱ·¢ÏÖ¶Ô·½²¢Ò»ÆðÆðÅÜ£¬ÍÃ×ÓÍùÕý±±60Ã×´¦µÄ³²Ñ¨ÅÜ£¬¶øÀÇÔÚ×·ÍÃ×Ó¡£ÒÑÖªÍÃ×Ó¡¢ÀÇÊÇÔÈËÙÅÜÇÒÀǵÄËÙ¶ÈÊÇÍÃ×ÓµÄÁ½±¶¡£ ÒªÇ󣺣¨1£©½¨Á¢ÀǵÄÔ˶¯¹ì¼£Î¢·ÖÄ£ÐÍ¡£ £¨2£©»³öÍÃ×ÓÓëÀǵÄÔ˶¯¹ì¼£Í¼ÐΡ£ £¨3£©ÓýâÎö·½·¨Çó½â£¬ÎÊÍÃ×ÓÄÜ·ñ°²È«»Øµ½³²Ñ¨£¿ £¨4£©ÓÃÊýÖµ·½·¨Çó½â£¬ÎÊÍÃ×ÓÄÜ·ñ°²È«»Øµ½³²Ñ¨£¿ ¡¾×¢¡¿³£Î¢·Ö·½³Ì¸ß½×³õÖµÎÊÌâµÄMATLAB¿âº¯ÊýΪ£ºode45¡£ Ó﷨Ϊ£º[t,Y] =ode45(odefun,tspan,y0) ÀýÈ纯Êý£º function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); ÉèÖÃÑ¡Ï options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); Çó½âµÃ£º [t,Y] = ode45(@rigid,[0 12],[0 1 1],options); »³ö½âº¯ÊýÇúÏßͼÐΣº plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') ÖÔÐÄлл¸÷λÓÐʶ֮ʿµÄ°ïÖú£¡~£¡£¡ [ Last edited by csfn on 2008-11-30 at 12:19 ] |
» ²ÂÄãϲ»¶
Ò»Ö¾Ô¸»ª¶«Ê¦·¶´óѧÓлú»¯Ñ§×¨Òµ£¬³õÊÔ351·Ö£¬¸´ÊÔ±»Ë¢Çóµ÷¼Á!
ÒѾÓÐ4È˻ظ´
298Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
²ÄÁÏѧ˶333Çóµ÷¼Á
ÒѾÓÐ11È˻ظ´
0856Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
356Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
332Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
348Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
305Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
Ò»Ö¾Ô¸±±¾©»¯¹¤´óѧ²ÄÁÏÓ뻯¹¤£¨085600£©296Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
µ÷¼Á310
ÒѾÓÐ9È˻ظ´
Ѿ³æ
½ð³æ (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 417
- Ìû×Ó: 58
- ÔÚÏß: 38.9Сʱ
- ³æºÅ: 557982
- ×¢²á: 2008-05-13
- רҵ: ΢ÉúÎïÒÅ´«Ñ§
2Â¥2008-09-19 16:32:30














»Ø¸´´ËÂ¥