今天咱们聊聊这个高斯积分，这家伙可不是省油的灯。我当初刚碰上它的时候，也是一头雾水，感觉无从下手。但实践出真知嘛多捣鼓捣鼓，总能找到点门道。
说起来，我琢磨这玩意儿，也不是一天两天了。最早是啥时候？好像是之前捣鼓一个啥统计相关的玩意儿，老是碰到它，像个拦路虎似的，就是那个形如 ∫ e^(-x²) dx 的积分，积分限从负无穷到正无穷。看着简单，但你用常规的那些换元法、分部积分法去试，根本行不通！
我的硬磕之路
我寻思着，不就一积分嘛硬积呗！结果？直接给我干懵了。常规的换元，分部积分，试了个遍，愣是积不出来。那叫一个头大。我就不信这个邪，还专门找了本微积分的书，翻来覆去地看，希望能找到点蛛丝马迹。
后来没办法了，就到处翻资料，看人家大佬们是怎么收拾它的。这一看，还真发现了个门道儿——就是先算它的平方！这个思路真是绝了，我当时就想，这谁想出来的，太聪明了。
具体咋操作？我这就把我当时捣鼓的过程给大伙儿捋一捋：
咱们设这个积不出来的玩意儿叫 I，也就是：
I = ∫ e^(-x²) dx (积分限从 -∞ 到 +∞)
然后，关键的一步来了，咱们算 I²。既然是平方，那就是自己乘自己嘛
I² = (∫ e^(-x²) dx) (∫ e^(-y²) dy)
这里我把第二个积分里的自变量换成了 y，这不影响结果，对？毕竟定积分的值跟积分变量叫啥名没关系。这么一弄，我们就可以把这两个积分合成一个二重积分：
I² = ∫∫ e^(-(x² + y²)) dx dy (积分区域是整个 x-y 平面)
看到 x² + y² 这个形式，是不是有点眼熟？没错！极坐标！这玩意儿简直是为极坐标量身定做的。咱就赶紧换上极坐标：
- 令 x = r cosθ
- 令 y = r sinθ
那么，x² + y² = r²，还有那个面积微元 dx dy = r dr dθ。积分限也得跟着换：r 是从 0 到 +∞，θ 是从 0 到 2π (转一圈嘛)。
换完之后，I² 就变成了这个样子：
I² = ∫ (从0到2π) [ ∫ (从0到+∞) e^(-r²) r dr ] dθ
这下看着就清爽多了！咱们先看里面那个对 r 的积分：∫ e^(-r²) r dr。这个好办，凑个微分就行了。令 u = r²，那么 du = 2r dr，所以 r dr = du/2。积分就变成了：
∫ e^(-u) (1/2) du = (-1/2) e^(-u) = (-1/2) e^(-r²)
把 r 的积分限 0 和 +∞ 代进去：
[(-1/2) e^(-r²)] 从 0 到 +∞ = (-1/2) (0 - 1) = 1/2
里面这块算出来了，是 1/2！那 I² 就只剩下外面那个对 θ 的积分了：
I² = ∫ (从0到2π) (1/2) dθ
这个就更简单了，积出来就是 (1/2) [θ] (从0到2π) = (1/2) (2π - 0) = π。
咱们辛辛苦苦搞了半天，得到了 I² = π。
那 I 是啥？开方呗！因为原来的被积函数 e^(-x²) 始终是正的，所以积分结果也肯定是正的。
I = √π
搞定！当时算出这个结果，我真是长舒一口气。这个过程虽然看着绕了点，但每一步都是实打实的推导，感觉特别踏实。
扩展一下思路
后来我还琢磨了一下，要是积分里面是 e^(-ax²) 这种形式？也差不多，稍微换元一下，套用咱们上面那个 √π 的结果，就能搞定。比如 ∫ e^(-ax²) dx (a > 0)，做个替换 u = √a x，那么 dx = du/√a。代进去之后，就能捣鼓出结果是 √(π/a)。这个也挺常用的，顺手记一下很有用。
解决高斯积分这个过程，对我来说算是一次挺有意思的探索。关键就是那个平方再换极坐标的思路，真是神来之笔。有时候解决问题，就得换个角度，不能死脑筋。这回的实践记录就到这儿，希望能给同样对这玩意儿犯迷糊的朋友一点启发。
还没有评论，来说两句吧...