|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。
: f" _, T" r1 ?" D' P3 @, ?这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。2 w2 C' K4 W ^8 A- `2 {$ B
但是我发现dsolve语句解微分方程组,好像不能含有变量。: `6 O, K9 l) U! d H2 n4 G
所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。9 G' z! _+ e- S2 ?. h% Z: {1 U
我编写的程序如下:( e* E3 g" ^" s
c(1)=0.1;
" {# J3 S4 W$ ~' E& |aa=17.26873*c(1)-3.04292;
* k, r; g$ m: C2 v% F- mbb=1.08304*c(1)-0.12098;
( l v. _, D; L' Z U4 {7 D2 Zcc=0.09918*c(1)^(-0.43023);
( i: C& }3 A8 L' z3 [3 Bz=0.05;/ X/ s+ k) t f
i=aa+bb*exp(cc*z);
) B. \) H3 S" P6 C! U[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');. J( S* d5 Y: g x
|
|