|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab解如下图的微分方程组: 我用的是dsolve语句。$ z4 x* t- [( V/ H
这个微分方程组中除了I,都是常数。而I则是一个变量,随着循环会不断改变。
& E- e \$ w! G0 [5 ?' m但是我发现dsolve语句解微分方程组,好像不能含有变量。
4 w2 n0 P1 H! I2 \所以想请问各路大神,我该如何使用Matlab解上图的微分方程组。
) Z% \0 E" [- N1 t' _我编写的程序如下:
. v/ i E" l1 D/ \c(1)=0.1;
+ {* t5 K! |" C% |" z4 vaa=17.26873*c(1)-3.04292;1 `+ H3 I7 m0 M# `1 S/ e# v4 m
bb=1.08304*c(1)-0.12098;( c9 ~0 a; t9 n8 `( L( Z6 y7 z
cc=0.09918*c(1)^(-0.43023);8 m9 T6 {& P1 k
z=0.05;
; `; t6 ]4 l n0 C( oi=aa+bb*exp(cc*z);& g/ I5 x- I) k2 Y
[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');
( C+ E+ t8 c. X) ]% I |
|