Recurrence relation with a matrix?

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
0
down vote

favorite













Let $ mathbfM = beginbmatrix 2 & 1 \ 1 & 2 endbmatrix$ and $mathbfX_n = beginbmatrixx_n \ y_n endbmatrix$. Solve $displaystyle mathbfX_n+1 = mathbfMX_n$ for $n ge 0$ with the i.c. $x_0 =2, y_0 = 0$.




I'm not sure how to solve this. My idea is to write this as $mathbfX_n+k =mathbfM^k mathbfX_n$ for $k ge 0$. But it's not clear how I'd get the solution from that, after I've gone through the calculation?







share|cite|improve this question





















  • You need to diagonalize $mathbf M$.
    – joriki
    Jul 29 at 8:48














up vote
0
down vote

favorite













Let $ mathbfM = beginbmatrix 2 & 1 \ 1 & 2 endbmatrix$ and $mathbfX_n = beginbmatrixx_n \ y_n endbmatrix$. Solve $displaystyle mathbfX_n+1 = mathbfMX_n$ for $n ge 0$ with the i.c. $x_0 =2, y_0 = 0$.




I'm not sure how to solve this. My idea is to write this as $mathbfX_n+k =mathbfM^k mathbfX_n$ for $k ge 0$. But it's not clear how I'd get the solution from that, after I've gone through the calculation?







share|cite|improve this question





















  • You need to diagonalize $mathbf M$.
    – joriki
    Jul 29 at 8:48












up vote
0
down vote

favorite









up vote
0
down vote

favorite












Let $ mathbfM = beginbmatrix 2 & 1 \ 1 & 2 endbmatrix$ and $mathbfX_n = beginbmatrixx_n \ y_n endbmatrix$. Solve $displaystyle mathbfX_n+1 = mathbfMX_n$ for $n ge 0$ with the i.c. $x_0 =2, y_0 = 0$.




I'm not sure how to solve this. My idea is to write this as $mathbfX_n+k =mathbfM^k mathbfX_n$ for $k ge 0$. But it's not clear how I'd get the solution from that, after I've gone through the calculation?







share|cite|improve this question














Let $ mathbfM = beginbmatrix 2 & 1 \ 1 & 2 endbmatrix$ and $mathbfX_n = beginbmatrixx_n \ y_n endbmatrix$. Solve $displaystyle mathbfX_n+1 = mathbfMX_n$ for $n ge 0$ with the i.c. $x_0 =2, y_0 = 0$.




I'm not sure how to solve this. My idea is to write this as $mathbfX_n+k =mathbfM^k mathbfX_n$ for $k ge 0$. But it's not clear how I'd get the solution from that, after I've gone through the calculation?









share|cite|improve this question












share|cite|improve this question




share|cite|improve this question








edited Jul 29 at 9:04
























asked Jul 29 at 8:46









Nebulae

505




505











  • You need to diagonalize $mathbf M$.
    – joriki
    Jul 29 at 8:48
















  • You need to diagonalize $mathbf M$.
    – joriki
    Jul 29 at 8:48















You need to diagonalize $mathbf M$.
– joriki
Jul 29 at 8:48




You need to diagonalize $mathbf M$.
– joriki
Jul 29 at 8:48










3 Answers
3






active

oldest

votes

















up vote
2
down vote



accepted










If you diagonalize the matrix,



$$M=PDP^-1$$ and



$$M^n=PDP^-1PDP^-1cdots PDP^-1=PD^nP^-1$$



so that the $n^th $ power of a diagonal matrix is the diagonal made of the $n^th$ powers.



Depending on the values of the Eigenvalues, the diagonal elements converge to zero or diverge to infinity (and $1$ remains and $-1$ alternates).



But in this exercise, all of this is useless as the initial conditions are null.






