|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。, z$ i A P) t# M) Y
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。0 R3 `0 w( R' u! A+ y2 _# f
但是我发现dsolve语句解微分方程组,好像不能含有变量。9 P% @1 g8 a% V+ T; c
所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。6 T- a% F, q5 {4 ^8 t- m
我编写的程序如下:4 B/ s) O3 r5 g
c(1)=0.1;3 B* g% r0 u" w
aa=17.26873*c(1)-3.04292;
* O1 n' ^ T7 ?bb=1.08304*c(1)-0.12098;
- o d* ^( ?8 ?4 Fcc=0.09918*c(1)^(-0.43023);
" ]' N3 g! W7 {& m. a$ s3 qz=0.05;- T6 L i, e& H- l0 L
i=aa+bb*exp(cc*z);7 g, h, N6 `% I4 R) d/ S4 m' H; T! ]
[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');
6 |& k' _! P s& R; N1 i- @ |
|