three order complex coefficient polynomial root matlab

Another classic solution is to find the matrix that has the same eigenvalues as your polynomial has roots. Then use eig to compute the eigenvalues of this "companion matrix". This is in fact what roots does.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 3

Walter Roberson 2023-5-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1974274-three-order-complex-coefficient-polynomial-root-matlab#answer_1245614

在 MATLAB Online 中打开

syms x A B C

eqn = 28*x^3 + A*x^2 + B*x - C == 0;

solutions = solve(eqn, x, 'MaxDegree', 3)

solutions =

char(solutions(1))

ans = '(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3) - (B/84 - A^2/7056)/(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3) - A/84'

char(solutions(2))

ans = '(B/84 - A^2/7056)/(2*(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3)) - A/84 - (C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3)/2 - (3^(1/2)*((B/84 - A^2/7056)/(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3) + (C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3))*1i)/2'

char(solutions(3))

ans = '(B/84 - A^2/7056)/(2*(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3)) - A/84 - (C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3)/2 + (3^(1/2)*((B/84 - A^2/7056)/(C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3) + (C/56 + ((C/56 - A^3/592704 + (A*B)/4704)^2 + (B/84 - A^2/7056)^3)^(1/2) - A^3/592704 + (A*B)/4704)^(1/3))*1i)/2'