share|cite|improve this answer





















  • I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
    – Nebulae
    Jul 29 at 9:05







  • 1




    @Farthing: no, $X_n=M^nX_0$.
    – Yves Daoust
    Jul 29 at 9:10











  • In general does $PDP^-1 = P^-1DP$?
    – Nebulae
    Jul 29 at 9:19










  • @Farthing: consider $Q=P^-1$.
    – Yves Daoust
    Jul 29 at 9:21

















up vote
2
down vote













Hint: Take a good look at that initial condition. What is $mathbf X_1$?






share|cite|improve this answer





















  • I'm sorry I typed the initial condition wrong.
    – Nebulae
    Jul 29 at 9:04











  • @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
    – Arthur
    Jul 29 at 9:10

















up vote
0
down vote













Alt. hint (without matrix powers):   the recurrence relations can be written as:



$$
beginalign
begincases
x_n+1=2 x_n + y_n \
y_n+1=x_n + 2y_n
endcases
endalign
$$



Subtracting the two equalities gives $,x_n+1-y_n+1=x_n-y_n,$, which telescopes to:



$$x_n-y_n=x_n-1-y_n-1=ldots=x_0-y_0=C$$



Thus $,y_n=x_n-C,$, and substituting back into the first relation gives the simple recurrence in $,x_n,$:



