It is currently Sun Jun 16, 2019 12:26 pm

 All times are UTC [ DST ]

 Page 1 of 1 [ 2 posts ]
 Print view Previous topic | Next topic
Author Message
 Post subject: Multiple integralPosted: Mon Jun 19, 2017 8:59 pm

Joined: Sat Nov 07, 2015 6:12 pm
Posts: 841
Location: Larisa
Let $\langle \cdot, \cdot \rangle$ denote the usual inner product of $\mathbb{R}^m$. Evaluate the integral

$$\mathcal{M} = \int \limits_{\mathbb{R}^m} \exp \left( - ( \langle x, \mathcal{S} ^{-1} x \rangle )^a \right) \, {\rm d}x$$

where $\mathcal{S}$ is a positive symmetric $m \times m$ matrix and $a>0$.

_________________
Imagination is much more important than knowledge.

Top

 Post subject: Re: Multiple integralPosted: Thu Nov 30, 2017 10:00 pm

Joined: Sat Nov 14, 2015 6:32 am
Posts: 159
Location: Melbourne, Australia
Since $\mathcal{S}$ is a positive symmetric matrix , so is $\mathcal{S}^{-1}$. For a positive symmetric matrix $\mathcal{A}$ there exists an $\mathcal{R}$ positive symmetric matrix such that $\mathcal{A} = \mathcal{R}^2$. Applying this to $\mathcal{S}^{-1}$ our integral becomes

$\mathcal{M} = \int \limits_{\mathbb{R}^m} \exp \left ( - \left \| \mathcal{R} x \right \|^{2a} \right ) \, {\rm d}x$

where $\left \| \cdot \right \|$ is the Euclidean norm. Applying a change of variables we have that

$\mathcal{M} = \det \left ( \mathcal{R}^{-1} \right ) \int \limits_{\mathbb{R}^m} e^{-\left \| y \right \|^{2a}} \, {\rm d}y$

Since $\det \left ( \mathcal{R}^{-1} \right ) = \sqrt{\det \left ( \mathcal{S} \right )}$ then by converting to polar coordinates we have that

\begin{align*} \mathcal{M} &= \omega_m \sqrt{\det \left ( \mathcal{S} \right )} \int_{0}^{\infty} r^{m-1} e^{-r^{2a}} \, {\rm d}r \\ &= \frac{\omega_m}{m} \sqrt{\det \left ( \mathcal{S} \right )}\Gamma \left ( \frac{m}{2a} + 1 \right ) \end{align*}

Here $\omega_m$ denotes the surface area measure of the unit sphere and it is known to be

$\omega_m = \frac{2 \pi^{m/2}}{\Gamma \left ( \frac{m}{2} \right )}$

hence

$\mathcal{M}=\frac{\sqrt{\det (\mathcal{S}) }\pi^{m/2}\Gamma\left(\frac{m}{2a}\right)}{2^{a-1}\Gamma\left(\frac{m}{2}\right)}$

where $\Gamma$ denotes the Gamma Euler function for which it holds that

$\Gamma(x+1) = x \Gamma(x) \quad \text{forall} \quad x>0$

_________________
$\displaystyle \sum_{n=1}^{\infty}\frac{1}{n^s}= \prod_{p \; \text{prime}}\frac{1}{1-p^{-s}}$

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 2 posts ]

 All times are UTC [ DST ]

#### Mathimatikoi Online

Users browsing this forum: No registered users and 1 guest

 You cannot post new topics in this forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum

Search for:
 Jump to:  Select a forum ------------------ Algebra    Linear Algebra    Algebraic Structures    Homological Algebra Analysis    Real Analysis    Complex Analysis    Calculus    Multivariate Calculus    Functional Analysis    Measure and Integration Theory Geometry    Euclidean Geometry    Analytic Geometry    Projective Geometry, Solid Geometry    Differential Geometry Topology    General Topology    Algebraic Topology Category theory Algebraic Geometry Number theory Differential Equations    ODE    PDE Probability & Statistics Combinatorics General Mathematics Foundation Competitions Archives LaTeX    LaTeX & Mathjax    LaTeX code testings Meta