| ²é¿´: 804 | »Ø¸´: 4 | |||
maggie0211гæ (СÓÐÃûÆø)
|
[ÇóÖú]
matlabÖÐÒ»µãСÎÊÌâ
|
|
³ÌÐòÈçÏ£º µÚÒ»¸ö£º function dx=ax3(t,x) dx=[x(2) -48.6*x(1)-1.26*x(2)+48.6*x(3)+21.6*sin(x(2)) 10*x(4) 1.95*x(1)-1.95*x(3)-0.333*sin(x(1)) x(6)+5.1734*(x(1)-x(5))-5.0069*(x(2)-x(6)) -48.6*x(5)-1.26*x(6)+48.6*x(7)+21.6*sin(x(6))-4.3407*(x(1)-x(5)) -4.1015*(x(2)-x(6)) 10*x(8) 1.95*x(5)-1.95*x(7)-0.333*sin(x(5))+1.0759*(x(1)-x(5))+1.0166*(x(2)-x(6))] µÚ¶þ¸ö£º tspan=[0 10]; x0=[1 2 1 0 1 2 0 0.5]'; [t,x]=ode15s(@ax3,tspan,x0); e1=x(:,1)-x(:,5); e2=x(:,2)-x(:,6); e3=x(:,3)-x(:,7); e4=x(:,4)-x(:,8); plot(t,x(:,1),'r*',t,x(:,5),'o') xlabel('t') legend('״̬Á¿x_{1}','״̬Á¿x_{1}µÄ¹À¼Æ') figure plot(t,e1) xlabel('t'); ylabel('״̬Á¿x_{1}ÓëÆä¹À¼ÆÁ¿µÄÎó²îe_{1}') figure plot(t,x(:,2),'r*',t,x(:,6),'o') xlabel('t') legend('״̬Á¿x_{2}','״̬Á¿x_{2}µÄ¹À¼Æ') figure plot(t,e2) xlabel('t') ylabel('״̬Á¿x_{2}ÓëÆä¹À¼ÆÁ¿µÄÎó²îe_{2}') figure plot(t,x(:,3),'r*',t,x(:,7),'o') xlabel('t') legend('״̬Á¿x_{3}','״̬Á¿x_{3}µÄ¹À¼Æ') figure plot(t,e3) xlabel('t'); ylabel('״̬Á¿x_{3}ÓëÆä¹À¼ÆÁ¿µÄÎó²îe_{3}') figure plot(t,x(:,4),'r*',t,x(:,8),'o') xlabel('t') legend('״̬Á¿x_{4}','״̬Á¿x_{4}µÄ¹À¼Æ') figure plot(t,e4) xlabel('t'); ylabel('״̬Á¿x_{4}ÓëÆä¹À¼ÆÁ¿µÄÎó²îe_{4}') Ö´Ðеڶþ¸ö³ÌÐòºó£¬³öÏÖÈçÏ´íÎó£º ??? Error using ==> funfun\private\odearguments Solving AX3 requires an initial condition vector of length 4. Error in ==> ode15s at 228 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Ó¦¸ÃÔõô¸Ä£¿ |
» ²ÂÄãϲ»¶
271Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
312Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
352Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
085602»¯Ñ§¹¤³ÌÇóµ÷¼Á¡£
ÒѾÓÐ3È˻ظ´
085600²ÄÁÏÓ뻯¹¤306
ÒѾÓÐ6È˻ظ´
»úеѧ˶310·Ö£¬ÊýÒ»Ó¢Ò»£¬Ò»Ö¾Ô¸211±¾¿ÆË«·ÇÕÒµ÷¼ÁÐÅÏ¢
ÒѾÓÐ3È˻ظ´
Ò»Ö¾Ô¸Î÷ÄϽ»Í¨ ר˶ ²ÄÁÏ355 ±¾¿ÆË«·Ç Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
085602 289·ÖÇóµ÷¼Á
ÒѾÓÐ8È˻ظ´
085601Çóµ÷¼Á×Ü·Ö293Ó¢Ò»Êý¶þ
ÒѾÓÐ4È˻ظ´
292Çóµ÷¼Á
ÒѾÓÐ10È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
×ܵª²â¶¨ÖУ¬±ê×¼ÇúÏßÎü¹âÖµÖÐÒ»¸öÊýֵƫСµÄÎÊÌâ¼°ÆäËûÎÊÌâ~
ÒѾÓÐ4È˻ظ´
·¢ÏÖMatlabÖÐÒ»¸öÓÐȤµÄÏÖÏó
ÒѾÓÐ4È˻ظ´
¡¾ÇóÖú¡¿ÇóÖúfortranÖÐÒ»¸öµ÷ÓÃÎļþµÄÎÊÌâ
ÒѾÓÐ3È˻ظ´
ÔÂÖ»À¶
Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)
-

ר¼Ò¾Ñé: +1059 - Ó¦Öú: 1712 (½²Ê¦)
- ¹ó±ö: 8.888
- ½ð±Ò: 68127.7
- É¢½ð: 1938
- ºì»¨: 443
- ɳ·¢: 4
- Ìû×Ó: 4373
- ÔÚÏß: 3291.4Сʱ
- ³æºÅ: 1122189
- ×¢²á: 2010-10-14
- רҵ: ÓîÖæÑ§
- ¹ÜϽ: ¼ÆËãÄ£ÄâÇø

2Â¥2013-07-03 11:34:34
maggie0211
гæ (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 57
- ºì»¨: 1
- Ìû×Ó: 53
- ÔÚÏß: 9.2Сʱ
- ³æºÅ: 2361581
- ×¢²á: 2013-03-20
3Â¥2013-07-03 12:16:32
Nonebull
ľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 13 (СѧÉú)
- ½ð±Ò: 5663.1
- ºì»¨: 6
- Ìû×Ó: 463
- ÔÚÏß: 306.9Сʱ
- ³æºÅ: 571775
- ×¢²á: 2008-06-10
- רҵ: »¯Ñ§¹¤³Ì¼°¹¤Òµ»¯Ñ§
4Â¥2013-07-08 03:41:21
Autoyang
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 0.5
- Ìû×Ó: 2
- ÔÚÏß: 1.1Сʱ
- ³æºÅ: 10278237
- ×¢²á: 2018-09-13
|
function xprim= xprim1(t,x) global a b c d e f q w r v s g xprim =[(e*s+a)*x(1)+(e*g+b)*x(2)+q*0.3*sin(t)+w*0.1*cos(2*t);(f*s+c)*x(1)+(f*g+d)*x(2)+r*0.3*sin(t)+v*0.1*cos(2*t)]; clear;clc; global a b c d e f q w r v s g a = 0; b = 1; c = 4; d = 0; e = 0; f = 1; q = 1; w = 0; r = 0; v = 1; s = 1; g = 3; [t,x] = ode45('xprim1',[0 9],[0.4;0.2]); plot(t,x) ³ö´í ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... ³ö´í event_trigger_ftb (line 10) [t,x] = ode45('xprim1',[0 9],[0.4;0.2]); ÄúºÃ£¬ÇëÎÊÕâ¸öÊÇÔõô»ØÊÂѽ? |
5Â¥2018-09-13 17:24:44













»Ø¸´´ËÂ¥