Computing the Integral of a Gaussian with a Cosine Factor in DTI Analysis

  • Thread starter Thread starter iamalexalright
  • Start date Start date
  • Tags Tags
    Gaussian Integral
Click For Summary

Homework Help Overview

The discussion revolves around computing a specific integral involving a cosine factor and a Gaussian function, relevant to diffusion tensor imaging (DTI). The integral in question is presented as part of a proof in a textbook, and the original poster expresses uncertainty about how to manipulate the cosine term within the integral.

Discussion Character

  • Exploratory, Mathematical reasoning, Problem interpretation

Approaches and Questions Raised

  • Participants explore the idea of rewriting the cosine as a complex exponential and completing the square in the exponent. There is also a discussion about using Cauchy's theorem to adjust the integration path.

Discussion Status

Some participants have shared their work and reasoning, indicating progress in understanding the manipulation of the integral. However, there is still uncertainty regarding the application of Cauchy's theorem, and the original poster has not confirmed a complete understanding or resolution of the problem.

Contextual Notes

The original poster is preparing for a presentation and notes that the textbook does not provide a detailed derivation of the integral, which adds to the challenge of their inquiry.

iamalexalright
Messages
157
Reaction score
0

Homework Statement


Looking through this proof in my DTI(diffusion tensor imaging) book and they don't show the work on how to compute this integral (and I'll need to know it for a presentation):


\int cos(\gamma G\delta x) * e^{-x^{2}/4D\Delta}dx


I haven't really gotten anywhere yet. I've looked at the proof of the integral of the gaussian but I don't know where/how to manipulate it with the cosine there.

I know the answer should be:


e^{-\gamma^{2}G^{2}\delta^{2}D\Delta}
 
Physics news on Phys.org
You can write the cosine as a complex exponential. You then have to complete the square in the exponential, the integral now is a Gaussian, but is along a path that is off the real axis. You can then use Cauchy's theorem to move the integration path back on the real axis.
 
I'll write out (most of) my work because I think I might have done something wrong:

Integral = \int (1/2)(e^{i\gamma G\delta x} + e^{-i\gamma G\delta x})e^{-x^{2}/4D\Delta}dx =

(1/2) \int (e^{-x^{2}/4D\Delta + i\gamma G\delta x} + (e^{-x^{2}/4D\Delta - i\gamma G\delta x}) dx

first term:
e^{(-1/4D\Delta)(x^{2} - 2iD\Delta \gamma G\delta x + 4i^{2}D^{2}\Delta^{2}\gamma^{2}G^{2}\delta^{2} - 4i^{2}D^{2}\Delta^{2}\gamma^{2}G^{2}\delta^{2})} =

e^{(-1/4D\Delta)((x - 2\Delta Di\gamma G\delta)^{2} + 4i^{2}D^{2}\Delta^{2}\gamma^{2}G^{2}\delta^{2})} =

e^{(-1/4D\Delta)(x - 2\Delta Di\gamma G\delta)^{2}}e^{D\Delta \gamma^{2}G^{2}\delta^{2}}

Okay, that's the first, second is similar and we should end up with:

(1/2)e^{D\Delta \gamma^{2}G^{2}\delta^{2}} \int e^{-(x - 2\Delta Di\gamma G\delta)^{2}} + e^{-(x + 2\Delta Di\gamma G\delta)^{2}}

Now, you mention Cauchy's theorem but I don't exactly know how to use it.
 
Nevermind, I got it.
 

Similar threads

  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 47 ·
2
Replies
47
Views
5K
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 105 ·
4
Replies
105
Views
11K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 13 ·
Replies
13
Views
6K
  • · Replies 4 ·
Replies
4
Views
2K