前阵子整理复变函数习题,重新做了 2.3 节的一道老题:习题2.3-5,标题写着“柯西积分计算-第一类零阶修正贝塞尔函数的特殊值”。刚看到这种表述,我以为是那种一眼套公式就能出结果的题,真正动手才发现,它把柯西积分公式和特殊函数串在了一起,绕了一个很有意思的弯。
先说结论:这道题的核心对象是积分
[
\int_0^{2\pi}e^{x\cos\theta},d\theta,
]
它算出来以后正好是 (2\pi I_0(x)),其中 (I_0) 是第一类零阶修正贝塞尔函数。以前只会在数学物理方法里见到这种积分,但用柯西积分公式推一遍,会发现比死记“特殊函数积分公式表”有用得多。这篇文章就当作一次习题复盘来写,把每一步怎么想到的、哪里容易踩坑、最后又能推出哪些特殊值,都展开聊聊,适合正在学复变函数、或者工作中需要查特殊函数的人。
1. 这道题到底在算什么?先看清题目背后的逻辑
1.1 柯西积分公式在这里不是“直接代公式”
柯西积分公式的常见形式是
[
f(a)=\frac{1}{2\pi i}\oint_{C}\frac{f(z)}{z-a},dz,
]
它说明一个区域内的解析函数值,可以由边界上的值完全确定。这是复分析里最反直觉也最漂亮的性质之一:局部信息被整体信息锁死。
但习题2.3-5不是那种让你算 (f(0)) 的送分题。它真正想让你体会的是:怎么把一个看起来完全没有“除法结构”的实积分,改写成具有 (\frac{\text{解析函数}}{z-a}) 形式的围道积分,然后再通过展开、逐项取积分把这个“潜在的 (f(a))”找出来。
很多初学者看到柯西积分公式,只会找“分母一次幂、分子解析”的固定模式。一旦被积函数变成 (e^{x\cos\theta}),就觉得这不是柯西积分公式能染指的东西。实际上,处理这类问题的常规套路是变量替换到单位圆上:令 (z=e^{i\theta}),把 (\cos\theta) 改写成 (\frac{z+z^{-1}}2),把 (d\theta) 改写成 (\frac{dz}{iz})。这么一换,被积函数就从三角形式变成一个包含 (z^{-1}) 的围道积分。
这里已经埋了一个隐患:换完以后,被积函数在单位圆内部会出现 (z^{-1}),也就是原点处不再是“解析点”,而是本性奇点或者至少是一个高阶极点。这也是很多人推到最后卡住的地方:柯西积分公式要求被积函数在围道内部解析,可这里原点明明有问题。怎么办?答案是不要对整体套公式,而是先把函数展开,再逐项用柯西积分公式的高阶形式。
1.2 修正贝塞尔函数为什么会突然出现
说一下修正贝塞尔函数。第一类修正贝塞尔函数 (I_n(x)) 的常见级数定义是
[
I_n(x)=\sum_{k=0}^{\infty}\frac{(x/2)^{n+2k}}{k!(n+k)!},
]
当 (n=0) 时就是
[
I_0(x)=\sum_{k=0}^{\infty}\frac{(x/2)^{2k}}{(k!)^2}.
]
从级数能看到它和普通贝塞尔函数 (J_0(x)) 很相似,差别只在有没有 ((-1)^k) 这一项。正因为缺了正负交替,(I_0(x)) 在实轴上不振荡,而是从 (I_0(0)=1) 开始单调增长。
这个函数在现实中的出场率很高。比如柱坐标系下的热传导方程、扩散方程,分离变量后常出现修正贝塞尔方程;又比如方向统计里的 von Mises 分布,归一化常数分母就是 (2\pi I_0(\kappa))。积分 (\in
