调频波的贝塞尔函数展开

---

前言

《锁相环（PLL）电路设计与应用》一书中在P12页用到了两个单音调频波的贝塞尔展开公式，我在想为什么叫做贝塞尔展开，开启来不是一个像傅里叶展开吗公式吗？因此这里绕过贝塞尔展开，我们直接尝试从傅里叶级数的角度求解两个展开问题，不需要贝塞尔函数的前置知识，并且可以从这个具体的问题来认识整数阶的贝塞尔函数：

---

一、要解决的问题

第一：$cos(\beta sin\theta )$的傅里叶级数
第二：$sin(\beta sin\theta )$的傅里叶级数

二、要用到的积分公式

由于后面有的积分需要计算，所以先引出一些公式，这考虑当$k\le n$时，积分：
$$\begin{array}{rl}& {\int }_0^{2\pi }si{n}^{2k}\theta [cos(2n\theta )+isin(2n\theta )]d\theta ={\int }_0^{2\pi }{\left(\frac{e^{i\theta }-e^{-i\theta }}{2i}\right)}^{2k}{e}^{i2n\theta }d\theta =\frac{-i(-1{)}^k}{2^{2k}}{\int }_0^{2\pi }\frac{(e^{2i\theta }-1{)}^{2k}}{e^{(2k-2n+1)i\theta }}d{e}^{i\theta }\end{array}$$
将$e^{i\theta }$换成复数$z$,则可以通过柯西积分得到它的值，可以注意到：
$$\begin{array}{rl}& {\int }_0^{2\pi }si{n}^{2k}\theta [cos(2n\theta )+isin(2n\theta )]d\theta \\ =& \frac{-i(-1{)}^k}{2^{2k}}\oint \frac{(z^2-1{)}^{2k}}{z^{2k-2n+1}}dz=\frac{2\pi (-1{)}^k}{2^{2k}}\frac{1}{2\pi i}\oint \frac{(z^2-1{)}^{2k}}{z^{2k-2n+1}}dz=\frac{2\pi (-1{)}^k}{2^{2k}}\frac{[(z^2+1{)}^{2k}{]}^{(2k-2n)}{|}_{z=0}}{(2k-2n)!}\\ =& \frac{2\pi (-1{)}^k}{2^{2k}}\frac{(-1{)}^{n+k}{C}_{2k}^{k-n}(2k-2n)!}{(2k-2n)!}\\ =& \frac{2\pi (-1{)}^n}{2^{2k}}\frac{(2k)!}{(k-n)!(k+n)!}\end{array}$$
当$k<n$时，被积分的函数时全纯的，积分为零，这个从上面的公式也可以看出来。

上面的积分两边同时取实部可以看出我们有积分公式：
$$\begin{array}{rl}& {\int }_0^{2\pi }si{n}^{2k}\theta cos(2n\theta )d\theta =\frac{2\pi (-1{)}^n}{2^{2k}}\frac{(2k)!}{(k-n)!(k+n)!}\cdots (1)\end{array}$$

仿照上面的办法，我们再考虑积分：
$$\begin{array}{rl}& {\int }_0^{2\pi }si{n}^{2k+1}\theta [cos[(2n+1)\theta ]+isin[(2n+1)\theta ]d\theta ={\int }_0^{2\pi }{\left(\frac{e^{i\theta }-e^{-i\theta }}{2i}\right)}^{2k+1}{e}^{i(2n+1)\theta }d\theta \\ =& \frac{(-1{)}^k}{2^{2k+1}}{\int }_0^{2\pi }\frac{(e^{2i\theta }-1{)}^{2k+1}}{e^{(2k-2n+1)i\theta }}d{e}^{i\theta }\\ =& \frac{(-1{)}^{k+1}2\pi i}{2^{2k+1}}\frac{1}{2\pi i}\oint \frac{(z^2-1{)}^{2k+1}}{z^{2k-2n+1}}dz=\frac{(-1{)}^{k+1}2\pi i}{2^{2k+1}}\frac{(-1{)}^{n+k+1}{C}_{2k+1}^{k-n}(2k-2n)!}{(2k-2n)!}\\ =& \frac{(-1{)}^{n}2\pi i}{2^{2k+1}}\frac{(2k+1)!}{(k-n)!(k+n+1)!}\end{array}$$
两边取虚部得到公式：
$$\begin{array}{rl}& {\int }_0^{2\pi }si{n}^{2k+1}\theta sin[(2n+1)\theta d\theta =\frac{(-1{)}^{n}2\pi }{2^{2k+1}}\frac{(2k+1)!}{(k-n)!(k+n+1)!}\cdots (2)\end{array}$$

三、$cos(\beta sin\theta )$的傅里叶系数

这个函数是一个偶函数，因此只存在偶数项，我们来使用系数求解公式：
$$\begin{array}{rl}{A}_{2n}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }cos(\beta sin\theta )cos(2n\theta )d\theta \end{array}$$
注意泰勒级数：
$$\begin{array}{r}cos(\beta sin\theta )=\sum _{k=0}^{\infty }(-1{)}^k\frac{si{n}^{2k}\theta }{(2k)!}{\beta }^{2k}\end{array}$$

于是我们得到傅里叶系数的积分公式：
$$\begin{array}{rl}{A}_{2n}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }cos(\beta sin\theta )cos(2n\theta )d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }(-1{)}^k\frac{{\beta }^{2k}}{(2k)!}{\int }_0^{2\pi }si{n}^{2k}\theta cos(2n\theta )d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }(-1{)}^k\frac{{\beta }^{2k}}{(2k)!}\frac{2\pi ((-1{)}^n)}{2^{2k}}\frac{(2k)!}{(k-n)!(k+n)!}\\ =& 2\sum _{k=0}^{\infty }\frac{(-1{)}^{k+n}}{(k-n)!(k+n)!}{\left(\frac{\beta }{2}\right)}^{2k}\\ =& 2\sum _{k=0}^{\infty }\frac{(-1{)}^k}{(k)!(k+2n)!}{\left(\frac{\beta }{2}\right)}^{2k+2n}\end{array}$$
上面的最后一步利用了$k<n$时$(k-n)!=\infty$的事实，并做了下标变换。

