|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。
; z3 u, E' Y. `/ U! q这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。5 v9 L" E8 V/ U6 z$ F. r( h* r! j& U
但是我发现dsolve语句解微分方程组,好像不能含有变量。
: t- }# D5 ?) _ J" K所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。; v n8 a- _0 ]- f- @
我编写的程序如下:
# ]1 Q& U& T& _5 cc(1)=0.1;0 d8 \* s7 ~8 A0 Y3 j$ E$ z
aa=17.26873*c(1)-3.04292;; G' Z e8 D0 }/ ~6 g; \
bb=1.08304*c(1)-0.12098;
& c ?/ O' K& k$ G7 U- g6 E9 hcc=0.09918*c(1)^(-0.43023);
9 B7 b* h/ G2 w* Mz=0.05;& [2 G1 s) x3 C' `- ? c
i=aa+bb*exp(cc*z);
; p) L$ ]7 c0 I2 n/ N7 C& e[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');
: L8 k, F0 i$ W* G: Z+ U! W |
|