Particular solution of second order differential equation
Given the ode:
$$
y''-2y'+y=e^t,
$$
how can I find the form of the particular solution?
At first, I tried the form $y=Ae^t$ but
$$
begin{split}
&frac{d^2y}{dt^2}Ae^t=frac{dy}{dt}Ae^t=Ae^t\
\
&Ae^t-2Ae^t+Ae^t=e^t\
\
&0=e^t\
end{split}.
$$
So this doesn't work.
I also tried the form $y=Ate^t$, but again
$$
begin{split}
&frac{d^2y}{dt^2}=A(2e^t+te^t)\
\
&frac{dy}{dt}=A(e^t+te^t)\
\
&A(2e^t+te^t)-2A(e^t+te^t)+Ate^t=e^t\
\
&2A+At-2A-2At+At=1\
end{split}.
$$
and again this doesn't work
Generally, what is the best way to guess the form of the solution?
differential-equations
add a comment |
Given the ode:
$$
y''-2y'+y=e^t,
$$
how can I find the form of the particular solution?
At first, I tried the form $y=Ae^t$ but
$$
begin{split}
&frac{d^2y}{dt^2}Ae^t=frac{dy}{dt}Ae^t=Ae^t\
\
&Ae^t-2Ae^t+Ae^t=e^t\
\
&0=e^t\
end{split}.
$$
So this doesn't work.
I also tried the form $y=Ate^t$, but again
$$
begin{split}
&frac{d^2y}{dt^2}=A(2e^t+te^t)\
\
&frac{dy}{dt}=A(e^t+te^t)\
\
&A(2e^t+te^t)-2A(e^t+te^t)+Ate^t=e^t\
\
&2A+At-2A-2At+At=1\
end{split}.
$$
and again this doesn't work
Generally, what is the best way to guess the form of the solution?
differential-equations
2
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
2
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17
add a comment |
Given the ode:
$$
y''-2y'+y=e^t,
$$
how can I find the form of the particular solution?
At first, I tried the form $y=Ae^t$ but
$$
begin{split}
&frac{d^2y}{dt^2}Ae^t=frac{dy}{dt}Ae^t=Ae^t\
\
&Ae^t-2Ae^t+Ae^t=e^t\
\
&0=e^t\
end{split}.
$$
So this doesn't work.
I also tried the form $y=Ate^t$, but again
$$
begin{split}
&frac{d^2y}{dt^2}=A(2e^t+te^t)\
\
&frac{dy}{dt}=A(e^t+te^t)\
\
&A(2e^t+te^t)-2A(e^t+te^t)+Ate^t=e^t\
\
&2A+At-2A-2At+At=1\
end{split}.
$$
and again this doesn't work
Generally, what is the best way to guess the form of the solution?
differential-equations
Given the ode:
$$
y''-2y'+y=e^t,
$$
how can I find the form of the particular solution?
At first, I tried the form $y=Ae^t$ but
$$
begin{split}
&frac{d^2y}{dt^2}Ae^t=frac{dy}{dt}Ae^t=Ae^t\
\
&Ae^t-2Ae^t+Ae^t=e^t\
\
&0=e^t\
end{split}.
$$
So this doesn't work.
I also tried the form $y=Ate^t$, but again
$$
begin{split}
&frac{d^2y}{dt^2}=A(2e^t+te^t)\
\
&frac{dy}{dt}=A(e^t+te^t)\
\
&A(2e^t+te^t)-2A(e^t+te^t)+Ate^t=e^t\
\
&2A+At-2A-2At+At=1\
end{split}.
$$
and again this doesn't work
Generally, what is the best way to guess the form of the solution?
differential-equations
differential-equations
edited Dec 2 at 10:50
amWhy
191k28224439
191k28224439
asked Dec 1 at 22:12
VakiPitsi
1637
1637
2
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
2
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17
add a comment |
2
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
2
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17
2
2
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
2
2
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17
add a comment |
3 Answers
3
active
oldest
votes
Hints/Guides on how to solve such differential equations :
$mathbf{1}$ - Method of Undetermined Coefficients :
Start of by solving the homogenous equation $y''-2y'+y= 0$ by assuming that a solution will be proportional to $e^{lambda t}$ for some $lambda$. Substitute in and calculate $lambda$. Notice the multiplicity of the solution for $lambda$ and adjust your general solution accordingly.
Then, use the method of undetermined coefficients to find a particular solution of the problem for $y''-2y'+y=e^t$.
The general solution of the initial differential equation, will then be the general solution of the homogenous plus the particular solution you found.
You can find more information and examples about that method, here.
$mathbf{2}$ - Laplace Transformation :
This is a very fast and straight forward way to approach the problem, but it needs some fluent handling of Laplace Transformation techniques. Note that you can apply the Laplace Transformation without even needing initial conditions, simply stating them as constants.
Start of by applying the Laplace Transformation
$$mathcal{L}_tbig[f(t)big](s) = int_0^infty f(t)e^{-st}mathrm{d}t$$
to both sides of the given differential equation :
$$mathcal{L}_tbig[y'' - 2y' + y'] = mathcal{L}_t[e^t]$$
$$Leftrightarrow$$
$$(s-1)^2big[mathcal{L}_t[y(t)](s)big] - (s-2)y(0) - y'(0) = frac{1}{s-1}$$
$$Leftrightarrow$$
$$mathcal{L}_tbig[y(t)big](s) = frac{y(0)(s^2-3s+2) + y'(0)(s-1) + 1}{(s-1)^3}$$
$$=$$
$$mathcal{L}_tbig[y(t)big](s) = frac{1}{(s-1)^3} - frac{y(0)}{(s-1)^2} + frac{y(0)}{s-1} + frac{y'(0)}{(s-1)^2} $$
$$implies$$
$$y(t) = frac{1}{2}e^t(t^2+2c_1 - 2c_1t + 2c_2t) = frac{e^tt^2}{2} + c_1e^t - c_1e^tt + c_2e^tt$$
$mathbf{3}$ - Variation of Parameters :
You must repeat the step of solving the homogenous equation by finding that $lambda$s mentioned. Then by listing the basis solution as $y_{b_1} = e^t$ and $_{b_2} = e^tt$ you can use variation of parameters to find the final general solution by computing the Wronskian and finding the integrals :
$$v_1(t) = - int frac{f(t)y_{b_2}(t)}{W(t)}mathrm{d}t quad text{and} quad v_2(t) = int frac{f(t)y_{b_1}(t)}{W(t)}mathrm{d}t$$
You can find more information and examples about that method, here.
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
add a comment |
Both your attempts are in fact right but fail because the fundamental set of solutions for your second order ODE is given by exactly your both guesses for the particular solution. It is not hard to show by using the characteristic equation that the fundamental set of solutions is given by
$$y(t)=c_1e^t+c_2te^t$$
Therefore in this case it is not possible to obtain a particular solution in the standard way hence the inhomogeneous term is in fact part of the solution.
Howsoever you second try was near the actual solution. Since both terms $e^t$ and $te^t$ did not worked out hence there are part of the solution you could further consider $t^2e^t$ as the next try. It turns out that this yields to the solution $($see here$)$.
This leads to the conjecture that at least for a inhomogeneous term of the form $Ae^{bt}$ you just have to try the specific exponential with different powers of the variable infront as particular solution until it works out.
add a comment |
I will answer to the question by considering the general $n$-th order constant coefficients linear ODE
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=flabel{1}tag{1}
$$
and the associated linear homogeneous equation
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=0label{a}tag{1'}
$$
where $fnotequiv 0$: the choice $a_n=1$ is only for formal simplification in the development and does not restrict the generality of the analysis.
There are basically two methods for finding a particular solution $y_p$ of the ODE eqref{1}:
by guessing, based on the solver's experience: this method will not analyzed here.
by choosing a solution $y_o$ of the associated homogeneous equation eqref{a} satisfying
$$
y_o(0)=y_o^{(1)}(0)=dots=y_o^{(n-2)}(0)quad y_o^{(n-1)}(0)=1label{2}tag{2}
$$
and forming the fundamental solution $mathscr{E}$ of eqref{1}
$$
mathscr{E}(t)=H(t)y_o(t),label{3}tag{3}
$$
where $H(t)$ is the Heaviside function. Then the sought for particular solution is
$$
y_p(t)=mathscr{E}ast f(t)=intlimits_{0}^ty_o(t-s)f(s)mathrm{d}s label{4}tag{4}
$$
Let's apply formula eqref{4} to the OP problem, before analyzing why it gives the sought for result. Since the characteristic equation is
$$
x^2-2x+1=0iff x=1 text{ with multiplicity 2}
$$
we have that a fundamental system of solutions of the homogeneous equation associated to the given one is
$$
y_1(t)=e^t,: y_2(t)=te^timplies y_o=y_2(t)
$$
since it is the only solution satisfying condition eqref{2}, i.e. $y_o(0)=0$ and $dfrac{mathrm{d}}{mathrm{d}t}y_o(0)=0$. Now we have that
$$
mathscr{E}(t)=H(t)y_o(t)=H(t)te^t
$$
and by applying eqref{4} we obtain
$$
begin{split}
y_p(t)=mathscr{E}ast exp(t)&=intlimits_{0}^{+infty}H(t-s)(t-s)e^{t-s}e^smathrm{d}s\
&=intlimits_{0}^ty_o(t-s)e^smathrm{d}s\
&=intlimits_{0}^t(t-s)e^{t}mathrm{d}s={t^2 over 2}e^t
end{split}
$$
which is the sought for particular solution of the ODE proposed as an example.
Why the methods works? Because of the properties of distributions ([1], §4.9, example 4.9.6 pp. 77-74 and §15.4, example 15.4.4): by using these properties and the condition eqref{2} we have
$$
begin{split}
frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)&=H(t)frac{mathrm{d}}{mathrm{d}t}y_o(t)\
frac{mathrm{d}^2}{mathrm{d}t^2}mathscr{E}(t)&=H(t)frac{mathrm{d}^2}{mathrm{d}t^2}y_o(t)\
&vdots\
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&=H(t)frac{mathrm{d}^n}{mathrm{d}t^n}y_o(t)+delta(t)
end{split}
$$
and thus
$$
begin{split}
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}mathscr{E}(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)+a_0mathscr{E}(t)\
=&H(t)Big[frac{mathrm{d}^n}{mathrm{d}t^n}y_0(t)+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}y_0(t)(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}y_0(t)+a_0y_0(t)Big]+delta(t)=delta(t)
end{split}
$$
The linearity of the differential operator and the properties of the Dirac $delta$-distribution do the rest.
Notes
The analyzed method is perhaps the simplest way to find a particular solution of the ODE proposed by the OP, and more generally of eqref{1}, because it requires only the knowledge of a complete systems of solutions of the associated equation eqref{2}. Indeed such a system of solutions is already required to solve any problem (Cauchy, boundary value, etc.) for eqref{1} and eqref{3}: the only further operation to do is calculating $y_o$ and the convolution integral eqref{3}.
The convolution integral used in eqref{3} is the standard one used in the operational calculus of one variable functions, i.e.
$$
mathscr{E}ast f(t)=intlimits_0^{+infty}mathscr{E}(t-s)f(s)mathrm{d}squad mathscr{E},fin L_mathrm{loc}^1(mathbb{R}_+)
$$
which can be deduced from the standard one by considering $H(t)f(t)$ instead of $f(t)$ as homogeneous term.The assumption $a_n=1$ does not restrict the generality of the above analysis because if we assume that we are dealing with an $n$-th order linear constant coefficient ordinary differential operator, we must assume $a_nneq 0$.
The distribution eqref{3} is called fundamental solution exactly because it can be used to construct the solution for every linear, constant coefficient non-homogeneous ODE.
[1] Vladimirov, V. S. (2002), Methods of the theory of generalized functions, Analytical Methods and Special Functions, 6, London–New York: Taylor & Francis, pp. XII+353, ISBN 0-415-27356-0, MR 2012831, Zbl 1078.46029.
add a comment |
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
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3021903%2fparticular-solution-of-second-order-differential-equation%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
3 Answers
3
active
oldest
votes
3 Answers
3
active
oldest
votes
active
oldest
votes
active
oldest
votes
Hints/Guides on how to solve such differential equations :
$mathbf{1}$ - Method of Undetermined Coefficients :
Start of by solving the homogenous equation $y''-2y'+y= 0$ by assuming that a solution will be proportional to $e^{lambda t}$ for some $lambda$. Substitute in and calculate $lambda$. Notice the multiplicity of the solution for $lambda$ and adjust your general solution accordingly.
Then, use the method of undetermined coefficients to find a particular solution of the problem for $y''-2y'+y=e^t$.
The general solution of the initial differential equation, will then be the general solution of the homogenous plus the particular solution you found.
You can find more information and examples about that method, here.
$mathbf{2}$ - Laplace Transformation :
This is a very fast and straight forward way to approach the problem, but it needs some fluent handling of Laplace Transformation techniques. Note that you can apply the Laplace Transformation without even needing initial conditions, simply stating them as constants.
Start of by applying the Laplace Transformation
$$mathcal{L}_tbig[f(t)big](s) = int_0^infty f(t)e^{-st}mathrm{d}t$$
to both sides of the given differential equation :
$$mathcal{L}_tbig[y'' - 2y' + y'] = mathcal{L}_t[e^t]$$
$$Leftrightarrow$$
$$(s-1)^2big[mathcal{L}_t[y(t)](s)big] - (s-2)y(0) - y'(0) = frac{1}{s-1}$$
$$Leftrightarrow$$
$$mathcal{L}_tbig[y(t)big](s) = frac{y(0)(s^2-3s+2) + y'(0)(s-1) + 1}{(s-1)^3}$$
$$=$$
$$mathcal{L}_tbig[y(t)big](s) = frac{1}{(s-1)^3} - frac{y(0)}{(s-1)^2} + frac{y(0)}{s-1} + frac{y'(0)}{(s-1)^2} $$
$$implies$$
$$y(t) = frac{1}{2}e^t(t^2+2c_1 - 2c_1t + 2c_2t) = frac{e^tt^2}{2} + c_1e^t - c_1e^tt + c_2e^tt$$
$mathbf{3}$ - Variation of Parameters :
You must repeat the step of solving the homogenous equation by finding that $lambda$s mentioned. Then by listing the basis solution as $y_{b_1} = e^t$ and $_{b_2} = e^tt$ you can use variation of parameters to find the final general solution by computing the Wronskian and finding the integrals :
$$v_1(t) = - int frac{f(t)y_{b_2}(t)}{W(t)}mathrm{d}t quad text{and} quad v_2(t) = int frac{f(t)y_{b_1}(t)}{W(t)}mathrm{d}t$$
You can find more information and examples about that method, here.
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
add a comment |
Hints/Guides on how to solve such differential equations :
$mathbf{1}$ - Method of Undetermined Coefficients :
Start of by solving the homogenous equation $y''-2y'+y= 0$ by assuming that a solution will be proportional to $e^{lambda t}$ for some $lambda$. Substitute in and calculate $lambda$. Notice the multiplicity of the solution for $lambda$ and adjust your general solution accordingly.
Then, use the method of undetermined coefficients to find a particular solution of the problem for $y''-2y'+y=e^t$.
The general solution of the initial differential equation, will then be the general solution of the homogenous plus the particular solution you found.
You can find more information and examples about that method, here.
$mathbf{2}$ - Laplace Transformation :
This is a very fast and straight forward way to approach the problem, but it needs some fluent handling of Laplace Transformation techniques. Note that you can apply the Laplace Transformation without even needing initial conditions, simply stating them as constants.
Start of by applying the Laplace Transformation
$$mathcal{L}_tbig[f(t)big](s) = int_0^infty f(t)e^{-st}mathrm{d}t$$
to both sides of the given differential equation :
$$mathcal{L}_tbig[y'' - 2y' + y'] = mathcal{L}_t[e^t]$$
$$Leftrightarrow$$
$$(s-1)^2big[mathcal{L}_t[y(t)](s)big] - (s-2)y(0) - y'(0) = frac{1}{s-1}$$
$$Leftrightarrow$$
$$mathcal{L}_tbig[y(t)big](s) = frac{y(0)(s^2-3s+2) + y'(0)(s-1) + 1}{(s-1)^3}$$
$$=$$
$$mathcal{L}_tbig[y(t)big](s) = frac{1}{(s-1)^3} - frac{y(0)}{(s-1)^2} + frac{y(0)}{s-1} + frac{y'(0)}{(s-1)^2} $$
$$implies$$
$$y(t) = frac{1}{2}e^t(t^2+2c_1 - 2c_1t + 2c_2t) = frac{e^tt^2}{2} + c_1e^t - c_1e^tt + c_2e^tt$$
$mathbf{3}$ - Variation of Parameters :
You must repeat the step of solving the homogenous equation by finding that $lambda$s mentioned. Then by listing the basis solution as $y_{b_1} = e^t$ and $_{b_2} = e^tt$ you can use variation of parameters to find the final general solution by computing the Wronskian and finding the integrals :
$$v_1(t) = - int frac{f(t)y_{b_2}(t)}{W(t)}mathrm{d}t quad text{and} quad v_2(t) = int frac{f(t)y_{b_1}(t)}{W(t)}mathrm{d}t$$
You can find more information and examples about that method, here.
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
add a comment |
Hints/Guides on how to solve such differential equations :
$mathbf{1}$ - Method of Undetermined Coefficients :
Start of by solving the homogenous equation $y''-2y'+y= 0$ by assuming that a solution will be proportional to $e^{lambda t}$ for some $lambda$. Substitute in and calculate $lambda$. Notice the multiplicity of the solution for $lambda$ and adjust your general solution accordingly.
Then, use the method of undetermined coefficients to find a particular solution of the problem for $y''-2y'+y=e^t$.
The general solution of the initial differential equation, will then be the general solution of the homogenous plus the particular solution you found.
You can find more information and examples about that method, here.
$mathbf{2}$ - Laplace Transformation :
This is a very fast and straight forward way to approach the problem, but it needs some fluent handling of Laplace Transformation techniques. Note that you can apply the Laplace Transformation without even needing initial conditions, simply stating them as constants.
Start of by applying the Laplace Transformation
$$mathcal{L}_tbig[f(t)big](s) = int_0^infty f(t)e^{-st}mathrm{d}t$$
to both sides of the given differential equation :
$$mathcal{L}_tbig[y'' - 2y' + y'] = mathcal{L}_t[e^t]$$
$$Leftrightarrow$$
$$(s-1)^2big[mathcal{L}_t[y(t)](s)big] - (s-2)y(0) - y'(0) = frac{1}{s-1}$$
$$Leftrightarrow$$
$$mathcal{L}_tbig[y(t)big](s) = frac{y(0)(s^2-3s+2) + y'(0)(s-1) + 1}{(s-1)^3}$$
$$=$$
$$mathcal{L}_tbig[y(t)big](s) = frac{1}{(s-1)^3} - frac{y(0)}{(s-1)^2} + frac{y(0)}{s-1} + frac{y'(0)}{(s-1)^2} $$
$$implies$$
$$y(t) = frac{1}{2}e^t(t^2+2c_1 - 2c_1t + 2c_2t) = frac{e^tt^2}{2} + c_1e^t - c_1e^tt + c_2e^tt$$
$mathbf{3}$ - Variation of Parameters :
You must repeat the step of solving the homogenous equation by finding that $lambda$s mentioned. Then by listing the basis solution as $y_{b_1} = e^t$ and $_{b_2} = e^tt$ you can use variation of parameters to find the final general solution by computing the Wronskian and finding the integrals :
$$v_1(t) = - int frac{f(t)y_{b_2}(t)}{W(t)}mathrm{d}t quad text{and} quad v_2(t) = int frac{f(t)y_{b_1}(t)}{W(t)}mathrm{d}t$$
You can find more information and examples about that method, here.
Hints/Guides on how to solve such differential equations :
$mathbf{1}$ - Method of Undetermined Coefficients :
Start of by solving the homogenous equation $y''-2y'+y= 0$ by assuming that a solution will be proportional to $e^{lambda t}$ for some $lambda$. Substitute in and calculate $lambda$. Notice the multiplicity of the solution for $lambda$ and adjust your general solution accordingly.
Then, use the method of undetermined coefficients to find a particular solution of the problem for $y''-2y'+y=e^t$.
The general solution of the initial differential equation, will then be the general solution of the homogenous plus the particular solution you found.
You can find more information and examples about that method, here.
$mathbf{2}$ - Laplace Transformation :
This is a very fast and straight forward way to approach the problem, but it needs some fluent handling of Laplace Transformation techniques. Note that you can apply the Laplace Transformation without even needing initial conditions, simply stating them as constants.
Start of by applying the Laplace Transformation
$$mathcal{L}_tbig[f(t)big](s) = int_0^infty f(t)e^{-st}mathrm{d}t$$
to both sides of the given differential equation :
$$mathcal{L}_tbig[y'' - 2y' + y'] = mathcal{L}_t[e^t]$$
$$Leftrightarrow$$
$$(s-1)^2big[mathcal{L}_t[y(t)](s)big] - (s-2)y(0) - y'(0) = frac{1}{s-1}$$
$$Leftrightarrow$$
$$mathcal{L}_tbig[y(t)big](s) = frac{y(0)(s^2-3s+2) + y'(0)(s-1) + 1}{(s-1)^3}$$
$$=$$
$$mathcal{L}_tbig[y(t)big](s) = frac{1}{(s-1)^3} - frac{y(0)}{(s-1)^2} + frac{y(0)}{s-1} + frac{y'(0)}{(s-1)^2} $$
$$implies$$
$$y(t) = frac{1}{2}e^t(t^2+2c_1 - 2c_1t + 2c_2t) = frac{e^tt^2}{2} + c_1e^t - c_1e^tt + c_2e^tt$$
$mathbf{3}$ - Variation of Parameters :
You must repeat the step of solving the homogenous equation by finding that $lambda$s mentioned. Then by listing the basis solution as $y_{b_1} = e^t$ and $_{b_2} = e^tt$ you can use variation of parameters to find the final general solution by computing the Wronskian and finding the integrals :
$$v_1(t) = - int frac{f(t)y_{b_2}(t)}{W(t)}mathrm{d}t quad text{and} quad v_2(t) = int frac{f(t)y_{b_1}(t)}{W(t)}mathrm{d}t$$
You can find more information and examples about that method, here.
edited Dec 1 at 22:47
answered Dec 1 at 22:19
Rebellos
14.3k31244
14.3k31244
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
add a comment |
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
Not "the" particular solution, but "a" particular solution.
– Jean Marie
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
@JeanMarie Thanks.
– Rebellos
Dec 1 at 22:22
add a comment |
Both your attempts are in fact right but fail because the fundamental set of solutions for your second order ODE is given by exactly your both guesses for the particular solution. It is not hard to show by using the characteristic equation that the fundamental set of solutions is given by
$$y(t)=c_1e^t+c_2te^t$$
Therefore in this case it is not possible to obtain a particular solution in the standard way hence the inhomogeneous term is in fact part of the solution.
Howsoever you second try was near the actual solution. Since both terms $e^t$ and $te^t$ did not worked out hence there are part of the solution you could further consider $t^2e^t$ as the next try. It turns out that this yields to the solution $($see here$)$.
This leads to the conjecture that at least for a inhomogeneous term of the form $Ae^{bt}$ you just have to try the specific exponential with different powers of the variable infront as particular solution until it works out.
add a comment |
Both your attempts are in fact right but fail because the fundamental set of solutions for your second order ODE is given by exactly your both guesses for the particular solution. It is not hard to show by using the characteristic equation that the fundamental set of solutions is given by
$$y(t)=c_1e^t+c_2te^t$$
Therefore in this case it is not possible to obtain a particular solution in the standard way hence the inhomogeneous term is in fact part of the solution.
Howsoever you second try was near the actual solution. Since both terms $e^t$ and $te^t$ did not worked out hence there are part of the solution you could further consider $t^2e^t$ as the next try. It turns out that this yields to the solution $($see here$)$.
This leads to the conjecture that at least for a inhomogeneous term of the form $Ae^{bt}$ you just have to try the specific exponential with different powers of the variable infront as particular solution until it works out.
add a comment |
Both your attempts are in fact right but fail because the fundamental set of solutions for your second order ODE is given by exactly your both guesses for the particular solution. It is not hard to show by using the characteristic equation that the fundamental set of solutions is given by
$$y(t)=c_1e^t+c_2te^t$$
Therefore in this case it is not possible to obtain a particular solution in the standard way hence the inhomogeneous term is in fact part of the solution.
Howsoever you second try was near the actual solution. Since both terms $e^t$ and $te^t$ did not worked out hence there are part of the solution you could further consider $t^2e^t$ as the next try. It turns out that this yields to the solution $($see here$)$.
This leads to the conjecture that at least for a inhomogeneous term of the form $Ae^{bt}$ you just have to try the specific exponential with different powers of the variable infront as particular solution until it works out.
Both your attempts are in fact right but fail because the fundamental set of solutions for your second order ODE is given by exactly your both guesses for the particular solution. It is not hard to show by using the characteristic equation that the fundamental set of solutions is given by
$$y(t)=c_1e^t+c_2te^t$$
Therefore in this case it is not possible to obtain a particular solution in the standard way hence the inhomogeneous term is in fact part of the solution.
Howsoever you second try was near the actual solution. Since both terms $e^t$ and $te^t$ did not worked out hence there are part of the solution you could further consider $t^2e^t$ as the next try. It turns out that this yields to the solution $($see here$)$.
This leads to the conjecture that at least for a inhomogeneous term of the form $Ae^{bt}$ you just have to try the specific exponential with different powers of the variable infront as particular solution until it works out.
answered Dec 1 at 22:23
mrtaurho
3,1682930
3,1682930
add a comment |
add a comment |
I will answer to the question by considering the general $n$-th order constant coefficients linear ODE
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=flabel{1}tag{1}
$$
and the associated linear homogeneous equation
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=0label{a}tag{1'}
$$
where $fnotequiv 0$: the choice $a_n=1$ is only for formal simplification in the development and does not restrict the generality of the analysis.
There are basically two methods for finding a particular solution $y_p$ of the ODE eqref{1}:
by guessing, based on the solver's experience: this method will not analyzed here.
by choosing a solution $y_o$ of the associated homogeneous equation eqref{a} satisfying
$$
y_o(0)=y_o^{(1)}(0)=dots=y_o^{(n-2)}(0)quad y_o^{(n-1)}(0)=1label{2}tag{2}
$$
and forming the fundamental solution $mathscr{E}$ of eqref{1}
$$
mathscr{E}(t)=H(t)y_o(t),label{3}tag{3}
$$
where $H(t)$ is the Heaviside function. Then the sought for particular solution is
$$
y_p(t)=mathscr{E}ast f(t)=intlimits_{0}^ty_o(t-s)f(s)mathrm{d}s label{4}tag{4}
$$
Let's apply formula eqref{4} to the OP problem, before analyzing why it gives the sought for result. Since the characteristic equation is
$$
x^2-2x+1=0iff x=1 text{ with multiplicity 2}
$$
we have that a fundamental system of solutions of the homogeneous equation associated to the given one is
$$
y_1(t)=e^t,: y_2(t)=te^timplies y_o=y_2(t)
$$
since it is the only solution satisfying condition eqref{2}, i.e. $y_o(0)=0$ and $dfrac{mathrm{d}}{mathrm{d}t}y_o(0)=0$. Now we have that
$$
mathscr{E}(t)=H(t)y_o(t)=H(t)te^t
$$
and by applying eqref{4} we obtain
$$
begin{split}
y_p(t)=mathscr{E}ast exp(t)&=intlimits_{0}^{+infty}H(t-s)(t-s)e^{t-s}e^smathrm{d}s\
&=intlimits_{0}^ty_o(t-s)e^smathrm{d}s\
&=intlimits_{0}^t(t-s)e^{t}mathrm{d}s={t^2 over 2}e^t
end{split}
$$
which is the sought for particular solution of the ODE proposed as an example.
Why the methods works? Because of the properties of distributions ([1], §4.9, example 4.9.6 pp. 77-74 and §15.4, example 15.4.4): by using these properties and the condition eqref{2} we have
$$
begin{split}
frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)&=H(t)frac{mathrm{d}}{mathrm{d}t}y_o(t)\
frac{mathrm{d}^2}{mathrm{d}t^2}mathscr{E}(t)&=H(t)frac{mathrm{d}^2}{mathrm{d}t^2}y_o(t)\
&vdots\
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&=H(t)frac{mathrm{d}^n}{mathrm{d}t^n}y_o(t)+delta(t)
end{split}
$$
and thus
$$
begin{split}
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}mathscr{E}(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)+a_0mathscr{E}(t)\
=&H(t)Big[frac{mathrm{d}^n}{mathrm{d}t^n}y_0(t)+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}y_0(t)(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}y_0(t)+a_0y_0(t)Big]+delta(t)=delta(t)
end{split}
$$
The linearity of the differential operator and the properties of the Dirac $delta$-distribution do the rest.
Notes
The analyzed method is perhaps the simplest way to find a particular solution of the ODE proposed by the OP, and more generally of eqref{1}, because it requires only the knowledge of a complete systems of solutions of the associated equation eqref{2}. Indeed such a system of solutions is already required to solve any problem (Cauchy, boundary value, etc.) for eqref{1} and eqref{3}: the only further operation to do is calculating $y_o$ and the convolution integral eqref{3}.
The convolution integral used in eqref{3} is the standard one used in the operational calculus of one variable functions, i.e.
$$
mathscr{E}ast f(t)=intlimits_0^{+infty}mathscr{E}(t-s)f(s)mathrm{d}squad mathscr{E},fin L_mathrm{loc}^1(mathbb{R}_+)
$$
which can be deduced from the standard one by considering $H(t)f(t)$ instead of $f(t)$ as homogeneous term.The assumption $a_n=1$ does not restrict the generality of the above analysis because if we assume that we are dealing with an $n$-th order linear constant coefficient ordinary differential operator, we must assume $a_nneq 0$.
The distribution eqref{3} is called fundamental solution exactly because it can be used to construct the solution for every linear, constant coefficient non-homogeneous ODE.
[1] Vladimirov, V. S. (2002), Methods of the theory of generalized functions, Analytical Methods and Special Functions, 6, London–New York: Taylor & Francis, pp. XII+353, ISBN 0-415-27356-0, MR 2012831, Zbl 1078.46029.
add a comment |
I will answer to the question by considering the general $n$-th order constant coefficients linear ODE
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=flabel{1}tag{1}
$$
and the associated linear homogeneous equation
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=0label{a}tag{1'}
$$
where $fnotequiv 0$: the choice $a_n=1$ is only for formal simplification in the development and does not restrict the generality of the analysis.
There are basically two methods for finding a particular solution $y_p$ of the ODE eqref{1}:
by guessing, based on the solver's experience: this method will not analyzed here.
by choosing a solution $y_o$ of the associated homogeneous equation eqref{a} satisfying
$$
y_o(0)=y_o^{(1)}(0)=dots=y_o^{(n-2)}(0)quad y_o^{(n-1)}(0)=1label{2}tag{2}
$$
and forming the fundamental solution $mathscr{E}$ of eqref{1}
$$
mathscr{E}(t)=H(t)y_o(t),label{3}tag{3}
$$
where $H(t)$ is the Heaviside function. Then the sought for particular solution is
$$
y_p(t)=mathscr{E}ast f(t)=intlimits_{0}^ty_o(t-s)f(s)mathrm{d}s label{4}tag{4}
$$
Let's apply formula eqref{4} to the OP problem, before analyzing why it gives the sought for result. Since the characteristic equation is
$$
x^2-2x+1=0iff x=1 text{ with multiplicity 2}
$$
we have that a fundamental system of solutions of the homogeneous equation associated to the given one is
$$
y_1(t)=e^t,: y_2(t)=te^timplies y_o=y_2(t)
$$
since it is the only solution satisfying condition eqref{2}, i.e. $y_o(0)=0$ and $dfrac{mathrm{d}}{mathrm{d}t}y_o(0)=0$. Now we have that
$$
mathscr{E}(t)=H(t)y_o(t)=H(t)te^t
$$
and by applying eqref{4} we obtain
$$
begin{split}
y_p(t)=mathscr{E}ast exp(t)&=intlimits_{0}^{+infty}H(t-s)(t-s)e^{t-s}e^smathrm{d}s\
&=intlimits_{0}^ty_o(t-s)e^smathrm{d}s\
&=intlimits_{0}^t(t-s)e^{t}mathrm{d}s={t^2 over 2}e^t
end{split}
$$
which is the sought for particular solution of the ODE proposed as an example.
Why the methods works? Because of the properties of distributions ([1], §4.9, example 4.9.6 pp. 77-74 and §15.4, example 15.4.4): by using these properties and the condition eqref{2} we have
$$
begin{split}
frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)&=H(t)frac{mathrm{d}}{mathrm{d}t}y_o(t)\
frac{mathrm{d}^2}{mathrm{d}t^2}mathscr{E}(t)&=H(t)frac{mathrm{d}^2}{mathrm{d}t^2}y_o(t)\
&vdots\
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&=H(t)frac{mathrm{d}^n}{mathrm{d}t^n}y_o(t)+delta(t)
end{split}
$$
and thus
$$
begin{split}
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}mathscr{E}(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)+a_0mathscr{E}(t)\
=&H(t)Big[frac{mathrm{d}^n}{mathrm{d}t^n}y_0(t)+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}y_0(t)(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}y_0(t)+a_0y_0(t)Big]+delta(t)=delta(t)
end{split}
$$
The linearity of the differential operator and the properties of the Dirac $delta$-distribution do the rest.
Notes
The analyzed method is perhaps the simplest way to find a particular solution of the ODE proposed by the OP, and more generally of eqref{1}, because it requires only the knowledge of a complete systems of solutions of the associated equation eqref{2}. Indeed such a system of solutions is already required to solve any problem (Cauchy, boundary value, etc.) for eqref{1} and eqref{3}: the only further operation to do is calculating $y_o$ and the convolution integral eqref{3}.
The convolution integral used in eqref{3} is the standard one used in the operational calculus of one variable functions, i.e.
$$
mathscr{E}ast f(t)=intlimits_0^{+infty}mathscr{E}(t-s)f(s)mathrm{d}squad mathscr{E},fin L_mathrm{loc}^1(mathbb{R}_+)
$$
which can be deduced from the standard one by considering $H(t)f(t)$ instead of $f(t)$ as homogeneous term.The assumption $a_n=1$ does not restrict the generality of the above analysis because if we assume that we are dealing with an $n$-th order linear constant coefficient ordinary differential operator, we must assume $a_nneq 0$.
The distribution eqref{3} is called fundamental solution exactly because it can be used to construct the solution for every linear, constant coefficient non-homogeneous ODE.
[1] Vladimirov, V. S. (2002), Methods of the theory of generalized functions, Analytical Methods and Special Functions, 6, London–New York: Taylor & Francis, pp. XII+353, ISBN 0-415-27356-0, MR 2012831, Zbl 1078.46029.
add a comment |
I will answer to the question by considering the general $n$-th order constant coefficients linear ODE
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=flabel{1}tag{1}
$$
and the associated linear homogeneous equation
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=0label{a}tag{1'}
$$
where $fnotequiv 0$: the choice $a_n=1$ is only for formal simplification in the development and does not restrict the generality of the analysis.
There are basically two methods for finding a particular solution $y_p$ of the ODE eqref{1}:
by guessing, based on the solver's experience: this method will not analyzed here.
by choosing a solution $y_o$ of the associated homogeneous equation eqref{a} satisfying
$$
y_o(0)=y_o^{(1)}(0)=dots=y_o^{(n-2)}(0)quad y_o^{(n-1)}(0)=1label{2}tag{2}
$$
and forming the fundamental solution $mathscr{E}$ of eqref{1}
$$
mathscr{E}(t)=H(t)y_o(t),label{3}tag{3}
$$
where $H(t)$ is the Heaviside function. Then the sought for particular solution is
$$
y_p(t)=mathscr{E}ast f(t)=intlimits_{0}^ty_o(t-s)f(s)mathrm{d}s label{4}tag{4}
$$
Let's apply formula eqref{4} to the OP problem, before analyzing why it gives the sought for result. Since the characteristic equation is
$$
x^2-2x+1=0iff x=1 text{ with multiplicity 2}
$$
we have that a fundamental system of solutions of the homogeneous equation associated to the given one is
$$
y_1(t)=e^t,: y_2(t)=te^timplies y_o=y_2(t)
$$
since it is the only solution satisfying condition eqref{2}, i.e. $y_o(0)=0$ and $dfrac{mathrm{d}}{mathrm{d}t}y_o(0)=0$. Now we have that
$$
mathscr{E}(t)=H(t)y_o(t)=H(t)te^t
$$
and by applying eqref{4} we obtain
$$
begin{split}
y_p(t)=mathscr{E}ast exp(t)&=intlimits_{0}^{+infty}H(t-s)(t-s)e^{t-s}e^smathrm{d}s\
&=intlimits_{0}^ty_o(t-s)e^smathrm{d}s\
&=intlimits_{0}^t(t-s)e^{t}mathrm{d}s={t^2 over 2}e^t
end{split}
$$
which is the sought for particular solution of the ODE proposed as an example.
Why the methods works? Because of the properties of distributions ([1], §4.9, example 4.9.6 pp. 77-74 and §15.4, example 15.4.4): by using these properties and the condition eqref{2} we have
$$
begin{split}
frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)&=H(t)frac{mathrm{d}}{mathrm{d}t}y_o(t)\
frac{mathrm{d}^2}{mathrm{d}t^2}mathscr{E}(t)&=H(t)frac{mathrm{d}^2}{mathrm{d}t^2}y_o(t)\
&vdots\
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&=H(t)frac{mathrm{d}^n}{mathrm{d}t^n}y_o(t)+delta(t)
end{split}
$$
and thus
$$
begin{split}
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}mathscr{E}(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)+a_0mathscr{E}(t)\
=&H(t)Big[frac{mathrm{d}^n}{mathrm{d}t^n}y_0(t)+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}y_0(t)(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}y_0(t)+a_0y_0(t)Big]+delta(t)=delta(t)
end{split}
$$
The linearity of the differential operator and the properties of the Dirac $delta$-distribution do the rest.
Notes
The analyzed method is perhaps the simplest way to find a particular solution of the ODE proposed by the OP, and more generally of eqref{1}, because it requires only the knowledge of a complete systems of solutions of the associated equation eqref{2}. Indeed such a system of solutions is already required to solve any problem (Cauchy, boundary value, etc.) for eqref{1} and eqref{3}: the only further operation to do is calculating $y_o$ and the convolution integral eqref{3}.
The convolution integral used in eqref{3} is the standard one used in the operational calculus of one variable functions, i.e.
$$
mathscr{E}ast f(t)=intlimits_0^{+infty}mathscr{E}(t-s)f(s)mathrm{d}squad mathscr{E},fin L_mathrm{loc}^1(mathbb{R}_+)
$$
which can be deduced from the standard one by considering $H(t)f(t)$ instead of $f(t)$ as homogeneous term.The assumption $a_n=1$ does not restrict the generality of the above analysis because if we assume that we are dealing with an $n$-th order linear constant coefficient ordinary differential operator, we must assume $a_nneq 0$.
The distribution eqref{3} is called fundamental solution exactly because it can be used to construct the solution for every linear, constant coefficient non-homogeneous ODE.
[1] Vladimirov, V. S. (2002), Methods of the theory of generalized functions, Analytical Methods and Special Functions, 6, London–New York: Taylor & Francis, pp. XII+353, ISBN 0-415-27356-0, MR 2012831, Zbl 1078.46029.
I will answer to the question by considering the general $n$-th order constant coefficients linear ODE
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=flabel{1}tag{1}
$$
and the associated linear homogeneous equation
$$
frac{mathrm{d}^{n}y}{mathrm{d}t^{n}}+a_{n-1}frac{mathrm{d}^{n-1}y}{mathrm{d}t^{n-1}}+dots+a_1frac{mathrm{d}y}{mathrm{d}t}+a_0y=0label{a}tag{1'}
$$
where $fnotequiv 0$: the choice $a_n=1$ is only for formal simplification in the development and does not restrict the generality of the analysis.
There are basically two methods for finding a particular solution $y_p$ of the ODE eqref{1}:
by guessing, based on the solver's experience: this method will not analyzed here.
by choosing a solution $y_o$ of the associated homogeneous equation eqref{a} satisfying
$$
y_o(0)=y_o^{(1)}(0)=dots=y_o^{(n-2)}(0)quad y_o^{(n-1)}(0)=1label{2}tag{2}
$$
and forming the fundamental solution $mathscr{E}$ of eqref{1}
$$
mathscr{E}(t)=H(t)y_o(t),label{3}tag{3}
$$
where $H(t)$ is the Heaviside function. Then the sought for particular solution is
$$
y_p(t)=mathscr{E}ast f(t)=intlimits_{0}^ty_o(t-s)f(s)mathrm{d}s label{4}tag{4}
$$
Let's apply formula eqref{4} to the OP problem, before analyzing why it gives the sought for result. Since the characteristic equation is
$$
x^2-2x+1=0iff x=1 text{ with multiplicity 2}
$$
we have that a fundamental system of solutions of the homogeneous equation associated to the given one is
$$
y_1(t)=e^t,: y_2(t)=te^timplies y_o=y_2(t)
$$
since it is the only solution satisfying condition eqref{2}, i.e. $y_o(0)=0$ and $dfrac{mathrm{d}}{mathrm{d}t}y_o(0)=0$. Now we have that
$$
mathscr{E}(t)=H(t)y_o(t)=H(t)te^t
$$
and by applying eqref{4} we obtain
$$
begin{split}
y_p(t)=mathscr{E}ast exp(t)&=intlimits_{0}^{+infty}H(t-s)(t-s)e^{t-s}e^smathrm{d}s\
&=intlimits_{0}^ty_o(t-s)e^smathrm{d}s\
&=intlimits_{0}^t(t-s)e^{t}mathrm{d}s={t^2 over 2}e^t
end{split}
$$
which is the sought for particular solution of the ODE proposed as an example.
Why the methods works? Because of the properties of distributions ([1], §4.9, example 4.9.6 pp. 77-74 and §15.4, example 15.4.4): by using these properties and the condition eqref{2} we have
$$
begin{split}
frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)&=H(t)frac{mathrm{d}}{mathrm{d}t}y_o(t)\
frac{mathrm{d}^2}{mathrm{d}t^2}mathscr{E}(t)&=H(t)frac{mathrm{d}^2}{mathrm{d}t^2}y_o(t)\
&vdots\
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&=H(t)frac{mathrm{d}^n}{mathrm{d}t^n}y_o(t)+delta(t)
end{split}
$$
and thus
$$
begin{split}
frac{mathrm{d}^n}{mathrm{d}t^n}mathscr{E}(t)&+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}mathscr{E}(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}mathscr{E}(t)+a_0mathscr{E}(t)\
=&H(t)Big[frac{mathrm{d}^n}{mathrm{d}t^n}y_0(t)+a_{n-1}frac{mathrm{d}^{n-1}}{mathrm{d}t^{n-1}}y_0(t)(t)+dots+a_{1}frac{mathrm{d}}{mathrm{d}t}y_0(t)+a_0y_0(t)Big]+delta(t)=delta(t)
end{split}
$$
The linearity of the differential operator and the properties of the Dirac $delta$-distribution do the rest.
Notes
The analyzed method is perhaps the simplest way to find a particular solution of the ODE proposed by the OP, and more generally of eqref{1}, because it requires only the knowledge of a complete systems of solutions of the associated equation eqref{2}. Indeed such a system of solutions is already required to solve any problem (Cauchy, boundary value, etc.) for eqref{1} and eqref{3}: the only further operation to do is calculating $y_o$ and the convolution integral eqref{3}.
The convolution integral used in eqref{3} is the standard one used in the operational calculus of one variable functions, i.e.
$$
mathscr{E}ast f(t)=intlimits_0^{+infty}mathscr{E}(t-s)f(s)mathrm{d}squad mathscr{E},fin L_mathrm{loc}^1(mathbb{R}_+)
$$
which can be deduced from the standard one by considering $H(t)f(t)$ instead of $f(t)$ as homogeneous term.The assumption $a_n=1$ does not restrict the generality of the above analysis because if we assume that we are dealing with an $n$-th order linear constant coefficient ordinary differential operator, we must assume $a_nneq 0$.
The distribution eqref{3} is called fundamental solution exactly because it can be used to construct the solution for every linear, constant coefficient non-homogeneous ODE.
[1] Vladimirov, V. S. (2002), Methods of the theory of generalized functions, Analytical Methods and Special Functions, 6, London–New York: Taylor & Francis, pp. XII+353, ISBN 0-415-27356-0, MR 2012831, Zbl 1078.46029.
edited Dec 3 at 10:46
answered Dec 2 at 21:33
Daniele Tampieri
1,6331619
1,6331619
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3021903%2fparticular-solution-of-second-order-differential-equation%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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
2
I am surprised you have no initial conditions ; having them I would advise you to use Laplace Transform.
– Jean Marie
Dec 1 at 22:24
2
Since both $e^t$ and $te^t$ are solutions of the homogeneous equation, for the particular solution you must try $At^2e^t$. See, for example, DE texts like Coddington & Levinson.
– GEdgar
Dec 1 at 22:39
@JeanMarie - he could create generic initial conditions and obtain a solution with it.
– DavidG
Dec 2 at 11:59
@DavidG that's right, but I think this issue should be mentionned in such a question.
– Jean Marie
Dec 2 at 12:42
@JeanMarie - Agree 100%.
– DavidG
Dec 3 at 0:17