$a x^2 +b xy + c y^2 + d x + e y = left( begin{array}{cc} x & yend{array}right) left( begin{array}{cc} a & b/2\ b/2 & cend{array}right) left( begin{array}{c} x\y end{array}right) + left( begin{array}{cc} d & eend{array}right) left( begin{array}{c} x\y end{array}right) $
上面的表达方式给我们一个启发,可以把这样的二次多项式写成如下形式:
$x A x + b cdot x$
其中 x 是向量 $x_i$, A 是一个对称矩阵,b 是一次项的系数向量。将A对角化,就可以把 $x_1 x_2$ 一类的交叉项去掉,将第一项变成独立的平方之和。但是怎么把第二个线性项 b*x 去掉呢?
借用我们一元二次方程的配平方方法,我们考虑如下的平移,
$x to x + frac{b}{2A}$
于是
$(x +frac{b}{2A}) A (x+frac{b}{2A}) = ??$
贸然进行代数运算之前,我们注意 A 是矩阵,左右顺序不能随便。为谨慎起见,我们明确写出矩阵计算,并利用 A 的对称性 $A^{top} = A$:
$(x +frac{b}{2A}) A (x+frac{b}{2A}) = (x + frac{1}{2} A^{-1} b)^top A (x + frac{1}{2} A^{-1} b) \=(x^{top} + frac{1}{2} b^top A^{-1}) A (x + frac{1}{2} A^{-1} b)\=x^{top} A x + frac{1}{2} b^{top} x + frac{1}{2} x^{top} b + frac{1}{4} b^{top} A^{-1} b\=x A x + b x + frac{1}{4}b A^{-1} b$
可见:
$x A x + b x = (x +frac{b}{2A}) A (x +frac{b}{2A}) - frac{1}{4} bA^{-1} b$
在物理中,我们经常要用到所谓高斯积分:
$int_{-infty}^{+infty} e^{- a x^2} dx = sqrt{pi/a}$
对 n 个变量的二次多项式 $x A x $, A 对角化之后,变成了独立的积分相乘,对角系数相乘就是矩阵 A 的行列式,利用一维的高斯积分公式,我们得出:
$int e^{- x A x } d^{n}x = sqrt{frac{pi^n}{det A}}$
如果增加线性项 bx ,则相应的积分是:
$int e^{- x A x + bx } d^{n}x \= int exp left[ - (x +frac{b}{2A} ) A (x+frac {b}{2A} ) + frac{1}{4} b A^{-1} b right] d^n x \= frac{1}{4} b A^{-1} b int exp left[ - (x +frac{b}{2A} ) A (x+frac {b}{2A} ) right] d^n x\ = frac{1}{4} b A^{-1} b sqrt{frac{pi^n}{det A}}$