主题: [MATLAB基础] 符号解三次方程有误?
查看单个帖子
旧 2012-11-10, 23:51   #3
zsy312
普通会员
 
注册日期: 2008-08-17
年龄: 43
帖子: 77
声望力: 18
zsy312 正向着好的方向发展
默认 回复: 符号解三次方程有误?

引用:
作者: 旁听生 查看帖子
谢谢回复,但你计算的结果不对。
这个三次方程判别式小于零,有三个不相等的实数根,请参见我主题贴里的附件图片。
我把x^2错输成x了。这个符号解如果引用maple函数就能得到化简后的结果:
代码:
maple('s:=evalc([solve(x^3+3*x^2-1,x)])')
maple('simplify(s)')
代码:
ans =

s := [1/2*8^(1/3)*cos(2/9*pi)+1/4*8^(2/3)*cos(2/9*pi)-1+i*(1/2*8^(1/3)*sin(2/9*pi)-1/4*8^(2/3)*sin(2/9*pi)), -1/4*8^(1/3)*cos(2/9*pi)-1/8*8^(2/3)*cos(2/9*pi)-1-1/2*3^(1/2)*(1/2*8^(1/3)*sin(2/9*pi)+1/4*8^(2/3)*sin(2/9*pi))+i*(-1/4*8^(1/3)*sin(2/9*pi)+1/8*8^(2/3)*sin(2/9*pi)+1/2*3^(1/2)*(1/2*8^(1/3)*cos(2/9*pi)-1/4*8^(2/3)*cos(2/9*pi))), -1/4*8^(1/3)*cos(2/9*pi)-1/8*8^(2/3)*cos(2/9*pi)-1+1/2*3^(1/2)*(1/2*8^(1/3)*sin(2/9*pi)+1/4*8^(2/3)*sin(2/9*pi))+i*(-1/4*8^(1/3)*sin(2/9*pi)+1/8*8^(2/3)*sin(2/9*pi)-1/2*3^(1/2)*(1/2*8^(1/3)*cos(2/9*pi)-1/4*8^(2/3)*cos(2/9*pi)))]


ans =

[2*cos(2/9*pi)-1, -cos(2/9*pi)-1-3^(1/2)*sin(2/9*pi), -cos(2/9*pi)-1+3^(1/2)*sin(2/9*pi)]
第二个ans还是没有虚部的。说明化简的还行。
zsy312 当前离线   回复时引用此帖