登录论坛

查看完整版本 : [MATLAB基础] 求一个微分方程组的解


zoom
2010-10-09, 18:15
遇到一个微分方程组的求解问题:
我编的程序如下:
syms A B C D E F G y z t
eq1='Dy=A+B*y-C*y^2';
eq2='Dz=z*(-B+2*C*y-2*D)+z*(((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^0.5+((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^(-0.5))-exp(-0.05*t)';
[y,z]=dsolve('eq1','eq2','y(0)=2000','z(0)=0','t')

可是运行以后,程序报错
??? Error using ==> maple
at offset 2, `}` unexpected

Error in ==> dsolve at 145
var_set = maple([var_set ' intersect ' var_set]);

不知道啥意思!请教论坛的高手!谢谢!

conquerking
2010-10-11, 10:50
帮你顶一下,我也很关心这个问题!

reader1
2010-10-13, 22:00
遇到一个微分方程组的求解问题:
我编的程序如下:
syms A B C D E F G y z t
eq1='Dy=A+B*y-C*y^2';
eq2='Dz=z*(-B+2*C*y-2*D)+z*(((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^0.5+((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^(-0.5))-exp(-0.05*t)';
[y,z]=dsolve('eq1','eq2','y(0)=2000','z(0)=0','t')

可是运行以后,程序报错
??? Error using ==> maple
at offset 2, `}` unexpected

Error in ==> dsolve at 145
var_set = maple([var_set ' intersect ' var_set]);

不知道啥意思!请教论坛的高手!谢谢!
可能是eq2中少了一个括号。
太长了,不容易发现。

zoom
2010-10-14, 00:34
可能是eq2中少了一个括号。
太长了,不容易发现。

我检查过了,应该没有少括号!但还是谢谢你!