Finding the 10th root of a matrix












13












$begingroup$


I want to find a $2 times 2$ matrix, named $A$ in this situation, such that:



$$A^{10}=begin {bmatrix} 1 & 1 \ 0 & 1 end {bmatrix} $$



How can I get started? I was thinking about filling $A$ with arbitrary values $a, b, c, d$ and then multiplying it by itself ten times, then setting those values equal to the given values but I quickly realized that would take too long. Is there a more efficient way?










share|cite|improve this question











$endgroup$












  • $begingroup$
    Generalized, here as well.
    $endgroup$
    – Jyrki Lahtonen
    Jan 29 '18 at 11:23












  • $begingroup$
    Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
    $endgroup$
    – mathreadler
    Jan 29 '18 at 20:20


















13












$begingroup$


I want to find a $2 times 2$ matrix, named $A$ in this situation, such that:



$$A^{10}=begin {bmatrix} 1 & 1 \ 0 & 1 end {bmatrix} $$



How can I get started? I was thinking about filling $A$ with arbitrary values $a, b, c, d$ and then multiplying it by itself ten times, then setting those values equal to the given values but I quickly realized that would take too long. Is there a more efficient way?










share|cite|improve this question











$endgroup$












  • $begingroup$
    Generalized, here as well.
    $endgroup$
    – Jyrki Lahtonen
    Jan 29 '18 at 11:23












  • $begingroup$
    Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
    $endgroup$
    – mathreadler
    Jan 29 '18 at 20:20
















13












13








13


3



$begingroup$


I want to find a $2 times 2$ matrix, named $A$ in this situation, such that:



$$A^{10}=begin {bmatrix} 1 & 1 \ 0 & 1 end {bmatrix} $$



How can I get started? I was thinking about filling $A$ with arbitrary values $a, b, c, d$ and then multiplying it by itself ten times, then setting those values equal to the given values but I quickly realized that would take too long. Is there a more efficient way?










share|cite|improve this question











$endgroup$




I want to find a $2 times 2$ matrix, named $A$ in this situation, such that:



$$A^{10}=begin {bmatrix} 1 & 1 \ 0 & 1 end {bmatrix} $$



How can I get started? I was thinking about filling $A$ with arbitrary values $a, b, c, d$ and then multiplying it by itself ten times, then setting those values equal to the given values but I quickly realized that would take too long. Is there a more efficient way?







linear-algebra matrices exponentiation matrix-equations






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Jan 29 '18 at 18:22









A Simmons

354111




354111










asked Jan 29 '18 at 8:33









fxcdfxcd

9414




9414












  • $begingroup$
    Generalized, here as well.
    $endgroup$
    – Jyrki Lahtonen
    Jan 29 '18 at 11:23












  • $begingroup$
    Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
    $endgroup$
    – mathreadler
    Jan 29 '18 at 20:20




















  • $begingroup$
    Generalized, here as well.
    $endgroup$
    – Jyrki Lahtonen
    Jan 29 '18 at 11:23












  • $begingroup$
    Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
    $endgroup$
    – mathreadler
    Jan 29 '18 at 20:20


















$begingroup$
Generalized, here as well.
$endgroup$
– Jyrki Lahtonen
Jan 29 '18 at 11:23






$begingroup$
Generalized, here as well.
$endgroup$
– Jyrki Lahtonen
Jan 29 '18 at 11:23














$begingroup$
Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
$endgroup$
– mathreadler
Jan 29 '18 at 20:20






$begingroup$
Here is another more general approach: math.stackexchange.com/questions/2576040/… which is proven to work here : math.stackexchange.com/questions/2601508/…
$endgroup$
– mathreadler
Jan 29 '18 at 20:20












6 Answers
6






active

oldest

votes


















34












$begingroup$

Take$$A=begin{bmatrix}1&x\0&1end{bmatrix}.$$Now, compute $A^2,A^3,ldots$ You'll find quickly which $x$ you should choose.






share|cite|improve this answer











