MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: DSolve with DiracDelta


It appears that you need a boundary condition at a time other than zero.

Clear[soln];
soln[t_]=y[t]/.
    DSolve[{y'[t] + a y[t] == DiracDelta[t], 
          y[1] == Exp[-a]}, y[t], t][[1]]

UnitStep[t]/E^(a*t)


Bob Hanlon

---- Steffen Paul <steffen.paul at item.uni-bremen.de> wrote: 
> Hi
> I tried to solve
> DSolve[{y'[t] + \[Alpha] y[t] == DiracDelta[t], y[0] == 0}, y, t]
> 
> and got
> 
> -\[ExponentialE]^(-t \[Alpha]) (1 - HeavisideTheta[t])
> 
> which is zero for t >0.
> 
> The solution is correct but I expected somthing else:
> 
> exp( - alpha t) UnitStep(t)
> 
> which is zero for t <0 and which is also a solution.
> 
> In engineering, these solutions are called impulse responses.
> 
> The last solution is physically more usefull , because the system responds 
> after the excitation (DiracDelta).
> 
> 
> 
> How can I force Mathematica to give only solutions with nonzero values for t 
>  >0 ?
> 
> 
> 
> Regards,
> 
> Steffen
> 
> 
> 
> 
> 



  • Prev by Date: Re: Re: How to get sqrt(Year^2)===Year?
  • Next by Date: reading non-adjacent cells
  • Previous by thread: Re: DSolve with DiracDelta
  • Next by thread: Re: Simplify by Recurrence Relations 2