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

Homework Help: Convolution with exponentials

  1. Apr 7, 2008 #1
    How do you do the convolution of

    exp (x(n))*u(x(n)) and exp(x(n-1))*u(x(n-1))

    where u(x) is the unit step function. Thanks.
  2. jcsd
  3. Apr 7, 2008 #2
    Because convolution is commutative,
    [tex] y(t) = u(t) * e^{t} = e^{t} * u(t)[/tex]
    therefore, using the general result
    [tex] h(t)*u(t) = \int_{-\infty} ^{t} *e^{\tau} d\tau= e^{t} - e^{-\infty } = e^{t}[/tex]

    So, as an example given that

    [tex] h(t) = e^{-t}*u(t)[/tex]
    [tex] f(t) = e^{-2t}*u(t) [/tex]

    and [tex] y(t) = h(t)* f(t)[/tex] , determine [tex] y(1). [/tex]


    [tex]y(t) = \int_{-\infty} ^{\infty} h(\tau)f(t-\tau) d\tau = \int_{-\infty} ^{\infty} [e^{-\tau}u(\tau)] [e^{-2(t-\tau)} u(t - \tau)] d\tau [/tex]

    for any value of t, it fallows that
    [tex] y(1) = \int_{-\infty} ^{\infty} [e^{-\tau}u(\tau)] [e^{-2(1-\tau)} u(1 - \tau)] d\tau = \int_{0} ^{1} e^{-\tau}e^{-2(1-\tau) }d\tau = e^{-2}(e^1 - 1)=e^{-1}-e^{-2} \approx 0.233.[/tex]

    I hope this helps a little bit.

    [tex] Reference. Convolution examples from Kudeki and Munson [/tex]
  4. Apr 7, 2008 #3
    Wait, the [itex]x(n)[/itex] is inside the [itex]u(x)[/itex]? That makes it difficult to do, unless you tell us something more about what [itex]x(n)[/itex] looks like.

    Are you sure it's not supposed to be just [itex]u(n)[/itex]?
  5. Apr 7, 2008 #4
    You are right, the other approach will be maniopulate the equation with the properties of convolution.
  6. Feb 11, 2009 #5

    Nice Job!

    I am new here and I bring a similar question for my first post.....

    I want to know how to do a convolution where the two functions are:

    Ca(t) - arbitrary Input function
    (it actually represents the time activity curve of a CT contrast bolus injection in the blood)

    R(t) - a piecewise function defined as follows:

    R(t) = 1, 0<t<Tm
    and E*(exp)^(kt), t>Tm

    (this R(t) is called the Impulse Residue Function for the Johnson WIlson model for capillary tracer exchange)

    so therefore:

    Ca(t)*R(t) = (from 0 to Tm){Ca(t) convolved with 1} + (from Tm to t)E*{Ca(t) convolved with (*exp)^(-kt)}

    Can anyone shed some light on this please?!

    If anyone is curious the context of this convolution is for determining the representation of CT tissue attenuation in tracer kinetics modelling, considering a distributive parameter model. A background link for those interested is below.


  7. Feb 17, 2009 #6
    what is matlab code for convolution in z domain?
  8. Feb 17, 2009 #7
    Nevermind, I worked this thing out

    I don't want Matlab code I just wanted to do it by hand, an analytical solution

    Matlab is worthless for giving insight into a solution, but it will provide you with a numerical solution...if you don't understand how it works and why it works the solution is useless for understanding your problem
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook