MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   数学建模专区 (https://www.labfans.com/bbs/forumdisplay.php?f=62)
-   -   [问题] 求助 matlab解常微分方程组答案有些不对怎么办 (https://www.labfans.com/bbs/showthread.php?t=13731)

leej 2011-07-26 08:23

求助 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??这不应该有的呀,求指点

lusunnycn 2011-07-26 12:26

回复: 求助 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)

leej 2011-07-27 09:02

回复: 求助 matlab解常微分方程组答案有些不对怎么办
 
是这样的啊,明明结果不对,找不到原因,急,求教!


所有时间均为北京时间。现在的时间是 19:44

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