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

Exponentially Modified Gaussian

  1. Sep 15, 2010 #1
    I am fitting a detector response using Matlab. I have been asked to fit the spectrum using a Exponentially Modified Gaussian. This is as follows.

    [ tex ] a^x_n [ /tex ]

    h(t) = [ tex ] \frac{A}{\sqrt{2 \pi} \sigma} e^{-\frac{(t-t_{R})^{2}}{2\sigma^{2}} [ /tex ]

    f(t) = [ tex ] \frac{1}{\tau}e^{-\frac{t}{tau}}[ /tex ]

    Using the convolution y(t) = [ tex ] \int^{\infinity}_{0} h(t^{\prime})f(t-t^{prime})dt^{\prime}[ /tex ]

    This gives y(t) = [ tex ] \frac{A}{\tau \sigma \sqrt{2 \pi}} \int^{\infinity}_{0}} e^{-{\frac{(t-t_{R} - t^{\prime})^{2}}{2\sigma^{2}}e^{\frac{-t^{\prime}}{tau}} [ /tex ]

    This is the convolution of h(t)*f(t) to give y(t)

    The answer is

    y(t) = [ tex ] \frac{A}{2 \tau}[1-erf(\frac{\sigma}{\sqrt{2}\tau} - \frac{t-t_{R}}{\sqrt{2}\sigma]e^{\frac{\sigma^{2)}{\sqrt{2}\tau} - \frac{t-t_{R}}{\tau}}[ /tex ]

    I have been unable to get the answer. Has anyone ever worked this out to get y(t)? Also this is my first post and I am not sure how to get the latex to appear?
    Last edited: Sep 15, 2010
  2. jcsd
  3. Sep 15, 2010 #2
    You must write your equations between [tex][tex][/tex][/tex]
    Last edited: Sep 15, 2010
  4. Sep 15, 2010 #3


    User Avatar
    Science Advisor

    You need to prescribe the limits of integration on t'.
  5. Sep 15, 2010 #4
    Yes...but I believe one must complete the square for the terms in the exponent. This gives you the z for the erf(z). However, I am unable to complete the square correctly. I am close and still working on it.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook