|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。4 l: P4 l$ T' h* e& \
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。
B J7 h0 t8 j" O- f0 V# X4 W# J- @但是我发现dsolve语句解微分方程组,好像不能含有变量。
' {. N$ Y% `' l: e9 |! J所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。
$ N% G. g _, x我编写的程序如下:6 i. Q3 c. ?4 ?+ g; y6 }- W
c(1)=0.1;# d! O& f, Y9 }9 c6 Z
aa=17.26873*c(1)-3.04292;. }0 p6 O+ @2 n2 K
bb=1.08304*c(1)-0.12098;
& @6 @+ m9 U1 E) D, v$ bcc=0.09918*c(1)^(-0.43023);, h; N0 W5 E) I1 ? \$ _
z=0.05;
. t. u: O I4 G* r o9 i; ii=aa+bb*exp(cc*z);" p* _4 @; w- H% |0 g2 Z/ h
[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');* f; p! e, ?% h g
|
|