Factorize reciprocal polynomial 4th-order










3














I try to factorize any polynomial like :



$x^4 + a.x^3 + b.x^2 + a.x + 1$ with $ a, b inBbbR$



into :



$(x^2 + c.x + d)(x^2 + e.x + f)$ with $ c, d, e, f inBbbR$



I also want $c(a, b)$, $d(a, b)$, $e(a, b)$, $f(a, b)$ to be continous, so I can be smooth at runtime when changing $(a, b)$.



The only way I succeeded for now is to compute all the roots and regroup them by conjugate pairs. But it's tricky because roots can be paired in many way, there are also cases where there are many solutions (when all roots are real), and I noticed roots may swap for specific values of $(a, b)$.



I would like to now if there is a simpler known method for this typical polynom.



Many thanks !










share|cite|improve this question























  • Please use MathJax to type the equations. This way they will be much more understandable.
    – Matti P.
    Nov 12 at 10:20










  • Thanks, edited.
    – user2443456
    Nov 12 at 10:35










  • Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
    – Matti P.
    Nov 12 at 10:42
















3














I try to factorize any polynomial like :



$x^4 + a.x^3 + b.x^2 + a.x + 1$ with $ a, b inBbbR$



into :



$(x^2 + c.x + d)(x^2 + e.x + f)$ with $ c, d, e, f inBbbR$



I also want $c(a, b)$, $d(a, b)$, $e(a, b)$, $f(a, b)$ to be continous, so I can be smooth at runtime when changing $(a, b)$.



The only way I succeeded for now is to compute all the roots and regroup them by conjugate pairs. But it's tricky because roots can be paired in many way, there are also cases where there are many solutions (when all roots are real), and I noticed roots may swap for specific values of $(a, b)$.



I would like to now if there is a simpler known method for this typical polynom.



Many thanks !










share|cite|improve this question























  • Please use MathJax to type the equations. This way they will be much more understandable.
    – Matti P.
    Nov 12 at 10:20










  • Thanks, edited.
    – user2443456
    Nov 12 at 10:35










  • Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
    – Matti P.
    Nov 12 at 10:42














3












3








3


1





I try to factorize any polynomial like :



$x^4 + a.x^3 + b.x^2 + a.x + 1$ with $ a, b inBbbR$



into :



$(x^2 + c.x + d)(x^2 + e.x + f)$ with $ c, d, e, f inBbbR$



I also want $c(a, b)$, $d(a, b)$, $e(a, b)$, $f(a, b)$ to be continous, so I can be smooth at runtime when changing $(a, b)$.



The only way I succeeded for now is to compute all the roots and regroup them by conjugate pairs. But it's tricky because roots can be paired in many way, there are also cases where there are many solutions (when all roots are real), and I noticed roots may swap for specific values of $(a, b)$.



I would like to now if there is a simpler known method for this typical polynom.



Many thanks !










share|cite|improve this question















I try to factorize any polynomial like :



$x^4 + a.x^3 + b.x^2 + a.x + 1$ with $ a, b inBbbR$



into :



$(x^2 + c.x + d)(x^2 + e.x + f)$ with $ c, d, e, f inBbbR$



I also want $c(a, b)$, $d(a, b)$, $e(a, b)$, $f(a, b)$ to be continous, so I can be smooth at runtime when changing $(a, b)$.



The only way I succeeded for now is to compute all the roots and regroup them by conjugate pairs. But it's tricky because roots can be paired in many way, there are also cases where there are many solutions (when all roots are real), and I noticed roots may swap for specific values of $(a, b)$.



I would like to now if there is a simpler known method for this typical polynom.



Many thanks !







polynomials






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Nov 12 at 10:32

























asked Nov 12 at 10:15









user2443456

424




