System of linear differential equation

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











up vote
3
down vote

favorite
1












How to solve the system of linear differential equation of the form $$x' = Ax + b$$ I can solve the homogeneous form by finding the eigenvalues and respective eigenvectors, but how to find the particular solution part. Also is there any limitation from getting eigenvalues positive, negative or complex. Any other different method involving matrix algebra is also welcome. Thank you.







share|cite|improve this question





















  • Thank you very much. All the answers are really helpful.
    – Mridul Dey
    Jul 21 at 18:08










  • You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
    – Davide Morgante
    Jul 21 at 19:36















up vote
3
down vote

favorite
1












How to solve the system of linear differential equation of the form $$x' = Ax + b$$ I can solve the homogeneous form by finding the eigenvalues and respective eigenvectors, but how to find the particular solution part. Also is there any limitation from getting eigenvalues positive, negative or complex. Any other different method involving matrix algebra is also welcome. Thank you.







share|cite|improve this question





















  • Thank you very much. All the answers are really helpful.
    – Mridul Dey
    Jul 21 at 18:08










  • You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
    – Davide Morgante
    Jul 21 at 19:36













up vote
3
down vote

favorite
1









up vote
3
down vote

favorite
1






1





How to solve the system of linear differential equation of the form $$x' = Ax + b$$ I can solve the homogeneous form by finding the eigenvalues and respective eigenvectors, but how to find the particular solution part. Also is there any limitation from getting eigenvalues positive, negative or complex. Any other different method involving matrix algebra is also welcome. Thank you.







share|cite|improve this question













How to solve the system of linear differential equation of the form $$x' = Ax + b$$ I can solve the homogeneous form by finding the eigenvalues and respective eigenvectors, but how to find the particular solution part. Also is there any limitation from getting eigenvalues positive, negative or complex. Any other different method involving matrix algebra is also welcome. Thank you.









share|cite|improve this question












share|cite|improve this question




share|cite|improve this question








edited Jul 21 at 17:39









Rodrigo de Azevedo

12.6k41751




12.6k41751









asked Jul 21 at 17:36









Mridul Dey

765




765











  • Thank you very much. All the answers are really helpful.
    – Mridul Dey
    Jul 21 at 18:08










  • You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
    – Davide Morgante
    Jul 21 at 19:36

















  • Thank you very much. All the answers are really helpful.
    – Mridul Dey
    Jul 21 at 18:08










  • You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
    – Davide Morgante
    Jul 21 at 19:36
















Thank you very much. All the answers are really helpful.
– Mridul Dey
Jul 21 at 18:08




Thank you very much. All the answers are really helpful.
– Mridul Dey
Jul 21 at 18:08












You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
– Davide Morgante
Jul 21 at 19:36





You're welcome! If you feel like one of the answers is what you were looking for, you should mark it as accepted!
– Davide Morgante
Jul 21 at 19:36











3 Answers
3






active

oldest

votes

















up vote
4
down vote













Being $b$ a fixed constant (suppose $det Aneq0$), we can choose a new variable $x_*$ that solves $$Ax_* + b = 0 Rightarrow x_* = -A^-1b$$ The now we can solve the equation in the fixed point $z = x-x_* = x+A^-1b$. This variable solves the following equation $$dotz = Az$$
In fact, form the definition of $z$ $$dotz=dotx\ Az = A(x+A^-1b) = Ax+b$$
So we're saying that, if $x$ is a solution to $dotx=Ax+b$ then the variable $z = x+A^-1b$ is a solution to the simpler equation $dotz=Az$.



The solution to this simpler equation is just $$z(t)= [e^At]z(0)$$ where $[e^At]$ is the exponential of the matrix $A$. At this point you can go back to find $x(t)$ simply by using the definition of $z$: $$x(t) = z(t)-A^-1b = [e^At](x(0)+A^-1b)-A^-1b = [e^At]x(0)+A^-1b([e^At]-1)$$ where we have used the fact that $$z(0) = x(0) + A^-1b$$



Easier way



You can even use this simpler method: given the system of equations $$dotx = Ax+f(t)$$ the solution to this is as follows $$x(t) = [e^At]x(0)+int_0^t e^A(t-t')f(t')dt'$$
in your case the function is just a constant function $f(t)=b$ for which you can see that gives the same solution as before! These methods are good because they let you solve every system of equations of this type because all the informations about the system is embedded into the exponential of the matrix






