±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1633  |  »Ø¸´: 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

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

·½ÕÉ´óʦ

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

%½¨Á¢·½³Ì
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µÄ»ØÌû

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µÄ»ØÌû

¿Èâ»ðÉÕ

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

ÒýÓûØÌû:
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µÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ¿Èâ»ðÉÕ µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 291Çóµ÷¼Á +6 ‹üÈA 2026-03-22 6/300 2026-03-27 19:56 by WYUMater
[¿¼ÑÐ] »¯Ñ§µ÷¼Á +4 °®³Ô·¬ÇѵÄÐñ 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[²ÄÁϹ¤³Ì] Ò»Ö¾Ô¸C9²ÄÁÏÓ뻯¹¤×¨Òµ×Ü·Ö300Çóµ÷¼Á +7 Âü111 2026-03-24 8/400 2026-03-27 13:53 by ÑîÑîÑî×Ï
[¿¼ÑÐ] Ò»Ö¾Ô¸»ª¶«Àí¹¤´óѧ081700£¬³õÊÔ·ÖÊý271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by »ÝÖÝÅíÓÚêÌ
[¿¼ÑÐ] 333Çóµ÷¼Á +3 questionÍì·ç 2026-03-23 3/150 2026-03-27 11:29 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 292Çóµ÷¼Á +4 ÇóÇóÁËÊÕÏÂÎÒ°É£ 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[¿¼ÑÐ] 324Çóµ÷¼Á +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[¿¼ÑÐ] 071000ÉúÎïѧÇóµ÷¼Á£¬³õÊԳɼ¨343 +6 ССÌðÃæÍÅ 2026-03-25 6/300 2026-03-26 23:01 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] Çóµ÷¼Á +6 °×QF 2026-03-21 6/300 2026-03-26 20:37 by fmesaito
[¿¼ÑÐ] ²ÄÁϵ÷¼Á 5+4 ÏëÒªÒ»ºøÌÒ»¨Ë® 2026-03-25 10/500 2026-03-26 19:56 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] Öйú¿ÆÑ§ÔºÉîÛÚÏȽø¼¼ÊõÑо¿Ôº-¹âÏË´«¸Ð¿ÎÌâ×éÕÐÉú-Öйú¿ÆÑ§Ôº´óѧ¡¢ÉîÛÚÀí¹¤´óѧÁªÅà +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by èßäèßäѽ
[¿¼ÑÐ] 081700 µ÷¼Á 267·Ö +11 ÃÔÈ˵Ĺþ¹þ 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤328·Öµ÷¼Á +6 ¡££¬¡££¬¡££¬¡£i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[¿¼ÑÐ] 26¿¼ÑÐ-291·Ö-ÏÃÃÅ´óѧ£¨085601£©-ÈáÐÔµç×ÓѧԺ²ÄÁϹ¤³ÌרҵÇóµ÷¼Á +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[¿¼ÑÐ] 0854AI CV·½ÏòÕÐÊÕµ÷¼Á +4 ÕÂСÓã567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[¿¼ÑÐ] »¯Ñ§µ÷¼Á +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[¿¼ÑÐ] 291Çóµ÷¼Á +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[ÂÛÎÄͶ¸å] ¼±·¢ºËÐÄÆÚ¿¯ÂÛÎÄ +3 ÏÍ´ïÎʽò 2026-03-23 5/250 2026-03-23 17:13 by ÃÃ×Ó²»ºÃÈÇ
[¿¼ÑÐ] 275Çóµ÷¼Á +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[¿¼ÑÐ] 336Çóµ÷¼Á +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by ѧԱ8dgXkO
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û