424











  • Please use MathJax to type the equations. This way they will be much more understandable.
    – Matti P.
    Nov 12 at 10:20










  • Thanks, edited.
    – user2443456
    Nov 12 at 10:35










  • Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
    – Matti P.
    Nov 12 at 10:42

















  • Please use MathJax to type the equations. This way they will be much more understandable.
    – Matti P.
    Nov 12 at 10:20










  • Thanks, edited.
    – user2443456
    Nov 12 at 10:35










  • Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
    – Matti P.
    Nov 12 at 10:42
















Please use MathJax to type the equations. This way they will be much more understandable.
– Matti P.
Nov 12 at 10:20




Please use MathJax to type the equations. This way they will be much more understandable.
– Matti P.
Nov 12 at 10:20












Thanks, edited.
– user2443456
Nov 12 at 10:35




Thanks, edited.
– user2443456
Nov 12 at 10:35












Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
– Matti P.
Nov 12 at 10:42





Well, just opening the parentheses in the second equation, we get $$ beginarrayrcc e + c &=& a \ f+d+ce &=& b \ cf + de &=& a \ df &=& 1 endarray $$
– Matti P.
Nov 12 at 10:42











2 Answers
2






active

oldest

votes


















5














As you pointed out yourself, the polynomial is palindromic (=equal to its own reciprocal). This means that $1/alpha$ is a zero whenever $alpha$ is.



In the case of four real roots this suggests using factors like $(x+alpha)(x+1/alpha)$. Such a factor has constant term $=1$, so you will have $d=f=1$.
This reduces your system of equations to $c+e=a$, $2+ce=b$. Meaning that $c$ and $e$ are the roots of the quadratic
$$
(T-c)(T-e)=T^2-(c+e)T+ce=T^2-aT+(b-2).
$$

The roots of this quadratic are real whenever we have the inequality
$$
a^2>4(b-2),
$$

when
$$
c,e=fracapmsqrta^2-4b+82.
$$




If this equality is not satisfied then you can proceed as follows.



  • Solve for $c$ and $e$ anyway (quadratic formula).

  • Find the zeros of $x^2+cx+1$ and $x^2+ex+1$ (quadratic formula again).

  • Match the complex zeros into conjugate pairs (one from each of the above pairs of solutions).

  • Build the factors.


Of course, when you have four real zeros you have many options with the factorization. The one I suggested is just easy to find, and the coefficients $c,e(,d,f)$ are easy to write down as functions of $a$ and $b$.






share|cite|improve this answer


















  • 1




    Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
    – user2443456
    Nov 13 at 11:11



















3














Write $$x^4+a x^3+b x^2+a x+1-left(x^2+c x+dright) left(x^2+e x+fright)=0$$ Expand and group terms to get
$$(1-d f)+x (a-c f-d e)+x^2 (b-c e-d-f)+x^3 (a-c-e)=0$$ Set each coefficient equal to $0$ and you should get four nasty equations in $c,d,e,f$