share|cite|improve this answer






























    up vote
    3
    down vote













    If $b$ is the columnspace of $A$, define $y$ so that $Ay = -b$. Then note that $x(t) = e^Atx(0) + y$ solves the differential equation since $x' = A e^Atx(0)$ and $Ax = A e^Atx(0) - b$.






    share|cite|improve this answer




























      up vote
      2
      down vote













      Take Laplace transform therefore $$sX(s)-x(0)=AX(s)+b\(sI-A)X(s)=b+x(0)$$or$$X(s)=(sI-A)^-1(b+x(0))$$but we can't proceed any further. We need to know the characteristics of the system.






      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%2f2858713%2fsystem-of-linear-differential-equation%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
        4
        down vote













        Being $b$ a fixed constant (suppose $det Aneq0$), we can choose a new variable $x_*$ that solves $$Ax_* + b = 0 Rightarrow x_* = -A^-1b$$ The now we can solve the equation in the fixed point $z = x-x_* = x+A^-1b$. This variable solves the following equation $$dotz = Az$$
        In fact, form the definition of $z$ $$dotz=dotx\ Az = A(x+A^-1b) = Ax+b$$
        So we're saying that, if $x$ is a solution to $dotx=Ax+b$ then the variable $z = x+A^-1b$ is a solution to the simpler equation $dotz=Az$.



        The solution to this simpler equation is just $$z(t)= [e^At]z(0)$$ where $[e^At]$ is the exponential of the matrix $A$. At this point you can go back to find $x(t)$ simply by using the definition of $z$: $$x(t) = z(t)-A^-1b = [e^At](x(0)+A^-1b)-A^-1b = [e^At]x(0)+A^-1b([e^At]-1)$$ where we have used the fact that $$z(0) = x(0) + A^-1b$$



        Easier way



        You can even use this simpler method: given the system of equations $$dotx = Ax+f(t)$$ the solution to this is as follows $$x(t) = [e^At]x(0)+int_0^t e^A(t-t')f(t')dt'$$
        in your case the function is just a constant function $f(t)=b$ for which you can see that gives the same solution as before! These methods are good because they let you solve every system of equations of this type because all the informations about the system is embedded into the exponential of the matrix






        share|cite|improve this answer



























          up vote
          4
          down vote













          Being $b$ a fixed constant (suppose $det Aneq0$), we can choose a new variable $x_*$ that solves $$Ax_* + b = 0 Rightarrow x_* = -A^-1b$$ The now we can solve the equation in the fixed point $z = x-x_* = x+A^-1b$. This variable solves the following equation $$dotz = Az$$
          In fact, form the definition of $z$ $$dotz=dotx\ Az = A(x+A^-1b) = Ax+b$$
          So we're saying that, if $x$ is a solution to $dotx=Ax+b$ then the variable $z = x+A^-1b$ is a solution to the simpler equation $dotz=Az$.



          The solution to this simpler equation is just $$z(t)= [e^At]z(0)$$ where $[e^At]$ is the exponential of the matrix $A$. At this point you can go back to find $x(t)$ simply by using the definition of $z$: $$x(t) = z(t)-A^-1b = [e^At](x(0)+A^-1b)-A^-1b = [e^At]x(0)+A^-1b([e^At]-1)$$ where we have used the fact that $$z(0) = x(0) + A^-1b$$



          Easier way



          You can even use this simpler method: given the system of equations $$dotx = Ax+f(t)$$ the solution to this is as follows $$x(t) = [e^At]x(0)+int_0^t e^A(t-t')f(t')dt'$$
          in your case the function is just a constant function $f(t)=b$ for which you can see that gives the same solution as before! These methods are good because they let you solve every system of equations of this type because all the informations about the system is embedded into the exponential of the matrix






          share|cite|improve this answer

























            up vote
            4
            down vote










            up vote
            4
            down vote









            Being $b$ a fixed constant (suppose $det Aneq0$), we can choose a new variable $x_*$ that solves $$Ax_* + b = 0 Rightarrow x_* = -A^-1b$$ The now we can solve the equation in the fixed point $z = x-x_* = x+A^-1b$. This variable solves the following equation $$dotz = Az$$
            In fact, form the definition of $z$ $$dotz=dotx\ Az = A(x+A^-1b) = Ax+b$$
            So we're saying that, if $x$ is a solution to $dotx=Ax+b$ then the variable $z = x+A^-1b$ is a solution to the simpler equation $dotz=Az$.



            The solution to this simpler equation is just $$z(t)= [e^At]z(0)$$ where $[e^At]$ is the exponential of the matrix $A$. At this point you can go back to find $x(t)$ simply by using the definition of $z$: $$x(t) = z(t)-A^-1b = [e^At](x(0)+A^-1b)-A^-1b = [e^At]x(0)+A^-1b([e^At]-1)$$ where we have used the fact that $$z(0) = x(0) + A^-1b$$



            Easier way



            You can even use this simpler method: given the system of equations $$dotx = Ax+f(t)$$ the solution to this is as follows $$x(t) = [e^At]x(0)+int_0^t e^A(t-t')f(t')dt'$$
            in your case the function is just a constant function $f(t)=b$ for which you can see that gives the same solution as before! These methods are good because they let you solve every system of equations of this type because all the informations about the system is embedded into the exponential of the matrix






            share|cite|improve this answer















            Being $b$ a fixed constant (suppose $det Aneq0$), we can choose a new variable $x_*$ that solves $$Ax_* + b = 0 Rightarrow x_* = -A^-1b$$ The now we can solve the equation in the fixed point $z = x-x_* = x+A^-1b$. This variable solves the following equation $$dotz = Az$$
            In fact, form the definition of $z$ $$dotz=dotx\ Az = A(x+A^-1b) = Ax+b$$
            So we're saying that, if $x$ is a solution to $dotx=Ax+b$ then the variable $z = x+A^-1b$ is a solution to the simpler equation $dotz=Az$.



            The solution to this simpler equation is just $$z(t)= [e^At]z(0)$$ where $[e^At]$ is the exponential of the matrix $A$. At this point you can go back to find $x(t)$ simply by using the definition of $z$: $$x(t) = z(t)-A^-1b = [e^At](x(0)+A^-1b)-A^-1b = [e^At]x(0)+A^-1b([e^At]-1)$$ where we have used the fact that $$z(0) = x(0) + A^-1b$$



            Easier way



            You can even use this simpler method: given the system of equations $$dotx = Ax+f(t)$$ the solution to this is as follows $$x(t) = [e^At]x(0)+int_0^t e^A(t-t')f(t')dt'$$
            in your case the function is just a constant function $f(t)=b$ for which you can see that gives the same solution as before! These methods are good because they let you solve every system of equations of this type because all the informations about the system is embedded into the exponential of the matrix







            share|cite|improve this answer















            share|cite|improve this answer



            share|cite|improve this answer








            edited Jul 21 at 17:55


























            answered Jul 21 at 17:49









            Davide Morgante

            1,812220




            1,812220




















                up vote
                3
                down vote













                If $b$ is the columnspace of $A$, define $y$ so that $Ay = -b$. Then note that $x(t) = e^Atx(0) + y$ solves the differential equation since $x' = A e^Atx(0)$ and $Ax = A e^Atx(0) - b$.






                share|cite|improve this answer

























                  up vote
                  3
                  down vote













                  If $b$ is the columnspace of $A$, define $y$ so that $Ay = -b$. Then note that $x(t) = e^Atx(0) + y$ solves the differential equation since $x' = A e^Atx(0)$ and $Ax = A e^Atx(0) - b$.






                  share|cite|improve this answer























                    up vote
                    3
                    down vote










                    up vote
                    3
                    down vote









                    If $b$ is the columnspace of $A$, define $y$ so that $Ay = -b$. Then note that $x(t) = e^Atx(0) + y$ solves the differential equation since $x' = A e^Atx(0)$ and $Ax = A e^Atx(0) - b$.






                    share|cite|improve this answer













                    If $b$ is the columnspace of $A$, define $y$ so that $Ay = -b$. Then note that $x(t) = e^Atx(0) + y$ solves the differential equation since $x' = A e^Atx(0)$ and $Ax = A e^Atx(0) - b$.







                    share|cite|improve this answer













                    share|cite|improve this answer



                    share|cite|improve this answer











                    answered Jul 21 at 17:48









                    Marcus M

                    8,1731847




                    8,1731847




















                        up vote
                        2
                        down vote













                        Take Laplace transform therefore $$sX(s)-x(0)=AX(s)+b\(sI-A)X(s)=b+x(0)$$or$$X(s)=(sI-A)^-1(b+x(0))$$but we can't proceed any further. We need to know the characteristics of the system.






                        share|cite|improve this answer

























                          up vote
                          2
                          down vote













                          Take Laplace transform therefore $$sX(s)-x(0)=AX(s)+b\(sI-A)X(s)=b+x(0)$$or$$X(s)=(sI-A)^-1(b+x(0))$$but we can't proceed any further. We need to know the characteristics of the system.






                          share|cite|improve this answer























                            up vote
                            2
                            down vote










                            up vote
                            2
                            down vote









                            Take Laplace transform therefore $$sX(s)-x(0)=AX(s)+b\(sI-A)X(s)=b+x(0)$$or$$X(s)=(sI-A)^-1(b+x(0))$$but we can't proceed any further. We need to know the characteristics of the system.






                            share|cite|improve this answer













                            Take Laplace transform therefore $$sX(s)-x(0)=AX(s)+b\(sI-A)X(s)=b+x(0)$$or$$X(s)=(sI-A)^-1(b+x(0))$$but we can't proceed any further. We need to know the characteristics of the system.







                            share|cite|improve this answer













                            share|cite|improve this answer



                            share|cite|improve this answer











                            answered Jul 21 at 17:44









                            Mostafa Ayaz

                            8,5773630




                            8,5773630






















                                 

                                draft saved


                                draft discarded


























                                 


                                draft saved


                                draft discarded














                                StackExchange.ready(
                                function ()
                                StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2858713%2fsystem-of-linear-differential-equation%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?