Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2010-10-09, 18:15   #1
zoom
初级会员
 
注册日期: 2010-08-30
帖子: 3
声望力: 0
zoom 正向着好的方向发展
默认 求一个微分方程组的解

遇到一个微分方程组的求解问题:
我编的程序如下:
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]);

不知道啥意思!请教论坛的高手!谢谢!
zoom 当前离线   回复时引用此帖
旧 2010-10-11, 10:50   #2
conquerking
初级会员
 
conquerking 的头像
 
注册日期: 2010-05-31
住址: 北京
帖子: 2
声望力: 0
conquerking 正向着好的方向发展
微笑 回复: 求一个微分方程组的解

帮你顶一下,我也很关心这个问题!
conquerking 当前离线   回复时引用此帖
旧 2010-10-13, 22:00   #3
reader1
初级会员
 
注册日期: 2008-01-07
年龄: 70
帖子: 9
声望力: 0
reader1 正向着好的方向发展
默认 回复: 求一个微分方程组的解

引用:
作者: zoom 查看帖子
遇到一个微分方程组的求解问题:
我编的程序如下:
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中少了一个括号。
太长了,不容易发现。
reader1 当前离线   回复时引用此帖
旧 2010-10-14, 00:34   #4
zoom
初级会员
 
注册日期: 2010-08-30
帖子: 3
声望力: 0
zoom 正向着好的方向发展
默认 回复: 求一个微分方程组的解

引用:
作者: reader1 查看帖子
可能是eq2中少了一个括号。
太长了,不容易发现。
我检查过了,应该没有少括号!但还是谢谢你!
zoom 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 20:15


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.