share|cite|improve this answer




















    Your Answer





    StackExchange.ifUsing("editor", function ()
    return StackExchange.using("mathjaxEditing", function ()
    StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
    StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
    );
    );
    , "mathjax-editing");

    StackExchange.ready(function()
    var channelOptions =
    tags: "".split(" "),
    id: "69"
    ;
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function()
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled)
    StackExchange.using("snippets", function()
    createEditor();
    );

    else
    createEditor();

    );

    function createEditor()
    StackExchange.prepareEditor(
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: true,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: 10,
    bindNavPrevention: true,
    postfix: "",
    imageUploader:
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    ,
    noCode: true, onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    );



    );













    draft saved

    draft discarded


















    StackExchange.ready(
    function ()
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2995124%2ffactorize-reciprocal-polynomial-4th-order%23new-answer', 'question_page');

    );

    Post as a guest















    Required, but never shown

























    2 Answers
    2






    active

    oldest

    votes








    2 Answers
    2






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    5














    As you pointed out yourself, the polynomial is palindromic (=equal to its own reciprocal). This means that $1/alpha$ is a zero whenever $alpha$ is.



    In the case of four real roots this suggests using factors like $(x+alpha)(x+1/alpha)$. Such a factor has constant term $=1$, so you will have $d=f=1$.
    This reduces your system of equations to $c+e=a$, $2+ce=b$. Meaning that $c$ and $e$ are the roots of the quadratic
    $$
    (T-c)(T-e)=T^2-(c+e)T+ce=T^2-aT+(b-2).
    $$

    The roots of this quadratic are real whenever we have the inequality
    $$
    a^2>4(b-2),
    $$

    when
    $$
    c,e=fracapmsqrta^2-4b+82.
    $$




    If this equality is not satisfied then you can proceed as follows.



    • Solve for $c$ and $e$ anyway (quadratic formula).

    • Find the zeros of $x^2+cx+1$ and $x^2+ex+1$ (quadratic formula again).

    • Match the complex zeros into conjugate pairs (one from each of the above pairs of solutions).

    • Build the factors.


    Of course, when you have four real zeros you have many options with the factorization. The one I suggested is just easy to find, and the coefficients $c,e(,d,f)$ are easy to write down as functions of $a$ and $b$.






    share|cite|improve this answer


















    • 1




      Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
      – user2443456
      Nov 13 at 11:11
















    5














    As you pointed out yourself, the polynomial is palindromic (=equal to its own reciprocal). This means that $1/alpha$ is a zero whenever $alpha$ is.



    In the case of four real roots this suggests using factors like $(x+alpha)(x+1/alpha)$. Such a factor has constant term $=1$, so you will have $d=f=1$.
    This reduces your system of equations to $c+e=a$, $2+ce=b$. Meaning that $c$ and $e$ are the roots of the quadratic
    $$
    (T-c)(T-e)=T^2-(c+e)T+ce=T^2-aT+(b-2).
    $$

    The roots of this quadratic are real whenever we have the inequality
    $$
    a^2>4(b-2),
    $$

    when
    $$
    c,e=fracapmsqrta^2-4b+82.
    $$




    If this equality is not satisfied then you can proceed as follows.



    • Solve for $c$ and $e$ anyway (quadratic formula).

    • Find the zeros of $x^2+cx+1$ and $x^2+ex+1$ (quadratic formula again).

    • Match the complex zeros into conjugate pairs (one from each of the above pairs of solutions).

    • Build the factors.


    Of course, when you have four real zeros you have many options with the factorization. The one I suggested is just easy to find, and the coefficients $c,e(,d,f)$ are easy to write down as functions of $a$ and $b$.






    share|cite|improve this answer


















    • 1




      Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
      – user2443456
      Nov 13 at 11:11














    5












    5








    5






    As you pointed out yourself, the polynomial is palindromic (=equal to its own reciprocal). This means that $1/alpha$ is a zero whenever $alpha$ is.



    In the case of four real roots this suggests using factors like $(x+alpha)(x+1/alpha)$. Such a factor has constant term $=1$, so you will have $d=f=1$.
    This reduces your system of equations to $c+e=a$, $2+ce=b$. Meaning that $c$ and $e$ are the roots of the quadratic
    $$
    (T-c)(T-e)=T^2-(c+e)T+ce=T^2-aT+(b-2).
    $$

    The roots of this quadratic are real whenever we have the inequality
    $$
    a^2>4(b-2),
    $$

    when
    $$
    c,e=fracapmsqrta^2-4b+82.
    $$




    If this equality is not satisfied then you can proceed as follows.



    • Solve for $c$ and $e$ anyway (quadratic formula).

    • Find the zeros of $x^2+cx+1$ and $x^2+ex+1$ (quadratic formula again).

    • Match the complex zeros into conjugate pairs (one from each of the above pairs of solutions).

    • Build the factors.


    Of course, when you have four real zeros you have many options with the factorization. The one I suggested is just easy to find, and the coefficients $c,e(,d,f)$ are easy to write down as functions of $a$ and $b$.






    share|cite|improve this answer














    As you pointed out yourself, the polynomial is palindromic (=equal to its own reciprocal). This means that $1/alpha$ is a zero whenever $alpha$ is.



    In the case of four real roots this suggests using factors like $(x+alpha)(x+1/alpha)$. Such a factor has constant term $=1$, so you will have $d=f=1$.
    This reduces your system of equations to $c+e=a$, $2+ce=b$. Meaning that $c$ and $e$ are the roots of the quadratic
    $$
    (T-c)(T-e)=T^2-(c+e)T+ce=T^2-aT+(b-2).
    $$

    The roots of this quadratic are real whenever we have the inequality
    $$
    a^2>4(b-2),
    $$

    when
    $$
    c,e=fracapmsqrta^2-4b+82.
    $$




    If this equality is not satisfied then you can proceed as follows.



    • Solve for $c$ and $e$ anyway (quadratic formula).

    • Find the zeros of $x^2+cx+1$ and $x^2+ex+1$ (quadratic formula again).

    • Match the complex zeros into conjugate pairs (one from each of the above pairs of solutions).

    • Build the factors.


    Of course, when you have four real zeros you have many options with the factorization. The one I suggested is just easy to find, and the coefficients $c,e(,d,f)$ are easy to write down as functions of $a$ and $b$.







    share|cite|improve this answer














    share|cite|improve this answer



    share|cite|improve this answer








    edited Nov 15 at 5:16

























    answered Nov 12 at 11:03









    Jyrki Lahtonen

    108k12166367




    108k12166367







    • 1




      Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
      – user2443456
      Nov 13 at 11:11













    • 1




      Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
      – user2443456
      Nov 13 at 11:11








    1




    1




    Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
    – user2443456
    Nov 13 at 11:11





    Thanks, Also, for the case where the equality is not satisfied, you can just calculate one root and get the 3 others with conjugate/inversion of this one.
    – user2443456
    Nov 13 at 11:11












    3














    Write $$x^4+a x^3+b x^2+a x+1-left(x^2+c x+dright) left(x^2+e x+fright)=0$$ Expand and group terms to get
    $$(1-d f)+x (a-c f-d e)+x^2 (b-c e-d-f)+x^3 (a-c-e)=0$$ Set each coefficient equal to $0$ and you should get four nasty equations in $c,d,e,f$






    share|cite|improve this answer

























      3














      Write $$x^4+a x^3+b x^2+a x+1-left(x^2+c x+dright) left(x^2+e x+fright)=0$$ Expand and group terms to get
      $$(1-d f)+x (a-c f-d e)+x^2 (b-c e-d-f)+x^3 (a-c-e)=0$$ Set each coefficient equal to $0$ and you should get four nasty equations in $c,d,e,f$






      share|cite|improve this answer























        3












        3








        3






        Write $$x^4+a x^3+b x^2+a x+1-left(x^2+c x+dright) left(x^2+e x+fright)=0$$ Expand and group terms to get
        $$(1-d f)+x (a-c f-d e)+x^2 (b-c e-d-f)+x^3 (a-c-e)=0$$ Set each coefficient equal to $0$ and you should get four nasty equations in $c,d,e,f$






        share|cite|improve this answer












        Write $$x^4+a x^3+b x^2+a x+1-left(x^2+c x+dright) left(x^2+e x+fright)=0$$ Expand and group terms to get
        $$(1-d f)+x (a-c f-d e)+x^2 (b-c e-d-f)+x^3 (a-c-e)=0$$ Set each coefficient equal to $0$ and you should get four nasty equations in $c,d,e,f$







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered Nov 12 at 10:44









        Claude Leibovici

        118k1157132




        118k1157132



























            draft saved

            draft discarded
















































            Thanks for contributing an answer to Mathematics Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.





            Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


            Please pay close attention to the following guidance:


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2995124%2ffactorize-reciprocal-polynomial-4th-order%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            這個網誌中的熱門文章

            How to read a connectionString WITH PROVIDER in .NET Core?

            In R, how to develop a multiplot heatmap.2 figure showing key labels successfully

            Museum of Modern and Contemporary Art of Trento and Rovereto