Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 基础科学 > 数学 > 数学建模专区
数学建模专区 A discussion board for Mathematical Modelling.
回复
 
主题工具 显示模式
旧 2011-07-26, 08:23   #1
leej
初级会员
 
注册日期: 2011-06-30
帖子: 2
声望力: 0
leej 正向着好的方向发展
帖子 求助 matlab解常微分方程组答案有些不对怎么办

本人系新手用的win7系统,结果发现解方程时结果异于或复杂于计算值,不知是安装问题还是别的,举例见下,特此求解,哪位大侠指点下谢谢了
>> syms x0 y0
>> [x,y]=dsolve('Dx=y','Dy=-2*x','x(0)=x0','y(0)=y0')


x =

(exp(2^(1/2)*i*t)*(2*x0 - 2^(1/2)*i*y0))/4 + (2*x0 + 2^(1/2)*i*y0)/(4*exp(2^(1/2)*i*t))


y =

(2^(1/2)*i*exp(2^(1/2)*i*t)*(2*x0 - 2^(1/2)*i*y0))/4 - (2^(1/2)*i*(2*x0 + 2^(1/2)*i*y0))/(4*exp(2^(1/2)*i*t))
里边出现了i??这不应该有的呀,求指点
leej 当前离线   回复时引用此帖
旧 2011-07-26, 12:26   #2
lusunnycn
初级会员
 
lusunnycn 的头像
 
注册日期: 2011-07-25
帖子: 3
声望力: 0
lusunnycn 正向着好的方向发展
微笑 回复: 求助 matlab解常微分方程组答案有些不对怎么办

根据你的操作,我也进行了求解,但是结果跟你的完全不同
>> syms x0 y0
>> [x,y]=dsolve('Dx=y','Dy=-2*x','x(0)=x0','y(0)=y0')

x =

x0*cos(2^(1/2)*t) + (2^(1/2)*y0*sin(2^(1/2)*t))/2


y =

y0*cos(2^(1/2)*t) - 2^(1/2)*x0*sin(2^(1/2)*t)
lusunnycn 当前离线   回复时引用此帖
旧 2011-07-27, 09:02   #3
leej
初级会员
 
注册日期: 2011-06-30
帖子: 2
声望力: 0
leej 正向着好的方向发展
默认 回复: 求助 matlab解常微分方程组答案有些不对怎么办

是这样的啊,明明结果不对,找不到原因,急,求教!
leej 当前离线   回复时引用此帖
回复

主题工具
显示模式

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

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



所有时间均为北京时间。现在的时间是 03:10


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