Integration, Marginalization

  • Thread starter touqra
  • Start date
  • #1
287
0
I'll like to know the probability density function for one of the x or y axis, given that there is an exponential decay of a material in two-dimensional space. So, that means I have to marginalize, say y and keep x, but I couldn't solve the integration. I even tried with Mathematica and Matlab. Mathematica couldn't solve it. Matlab gives a Bessel function when x == 1, but when x != 1, it couldn't solve it. Please help.

[tex] PDF(x) = \int e^{-r} \, dy = \int e^{-\sqrt{x^2+y^2}} \, dy [/tex]
 

Answers and Replies

  • #2
DrDu
Science Advisor
6,107
807
Something like ## K_0(|x|)##, where K is the MacDonald function.
so $$
f(x_0)=\int \exp(-r) dy|_{x=x_0}=\int dx dy \exp(-r) \delta(r \cos \phi -x_0)=$$
$$=\int dr d\phi r \exp(-r) \frac{1}{2\pi} \int dk \exp(ik (r \cos \phi -x))$$
Now we first integrate over ##\phi## using ##\int d\phi \exp(ikr\cos\phi)=2\pi J_0(kr)##:
$$
f(x)=\frac{1}{2\pi}\int dk \int dr r \exp(-r)2 \pi J_0(kr)\exp(ikx_0)
$$
The integral over r I found in Magnus Oberhettinger, Formeln und Saetze fuer die speziellen Funktionen der Physik, p33
##\int_0^\infty \exp(-at)J_\nu(bt)t^\nu dt=\frac{(2b)^\nu \Gamma(\nu+1/2)}{(a^2+b^2)^{\nu+1/2} \sqrt{\pi}}##
so that (setting x_0=x)
$$ f(x)= \int dk (k^2+1)^{-1/2} \exp(ikx)$$
The Fourier transform is standard and yields the MacDonald function.
 

Related Threads on Integration, Marginalization

  • Last Post
Replies
1
Views
3K
  • Last Post
Replies
4
Views
694
  • Last Post
Replies
11
Views
7K
  • Last Post
Replies
11
Views
1K
  • Last Post
Replies
6
Views
5K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
13
Views
5K
  • Last Post
Replies
6
Views
9K
Replies
3
Views
5K
  • Last Post
Replies
10
Views
8K
Top