四、$sin(\beta sin\theta )$的傅里叶系数

这是一个奇函数，因此只存在奇数项，我们来使用系数求解公式：
$$\begin{array}{rl}{A}_{2n+1}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }sin(\beta sin\theta )sin[(2n+1)\theta ]d\theta \end{array}$$
注意泰勒级数：
$$\begin{array}{r}sin(\beta sin\theta )=\sum _{k=0}^{\infty }(-1{)}^k\frac{si{n}^{2k+1}\theta }{(2k+1)!}{\beta }^{2k+1}\end{array}$$
所以对应的傅里叶系数的积分公式如下：
$$\begin{array}{rl}{A}_{2n+1}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }sin(\beta sin\theta )sin[(2n+1)\theta ]d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }(-1{)}^k\frac{{\beta }^{2k+1}}{(2k+1)!}{\int }_0^{2\pi }si{n}^{2k+1}\theta sin[(2n+1)\theta d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }(-1{)}^k\frac{{\beta }^{2k+1}}{(2k+1)!}\frac{(-1{)}^{n}2\pi }{2^{2k+1}}\frac{(2k+1)!}{(k-n)!(k+n+1)!}\\ =& 2\sum _{k=0}^{\infty }\frac{(-1{)}^{k+n}}{(k-n)!(k+n+1)!}{\left(\frac{\beta }{2}\right)}^{2k+1}\\ =& 2\sum _{k=0}^{\infty }\frac{(-1{)}^k}{k!(k+2n+1)!}{\left(\frac{\beta }{2}\right)}^{2k+2n+1}\end{array}$$

五、认识整数阶贝塞尔函数

我们先把上述两个函数的傅里叶系数整理在下面：
$$\begin{array}{rl}{A}_{2n}(\beta )=& 2\sum _{k=0}^{\infty }\frac{(-1{)}^k}{(k)!(k+2n)!}{\left(\frac{\beta }{2}\right)}^{2k+2n}\\ A_{2n+1}(\beta )=& 2\sum _{k=0}^{\infty }\frac{(-1{)}^k}{k!(k+2n+1)!}{\left(\frac{\beta }{2}\right)}^{2k+2n+1}\end{array}$$
如果我们定义一个函数：
$$\begin{array}{rl}{J}_m(\beta )=& \sum _{k=0}^{\infty }\frac{(-1{)}^k}{(k)!(k+m)!}{\left(\frac{\beta }{2}\right)}^{2k+m}\end{array}$$
那么两个系数可以分别简记为：
$$\begin{array}{rl}{A}_{2n}(\beta )=& 2J_{2n}(\beta )\\ A_{2n+1}(\beta )=& 2J_{2n+1}(\beta )\end{array}$$
这里我们定义的函数$J_m(\beta )$就叫做贝塞尔函数了。

六、$cos(\beta sin\theta )$以及$sin(\beta sin\theta )$的傅里叶展开(贝塞尔展开)

现在我们已经有了傅里叶系数了，我们来按照常规做法把它们展开吧！
$$\begin{array}{rl}cos(\beta sin\theta )=& J_0(\beta )+2\sum _{n=0}^{\infty }{J}_{2n}(\beta )cos(2n\theta )\\ sin(\beta sin\theta )=& 2\sum _{n=0}^{\infty }{J}_{2n+1}(\beta )sin[(2n+1)\theta ]\end{array}$$

七、其他尝试方向的附加产物——两个组合恒等式

