Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2007
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2007

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

Search the Archive

Re: DSolve with DiracDelta

  • To: mathgroup at smc.vnet.net
  • Subject: [mg75803] Re: DSolve with DiracDelta
  • From: dimitris <dimmechan at yahoo.com>
  • Date: Fri, 11 May 2007 05:20:35 -0400 (EDT)
  • References: <f1unqt$7si$1@smc.vnet.net>

I guess you use the new mathematica 6?

How about 5.2?

In[109]:=
DSolve[{Derivative[1][y][t] + o*y[t] == DiracDelta[t], y[0] == 0}, y,
t]

% /. o -> 7/10

Show[Block[{$DisplayFunction = Identity}, (Plot[y[t] /. %[[1]], {t,
#1[[1]], #1[[2]]}, Axes -> False, Frame -> True] & ) /@
Partition[Range[-2, 2], 2, 1]]]

Out[109]=
{{y -> Function[{t}, (-1 + UnitStep[t])/E^(o*t)]}}

Out[110]=
{{y -> Function[{t}, (-1 + UnitStep[t])/E^((7*t)/10)]}}

Out[111]=
Graphics[]

Dimitris


=CF/=C7 Steffen Paul =DD=E3=F1=E1=F8=E5:
> 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: DSolve with DiracDelta
  • Next by Date: Re: Making sense of Directive[ ]
  • Previous by thread: Re: DSolve with DiracDelta
  • Next by thread: Re: DSolve with DiracDelta