|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。
6 o4 y: `" L/ ]6 S: L这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。2 t/ S) C* X4 y* j8 L( X( Z
但是我发现dsolve语句解微分方程组,好像不能含有变量。
9 J3 R0 ~* d) a! L' L0 y2 I. J所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。. L n% S1 E+ `5 c% c2 @! a' v
我编写的程序如下:. i2 b- \+ a+ w
c(1)=0.1;$ u2 g a( o; m
aa=17.26873*c(1)-3.04292;
5 U% y8 Y2 f5 G3 f1 c: Qbb=1.08304*c(1)-0.12098;$ ^! K; Y! @) N; |0 i
cc=0.09918*c(1)^(-0.43023);
) }( Y& V" V9 d: G% m' e% lz=0.05;- m" b5 ?3 L, V2 t- W1 f) a+ p
i=aa+bb*exp(cc*z);5 M( H% v. g2 X. v% `% n
[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');
i, |9 |, h( m |
|