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

glgjw

Òø³æ (СÓÐÃûÆø)

[ÇóÖú] matlabÇó´øÊ¼þµÄ΢·Ö·½³Ì×ܳö´í ÒÑÓÐ1È˲ÎÓë

>> opts=odeset('events',@lift)
>> [t,y,timeend]=ode45('mass1',[0 10],[0;0],opts)
ÉÏÃæµÄº¯ÊýliftºÍmass1±àдÈçÏ£º
º¯Êýlift£º
>> function[value,isterminal,direction]=lift(t,y)
value=160.5*y(1)+448*t-75
isterminal=1
direction=[]
º¯Êýmass1£º
function dy=mass1(t,y)
dy=[y(2);(4.48*10^6*8*t/15/(6-y(1))-(7.35*10^6+4.48*10^6/(6-y(1)))*y(1)-17.785*10^3*y(2))/4080]

ÔËÐкóµÄ³ö´íÐÅÏ¢£º
dy =

     0
     0

??? Function 'lower' is not defined for values of class 'function_handle'.

Error in ==> lower at 26
  [varargout{1:nargout}] = builtin('lower', varargin{:});

Error in ==> funfun\private\odeevents at 32
  switch lower(eventFcn)

Error in ==> ode45 at 209
[haveEventFcn,eventFcn,eventArgs,valt,teout,yeout,ieout] = ...

º¯Êýmass1¿Ï¶¨Ã»ÓдíÎ󣬴ó¼Ò¼ì²éÆäËûµÄ²¿·Ö¼´¿É¡£Ð»Ð»£¡
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

wurongjun

ר¼Ò¹ËÎÊ (Ö°Òµ×÷¼Ò)

ÒýÓûØÌû:
4Â¥: Originally posted by glgjw at 2016-02-18 11:04:42
ÇëÎÊÄú¾ÍÊǰÑ'mass1'¸ÄΪ@mass1ÁËÂð£¿ÎÒ֮ǰҲÕâô¸Ä¹ý£¬µ«»¹ÊdzöÏÖͬÑùµÄ´íÎó¡£Çë´Í½Ì£¡
...

ÊǵÄ!
ÉÆ¶ñµ½Í·ÖÕÓб¨,È˼äÕýµÀÊDz×É£.
5Â¥2016-02-18 11:16:38
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 9 ¸ö»Ø´ð
2Â¥2016-02-17 22:39:49
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wurongjun

ר¼Ò¹ËÎÊ (Ö°Òµ×÷¼Ò)

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
Ò»µãµãÎÊÌâ!
¼òµ¥ÐÞ¸ÄÈçÏÂ!
opts=odeset('events',@lift);
[t,y,timeend]=ode45(@mass1,[0 10],[0;0],opts)
function[value,isterminal,direction]=lift(t,y)
value=160.5*y(1)+448*t-75;
isterminal=1;
direction=[];
function dy=mass1(t,y)
dy=[y(2);(4.48*10^6*8*t/15/(6-y(1))-(7.35*10^6+4.48*10^6/(6-y(1)))*y(1)-17.785*10^3*y(2))/4080];
ÔËÐнá¹ûÈçÏÂ:
>>
t =

         0
    0.0040
    0.0080
    0.0120
    0.0160
    0.0200
    0.0240
    0.0280
    0.0320
    0.0367
    0.0415
    0.0463
    0.0510
    0.0562
    0.0614
    0.0666
    0.0718
    0.0773
    0.0828
    0.0884
    0.0939
    0.1005
    0.1070
    0.1136
    0.1202
    0.1246
    0.1291
    0.1336
    0.1380
    0.1425
    0.1470
    0.1514
    0.1559
    0.1581
    0.1603
    0.1625
    0.1648


y =

         0         0
    0.0000    0.0008
    0.0000    0.0030
    0.0000    0.0067
    0.0001    0.0117
    0.0001    0.0177
    0.0002    0.0247
    0.0003    0.0322
    0.0005    0.0402
    0.0007    0.0499
    0.0009    0.0593
    0.0012    0.0681
    0.0016    0.0759
    0.0020    0.0828
    0.0024    0.0879
    0.0029    0.0907
    0.0034    0.0913
    0.0039    0.0896
    0.0044    0.0854
    0.0048    0.0792
    0.0052    0.0713
    0.0057    0.0606
    0.0060    0.0492
    0.0063    0.0382
    0.0065    0.0283
    0.0067    0.0228
    0.0067    0.0184
    0.0068    0.0152
    0.0069    0.0135
    0.0069    0.0132
    0.0070    0.0143
    0.0071    0.0168
    0.0072    0.0205
    0.0072    0.0227
    0.0073    0.0252
    0.0073    0.0278
    0.0074    0.0307


