Is there an analytical solution for the following inverse Fourier transform?
$$f(x)=\frac{1}{\sqrt{2\pi}}\,\mathrm{j}^n \int_{-\infty}^\infty\frac{1}{\sqrt{k}\,(k^4-\lambda^4)}\mathrm{J}_{n+\frac{1}{2}}(k)\, \mathrm{e}^{-\mathrm{j}\,k\,x}\,\mathrm{d}k$$
with $n \in \mathbb{N_0}$, $\mathrm{Re}(\lambda)>0$, $\mathrm{Im}(\lambda)<0$, $x \in \mathbb{R}$ and $\mathrm{J}_\bullet$ is the Bessel function of the first kind.
Edit: I found out (through a completely different approach) that
$$f(x)=-\frac{1}{4\,\lambda^3}\left(\mathrm{e}^{-\lambda\,x}\,\mathrm{j}^n\,\sqrt{\frac{2\pi\mathrm{j}}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}\!\!\left(\frac{\lambda}{\mathrm{j}}\right)+\mathrm{e}^{-\mathrm{j}\,\lambda\,x}\,\mathrm{j}^{n+1}\,\sqrt{\frac{2\pi}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}\!\!\left(\lambda\right)\right)\quad\text{for}\quad x\geq1$$ and $$f(x)=-\frac{1}{4\,\lambda^3}\left(\mathrm{e}^{\lambda\,x}\,\mathrm{j}^n\,\sqrt{-\frac{2\pi\mathrm{j}}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}\!\!\left(-\frac{\lambda}{\mathrm{j}}\right)+\mathrm{e}^{\,\mathrm{j}\,\lambda\,x}\,\mathrm{j}^{n+1}\,\sqrt{-\frac{2\pi}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}\!\!\left(-\lambda\right)\right)\quad\text{for}\quad x\leq-1$$
Now the question remains, how to solve it for $-1\leq x \leq1$. I am accually surprised that the integral even changes. Any idea, how to get the solution for $x\geq1$ and $x\leq-1$ directly from the integral above?
Edit2: Thanks to Mathematica I got the solution for $x=0$
$$f(0)=\frac{\sqrt{\pi}\,(-1)^{n+1}}{8\,\lambda^{7/2}}\left((2+2\,\mathrm{j})\, \mathrm{I}_{n+\frac{1}{2}}\!\!\left(\mathrm{j}\,\lambda\right)+2\,\sqrt{2}\, \mathrm{I}_{n+\frac{1}{2}}\!\!\left(\lambda\right)+\lambda^{3/2}\left(_1\tilde{F}_2\left(1;\frac{3}{2}-\frac{n}{2},2+\frac{n}{2};-\frac{\lambda^2}{4}\right)-\,_1\tilde{F}_2\left(1;\frac{3}{2}-\frac{n}{2},2+\frac{n}{2};\frac{\lambda^2}{4}\right)\right)\right)$$
with $\mathrm{I}_\bullet$ is the modified Bessel function of the first kind and $_1\tilde{F}_2\left(\cdots;\cdots;\bullet\right)$ the regularized hypergeometric function. The solution for $x=0$ vanishes for odd $n$.
Not an answer
Residue theorem can be used to calculate the integral \begin{equation} f(x)=\frac{1}{\sqrt{2\pi}}\,\mathrm{j}^n \int_{-\infty}^\infty\frac{1}{\sqrt{k}\,(k^4-\lambda^4)}\mathrm{J}_{n+\frac{1}{2}}(k)\, \mathrm{e}^{-\mathrm{j}\,k\,x}\,\mathrm{d}k \end{equation} When $\Re x\le -1$, we close the contour by the large upper half-circle. This additionnal contribution vanishes owing to Jordan lemma ($\left|J(it)\right|\sim e^{\left|t\right|}/\sqrt t$), we have then \begin{equation} f(x)=\sqrt{2\pi}\mathrm{j}^{n+1}\sum_{k_p/\Im k_p>0}\operatorname{Res}\left( \phi(k),k=k_p \right) \end{equation} where $k_p$ are the 4 poles of \begin{equation} \phi(k)=\frac{\mathrm{J}_{n+\frac{1}{2}}(k)\, \mathrm{e}^{-\mathrm{j}\,k\,x}}{\sqrt{k}\,(k^4-\lambda^4)} \end{equation} defined on the complex plane cut along the negative real axis. The poles are $\lambda,\mathrm{j}\lambda,-\lambda,-\mathrm{j}\lambda$.
As $-\pi/2<\operatorname{Arg}(\lambda)<0$, only $\mathrm{j}\lambda$ and $-\lambda$ have to be taken into account in the summation. Their residues are \begin{align} \operatorname{Res}\left( \phi(k),k=\mathrm{j}\lambda\right)&=\frac{\mathrm{J}_{n+\frac{1}{2}}(\mathrm{j}\lambda)\, \mathrm{e}^{\lambda\,x}}{4\sqrt{\mathrm{j}\lambda}\,\mathrm{j}^3\lambda^3}\\ &=-\frac{1}{4\mathrm{j}\lambda^3}\sqrt{\frac{-\mathrm{j}}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}(\mathrm{j}\lambda)\, \mathrm{e}^{\lambda\,x}\\ \operatorname{Res}\left( \phi(k),k=-\lambda\right)&=\frac{\mathrm{J}_{n+\frac{1}{2}}(-\lambda)\, \mathrm{e}^{\mathrm{j}\lambda\,x}}{4\sqrt{-\lambda}\,(-\lambda)^3}\\ &=-\frac{1}{4\lambda^3}\sqrt{\frac{1}{-\lambda}}\mathrm{J}_{n+\frac{1}{2}}(-\lambda)\, \mathrm{e}^{\mathrm{j}\lambda\,x} \end{align} and thus \begin{equation} f(x)=-\frac{\sqrt{2\pi}\mathrm{j}^{n}}{4\lambda^3}\left[\sqrt{\frac{-\mathrm{j}}{\lambda}}\mathrm{J}_{n+\frac{1}{2}}(\mathrm{j}\lambda)\, \mathrm{e}^{\lambda\,x} +\mathrm{j}\sqrt{\frac{1}{-\lambda}}\mathrm{J}_{n+\frac{1}{2}}(-\lambda)\, \mathrm{e}^{\mathrm{j}\lambda\,x} \right] \end{equation} which is identical to the quoted result for $x\le-1$.
For $\Re x\ge1$, the same calculation can be made by closing the contour downwards. It can be easier to calculate the complex conjugate and use the previous result. Using the symmetry property \begin{equation} J_{n+1/2}\left(ze^{\mathrm{j}\pi }\right)=e^{\mathrm{j}(n+1/2)\pi }J_{n+1/2}\left(z\right) \end{equation} we obtain the proposed result.