$$
x_n+1=2x_n+(x_n-C) quadiffquad x_n+1=3x_n - C
$$






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',
    convertImagesToLinks: true,
    noModals: false,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: 10,
    bindNavPrevention: true,
    postfix: "",
    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%2f2865902%2frecurrence-relation-with-a-matrix%23new-answer', 'question_page');

    );

    Post as a guest






























    3 Answers
    3






    active

    oldest

    votes








    3 Answers
    3






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes








    up vote
    2
    down vote



    accepted










    If you diagonalize the matrix,



    $$M=PDP^-1$$ and



    $$M^n=PDP^-1PDP^-1cdots PDP^-1=PD^nP^-1$$



    so that the $n^th $ power of a diagonal matrix is the diagonal made of the $n^th$ powers.



    Depending on the values of the Eigenvalues, the diagonal elements converge to zero or diverge to infinity (and $1$ remains and $-1$ alternates).



    But in this exercise, all of this is useless as the initial conditions are null.






    share|cite|improve this answer





















    • I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
      – Nebulae
      Jul 29 at 9:05







    • 1




      @Farthing: no, $X_n=M^nX_0$.
      – Yves Daoust
      Jul 29 at 9:10











    • In general does $PDP^-1 = P^-1DP$?
      – Nebulae
      Jul 29 at 9:19










    • @Farthing: consider $Q=P^-1$.
      – Yves Daoust
      Jul 29 at 9:21














    up vote
    2
    down vote



    accepted










    If you diagonalize the matrix,



    $$M=PDP^-1$$ and



    $$M^n=PDP^-1PDP^-1cdots PDP^-1=PD^nP^-1$$



    so that the $n^th $ power of a diagonal matrix is the diagonal made of the $n^th$ powers.



    Depending on the values of the Eigenvalues, the diagonal elements converge to zero or diverge to infinity (and $1$ remains and $-1$ alternates).



    But in this exercise, all of this is useless as the initial conditions are null.






    share|cite|improve this answer





















    • I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
      – Nebulae
      Jul 29 at 9:05







    • 1




      @Farthing: no, $X_n=M^nX_0$.
      – Yves Daoust
      Jul 29 at 9:10











    • In general does $PDP^-1 = P^-1DP$?
      – Nebulae
      Jul 29 at 9:19










    • @Farthing: consider $Q=P^-1$.
      – Yves Daoust
      Jul 29 at 9:21












    up vote
    2
    down vote



    accepted







    up vote
    2
    down vote



    accepted






    If you diagonalize the matrix,



    $$M=PDP^-1$$ and



    $$M^n=PDP^-1PDP^-1cdots PDP^-1=PD^nP^-1$$



    so that the $n^th $ power of a diagonal matrix is the diagonal made of the $n^th$ powers.



    Depending on the values of the Eigenvalues, the diagonal elements converge to zero or diverge to infinity (and $1$ remains and $-1$ alternates).



    But in this exercise, all of this is useless as the initial conditions are null.






    share|cite|improve this answer













    If you diagonalize the matrix,



    $$M=PDP^-1$$ and



    $$M^n=PDP^-1PDP^-1cdots PDP^-1=PD^nP^-1$$



    so that the $n^th $ power of a diagonal matrix is the diagonal made of the $n^th$ powers.



    Depending on the values of the Eigenvalues, the diagonal elements converge to zero or diverge to infinity (and $1$ remains and $-1$ alternates).



    But in this exercise, all of this is useless as the initial conditions are null.







    share|cite|improve this answer













    share|cite|improve this answer



    share|cite|improve this answer











    answered Jul 29 at 9:02









    Yves Daoust

    110k665203




    110k665203











    • I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
      – Nebulae
      Jul 29 at 9:05







    • 1




      @Farthing: no, $X_n=M^nX_0$.
      – Yves Daoust
      Jul 29 at 9:10











    • In general does $PDP^-1 = P^-1DP$?
      – Nebulae
      Jul 29 at 9:19










    • @Farthing: consider $Q=P^-1$.
      – Yves Daoust
      Jul 29 at 9:21
















    • I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
      – Nebulae
      Jul 29 at 9:05







    • 1




      @Farthing: no, $X_n=M^nX_0$.
      – Yves Daoust
      Jul 29 at 9:10











    • In general does $PDP^-1 = P^-1DP$?
      – Nebulae
      Jul 29 at 9:19










    • @Farthing: consider $Q=P^-1$.
      – Yves Daoust
      Jul 29 at 9:21















    I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
    – Nebulae
    Jul 29 at 9:05





    I got the initial condition wrong. Also I thought it was $M = P^-1DP$? Thanks. So after I find $M^kX_n $, I let $k=0$?
    – Nebulae
    Jul 29 at 9:05





    1




    1




    @Farthing: no, $X_n=M^nX_0$.
    – Yves Daoust
    Jul 29 at 9:10





    @Farthing: no, $X_n=M^nX_0$.
    – Yves Daoust
    Jul 29 at 9:10













    In general does $PDP^-1 = P^-1DP$?
    – Nebulae
    Jul 29 at 9:19




    In general does $PDP^-1 = P^-1DP$?
    – Nebulae
    Jul 29 at 9:19












    @Farthing: consider $Q=P^-1$.
    – Yves Daoust
    Jul 29 at 9:21




    @Farthing: consider $Q=P^-1$.
    – Yves Daoust
    Jul 29 at 9:21










    up vote
    2
    down vote













    Hint: Take a good look at that initial condition. What is $mathbf X_1$?






    share|cite|improve this answer





















    • I'm sorry I typed the initial condition wrong.
      – Nebulae
      Jul 29 at 9:04











    • @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
      – Arthur
      Jul 29 at 9:10














    up vote
    2
    down vote













    Hint: Take a good look at that initial condition. What is $mathbf X_1$?






    share|cite|improve this answer





















    • I'm sorry I typed the initial condition wrong.
      – Nebulae
      Jul 29 at 9:04











    • @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
      – Arthur
      Jul 29 at 9:10












    up vote
    2
    down vote










    up vote
    2
    down vote









    Hint: Take a good look at that initial condition. What is $mathbf X_1$?






    share|cite|improve this answer













    Hint: Take a good look at that initial condition. What is $mathbf X_1$?







    share|cite|improve this answer













    share|cite|improve this answer



    share|cite|improve this answer











    answered Jul 29 at 8:48









    Arthur

    98.4k793174




    98.4k793174











    • I'm sorry I typed the initial condition wrong.
      – Nebulae
      Jul 29 at 9:04











    • @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
      – Arthur
      Jul 29 at 9:10
















    • I'm sorry I typed the initial condition wrong.
      – Nebulae
      Jul 29 at 9:04











    • @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
      – Arthur
      Jul 29 at 9:10















    I'm sorry I typed the initial condition wrong.
    – Nebulae
    Jul 29 at 9:04





    I'm sorry I typed the initial condition wrong.
    – Nebulae
    Jul 29 at 9:04













    @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
    – Arthur
    Jul 29 at 9:10




    @Farthing Well, then. Diagonalisation it is. Or you could try to first solve the initial value problem $M_0=M, M_n=MM_n-1$. Whichever you think is easiest.
    – Arthur
    Jul 29 at 9:10










    up vote
    0
    down vote













    Alt. hint (without matrix powers):   the recurrence relations can be written as:



    $$
    beginalign
    begincases
    x_n+1=2 x_n + y_n \
    y_n+1=x_n + 2y_n
    endcases
    endalign
    $$



    Subtracting the two equalities gives $,x_n+1-y_n+1=x_n-y_n,$, which telescopes to:



    $$x_n-y_n=x_n-1-y_n-1=ldots=x_0-y_0=C$$



    Thus $,y_n=x_n-C,$, and substituting back into the first relation gives the simple recurrence in $,x_n,$:



    $$
    x_n+1=2x_n+(x_n-C) quadiffquad x_n+1=3x_n - C
    $$






    share|cite|improve this answer

























      up vote
      0
      down vote













      Alt. hint (without matrix powers):   the recurrence relations can be written as:



      $$
      beginalign
      begincases
      x_n+1=2 x_n + y_n \
      y_n+1=x_n + 2y_n
      endcases
      endalign
      $$



      Subtracting the two equalities gives $,x_n+1-y_n+1=x_n-y_n,$, which telescopes to:



      $$x_n-y_n=x_n-1-y_n-1=ldots=x_0-y_0=C$$



      Thus $,y_n=x_n-C,$, and substituting back into the first relation gives the simple recurrence in $,x_n,$:



      $$
      x_n+1=2x_n+(x_n-C) quadiffquad x_n+1=3x_n - C
      $$






      share|cite|improve this answer























        up vote
        0
        down vote










        up vote
        0
        down vote









        Alt. hint (without matrix powers):   the recurrence relations can be written as:



        $$
        beginalign
        begincases
        x_n+1=2 x_n + y_n \
        y_n+1=x_n + 2y_n
        endcases
        endalign
        $$



        Subtracting the two equalities gives $,x_n+1-y_n+1=x_n-y_n,$, which telescopes to:



        $$x_n-y_n=x_n-1-y_n-1=ldots=x_0-y_0=C$$



        Thus $,y_n=x_n-C,$, and substituting back into the first relation gives the simple recurrence in $,x_n,$:



        $$
        x_n+1=2x_n+(x_n-C) quadiffquad x_n+1=3x_n - C
        $$






        share|cite|improve this answer













        Alt. hint (without matrix powers):   the recurrence relations can be written as:



        $$
        beginalign
        begincases
        x_n+1=2 x_n + y_n \
        y_n+1=x_n + 2y_n
        endcases
        endalign
        $$



        Subtracting the two equalities gives $,x_n+1-y_n+1=x_n-y_n,$, which telescopes to:



        $$x_n-y_n=x_n-1-y_n-1=ldots=x_0-y_0=C$$



        Thus $,y_n=x_n-C,$, and substituting back into the first relation gives the simple recurrence in $,x_n,$:



        $$
        x_n+1=2x_n+(x_n-C) quadiffquad x_n+1=3x_n - C
        $$







        share|cite|improve this answer













        share|cite|improve this answer



        share|cite|improve this answer











        answered Jul 29 at 18:11









        dxiv

        53.8k64796




        53.8k64796






















             

            draft saved


            draft discarded


























             


            draft saved


            draft discarded














            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2865902%2frecurrence-relation-with-a-matrix%23new-answer', 'question_page');

            );

            Post as a guest













































































            Comments

            Popular posts from this blog

            What is the equation of a 3D cone with generalised tilt?

            Color the edges and diagonals of a regular polygon

            Relationship between determinant of matrix and determinant of adjoint?