timeend =

    0.1648

>>
ÉÆ¶ñµ½Í·ÖÕÓб¨,È˼äÕýµÀÊDz×É£.
3Â¥2016-02-18 09:54:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

glgjw

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
3Â¥: Originally posted by wurongjun at 2016-02-18 09:54:32
Ò»µãµãÎÊÌâ!
¼òµ¥ÐÞ¸ÄÈçÏÂ!
opts=odeset('events',@lift);
=ode45(@mass1,,,opts)
function=lift(t,y)
value=160.5*y(1)+448*t-75;
isterminal=1;
direction=[];
function dy=mass1(t,y)
dy=;
ÔËÐнá¹ûÈç ...

ÇëÎÊÄú¾ÍÊǰÑ'mass1'¸ÄΪ@mass1ÁËÂð£¿ÎÒ֮ǰҲÕâô¸Ä¹ý£¬µ«»¹ÊdzöÏÖͬÑùµÄ´íÎó¡£Çë´Í½Ì£¡

·¢×ÔСľ³æAndroid¿Í»§¶Ë
4Â¥2016-02-18 11:04:42
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ÉúÎïѧ308Çóµ÷¼Á +4 ÏàÐűػá¹ââÍòÕ 2026-04-11 4/200 2026-04-12 11:48 by babysonlkd
[¿¼ÑÐ] 291Çóµ÷¼Á +11 ¹ØÒä±±. 2026-04-09 12/600 2026-04-12 10:32 by ÄæË®³Ë·ç
[¿¼ÑÐ] Çóµ÷¼Á +11 ÔÂ@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[¿¼ÑÐ] 303Çóµ÷¼Á +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 366Çóµ÷¼Á +7 ²»ÖªÃûµÄСئ 2026-04-11 7/350 2026-04-11 16:45 by PeggyPeng17
[¿¼ÑÐ] 300·ÖÇóµ÷¼Á £¨085501»úеר˶£¬±¾¿ÆÑï´ó£© +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[¿¼ÑÐ] ũҵ¹ÜÀí302·ÖÇóµ÷¼Á +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[¿¼ÑÐ] Çóµ÷¼Á288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[¿¼ÑÐ] 314Çóµ÷¼Á +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[¿¼ÑÐ] һ־Ը³¶«´óѧ071000ÉúÎïѧѧ˶³õÊÔ·ÖÊý276Çóµ÷¼Á +20 Ľ¾øcc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[¿¼ÑÐ] Òѵ÷¼Á +18 ²ñ¿¤Ã¨_ 2026-04-09 19/950 2026-04-09 22:10 by ²ñ¿¤Ã¨_
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤301·ÖÇóµ÷¼ÁԺУ +33 ´ÌÍ´jk 2026-04-06 34/1700 2026-04-09 18:31 by hy861222
[¿¼ÑÐ] 348Çóµ÷¼Á +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by ¶ÎΰÑÞ
[¿¼ÑÐ] Ò»Ö¾Ô¸ÏôóÉúÎïѧ332Çóµ÷¼Á +10 ³Ø³Ø³Ø³Ø³Ø³Ø 2026-04-08 10/500 2026-04-09 17:10 by ¶À×íÃι³Ç
[¿¼ÑÐ] 086004 Çóµ÷¼Á 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[¿¼ÑÐ] ²ÄÁϹ¤³Ì322 +18 ¹þ¹þ¹þºðºðºð¹þ 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[¿¼ÑÐ] 0860004 Çóµ÷¼Á 309·Ö +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by °¡Àî999
[¿¼ÑÐ] 285Çóµ÷¼Á +12 AZMK 2026-04-05 18/900 2026-04-08 20:43 by ÄæË®³Ë·ç
[¿¼ÑÐ] 307Çóµ÷¼Á +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by СºÚ²»ÅÂÄÑ
[¿¼ÑÐ] 285Çóµ÷¼Á +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by À¶ÔÆË¼Óê
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û