如果我们采用下面的办法则可以得到两个组合恒等式：
利用：
$$\begin{array}{rl}cos(2n\theta )=& \sum _{l=0}^n(-1{)}^l{C}_{2n}^{2l}co{s}^{2n-2l}\theta si{n}^{2l}\theta \end{array}$$
我们把傅里叶系数的积分公式写成：
$$\begin{array}{rl}{A}_{2n}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }cos(\beta sin\theta )cos(2n\theta )d\theta \\ =& \frac{1}{\pi }{\int }_0^{2\pi }\left[\sum _{k=0}^{\infty }(-1{)}^k\frac{si{n}^{2k}\theta }{(2k)!}{\beta }^{2k}\right]\left[\sum _{l=0}^n(-1{)}^l{C}_{2n}^{2l}co{s}^{2n-2l}\theta si{n}^{2l}\theta \right]d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }\sum _{l=0}^n(-1{)}^{k+l}\frac{C_{2n}^{2l}}{(2k)!}{\beta }^{2k}{\int }_0^{2\pi }co{s}^{2n-2l}\theta si{n}^{2l+2k}\theta d\theta \end{array}$$
利用计算公式见一个一个特别地参数积分I(p,q)：
$$\begin{array}{r}I(p,q)={\int }_0^{2\pi }co{s}^{2p}\theta si{n}^{2q}\theta d\theta =\frac{2\pi }{2^{2p+2q}}\frac{(2p)!(2q)!}{p!(p+q)!q!}\end{array}$$
$$\begin{array}{rl}{A}_{2n}(\beta )=& \frac{1}{\pi }\sum _{k=0}^{\infty }\sum _{l=0}^n(-1{)}^{k+l}\frac{C_{2n}^{2l}}{(2k)!}{\beta }^{2k}I(n-l,l+k)\\ =& 2\sum _{k=0}^{\infty }(-1{)}^k{\left(\frac{\beta }{2}\right)}^{2k}\frac{(2n)!}{2^{2n}(2k)!(n+k)!}\sum _{l=0}^n\frac{(-1{)}^l[2(l+k)]!}{(2l)!(n-l)!(l+k)!}\end{array}$$
要使$A_{2n}(\beta )=2J_{2n}(\beta )$，我们必定就得到了恒等式：
$$\begin{array}{r}\sum _{l=0}^n\frac{(-1{)}^l[2(l+k)]!}{(2l)!(n-l)!(l+k)!}=\frac{(-1{)}^n{2}^{2n}(2k)!}{(2n)!(k-n)!}\end{array}$$

利用：
$$\begin{array}{rl}sin[(2n+1)\theta ]=& \sum _{l=0}^n(-1{)}^l{C}_{2n+1}^{2l+1}co{s}^{2n-2l}\theta si{n}^{2l+1}\theta \end{array}$$
我们把傅里叶系数的积分公式写成：
$$\begin{array}{rl}{A}_{2n+1}(\beta )=& \frac{1}{\pi }{\int }_0^{2\pi }sin(\beta sin\theta )sin[(2n+1)\theta ]d\theta \\ =& \frac{1}{\pi }{\int }_0^{2\pi }\left[\sum _{k=0}^{\infty }(-1{)}^k\frac{si{n}^{2k+1}\theta }{(2k+1)!}{\beta }^{2k+1}\right]\left[\sum _{l=0}^n(-1{)}^l{C}_{2n+1}^{2l+1}co{s}^{2n-2l}\theta si{n}^{2l+1}\theta \right]d\theta \\ =& \frac{1}{\pi }\sum _{k=0}^{\infty }\sum _{l=0}^n(-1{)}^{k+l}\frac{C_{2n+1}^{2l+1}}{(2k+1)!}{\beta }^{2k+1}{\int }_0^{2\pi }co{s}^{2n-2l}\theta si{n}^{2l+2k+2}\theta d\theta \\ =& 2\sum _{k=0}^{\infty }\sum _{l=0}^n(-1{)}^{k+l}\frac{C_{2n+1}^{2l+1}}{(2k+1)!}{\beta }^{2k+1}\frac{1}{2^{2n+2k+2}}\frac{(2n-2l)!(2l+2k+2)!}{(n-l)!(n+k+1)!(l+k+1)!}\\ =& 2\sum _{k=0}^{\infty }(-1{)}^k{\left(\frac{\beta }{2}\right)}^{2k+1}\frac{(2n+1)!}{2^{2n+1}(2k+1)!(n+k+1)!}\sum _{l=0}^n\frac{(-1{)}^l(2l+2k+2)!}{(2l+1)!(n-l)!(l+k+1)!}\end{array}$$
要使$A_{2n+1}(\beta )=2J_{2n+1}(\beta )$，我们必定就得到了恒等式：
$$\begin{array}{r}\sum _{l=0}^n\frac{(-1{)}^l(2l+2k+2)!}{(2l+1)!(n-l)!(l+k+1)!}=\frac{(-1{)}^n{2}^{2n+1}(2k+1)!}{(2n+1)!(k-n)!}\end{array}$$

总结

通过努力，我们成功地从傅里叶级数分解的角度解释了调频波中使用的贝塞尔展开式，并且最后介绍了先前未能成功证明展开式的思路，反而我们还附带地取得了两个组合恒等式，收获不小。