Closed form solution for Double Integrators dynamics












0












$begingroup$


I have a relatively specific problem. Consider a linear system, $dot{x}(t) = Ax(t)+Bu(t)$, where $A in mathbb{R}^{n times n}$ and $B in mathbb{R}^{n times 1}$ are the constant matrix and constant vector that defines the system.



The closed-form for a linear system with a constant control input $u(t) := u_k, forall tin [0, T]$ has the following form:
$x(t) = e^{At}x(0) + int_{0}^{T} e^{A(t-tau)} dtau Bu$. If $A$ is non-singular (invertible), then $int_{0}^{T} e^{A(t-tau)} dtau $ can be rewrite to $A^{-1}(e^{AT}-I)$, where $I$ is an identity matrix with appropriate dimension. At this point, I can rewrite $x(t) = e^{At}x(0) + A^{-1}(e^{AT}-I) Bu$, which does not have the integral.



But a lot of systems in real-life have a singular matrix $A$ for its dynamics ($A$ is not invertible ), such as a double integrators with $A = [0,1;0,0]$. In this case, are there any alternatives for me to write its closed-form solution that does not contain integral?










share|cite|improve this question









$endgroup$

















    0












    $begingroup$


    I have a relatively specific problem. Consider a linear system, $dot{x}(t) = Ax(t)+Bu(t)$, where $A in mathbb{R}^{n times n}$ and $B in mathbb{R}^{n times 1}$ are the constant matrix and constant vector that defines the system.



    The closed-form for a linear system with a constant control input $u(t) := u_k, forall tin [0, T]$ has the following form:
    $x(t) = e^{At}x(0) + int_{0}^{T} e^{A(t-tau)} dtau Bu$. If $A$ is non-singular (invertible), then $int_{0}^{T} e^{A(t-tau)} dtau $ can be rewrite to $A^{-1}(e^{AT}-I)$, where $I$ is an identity matrix with appropriate dimension. At this point, I can rewrite $x(t) = e^{At}x(0) + A^{-1}(e^{AT}-I) Bu$, which does not have the integral.



    But a lot of systems in real-life have a singular matrix $A$ for its dynamics ($A$ is not invertible ), such as a double integrators with $A = [0,1;0,0]$. In this case, are there any alternatives for me to write its closed-form solution that does not contain integral?










    share|cite|improve this question









    $endgroup$















      0












      0








      0





      $begingroup$


      I have a relatively specific problem. Consider a linear system, $dot{x}(t) = Ax(t)+Bu(t)$, where $A in mathbb{R}^{n times n}$ and $B in mathbb{R}^{n times 1}$ are the constant matrix and constant vector that defines the system.



      The closed-form for a linear system with a constant control input $u(t) := u_k, forall tin [0, T]$ has the following form:
      $x(t) = e^{At}x(0) + int_{0}^{T} e^{A(t-tau)} dtau Bu$. If $A$ is non-singular (invertible), then $int_{0}^{T} e^{A(t-tau)} dtau $ can be rewrite to $A^{-1}(e^{AT}-I)$, where $I$ is an identity matrix with appropriate dimension. At this point, I can rewrite $x(t) = e^{At}x(0) + A^{-1}(e^{AT}-I) Bu$, which does not have the integral.



      But a lot of systems in real-life have a singular matrix $A$ for its dynamics ($A$ is not invertible ), such as a double integrators with $A = [0,1;0,0]$. In this case, are there any alternatives for me to write its closed-form solution that does not contain integral?










      share|cite|improve this question









      $endgroup$




      I have a relatively specific problem. Consider a linear system, $dot{x}(t) = Ax(t)+Bu(t)$, where $A in mathbb{R}^{n times n}$ and $B in mathbb{R}^{n times 1}$ are the constant matrix and constant vector that defines the system.



      The closed-form for a linear system with a constant control input $u(t) := u_k, forall tin [0, T]$ has the following form:
      $x(t) = e^{At}x(0) + int_{0}^{T} e^{A(t-tau)} dtau Bu$. If $A$ is non-singular (invertible), then $int_{0}^{T} e^{A(t-tau)} dtau $ can be rewrite to $A^{-1}(e^{AT}-I)$, where $I$ is an identity matrix with appropriate dimension. At this point, I can rewrite $x(t) = e^{At}x(0) + A^{-1}(e^{AT}-I) Bu$, which does not have the integral.



      But a lot of systems in real-life have a singular matrix $A$ for its dynamics ($A$ is not invertible ), such as a double integrators with $A = [0,1;0,0]$. In this case, are there any alternatives for me to write its closed-form solution that does not contain integral?







      linear-algebra control-theory






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Dec 3 '18 at 17:57









      axpilllaxpilll

      32




      32






















          1 Answer
          1






          active

          oldest

          votes


















          0












          $begingroup$

          You can use a trick that is also often used to obtain a zero-order-hold discetization. Namely



          $$
          e^{begin{bmatrix}A & B \ 0 & 0end{bmatrix}t} =
          begin{bmatrix}e^{A,t} & int_0^t e^{A(t-tau)} dtau,B \ 0 & Iend{bmatrix}.
          $$



          If you want to verify this for yourself you can try to expand the matrix exponential, which can be defined as



          $$
          e^M = I + M + tfrac{1}{2}M^2 + tfrac{1}{3!}M^3 + cdots.
          $$






          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%2f3024436%2fclosed-form-solution-for-double-integrators-dynamics%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown

























            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            0












            $begingroup$

            You can use a trick that is also often used to obtain a zero-order-hold discetization. Namely



            $$
            e^{begin{bmatrix}A & B \ 0 & 0end{bmatrix}t} =
            begin{bmatrix}e^{A,t} & int_0^t e^{A(t-tau)} dtau,B \ 0 & Iend{bmatrix}.
            $$



            If you want to verify this for yourself you can try to expand the matrix exponential, which can be defined as



            $$
            e^M = I + M + tfrac{1}{2}M^2 + tfrac{1}{3!}M^3 + cdots.
            $$






            share|cite|improve this answer









            $endgroup$


















              0












              $begingroup$

              You can use a trick that is also often used to obtain a zero-order-hold discetization. Namely



              $$
              e^{begin{bmatrix}A & B \ 0 & 0end{bmatrix}t} =
              begin{bmatrix}e^{A,t} & int_0^t e^{A(t-tau)} dtau,B \ 0 & Iend{bmatrix}.
              $$



              If you want to verify this for yourself you can try to expand the matrix exponential, which can be defined as



              $$
              e^M = I + M + tfrac{1}{2}M^2 + tfrac{1}{3!}M^3 + cdots.
              $$






              share|cite|improve this answer









              $endgroup$
















                0












                0








                0





                $begingroup$

                You can use a trick that is also often used to obtain a zero-order-hold discetization. Namely



                $$
                e^{begin{bmatrix}A & B \ 0 & 0end{bmatrix}t} =
                begin{bmatrix}e^{A,t} & int_0^t e^{A(t-tau)} dtau,B \ 0 & Iend{bmatrix}.
                $$



                If you want to verify this for yourself you can try to expand the matrix exponential, which can be defined as



                $$
                e^M = I + M + tfrac{1}{2}M^2 + tfrac{1}{3!}M^3 + cdots.
                $$






                share|cite|improve this answer









                $endgroup$



                You can use a trick that is also often used to obtain a zero-order-hold discetization. Namely



                $$
                e^{begin{bmatrix}A & B \ 0 & 0end{bmatrix}t} =
                begin{bmatrix}e^{A,t} & int_0^t e^{A(t-tau)} dtau,B \ 0 & Iend{bmatrix}.
                $$



                If you want to verify this for yourself you can try to expand the matrix exponential, which can be defined as



                $$
                e^M = I + M + tfrac{1}{2}M^2 + tfrac{1}{3!}M^3 + cdots.
                $$







                share|cite|improve this answer












                share|cite|improve this answer



                share|cite|improve this answer










                answered Dec 3 '18 at 23:41









                Kwin van der VeenKwin van der Veen

                5,5152828




                5,5152828






























                    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%2f3024436%2fclosed-form-solution-for-double-integrators-dynamics%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

                    How to change which sound is reproduced for terminal bell?

                    Title Spacing in Bjornstrup Chapter, Removing Chapter Number From Contents

                    Can I use Tabulator js library in my java Spring + Thymeleaf project?