필터 지우기
필터 지우기

'Solve' 10th order polynomial symbolically; answer comes in the form of 'z' although defining MaxDegree!

조회 수: 5 (최근 30일)
Hello,
I am trying to solve a 10th order polynomial symbolically using 'solve' command. I am using Matlab R2019a academic version. I am using the following the code:
syms zita lamda omega r Fe positive
DT=((2*zita*omega).^2)*r^2+((1-(1+lamda*r^2/2+lamda*r^4/4).*omega.^2).^2)*r^2-(Fe^2)==0;
E=expand(DT);
R=solve(E,r,'MaxDegree',10)
What I'm getting for answer is root (..) in terms of 'z'. Previously I've ran into such problem but when I provided the MaxDegree, Matlab was able to provide the solution. In those cases, maximum degree was 6. I was hoping if anyone could suggest a different approach to solving this problem. I'd really appreciate it.
Anika

채택된 답변

Walter Roberson
Walter Roberson 2020년 3월 10일
You have a quintic (degree 5). There are no general closed form solutions for degree 5 or higher.
The 'MaxDegree' option is only useful if the polynomial can be factored into expressions that are degree 4 or less. MaxDegree is not magic: it just tells MATLAB how aggressive to be in substituting in long closed form formulas when they are available . And degree 4 is the highest degree that you can be sure that the long formulas are available.
  댓글 수: 2
Anika Tabassum Sarkar
Anika Tabassum Sarkar 2020년 3월 10일
Thank you Walter! This explains it. I was wondering if you know of any other tool that is capable of solving higher degree closed form equation.
Walter Roberson
Walter Roberson 2020년 3월 10일
It has been proven theoretically that there are quintics (degree 5) that have no "algebraic root" . So for degree 5 or higher, finding a symbolic root is a matter of managing to factor the expression through various techniques such as "completing the square" in order to reduce it to the product of terms that are of lower degree. The standard symbolic packages are typically not bad at that -- though not fool-proof.
As best I can tell, the equation you show us does not have closed form roots in r.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Assumptions에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by