I would like to solve the following integral that is a variation of this one (Integral involving Modified Bessel Function of the First Kind).
Namely, I have:
$$\frac{1}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty} \, x \, e^{(-\frac{x}{\alpha}-\frac{1}{2w^2}(x-\hat{x})^2)} \, I_0\left(\frac{x}{\beta}\right)\,dx .$$
By using the series representation for $I_0(x)$, I obtain the Tricomi hypergeometric function and I don't know how to calculate the corresponding series.
Any suggestion is highly appreciated.
Best Answer
One can try to do the following derivations. $$\mathrm{Int}=\frac{1}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty} \, x \, e^{(-\frac{x}{\alpha}-\frac{1}{2w^2}(x-\hat{x})^2)} \, I_0\left(\frac{x}{\beta}\right)\,dx$$ You can simplify the power of the exponent: $$-\left(\frac{x}{\alpha}+\frac{1}{2w^2}(x-\hat{x})^2\right)=-\left(\frac{(x-\mu)^2}{2w^2}+\gamma \right),$$ where $\gamma=\left(\frac{\hat{x}}{\alpha}-\frac{w^2}{2\alpha^2}\right)$, $\mu=\hat{x}-\frac{w^2}{\alpha^2}$. So: $$\mathrm{Int}=\frac{1}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty} \, x \, e^{(-\frac{x}{\alpha}-\frac{1}{2w^2}(x-\hat{x})^2)} \, I_0\left(\frac{x}{\beta}\right)\,dx=\frac{e^{-\gamma}}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty} \, x \, e^{-\frac{(x-\mu)^2}{2w^2}} \, I_0\left(\frac{x}{\beta}\right)\,dx$$ Then you can change the variable of integration $y=x-\mu$: $$\mathrm{Int}=\frac{e^{-\gamma}}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty} \, (y+\mu) \, e^{-\frac{y^2}{2w^2}} \, I_0\left(\frac{y+\mu}{\beta}\right)\,dx$$ and make use of the Neumann’s addition theorem: $$\mathop{I_{{\nu}}}\nolimits\!\left(u\pm v\right)=\sum _{{k=-\infty}}^{\infty}(\pm 1)^{k}\mathop{I_{{\nu+k}}}\nolimits\!\left(u\right)\mathop{I_{{k}}}\nolimits\!\left(v\right)$$ setting $\nu=0$, assuming that $y,\mu,\beta\in\mathbb{R}$ and using connection formulas one will get: $$\mathop{I_{0}}\nolimits\!\left(\frac{y+\mu}{\beta}\right)=\sum _{{k=-\infty}}^{\infty}\mathop{I_{k}}\nolimits\!\left(\frac{y}{\beta}\right)\mathop{I_{{k}}}\nolimits\!\left(\frac{\mu}{\beta}\right)=\mathop{I_{0}}\nolimits\!\left(\frac{y}{\beta}\right)\mathop{I_{{0}}}\nolimits\!\left(\frac{\mu}{\beta}\right)+2\sum _{{k=1}}^{\infty}\mathop{I_{k}}\nolimits\!\left(\frac{y}{\beta}\right)\mathop{I_{{k}}}\nolimits\!\left(\frac{\mu}{\beta}\right)$$
$$\mathrm{Int}=\! \frac{\mathop{I_0}\nolimits\!\left(\frac{\mu}{\beta}\right)e^{-\gamma}}{\sqrt{2\pi w^2}}\!\int_{-\infty}^{+\infty} \, (y+\mu)\! e^{-\frac{y^2}{2w^2}}\! \mathop{I_{0}}\nolimits\!\left(\frac{y}{\beta}\right)\!dx\! +\!\frac{2 e^{-\gamma}}{\sqrt{2\pi w^2}}\int_{-\infty}^{+\infty}\! (y+\mu) \, e^{-\frac{y^2}{2w^2}} \! \sum _{k=1}^{\infty}\!\mathop{I_k}\!\!\left(\frac{y}{\beta}\right)\!\!\mathop{I_k}\nolimits\!\left(\frac{\mu}{\beta}\right) \,dx$$ Or reorganising terms: $\mathop{I_0}\nolimits\!\left(\frac{\mu}{\beta}\right)$
$$\mathrm{Int}=\! \frac{e^{-\gamma}}{\sqrt{2\pi w^2}}\!\left(\mathop{I_0}\nolimits\!\left(\frac{\mu}{\beta}\right)\int_{-\infty}^{+\infty} \, (y\!+\!\mu)\! e^{-\frac{y^2}{2w^2}}\! \mathop{I_{0}}\nolimits\!\left(\frac{y}{\beta}\right)\!dx\! +\!2\!\sum _{k=1}^{\infty}\!\mathop{I_k}\!\left(\frac{\mu}{\beta}\right)\!\int_{-\infty}^{+\infty}\! (y\!+\!\mu) \, e^{-\frac{y^2}{2w^2}} \! \!\mathop{I_k}\!\left(\frac{y}{\beta}\right)\!\!\,dx \right)$$ Let's set $$\mathrm{Int}_k(\mu,\alpha,\beta)\!=\!\int_{-\infty}^{+\infty} \, (y\!+\!\mu) e^{-\frac{y^2}{2w^2}}\! \mathop{I_{k}}\nolimits\!\left(\frac{y}{\beta}\right)\!dx$$ $$\mathrm{Int}_0(\mu,\alpha,\beta)\!=\!\int_{-\infty}^{+\infty} \, (y\!+\!\mu) e^{-\frac{y^2}{2w^2}}\! \mathop{I_{0}}\nolimits\!\left(\frac{y}{\beta}\right)\!dx$$ Here one can play with the evenness/unevenness of the functions under integrals (as I did at first), or just use CAS and get:
$$\mathrm{Int}_k(\mu,\alpha,\beta)\!=\!\sqrt{\frac{\pi \alpha ^2}{8 \beta ^2}}\! e^{\frac{\alpha ^2}{4 \beta ^2}} \left(\!\!2 \beta \left(\!\!(-1)^k+1\!\right)\! \mu I_{\frac{k}{2}}\!\left(\!\!\frac{\alpha ^2}{4 \beta ^2}\!\!\right)+\alpha ^2 \left(\!1-(-1)^k\!\right)\!\! \left(\!\!I_{\frac{k-1}{2}}\!\!\left(\!\!\frac{\alpha ^2}{4 \beta ^2}\!\!\right)+I_{\frac{k+1}{2}}\!\left(\!\!\frac{\alpha ^2}{4 \beta ^2}\!\!\right)\!\!\right)\!\!\right)$$ $$\mathrm{Int}_0(\mu,\alpha,\beta)\!=\! \sqrt{2 \pi } \alpha \mu e^{\frac{\alpha ^2}{4 \beta ^2}} I_0\left(\frac{\alpha ^2}{4 \beta ^2}\right)$$ So the original integral can be represented in the following way: $$\mathrm{Int}(\mu,\alpha,\beta)\!=\! \frac{e^{-\gamma}}{\sqrt{2\pi w^2}}\left(\mathrm{Int}_0(\mu,\alpha,\beta)+ 2\!\sum _{k=1}^{\infty}\mathop{I_k}\left(\frac{\mu}{\beta}\right) \mathrm{Int}_k(\mu,\alpha,\beta) \right)$$ And I am not shure that there is any chance to simplify it futher (meaning finding the sum of the series).