Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Taking the first derivative of a function with improper integral

  1. Mar 22, 2010 #1

    I have a function:

    u(t,\tau)=\frac{1}{\pi}\int_{0}^{\infty}\! G(\omega)\, d\omega

    G(\omega)=4\sqrt{\pi}\frac{\omega^{2}}{\omega_{0}^{3}}\mbox{exp}\left(-\frac{\omega^{2}}{\omega_{0}^{2}}\right)\mbox{cos\left(\omega t-\left(\frac{\omega}{\omega_{0}}\right)^{-\gamma}\omega\tau\right)\mbox{exp}\left(-\frac{1}{2Q}\left(\frac{\omega}{\omega_{0}}\right)^{-\gamma}\omega t\right)}

    Now what I would like to do is to take the first time derivative of [tex]u(t,\tau)[/tex] to obtain the roots of [tex]\partial u(t,\tau) / \partial t = 0[/tex], where [tex]\partial u(t,\tau) / \partial t[/tex] is the derivative with respect to [tex]t[/tex].

    How would I get started? I think that I need to somehow get rid of the improper integral so that I can then take the first derivative. I've noticed that for [tex]\omega \rightarrow \infty[/tex], [tex]G(\omega) \rightarrow 0[/tex].
  2. jcsd
  3. Mar 22, 2010 #2
    Once you've demonstrated convergence of the improper integral, you could probably differentiate under the integral sign, i.e., switch the order of the differentiation and integration, here. I don't remember what the precise conditions are, but I do know that for distributions, they are pretty general. If I recall correctly, it's OK to differentiate under the integral sign to determine if the resulting expression is a continuous linear functional because then you can conclude that the original expression was indeed a distribution and this would justify the interchange.

    Be warned that I was exposed to a little bit of distribution theory before I learned anything beyond basic analysis, so perhaps someone more knowledgeable on the topic could give better/correct suggestions.
  4. Mar 22, 2010 #3
    That sounds good, snipez90 - I'll have to investigate. Thank you very much for your response.

    The only thing that is slightly perplexing is being able to show convergence of [tex]G(\omega)[/tex]. How should I begin?
  5. Mar 22, 2010 #4
    Well, for a start I think that

    4\sqrt{\pi}\frac{\omega^{2}}{\omega_{0}^{3}} \rightarrow \infty

    as [tex]\omega \rightarrow \infty [/tex] and

    \mbox{cos}\left(\omega t-\left(\frac{\omega}{\omega_{0}}\right)^{-\gamma}\omega\tau\right)

    will oscillate as [tex]\omega \rightarrow \infty [/tex] since we are dealing with a cosine function, and

    \mbox{exp}\left(-\frac{1}{2Q}\left(\frac{\omega}{\omega_{0}}\right)^{-\gamma}\omega t\right) \rightarrow 0

    as [tex]\omega \rightarrow \infty [/tex], for choices of [tex]\gamma, Q, \omega_0, t[/tex] particular to my problem. This may suggest that [tex]G(\omega) \rightarrow 0[/tex] as [tex]\omega \rightarrow \infty [/tex].

    Now how do I analytically perform the integration required for [tex]G(\omega)[/tex] so that I can find [tex]u(t,\tau)[/tex]? I understand that the integration could be done before or after the differentiation with respect to time.
    Last edited: Mar 22, 2010
  6. Mar 23, 2010 #5
    Do you have the specific bounds on say, gamma? For instance, as [itex]\omega \rightarrow \infty[/itex] the argument within the exponential goes to zero if [itex]\gamma > 1.[/itex]
  7. Mar 23, 2010 #6
    Thank you for your response, snipez90.

    Yes, often [tex]\gamma > 1[/tex], but also I also think that there's a possibility of [tex]\gamma \approx \frac{1}{\pi Q}[/tex], where [tex]Q[/tex] is a real number.
  8. Mar 23, 2010 #7
    Moreover, [tex]25 \leq Q \leq 400[/tex].
  9. Mar 23, 2010 #8
    All right, well I guess if [itex]\gamma \leq 1[/itex] it doesn't really matter since you'll still get convergence in that last exponential product of the integrand. Moreover, since exponential growth dominates polynomial growth,

    [tex]\frac{\omega^{2}}{\omega_{0}^ {3}}\cdot\mbox{exp}\left(-\frac{\omega^{2}}{\omega_{0}^{2}}\right) \rightarrow 0 \hspace{2mm} \mbox{as} \hspace{2mm} \omega \rightarrow \infty,[/tex]

    so convergence of [itex]G(\omega)[/itex] certainly doesn't seem like a problem.
  10. Mar 23, 2010 #9
    Thank you once again for your response!

    I agree that convergence of [tex]G(\omega)[/tex] doesn't seem to be a problem. Now I wonder if it would be possible to eliminate the improper integral before or after the differentiation.

    Although the improper integral can be dealt with numerically, could I get rid of it analytically? Or perhaps my best bet would be to use a series expansion to approximate the improper integral?
  11. Mar 23, 2010 #10
    If that integral can be dealt with analytically, I don't really see where to begin. I mean omega appears everywhere and we have a product involving a quadratic, a cosine, and two exponential functions. Ignoring what I said earlier about distributions for now, I think differentiating under the integral sign (i.e. differentiating the integrand first) with respect to the variable of your choice (I guess here it's t) can be justified if you ensure that the partial derivative of the integrand with respect to t exists and if you can find a dominating function H such that

    [tex]\left|\frac{\partial G}{\partial t}\right| \leq H.[/tex]

    (I left out some conditions regarding measurability/lebesgue integrability of the functions appearing in the integrand because everything making up the integrand is continuous)
    These certainly aren't the most general conditions, but I have a feeling it suffices for this purpose. Out of curiosity, what exactly is the context of this calculation?
    Last edited: Mar 23, 2010
  12. Mar 23, 2010 #11
    Yeah, that kind of confirms my own suspicions that it's not possible to really deal with the improper integral using analysis. Like you say, the best way to proceed would be to find the partial derivative of the integrand.

    [tex]u(t, \tau)[/tex] models an acoustic wavelet which has propagated through a lossy material. The speed and attenuation of the wavelet are dependent on the angular frequency [tex]\omega[/tex] and the [tex]Q[/tex] factor, which governs the attenuation.

    By taking the first time derivative of [tex]u(t, \tau)[/tex], and finding the roots of [tex]\partial u(t,\tau) / \partial t = 0[/tex], I would like to find the width of the wavelet, which is the total time between the two side-lobes.

    Thanks snipez90.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook