大家可以跑一下
syms x;
A_pl =x^24 - 24*x^22 + 252*x^20 - 1522*x^18 + 5850*x^16 - 14958*x^14 + 25845*x^12 - 30042*x^10 + 22923*x^8 - 10894*x^6 + 2913*x^4 - 360*x^2 + 16; %待解的方程
A_ans = solve(A_pl);%求解
real(A_ans);%只取实部,可是输出的结果,虚部还在,好奇怪哦?找不到原因。
%对于简单一点的方程却可以把实部求出来,难道是bug?
ans =
2
1
1
1
2
1
-2
-2
-1
-1
-1
-1
1/((3^(1/2)*i)/2 + 1/2)^(1/3) + ((3^(1/2)*i)/2 + 1/2)^(1/3)
1/((3^(1/2)*i)/2 + 1/2)^(1/3) + ((3^(1/2)*i)/2 + 1/2)^(1/3)
1/((3^(1/2)*i)/2 - 1/2)^(1/3) + ((3^(1/2)*i)/2 - 1/2)^(1/3)
1/((3^(1/2)*i)/2 - 1/2)^(1/3) + ((3^(1/2)*i)/2 - 1/2)^(1/3)
(3^(1/2)*(1/((3^(1/2)*i)/2 + 1/2)^(1/3) - ((3^(1/2)*i)/2 + 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) + 1/2)^(1/3)) - ((3^(1/2)*i)/2 + 1/2)^(1/3)/2
(3^(1/2)*(1/((3^(1/2)*i)/2 + 1/2)^(1/3) - ((3^(1/2)*i)/2 + 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) + 1/2)^(1/3)) - ((3^(1/2)*i)/2 + 1/2)^(1/3)/2
- (3^(1/2)*(1/((3^(1/2)*i)/2 + 1/2)^(1/3) - ((3^(1/2)*i)/2 + 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) + 1/2)^(1/3)) - ((3^(1/2)*i)/2 + 1/2)^(1/3)/2
- (3^(1/2)*(1/((3^(1/2)*i)/2 + 1/2)^(1/3) - ((3^(1/2)*i)/2 + 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) + 1/2)^(1/3)) - ((3^(1/2)*i)/2 + 1/2)^(1/3)/2
(3^(1/2)*(1/((3^(1/2)*i)/2 - 1/2)^(1/3) - ((3^(1/2)*i)/2 - 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) - 1/2)^(1/3)) - ((3^(1/2)*i)/2 - 1/2)^(1/3)/2
(3^(1/2)*(1/((3^(1/2)*i)/2 - 1/2)^(1/3) - ((3^(1/2)*i)/2 - 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) - 1/2)^(1/3)) - ((3^(1/2)*i)/2 - 1/2)^(1/3)/2
- (3^(1/2)*(1/((3^(1/2)*i)/2 - 1/2)^(1/3) - ((3^(1/2)*i)/2 - 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) - 1/2)^(1/3)) - ((3^(1/2)*i)/2 - 1/2)^(1/3)/2
- (3^(1/2)*(1/((3^(1/2)*i)/2 - 1/2)^(1/3) - ((3^(1/2)*i)/2 - 1/2)^(1/3))*i)/2 - 1/(2*(3^(1/2)*(i/2) - 1/2)^(1/3)) - ((3^(1/2)*i)/2 - 1/2)^(1/3)/2