Solving ODE with Laplace transform and Existence and Uniqueness Theorem

  • #1
lys04
113
4
1729847743236.png

1729847806788.png

This is my answer to the ODE (I think it's correct) via Laplace transform, but I'm more concerned about whether or not my explanation for part b is correct or not?
Any help would be greatly appreciated!!
 
Physics news on Phys.org
  • #2
I observe your solution (a) is continuous at t=##\pi##.
 
  • #3
anuttarasammyak said:
I observe your solution (a) is continuous at t=π.
Oh yeah you're right, what might be a good explanation for (b) then?
 
  • #4
I graphed it and it seems like it’s twice differentiable as well
 
  • #5
I would expect [itex]y[/itex] and [itex]y'[/itex] to be continuous, and [itex]y''[/itex] to have a jump discontinuity.
 
  • #6
pasmith said:
I would expect [itex]y[/itex] and [itex]y'[/itex] to be continuous, and [itex]y''[/itex] to have a jump discontinuity.
Why?
 
  • #7
Because
If y has a jump, derivatives y' and y" would have "hyper" jump which does not coincide with RHS of the ODE.
If y' has a jump, derivative y" would have "hyper" jump which does not coincide with RHS of the ODE.
y" should have a jump as RHS of the ODE has. "hyper" jump here is actually delta function and its derivative.
 
Last edited:
  • #8
anuttarasammyak said:
If y has a jump, derivatives y' and y" would have "hyper" jump which does not coincide with RHS of the ODE.
If y' has a jump, derivative y" would have "hyper" jump which does not coincide with RHS of the ODE.
y" should have a jump as RHS of the ODE has. "hyper" jump here is actually delta function and its derivative.
By a jump do you mean a discontinuity? A heaviside function is discontinuous though?
 
  • #9
lys04 said:
By a jump do you mean a discontinuity? A heaviside function is discontinuous though?
Yes. Is there anything you feel uneasy ?
 
Last edited:
  • #10
anuttarasammyak said:
Yes. Is there anything you feel uneasy ?
Yeah. By y do you mean my solution found through Laplace transform? But then you speak of the RHS, is that referring to the Heaviside function?
So if you’re talking about the Heaviside function then I’m not sure what you mean by y’ and y’’ Would have hyper jumps which do not coincide with the RHS of the ODE because y (Heaviside does have a jump)?
Also what does a hyper jump mean?
 
  • #11
OK. First could you calculate second derivative of your solution and show it to us?
 
  • #12
anuttarasammyak said:
OK. First could you calculate second derivative of your solution and show it to us?
Yeah sure I could do that but the question says without doing any explicit calculations?
 
  • #13
Okay, the for $$t<=pi$$,

the first derivative is $$e^{t}sin(3t)+3e^{t}cos(3t)$$ and the second derivative is $$-8e^{t}sin(3t)+6e^{t}cos(3t)$$

and for $t>pi$,
the first derivative is $$e^{t}sin(3t)+3e^{t}cos(3t)+10e^{t-\pi}sin(3(t-\pi))$$ and the second derivative is $$-8e^{t}sin(3t)+6e^{t}cos(3t)+10e^{t-\pi}sin(3(t-\pi))+30e^{t-\pi}cos(3(t-\pi))$$
 
Last edited:
  • #14
Thanks for your effort. Say your sokution is y
[tex]y^"=-8e^{t}sin(3t)+6e^{t}cos(3t)=8e^\pi[/tex] at t=##\pi-0##
[tex]y^"=-8e^{t}sin(3t)+6e^{t}cos(3t)+10e^{t-\pi}sin(3(t-\pi))+30e^{t-\pi}cos(3(t-\pi))=8e^\pi+30[/tex] at t=##\pi+0##. We see discotinuity of gap 30 at at t=##\pi## in y".

You have already found that y is continuas at t=##\pi##. You can easily check that y' is continuous too there.

This is an example of post #7.
 
Last edited:
  • #15
anuttarasammyak said:
Thanks for your effort. Say your sokution is y
[tex]y^"=-8e^{t}sin(3t)+6e^{t}cos(3t)=8e^\pi[/tex] at t=##\pi-0##
[tex]y^"=-8e^{t}sin(3t)+6e^{t}cos(3t)+10e^{t-\pi}sin(3(t-\pi))+30e^{t-\pi}cos(3(t-\pi))=8e^\pi+30[/tex] at t=##\pi+0##. We see discotinuity of gap 30 at at t=##\pi## in y".

You have already found that y is continuas at t=##\pi##. You can easily check that y' is continuous too there.

This is an example of post #7.
Hm yeah, but how can I tell all of this without doing explicit calculations?
 
  • #16
lys04 said:
Why?

Because the highest derivative is the most badly behaved. This accords with our physical intuition, and guarantees that the lower order derivatives will actually exist.

The way of solving this equation without transforms is as a particular integral (which is zero for [itex]t \leq \pi[/itex]) plus a linear combination of complementary functions. The coefficients of the complementary functions in [itex][0,\pi][/itex] are fixed by the conditions at zero; the cofficients in [itex](\pi,\infty)[/itex] are fixed by requiring [itex]y[/itex] and [itex]y'[/itex] to be continuous at [itex]\pi[/itex].

This is also the solution obtained by the Laplace transform method, if you use the result [tex]
\mathcal{L}(y') = p\mathcal{L}(y) - y(0).[/tex] If [itex]y[/itex] were discontinuous at [itex]\pi[/itex], you would instead have [tex]
\mathcal{L}(y') = p\mathcal{L}(y) + (y(\pi^{-}) - y(\pi^{+}))e^{-p\pi} - y(0)[/tex] with a similar result for [itex]\mathcal{L}(y'')[/itex].
 
  • Like
Likes anuttarasammyak
  • #17
ODE is
[tex]y^"-2y'+10y=30H(t-\pi)[/tex]
H is in either y,y' or y". If H is in y or y',
[tex]H'(x)=\delta (x)[/tex]
would appear from y' or y" which is calculated from such y or y'. But we do not find it in ODE.
This is just another rephrasing of post #7.
 
Last edited:
  • #18
pasmith said:
This is also the solution obtained by the Laplace transform method, if you use the result L(y′)=pL(y)−y(0). If y were discontinuous at π, you would instead have L(y′)=pL(y)+(y(π−)−y(π+))e−pπ−y(0) with a similar result for L(y″).
I have been a lazy learner of Laplace transform so let me know it more.
Laplace transform of ODE in post #17 with the initial condition is
[tex](p^2-2p+10) \mathcal{L}(y)=3+30\frac{e^{-\pi p}}{p}[/tex]
##\mathcal{L}(y)## is written as sum of
[tex]\frac{1}{p-1-3i}, \frac{1}{p-1+3i}, \frac{1}{p}[/tex]
terms. Is it a right track to solve this ODE ?
 
  • #19
Once you have [tex]
\mathcal{L}(y) = \frac{3}{(p- 1)^2 + 3^2}\left(1 + \frac{10e^{-\pi p}}{p}\right)[/tex] you can look these up in transform tables.
 
  • Like
Likes anuttarasammyak
  • #20
Hello. It is a second order linear ordinary differential equation with constant coefficients and initial values.

[tex]y(t) = y_{h}(t) + y_{p}(t) [/tex]
In general one may ask themselves what happened to the transient (general, homogenous) and steady state response (particular, inhomogenous)?
[tex] \textbf{Transient Response: Solution Families/Categories Of Finite, Known Cases Of Solutions} [/tex]
The transient response is the homogenous solution. The homogenous solution occurs for:
- Distinct real roots,
- Repeated real roots,
- Complex roots
[tex] \textbf{Steady State response: inhomogenous solution/forcing function } [/tex]
The transient response is the homogenous solution. The homogenous solution occurs for:
The steady state response is the inhomogeneous solution.

The inhomogeneous solution requires a forcing function. Here, the forcing function is a unit step/Heaviside delayed by pi units.

[tex] h(t - \pi) [/tex]
In general the unit step function scaled and shifted is defined as:
[tex]
K h(t - T_{0}) =
\begin{cases}
0 & \text{for } t < T_{0} \\
K & \text{for } t \geq T_{0}
\end{cases}

[/tex]
Whose derivative is an equally scaled and shifted delta pulse:
[tex]
\frac{d}{dt} (K h(t - T_{0})) = K \delta(t - T_{0})

[/tex]
The derivative of the dirac delta (an term included within and of the second derivative of y(t)) is a distribution defined over a test function phi:
[tex]
\delta'(t) = \int_{-\infty}^{\infty} \delta'(t) \phi(t) \, dt = -\phi'(0)

[/tex]


For O.D.E, the transient response is the homogenous solution to the equation.
The steady state response is directly related to the forcing function.

it is fair to say that the transient response of a second-order linear differential equation is both differentiable everywhere, continuous in all cases of roots:
- real and distinct,
-real and repeated,
- and complex

Hence y(t) would be discontinuous (jump discontinuity)...and the nature of the y'(t) and y"(t) follows from the jump discontinuity of the step function, a term in y(t)...
....

Discontinuities in y(t) occur at the point where the unit step/heaviside function is applied

So for the solution for this family/category of these equations as the behavior of the homogenous parts are well defined across all finite and known cases, the discontinuity behavior arises from the forcing function, inhomogeneous solution.
 
Last edited:

Similar threads

Replies
17
Views
2K
Replies
1
Views
1K
Replies
5
Views
2K
Replies
10
Views
3K
Replies
1
Views
2K
Replies
1
Views
3K
Back
Top