$endgroup$









  • 2




    $begingroup$
    Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
    $endgroup$
    – Jeppe Stig Nielsen
    Jan 29 '18 at 18:08






  • 2




    $begingroup$
    @JeppeStigNielsen Indeed, that's a more difficult problem.
    $endgroup$
    – José Carlos Santos
    Jan 29 '18 at 19:10










  • $begingroup$
    (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
    $endgroup$
    – robjohn
    Jan 30 '18 at 18:10



















22












$begingroup$

You can also consider $$A=begin {bmatrix} 1 & 0 \ 0 & 1 end {bmatrix} + begin {bmatrix} 0 & a \ 0 & 0 end {bmatrix} =I+N $$ and try to use binomial formula for expansion of powers of a binomial.



Notice that $N^k=0$ for $k>1$.






share|cite|improve this answer









$endgroup$









  • 22




    $begingroup$
    Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
    $endgroup$
    – M.Herzkamp
    Jan 29 '18 at 11:45






  • 4




    $begingroup$
    @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
    $endgroup$
    – Widawensen
    Jan 29 '18 at 11:49





















10












$begingroup$

The answers given by José Carlos Santos and Widawensen are correct. There is a more general alternative method, even through it does not work for your particular example. But I thought I'd give it anyway in case it is useful to others.



If you want the $n^{text{th}}$ root of the matrix $M$, and $M$ is diagonalizable - i.e., there exists $Q$ such that $D = Q^{-1}MQ$ is a diagonal matrix, then $D^{1/n}$ is just the matrix whose diagonal elements are those of $D$ raised to the $1/n$ power. So you can set $$A = QD^{1/n}Q^{-1}$$ and find that $$A^n = left(QD^{1/n}Q^{-1}right)^n = Qleft(D^{1/n}right)^nQ^{-1} = QDQ^{-1} = M$$



To diagonalize $M$, you find it's eigenvalues and their corresponding eigenvectors. If the eigenvectors span the entire space, then $M$ is diagonalizable, the columns of $Q$ are eigenvectors of $M$, and the diagonal elements of $D$ are the eigenvalues.



Alas, it is not applicable here: $$begin{bmatrix}1 & 1\0&1end{bmatrix}$$
is not diagonalizable.






share|cite|improve this answer









$endgroup$













  • $begingroup$
    If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
    $endgroup$
    – Kyle Miller
    Jan 30 '18 at 5:52



















5












$begingroup$

Ok this may seem a bit overkill for this particular question, but please bear with me as it could help others.



You can use the Newton-Rhapson method, trying to solve the equation



$$f(x) = x^{10}-d=0$$
where
$$f'(x) = 10x^9$$
and $d$ is that matrix of yours and the iteration:



$${bf X_{n+1}} ={bf X_n}- f({bf X_n})f'({bf X_n})^{-1}$$
As long as you set initial $bf X_0$ to nothing too crazy, it will work. In fact, this is the simplest example of an discrete fractional integration solution, if you look at larger matrices $d$ filled with 1 on and below the diagonal and 0 otherwise (ok, transpose of the same matrix but same idea).



The fractional integral operators in fractional calculus found if guessing at the $bf X_0 = bf I$ matrix (described as linear convolutional filters, corresponding to a row in the solution matrix):



enter image description here



The two we recognize is the one in the middle (constant, corresponding to normal integral) and the one at the end (linear, corresponding to double integration).






share|cite|improve this answer











$endgroup$













  • $begingroup$
    Is that supposed to be $f'(x)=10x^9$?
    $endgroup$
    – ASKASK
    Jan 30 '18 at 11:52










  • $begingroup$
    @ASKASK yep must have gotten an elf in my computer who ate it away ;)
    $endgroup$
    – mathreadler
    Jan 30 '18 at 11:53





















3












$begingroup$

Such a matrix $A$ will commute with $pmatrix{1&1\0&1}$, and so be of the form $pmatrix{a&b\0&a}$. This means that $a^{10}=1$ and so
$A=apmatrix{1&c\0&1}$ for some $c$. We immediately get $c=1/10$
(as long as $10$ is invertible in your field).






share|cite|improve this answer









