|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。- J$ t9 q( F0 O# p, b, ~7 Z
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。3 g: T/ |4 s; L
但是我发现dsolve语句解微分方程组,好像不能含有变量。
' H2 v3 L3 o8 x% S+ ]" D4 o所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。5 D+ f+ P: B& G# r1 W4 H c
我编写的程序如下:( C9 v: h1 j$ D# q7 u% a1 t
c(1)=0.1;9 P, {# E4 Q: t% N* m7 k; K, E
aa=17.26873*c(1)-3.04292;
$ I6 x4 n Z' Wbb=1.08304*c(1)-0.12098;
: G8 i4 W$ P3 S( O" ucc=0.09918*c(1)^(-0.43023);+ w: J' I4 ]$ p/ f. @
z=0.05;
' p# E. r+ r; xi=aa+bb*exp(cc*z);6 [7 L: b l: Y; U6 S# P
[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');! C3 ~' f- `- t3 P' f ^
|
|