|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。8 B8 m$ k$ h8 }2 t$ D& f
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。4 c/ w: K, M# e" {& c
但是我发现dsolve语句解微分方程组,好像不能含有变量。" [( N! g0 ]6 ]: H3 _: L7 O% l- R- W
所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。! j* Y0 m+ f0 Q3 o* l
我编写的程序如下:
5 W; O3 |" m9 u6 e3 ~c(1)=0.1;
! \8 X4 e5 [' p9 N. z+ }9 taa=17.26873*c(1)-3.04292;
. v/ l$ v% I+ g0 z; cbb=1.08304*c(1)-0.12098;
+ J8 Z8 I3 z3 L4 I" zcc=0.09918*c(1)^(-0.43023);0 l1 l/ s$ U$ g+ N3 C) V- `. }- ? k8 B
z=0.05;
. I& b, J* g- }% g* gi=aa+bb*exp(cc*z);# @ a. E( E& V: S3 Z8 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');
' S' L/ n$ o; m! k, z |
|