|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。
1 h$ l! a5 U) E" e# _( I4 r这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。2 F( [ ?+ k; Z, V" Y
但是我发现dsolve语句解微分方程组,好像不能含有变量。
7 f& H6 I( O( V1 p$ T* A所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。. }. o K! r+ b6 g2 e; x2 @
我编写的程序如下:
3 p! F7 y2 a3 U$ `c(1)=0.1;
\" D" u! T7 D# g. gaa=17.26873*c(1)-3.04292;, G- g/ ~4 W! v7 J
bb=1.08304*c(1)-0.12098;( U( G% x. n6 J
cc=0.09918*c(1)^(-0.43023);
7 M) u4 `: j p5 F9 P! l6 f3 j: Fz=0.05;
( E& \4 ~7 r0 d3 u7 w* di=aa+bb*exp(cc*z);0 ^' C3 M7 d1 g* K
[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');
) @; D2 V2 Z' [. e; s# V4 X/ k |
|