|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。4 d% S, g3 T1 `
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。# K0 r% b7 A) I+ o
但是我发现dsolve语句解微分方程组,好像不能含有变量。9 M# i6 ?4 K* n8 d
所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。0 {; t2 c: T+ j9 i9 x5 @! S" g0 _
我编写的程序如下:
0 n- l6 T( h- sc(1)=0.1;
5 o6 |# U" E+ M/ Kaa=17.26873*c(1)-3.04292;+ _6 ]! E8 w' q
bb=1.08304*c(1)-0.12098;% U# @! Y+ Q$ E) v: _* U! C. Y
cc=0.09918*c(1)^(-0.43023);8 ^. j3 A# X. G* D8 w) }
z=0.05;* h2 p, L; r7 R' }3 L, w
i=aa+bb*exp(cc*z);# F3 n9 D/ B9 z; n, I$ Q$ ~* M
[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');0 A# g* }& g* E+ a
|
|