MathGroup Archive 2004

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

Search the Archive

Re: need help with integration

  • To: mathgroup at smc.vnet.net
  • Subject: [mg52045] Re: need help with integration
  • From: "symbio" <symbio at has.com>
  • Date: Mon, 8 Nov 2004 03:13:33 -0500 (EST)
  • References: <cmhv6l$pqv$1@smc.vnet.net> <cmklju$jk2$1@smc.vnet.net>
  • Reply-to: "symbio" <symbio at sha.com>
  • Sender: owner-wri-mathgroup at wolfram.com

Also, please note, that FullSimplify with assumptions also does not work, 
see below:
\!\(x[t_] =
3*\(E\^\(-\(t\/2\)\)\) UnitStep[t] + DiracDelta[t + 3]\[IndentingNewLine]
h[t_] = UnitStep[t] - UnitStep[t - 2]\[IndentingNewLine]
Integrate[
x[\[Tau]]
h[t - \[Tau]], {\[Tau], \(-\[Infinity]\), \[Infinity]}]\
\[IndentingNewLine]
eq1 = FullSimplify[%,
t \[Element] Reals\ && \ \[Tau]\ \[Element] Reals]\[IndentingNewLine]
Plot[eq1, {t, \(-10\), 10}, PlotRange -> All]\[IndentingNewLine]
Integrate[
x[\[Tau]] h[t - \[Tau]], {\[Tau], 0, \[Infinity]}]\[IndentingNewLine]
eq2 = FullSimplify[%,
t \[Element] Reals\ && \ \[Tau]\ \[Element] Reals]\[IndentingNewLine]
Plot[eq2, {t, \(-10\), 10}, PlotRange -> All]\[IndentingNewLine]
Integrate[
x[\[Tau]] h[t - \[Tau]], {\[Tau], \(-\[Infinity]\),
0}]\[IndentingNewLine]
eq3 = FullSimplify[%,
t \[Element] Reals\ && \ \[Tau]\ \[Element] Reals]\[IndentingNewLine]
Plot[eq3, {t, \(-10\), 0}, PlotRange -> All]\)



"David W. Cantrell" <DWCantrell at sigmaxi.org> wrote in message 
news:cmklju$jk2$1 at smc.vnet.net...
> "symbio" <symbio at sha.com> wrote:
>>  Below I define two functions x[t] and h[t], then in eq1 I integrate
>> the integrand ( x[tau] h[t-tau] ) from -inf to +inf and get one set of
>> results and plots, then I integrate the same integrand as before but this
>> time in two steps, once from -inf to 0 and once from 0 to +inf, but the
>> results and plots from the integration performed in two steps are NOT the
>> same as the results and plots from integration in one step.  Is this a
>> bug and if so what's the work around?
>
> This looks essentially like what you asked here about two weeks ago in
> "Integration of UnitStep has bugs!? help!"
>
> <http://groups.google.com/groups?threadm=clneor%24obq%241%40smc.vnet.net>
>
> I would have answered your question in the previous thread had you not
> already gotten two good answers. Look at them carefully. I suspect that 
> the
> trouble you're having now is the same as you were having previously; if 
> so,
> then your question has already been answered, twice.
>
> David
>
>> Try this in the input cell
>>
>> In[19]:=
>> x[t_] = (3Exp[-.5t]UnitStep[t]) + DiracDelta[t + 3]
>> h[t_] = UnitStep[t] - UnitStep[t - 2]
>> eq1 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], -\[Infinity],
>> \[Infinity]}] Plot[eq1, {t, -10, 10}, PlotRange -> All]
>> eq2 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], 0, \[Infinity]}]
>> Plot[eq2, {t, -10, 10}, PlotRange -> All]
>> eq3 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], -\[Infinity], 0}]
>> Plot[eq3, {t, -10, 0}, PlotRange -> All]
>>
>> You will get the following output expressions (I couldn't paste the plots
>> here, but if you run the above input cells you should get the plots too,
>> then it will be more obvious what the problem is)
>>
>> Out[21]=
>> \!\(1\/2\ \((1 + \(6.`\ \@\((\(-2\) + t)\)\^2\)\/\(\(\(2.`\)\(\
>> \[InvisibleSpace]\)\) - 1.`\ t\) + \(6.`\ \@t\^2\)\/t + \@\((3 +
>> t)\)\^2\/\(3 \
>> + t\) - \(1 + t + \@\((1 + t)\)\^2\)\/\(1 + t\))\)\)
>>
>> Out[23]=
>> \!\(\((\(-6.`\) +
>> 16.30969097075427`\ \[ExponentialE]\^\(\(-0.5`\)\ t\))\)\ UnitStep[\
>> \(-2\) + t] + \((\(\(6.`\)\(\[InvisibleSpace]\)\) -
>> 6.`\ \[ExponentialE]\^\(\(-0.5`\)\ t\))\)\ UnitStep[t]\)
>>
>> Out[25]=
>> \!\(1\/2\ \((\(-\(\(1 +
>> t\)\/\@\((1 + t)\)\^2\)\) + \(3 + t\)\/\@\((3 + t)\)\^2)\)\)
> 


  • Prev by Date: Re: Re: Adding Vectors -- Newbie help please
  • Next by Date: Re: using the prime gaps to make a convergent series
  • Previous by thread: Re: Re: need help with integration
  • Next by thread: Re: MathGroup /: Descriptive headings