Solving DE Homework: g'(t) + ag(t) = \delta(t-\xi)

  • Thread starter Thread starter squenshl
  • Start date Start date
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 2K views
squenshl
Messages
468
Reaction score
4

Homework Statement


Solve, sketch & check it's accuracy by operating on a test function of g'(t) + ag(t) = [tex]\delta[/tex](t-[tex]\xi[/tex]), g(t) = 0, t < [tex]\xi[/tex]

Homework Equations





The Attempt at a Solution


I have solved it getting g(t) = c1e-at + ea([tex]\xi[/tex]-t)H(t-[tex]\xi[/tex]) where H(x) is the Heaviside function. How do I sketch this & operate it on a test function.
 
Physics news on Phys.org
squenshl said:

Homework Statement


Solve, sketch & check it's accuracy by operating on a test function of g'(t) + ag(t) = [tex]\delta[/tex](t-[tex]\xi[/tex]), g(t) = 0, t < [tex]\xi[/tex]

Homework Equations



The Attempt at a Solution


I have solved it getting g(t) = c1e-at + ea([tex]\xi[/tex]-t)H(t-[tex]\xi[/tex]) where H(x) is the Heaviside function. How do I sketch this & operate it on a test function.

You can plot it manually if you wish. Just choose some values for a and e and c. But why not, if you're studying DEs, learn how to work with them in Mathematica. Here's the code I would use to "sketch" it and use a "test" function with the solution. Try and figure out what I'm doing if you're interested:

Code:
e = 5;
a = 1;
mysol = y /. DSolve[y'[t] + a y[t] == DiracDelta[t - e], y, t] // First
Plot[mysol[t] /. C[1] -> 1, {t, 0, 5}]
 
I got the sketch, so do I investigate a test function in the usual way.
 
squenshl said:
I got the sketch, so do I investigate a test function in the usual way.

Hi. The act of solving it numerically implicitly generates a particular (test) solution to the DE. I mean you choose some values for a, e and the arbitrary integration constant, C[1], in my code, then run the numeric integrator NDSolve. That "creates" a numeric function which satisfies the DE.