|  | 
 
| 
用Matlab解如下图的微分方程组:
x
EDA365欢迎您登录!您需要 登录 才可以下载或查看,没有帐号?注册  ![]() 我用的是dsolve语句。+ ~8 R/ c' j+ O2 o  L7 N/ N- K 这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。
 : Z7 a8 m8 H: G但是我发现dsolve语句解微分方程组,好像不能含有变量。
 ) j8 N3 w& `8 F. p/ I& r& m所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。- G; b; R7 R+ S$ q9 W
 我编写的程序如下:
 4 v: P$ E! Q, o1 m7 Uc(1)=0.1;
 , L+ O! `2 A8 x8 @aa=17.26873*c(1)-3.04292;
 6 a# V6 z4 W# g: p* Ybb=1.08304*c(1)-0.12098;
 / @& d; [7 k- Y. o9 |  G! w6 }; kcc=0.09918*c(1)^(-0.43023);( p( M4 b" b& y1 r; S, ?
 z=0.05;' u3 S' v; a/ D
 i=aa+bb*exp(cc*z);
 0 [, `$ r& v! K+ B[x y]=dsolve('Dx=-(7.253*10^(-4)*i+9.968*10^(-6))*x+(0.009485-9.968*10^(-6))*y+9.968*10^(-6),Dy=7.253*10^(-4)*i*x-(0.009485+1.703*10^(-8)*i)*y','x(0)=1,y(0)=0','t');
 . z' q# s' h0 @3 Q& [9 `
 | 
 |