24СʱÈÈÃŰæ¿éÅÅÐаñ    

²é¿´: 1807  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

¿Èâ»ðÉÕ

½ð³æ (СÓÐÃûÆø)

[ÇóÖú] matlab½¨Ä£ÖÐsolverµÄÎÊÌâ

%½¨Á¢·½³Ì
function dy=fai(x,y)
dy=2*(1/2*1.5*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*0.7*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500)
%ÔËÐÐmatlab³ÌÐòÇó½â
[X,Y]=ode45('fai',[0 0.008],[0 1.5],y(0)=0)
plot(x,y)

ÔËÐнá¹û
??? [X,Y]=ode45('fai',[0 0.008],[0 1.5],y(0)=0)
                                            |
Error: The expression to the left of the equals sign is not a valid target for an assignment.

Ïò¸÷λ´óÏÀÇóÖúÔõô¸Ä£¬Ð»Ð»£¡
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

change0618

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

·½ÕÉ´óʦ

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

¿Èâ»ðÉÕ(½ð±Ò+1): Ê®·Ö¸Ðл£¡£¡£¡ 2011-10-06 09:48:36
%½¨Á¢·½³Ì
function dy=fai(t,y)
global k1 k2
dy=2*(1/2*1.5*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*0.7*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500);
%ÔËÐÐmatlab³ÌÐòÇó½â
>> global k1 k2
>> k1=0;
>> k2=0.008;
>> [x,y]=ode45('fai',[0 15],0);
>> plot(x,y)
3Â¥2011-10-06 09:34:19
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

change0618

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

·½ÕÉ´óʦ

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

%½¨Á¢·½³Ì
function dy=fai(t,y,k1,k2)
dy=2*(1/2*1.5*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*0.7*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500);
%ÔËÐÐmatlab³ÌÐòÇó½â
[x,y]=ode45('fai',[0 15],0,[],0,0.008)
plot(x,y)
2Â¥2011-10-06 08:22:25
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

¿Èâ»ðÉÕ

½ð³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by change0618 at 2011-10-06 08:22:25:
%½¨Á¢·½³Ì
function dy=fai(t,y,k1,k2)
dy=2*(1/2*1.5*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*0.7*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*c ...

%ÔËÐÐmatlab³ÌÐòÇó½â
[x,y]=ode45(@fai,[0 15],0,[],0,0.008);
plot(x,y)


%½¨Á¢·½³Ì
function dy=fai(x,y,k1,k2)
dy=2*(1/2*1.5*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*0.7*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500);
ÒÔÉÏÊÇchange0618Õâ¸öÅóÓÑͨ¹ýÕ¾ÄÚÏûÏ¢°ïæ¸ÄÕýºóµÄÕýÈ·´ð°¸£¬Ð»Ð»£¡
4Â¥2011-10-06 09:41:55
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] 2026Äê8ÔÂ25ÈÕ¹ú×ÔÈ»·Å°ñǰͻȻÊÕµ½ÁÐÈëÆÀÉóר¼ÒÓʼþ£¬ÓйØÏµÂ𣿠+15 ľˮ˼¶¹ 2026-08-25 18/900 2026-08-25 23:28 by wj910103
[¿¼ÑÐ] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.55.1.O.5.4,¿ÆÄ¿È«,¿É+¼± +3 7K1CJE38xLG4 2026-08-25 3/150 2026-08-25 23:20 by cNXvBfCpiZOM
[»ù½ðÉêÇë] ½ñÌìÎñί»á¿ªÍêÁË£¬Ã÷Ìì³ö½á¹ûÂð +17 angus9576 2026-08-25 21/1050 2026-08-25 22:39 by ÓÆÓÆÃÀÅ®
[»ù½ðÉêÇë] ÓÐûÓдóÉñ°ïÎÒ¿´¿´»ù½ð´úÂë 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[»ù½ðÉêÇë] ·Å°ñǰµÄ²»µ­¶¨ 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[»ù½ðÉêÇë] 2026¹ú×ÔÈ»º¯ÆÀ·Ñµ½ÕË +20 ÑòÑü°å 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[»ù½ðÉêÇë] ½ñÈÕ²»·Å°ñ£¿Íø´«¹ú×ÔȻԤ¼Æ 8 Ô 27 Èտɲé½á¹û +17 ҽѧÀÏÄк¢ 2026-08-20 22/1100 2026-08-25 15:36 by ҽѧÀÏÄк¢
[½Ìʦ֮¼Ò] µ¼Ê¦Í²ۣºÎÒÔõô̯ÉÏÁËÕâô¸ö¼«Æ·Ñо¿Éú£¡ +3 ËÕ¶«ÆÂ¶þÊÀ 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[»ù½ðÉêÇë] ½ñÌì»ù½ð»á³ö½á¹ûÂð£¿20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[»ù½ðÉêÇë] ÎÒÃæÉÏÍêµ°ÁË +13 ÇÒÌý»¢Ð¥ 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[»ù½ðÉêÇë] Ö»ÓÐÿÄêÕâÖÖʱºòÀ´¹ä¹äСľ³æ +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[»ù½ðÉêÇë] ÄÜ·ñÍ˳ö²ÎÓëµÄÃæÉÏÏîÄ¿½â³ýÏÞÏî +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by ¼ÒÓëÔ¶·½
[»ù½ðÉêÇë] filecode£¬4¸öjtjcÁË +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by ¹þ¹þ¸ò£¿
[»ù½ðÉêÇë] ÈÃÎÒÖÐÒ»¸öÃæÉϰɣ¡ +13 ´óƼ1987 2026-08-20 16/800 2026-08-24 10:23 by ̫ɵÁË
[»ù½ðÉêÇë] ¿ÆÑй¶ùÌ«ÄÑÁË +18 ÎÒ4´ó°×²Ë 2026-08-20 19/950 2026-08-24 09:47 by ¿­¶÷¹ãÊ¢´ß»¯·ÖÎ
[½Ìʦ֮¼Ò] Ìø²ÛºóÔÚÑÐÏîÄ¿Ôõô°ì£¿ +5 ¼òµ¥»¯xn 2026-08-22 10/500 2026-08-23 12:38 by ¼òµ¥»¯xn
[»ù½ðÉêÇë] ½ñÌì·Å°ñÂ𣿠+15 ²¼²¼ºÍÒ»¶þ 2026-08-19 16/800 2026-08-23 09:55 by ÕÅ´ºÉú
[»ù½ðÉêÇë] ʱ¼ä´ÁÓÖ±äÁË +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by ×Ïɼ´¼
[»ù½ðÉêÇë] Ó¦¸ÃÊÇÏÂÖÜÈý26ÈÕ¹«²¼Á˰ɣ¿ +4 ¹þ¹þ¸ò£¿ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[»ù½ðÉêÇë] ½ñÌì·Å°ñûϷÁ衃 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û