| ²é¿´: 177 | »Ø¸´: 0 | |||
| µ±Ç°Ö÷ÌâÒѾ´æµµ¡£ | |||
sspxiaojiÒø³æ (ÕýʽдÊÖ)
|
[½»Á÷]
¡¾ÇóÖú¡¿matlab½â΢·Ö·½³Ì×é
|
||
|
ÓÃ[t,y]=ode45('nov',[0 1],[1 1])½â³£Î¢·Ö·½³Ì×éʱ³öÏÖ£º ??? Inputs to odearguments must be floats, namely single or double ÒªÈçºÎÐ޸ģ¿ MÎļþΪ£º function dydt=nov7(t,y) global A B n digits(7); beta=0.5; T=sym('[y(1) y(2) y(3) y(4) y(5) y(6) y(7) y(8) y(9) y(10) y(11) y(12) y(13) y(14) y(15) y(16) y(17) y(18) y(19) y(20) y(21)]'); T(n+1)=0.05; for i=1:n dy(i)=T(i)^beta*(sum(B(i,1:n).*T(1:n))+B(i,n+1)*T(n+1))+beta*T(i)^(-beta)*(sum(A(i,1:n).*T(1:n))-A(i,n+1)*T(n+1))^2; end dydt=vpa(transpose(dy)); ÎÊÌâºÃÏñ³öÔÚÊý×ÖºÍy£¨1£©¡£¡£¡£y£¨2£©µÄÀàÐͲ»Í¬£¬µ«ÎÒ²»ÖªµÀÈçºÎ°ÑËûÃÇÀàÐ͸ijÉÒ»Ö Çë¸÷λ³æÓÑÖ¸½Ì [ Last edited by cuplgz on 2008-1-18 at 14:18 ] |
» ²ÂÄãϲ»¶
ҩѧ105500ר˶289·ÖÇóµ÷¼Á
ÒѾÓÐ3È˻ظ´
0703×Ü·Ö331Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
Ò»Ö¾Ô¸±±¾©2£¬²ÄÁÏÓ뻯¹¤308Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
318Çóµ÷¼Á
ÒѾÓÐ14È˻ظ´
285Çóµ÷¼Á
ÒѾÓÐ10È˻ظ´
312Çóµ÷¼Á
ÒѾÓÐ11È˻ظ´
273Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£©
ÒѾÓÐ8È˻ظ´
ÉúÎïѧÇóµ÷¼Á Ò»Ö¾Ô¸»¦9£¬326·Ö
ÒѾÓÐ3È˻ظ´
080500Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´














»Ø¸´´ËÂ¥