$endgroup$





















    2












    $begingroup$

    Hint: Another approach is to note that
    $$
    expleft(begin{bmatrix}0&x\0&0end{bmatrix}right)=begin{bmatrix}1&x\0&1end{bmatrix}
    $$






    share|cite|improve this answer









    $endgroup$













      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%2f2626206%2ffinding-the-10th-root-of-a-matrix%23new-answer', 'question_page');
      }
      );

      Post as a guest















      Required, but never shown

























      6 Answers
      6






      active

      oldest

      votes








      6 Answers
      6






      active

      oldest

      votes









      active

      oldest

      votes






      active

      oldest

      votes









      34












      $begingroup$

      Take$$A=begin{bmatrix}1&x\0&1end{bmatrix}.$$Now, compute $A^2,A^3,ldots$ You'll find quickly which $x$ you should choose.






      share|cite|improve this answer











      $endgroup$









      • 2




        $begingroup$
        Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
        $endgroup$
        – Jeppe Stig Nielsen
        Jan 29 '18 at 18:08






      • 2




        $begingroup$
        @JeppeStigNielsen Indeed, that's a more difficult problem.
        $endgroup$
        – José Carlos Santos
        Jan 29 '18 at 19:10










      • $begingroup$
        (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
        $endgroup$
        – robjohn
        Jan 30 '18 at 18:10
















      34












      $begingroup$

      Take$$A=begin{bmatrix}1&x\0&1end{bmatrix}.$$Now, compute $A^2,A^3,ldots$ You'll find quickly which $x$ you should choose.






      share|cite|improve this answer











      $endgroup$









      • 2




        $begingroup$
        Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
        $endgroup$
        – Jeppe Stig Nielsen
        Jan 29 '18 at 18:08






      • 2




        $begingroup$
        @JeppeStigNielsen Indeed, that's a more difficult problem.
        $endgroup$
        – José Carlos Santos
        Jan 29 '18 at 19:10










      • $begingroup$
        (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
        $endgroup$
        – robjohn
        Jan 30 '18 at 18:10














      34












      34








      34





      $begingroup$

      Take$$A=begin{bmatrix}1&x\0&1end{bmatrix}.$$Now, compute $A^2,A^3,ldots$ You'll find quickly which $x$ you should choose.






      share|cite|improve this answer











      $endgroup$



      Take$$A=begin{bmatrix}1&x\0&1end{bmatrix}.$$Now, compute $A^2,A^3,ldots$ You'll find quickly which $x$ you should choose.







      share|cite|improve this answer














      share|cite|improve this answer



      share|cite|improve this answer








      edited Dec 26 '18 at 23:29

























      answered Jan 29 '18 at 8:35









      José Carlos SantosJosé Carlos Santos

      167k22132235




      167k22132235








      • 2




        $begingroup$
        Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
        $endgroup$
        – Jeppe Stig Nielsen
        Jan 29 '18 at 18:08






      • 2




        $begingroup$
        @JeppeStigNielsen Indeed, that's a more difficult problem.
        $endgroup$
        – José Carlos Santos
        Jan 29 '18 at 19:10










      • $begingroup$
        (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
        $endgroup$
        – robjohn
        Jan 30 '18 at 18:10














      • 2




        $begingroup$
        Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
        $endgroup$
        – Jeppe Stig Nielsen
        Jan 29 '18 at 18:08






      • 2




        $begingroup$
        @JeppeStigNielsen Indeed, that's a more difficult problem.
        $endgroup$
        – José Carlos Santos
        Jan 29 '18 at 19:10










      • $begingroup$
        (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
        $endgroup$
        – robjohn
        Jan 30 '18 at 18:10








      2




      2




      $begingroup$
      Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
      $endgroup$
      – Jeppe Stig Nielsen
      Jan 29 '18 at 18:08




      $begingroup$
      Yes. This idea can be conceived by noting that the given matrix is upper triangular, then maybe an upper triangular solution $A$ is possible, and the eigenvalues of the original matrix are $1$ and $1$, then maybe the eigenvalues of $A$ can be taken to be also $1$ and $1$ (a particular tenth root of $1$ is $1$). So if you try it out, you will succeed. A more difficult problem is to find all solutions $A$ with entries in some field (such as $mathbb{Q}$ or $mathbb{C}$).
      $endgroup$
      – Jeppe Stig Nielsen
      Jan 29 '18 at 18:08




      2




      2




      $begingroup$
      @JeppeStigNielsen Indeed, that's a more difficult problem.
      $endgroup$
      – José Carlos Santos
      Jan 29 '18 at 19:10




      $begingroup$
      @JeppeStigNielsen Indeed, that's a more difficult problem.
      $endgroup$
      – José Carlos Santos
      Jan 29 '18 at 19:10












      $begingroup$
      (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
      $endgroup$
      – robjohn
      Jan 30 '18 at 18:10




      $begingroup$
      (+1) This is the approach that first came to mind. Now I need to approach this slightly differently.
      $endgroup$
      – robjohn
      Jan 30 '18 at 18:10











      22












      $begingroup$

      You can also consider $$A=begin {bmatrix} 1 & 0 \ 0 & 1 end {bmatrix} + begin {bmatrix} 0 & a \ 0 & 0 end {bmatrix} =I+N $$ and try to use binomial formula for expansion of powers of a binomial.



      Notice that $N^k=0$ for $k>1$.






      share|cite|improve this answer









      $endgroup$









      • 22




        $begingroup$
        Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
        $endgroup$
        – M.Herzkamp
        Jan 29 '18 at 11:45






      • 4




        $begingroup$
        @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
        $endgroup$
        – Widawensen
        Jan 29 '18 at 11:49


















      22












      $begingroup$

      You can also consider $$A=begin {bmatrix} 1 & 0 \ 0 & 1 end {bmatrix} + begin {bmatrix} 0 & a \ 0 & 0 end {bmatrix} =I+N $$ and try to use binomial formula for expansion of powers of a binomial.



      Notice that $N^k=0$ for $k>1$.






      share|cite|improve this answer









      $endgroup$









      • 22




        $begingroup$
        Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
        $endgroup$
        – M.Herzkamp
        Jan 29 '18 at 11:45






      • 4




        $begingroup$
        @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
        $endgroup$
        – Widawensen
        Jan 29 '18 at 11:49
















      22












      22








      22





      $begingroup$

      You can also consider $$A=begin {bmatrix} 1 & 0 \ 0 & 1 end {bmatrix} + begin {bmatrix} 0 & a \ 0 & 0 end {bmatrix} =I+N $$ and try to use binomial formula for expansion of powers of a binomial.



      Notice that $N^k=0$ for $k>1$.






      share|cite|improve this answer









      $endgroup$



      You can also consider $$A=begin {bmatrix} 1 & 0 \ 0 & 1 end {bmatrix} + begin {bmatrix} 0 & a \ 0 & 0 end {bmatrix} =I+N $$ and try to use binomial formula for expansion of powers of a binomial.



      Notice that $N^k=0$ for $k>1$.







      share|cite|improve this answer












      share|cite|improve this answer



      share|cite|improve this answer










      answered Jan 29 '18 at 11:13









      WidawensenWidawensen

      4,60821446




      4,60821446








      • 22




        $begingroup$
        Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
        $endgroup$
        – M.Herzkamp
        Jan 29 '18 at 11:45






      • 4




        $begingroup$
        @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
        $endgroup$
        – Widawensen
        Jan 29 '18 at 11:49
















      • 22




        $begingroup$
        Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
        $endgroup$
        – M.Herzkamp
        Jan 29 '18 at 11:45






      • 4




        $begingroup$
        @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
        $endgroup$
        – Widawensen
        Jan 29 '18 at 11:49










      22




      22




      $begingroup$
      Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
      $endgroup$
      – M.Herzkamp
      Jan 29 '18 at 11:45




      $begingroup$
      Note that it is usually not possible to use the binomial formula for matrices, because the product is not commutative. In this case it is fine, because one of the matrices is the identity.
      $endgroup$
      – M.Herzkamp
      Jan 29 '18 at 11:45




      4




      4




      $begingroup$
      @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
      $endgroup$
      – Widawensen
      Jan 29 '18 at 11:49






      $begingroup$
      @M.Herzkamp Yes, we can use the formula because matrices commute. Thank you for valuable remark..
      $endgroup$
      – Widawensen
      Jan 29 '18 at 11:49













      10












      $begingroup$

      The answers given by José Carlos Santos and Widawensen are correct. There is a more general alternative method, even through it does not work for your particular example. But I thought I'd give it anyway in case it is useful to others.



      If you want the $n^{text{th}}$ root of the matrix $M$, and $M$ is diagonalizable - i.e., there exists $Q$ such that $D = Q^{-1}MQ$ is a diagonal matrix, then $D^{1/n}$ is just the matrix whose diagonal elements are those of $D$ raised to the $1/n$ power. So you can set $$A = QD^{1/n}Q^{-1}$$ and find that $$A^n = left(QD^{1/n}Q^{-1}right)^n = Qleft(D^{1/n}right)^nQ^{-1} = QDQ^{-1} = M$$



      To diagonalize $M$, you find it's eigenvalues and their corresponding eigenvectors. If the eigenvectors span the entire space, then $M$ is diagonalizable, the columns of $Q$ are eigenvectors of $M$, and the diagonal elements of $D$ are the eigenvalues.



      Alas, it is not applicable here: $$begin{bmatrix}1 & 1\0&1end{bmatrix}$$
      is not diagonalizable.






      share|cite|improve this answer









      $endgroup$













      • $begingroup$
        If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
        $endgroup$
        – Kyle Miller
        Jan 30 '18 at 5:52
















      10












      $begingroup$

      The answers given by José Carlos Santos and Widawensen are correct. There is a more general alternative method, even through it does not work for your particular example. But I thought I'd give it anyway in case it is useful to others.



      If you want the $n^{text{th}}$ root of the matrix $M$, and $M$ is diagonalizable - i.e., there exists $Q$ such that $D = Q^{-1}MQ$ is a diagonal matrix, then $D^{1/n}$ is just the matrix whose diagonal elements are those of $D$ raised to the $1/n$ power. So you can set $$A = QD^{1/n}Q^{-1}$$ and find that $$A^n = left(QD^{1/n}Q^{-1}right)^n = Qleft(D^{1/n}right)^nQ^{-1} = QDQ^{-1} = M$$



      To diagonalize $M$, you find it's eigenvalues and their corresponding eigenvectors. If the eigenvectors span the entire space, then $M$ is diagonalizable, the columns of $Q$ are eigenvectors of $M$, and the diagonal elements of $D$ are the eigenvalues.



      Alas, it is not applicable here: $$begin{bmatrix}1 & 1\0&1end{bmatrix}$$
      is not diagonalizable.






      share|cite|improve this answer









      $endgroup$













      • $begingroup$
        If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
        $endgroup$
        – Kyle Miller
        Jan 30 '18 at 5:52














      10












      10








      10





      $begingroup$

      The answers given by José Carlos Santos and Widawensen are correct. There is a more general alternative method, even through it does not work for your particular example. But I thought I'd give it anyway in case it is useful to others.



      If you want the $n^{text{th}}$ root of the matrix $M$, and $M$ is diagonalizable - i.e., there exists $Q$ such that $D = Q^{-1}MQ$ is a diagonal matrix, then $D^{1/n}$ is just the matrix whose diagonal elements are those of $D$ raised to the $1/n$ power. So you can set $$A = QD^{1/n}Q^{-1}$$ and find that $$A^n = left(QD^{1/n}Q^{-1}right)^n = Qleft(D^{1/n}right)^nQ^{-1} = QDQ^{-1} = M$$



      To diagonalize $M$, you find it's eigenvalues and their corresponding eigenvectors. If the eigenvectors span the entire space, then $M$ is diagonalizable, the columns of $Q$ are eigenvectors of $M$, and the diagonal elements of $D$ are the eigenvalues.



      Alas, it is not applicable here: $$begin{bmatrix}1 & 1\0&1end{bmatrix}$$
      is not diagonalizable.






      share|cite|improve this answer









      $endgroup$



      The answers given by José Carlos Santos and Widawensen are correct. There is a more general alternative method, even through it does not work for your particular example. But I thought I'd give it anyway in case it is useful to others.



      If you want the $n^{text{th}}$ root of the matrix $M$, and $M$ is diagonalizable - i.e., there exists $Q$ such that $D = Q^{-1}MQ$ is a diagonal matrix, then $D^{1/n}$ is just the matrix whose diagonal elements are those of $D$ raised to the $1/n$ power. So you can set $$A = QD^{1/n}Q^{-1}$$ and find that $$A^n = left(QD^{1/n}Q^{-1}right)^n = Qleft(D^{1/n}right)^nQ^{-1} = QDQ^{-1} = M$$



      To diagonalize $M$, you find it's eigenvalues and their corresponding eigenvectors. If the eigenvectors span the entire space, then $M$ is diagonalizable, the columns of $Q$ are eigenvectors of $M$, and the diagonal elements of $D$ are the eigenvalues.



      Alas, it is not applicable here: $$begin{bmatrix}1 & 1\0&1end{bmatrix}$$
      is not diagonalizable.







      share|cite|improve this answer












      share|cite|improve this answer



      share|cite|improve this answer










      answered Jan 29 '18 at 18:06









      Paul SinclairPaul Sinclair

      20.4k21443




      20.4k21443












      • $begingroup$
        If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
        $endgroup$
        – Kyle Miller
        Jan 30 '18 at 5:52


















      • $begingroup$
        If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
        $endgroup$
        – Kyle Miller
        Jan 30 '18 at 5:52
















      $begingroup$
      If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
      $endgroup$
      – Kyle Miller
      Jan 30 '18 at 5:52




      $begingroup$
      If you can figure out how exponentials of Jordan blocks work, you can use Jordan normal form instead of just diagonalization.
      $endgroup$
      – Kyle Miller
      Jan 30 '18 at 5:52











      5












      $begingroup$

      Ok this may seem a bit overkill for this particular question, but please bear with me as it could help others.



      You can use the Newton-Rhapson method, trying to solve the equation



      $$f(x) = x^{10}-d=0$$
      where
      $$f'(x) = 10x^9$$
      and $d$ is that matrix of yours and the iteration:



      $${bf X_{n+1}} ={bf X_n}- f({bf X_n})f'({bf X_n})^{-1}$$
      As long as you set initial $bf X_0$ to nothing too crazy, it will work. In fact, this is the simplest example of an discrete fractional integration solution, if you look at larger matrices $d$ filled with 1 on and below the diagonal and 0 otherwise (ok, transpose of the same matrix but same idea).



      The fractional integral operators in fractional calculus found if guessing at the $bf X_0 = bf I$ matrix (described as linear convolutional filters, corresponding to a row in the solution matrix):



      enter image description here



      The two we recognize is the one in the middle (constant, corresponding to normal integral) and the one at the end (linear, corresponding to double integration).






      share|cite|improve this answer











      $endgroup$













      • $begingroup$
        Is that supposed to be $f'(x)=10x^9$?
        $endgroup$
        – ASKASK
        Jan 30 '18 at 11:52










      • $begingroup$
        @ASKASK yep must have gotten an elf in my computer who ate it away ;)
        $endgroup$
        – mathreadler
        Jan 30 '18 at 11:53


















      5












      $begingroup$

      Ok this may seem a bit overkill for this particular question, but please bear with me as it could help others.



      You can use the Newton-Rhapson method, trying to solve the equation



      $$f(x) = x^{10}-d=0$$
      where
      $$f'(x) = 10x^9$$
      and $d$ is that matrix of yours and the iteration:



      $${bf X_{n+1}} ={bf X_n}- f({bf X_n})f'({bf X_n})^{-1}$$
      As long as you set initial $bf X_0$ to nothing too crazy, it will work. In fact, this is the simplest example of an discrete fractional integration solution, if you look at larger matrices $d$ filled with 1 on and below the diagonal and 0 otherwise (ok, transpose of the same matrix but same idea).



      The fractional integral operators in fractional calculus found if guessing at the $bf X_0 = bf I$ matrix (described as linear convolutional filters, corresponding to a row in the solution matrix):



      enter image description here



      The two we recognize is the one in the middle (constant, corresponding to normal integral) and the one at the end (linear, corresponding to double integration).






      share|cite|improve this answer











      $endgroup$













      • $begingroup$
        Is that supposed to be $f'(x)=10x^9$?
        $endgroup$
        – ASKASK
        Jan 30 '18 at 11:52










      • $begingroup$
        @ASKASK yep must have gotten an elf in my computer who ate it away ;)
        $endgroup$
        – mathreadler
        Jan 30 '18 at 11:53
















      5












      5








      5





      $begingroup$

      Ok this may seem a bit overkill for this particular question, but please bear with me as it could help others.



      You can use the Newton-Rhapson method, trying to solve the equation



      $$f(x) = x^{10}-d=0$$
      where
      $$f'(x) = 10x^9$$
      and $d$ is that matrix of yours and the iteration:



      $${bf X_{n+1}} ={bf X_n}- f({bf X_n})f'({bf X_n})^{-1}$$
      As long as you set initial $bf X_0$ to nothing too crazy, it will work. In fact, this is the simplest example of an discrete fractional integration solution, if you look at larger matrices $d$ filled with 1 on and below the diagonal and 0 otherwise (ok, transpose of the same matrix but same idea).



      The fractional integral operators in fractional calculus found if guessing at the $bf X_0 = bf I$ matrix (described as linear convolutional filters, corresponding to a row in the solution matrix):



      enter image description here



      The two we recognize is the one in the middle (constant, corresponding to normal integral) and the one at the end (linear, corresponding to double integration).






      share|cite|improve this answer











      $endgroup$



      Ok this may seem a bit overkill for this particular question, but please bear with me as it could help others.



      You can use the Newton-Rhapson method, trying to solve the equation



      $$f(x) = x^{10}-d=0$$
      where
      $$f'(x) = 10x^9$$
      and $d$ is that matrix of yours and the iteration:



      $${bf X_{n+1}} ={bf X_n}- f({bf X_n})f'({bf X_n})^{-1}$$
      As long as you set initial $bf X_0$ to nothing too crazy, it will work. In fact, this is the simplest example of an discrete fractional integration solution, if you look at larger matrices $d$ filled with 1 on and below the diagonal and 0 otherwise (ok, transpose of the same matrix but same idea).



      The fractional integral operators in fractional calculus found if guessing at the $bf X_0 = bf I$ matrix (described as linear convolutional filters, corresponding to a row in the solution matrix):



      enter image description here



      The two we recognize is the one in the middle (constant, corresponding to normal integral) and the one at the end (linear, corresponding to double integration).







      share|cite|improve this answer














      share|cite|improve this answer



      share|cite|improve this answer








      edited Jan 30 '18 at 17:40

























      answered Jan 29 '18 at 21:01









      mathreadlermathreadler

      15.1k72263




      15.1k72263












      • $begingroup$
        Is that supposed to be $f'(x)=10x^9$?
        $endgroup$
        – ASKASK
        Jan 30 '18 at 11:52










      • $begingroup$
        @ASKASK yep must have gotten an elf in my computer who ate it away ;)
        $endgroup$
        – mathreadler
        Jan 30 '18 at 11:53




















      • $begingroup$
        Is that supposed to be $f'(x)=10x^9$?
        $endgroup$
        – ASKASK
        Jan 30 '18 at 11:52










      • $begingroup$
        @ASKASK yep must have gotten an elf in my computer who ate it away ;)
        $endgroup$
        – mathreadler
        Jan 30 '18 at 11:53


















      $begingroup$
      Is that supposed to be $f'(x)=10x^9$?
      $endgroup$
      – ASKASK
      Jan 30 '18 at 11:52




      $begingroup$
      Is that supposed to be $f'(x)=10x^9$?
      $endgroup$
      – ASKASK
      Jan 30 '18 at 11:52












      $begingroup$
      @ASKASK yep must have gotten an elf in my computer who ate it away ;)
      $endgroup$
      – mathreadler
      Jan 30 '18 at 11:53






      $begingroup$
      @ASKASK yep must have gotten an elf in my computer who ate it away ;)
      $endgroup$
      – mathreadler
      Jan 30 '18 at 11:53













      3












      $begingroup$

      Such a matrix $A$ will commute with $pmatrix{1&1\0&1}$, and so be of the form $pmatrix{a&b\0&a}$. This means that $a^{10}=1$ and so
      $A=apmatrix{1&c\0&1}$ for some $c$. We immediately get $c=1/10$
      (as long as $10$ is invertible in your field).






      share|cite|improve this answer









      $endgroup$


















        3












        $begingroup$

        Such a matrix $A$ will commute with $pmatrix{1&1\0&1}$, and so be of the form $pmatrix{a&b\0&a}$. This means that $a^{10}=1$ and so
        $A=apmatrix{1&c\0&1}$ for some $c$. We immediately get $c=1/10$
        (as long as $10$ is invertible in your field).






        share|cite|improve this answer









        $endgroup$
















          3












          3








          3





          $begingroup$

          Such a matrix $A$ will commute with $pmatrix{1&1\0&1}$, and so be of the form $pmatrix{a&b\0&a}$. This means that $a^{10}=1$ and so
          $A=apmatrix{1&c\0&1}$ for some $c$. We immediately get $c=1/10$
          (as long as $10$ is invertible in your field).






          share|cite|improve this answer









          $endgroup$



          Such a matrix $A$ will commute with $pmatrix{1&1\0&1}$, and so be of the form $pmatrix{a&b\0&a}$. This means that $a^{10}=1$ and so
          $A=apmatrix{1&c\0&1}$ for some $c$. We immediately get $c=1/10$
          (as long as $10$ is invertible in your field).







          share|cite|improve this answer












          share|cite|improve this answer



          share|cite|improve this answer










          answered Jan 30 '18 at 6:24









          Lord Shark the UnknownLord Shark the Unknown

          106k1161133




          106k1161133























              2












              $begingroup$

              Hint: Another approach is to note that
              $$
              expleft(begin{bmatrix}0&x\0&0end{bmatrix}right)=begin{bmatrix}1&x\0&1end{bmatrix}
              $$






              share|cite|improve this answer









              $endgroup$


















                2












                $begingroup$

                Hint: Another approach is to note that
                $$
                expleft(begin{bmatrix}0&x\0&0end{bmatrix}right)=begin{bmatrix}1&x\0&1end{bmatrix}
                $$






                share|cite|improve this answer









                $endgroup$
















                  2












                  2








                  2





                  $begingroup$

                  Hint: Another approach is to note that
                  $$
                  expleft(begin{bmatrix}0&x\0&0end{bmatrix}right)=begin{bmatrix}1&x\0&1end{bmatrix}
                  $$






                  share|cite|improve this answer









                  $endgroup$



                  Hint: Another approach is to note that
                  $$
                  expleft(begin{bmatrix}0&x\0&0end{bmatrix}right)=begin{bmatrix}1&x\0&1end{bmatrix}
                  $$







                  share|cite|improve this answer












                  share|cite|improve this answer



                  share|cite|improve this answer










                  answered Jan 30 '18 at 18:09









                  robjohnrobjohn

                  269k27311638




                  269k27311638






























                      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.




                      draft saved


                      draft discarded














                      StackExchange.ready(
                      function () {
                      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2626206%2ffinding-the-10th-root-of-a-matrix%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







                      Popular posts from this blog

                      Tonle Sap (See)

                      I get strange results when I access the Sqlitedatabase with Unity C# via XAMPP

                      Guatemaltekische Davis-Cup-Mannschaft