平方完成

中学校で習う二次方程式の解法. 因数分解が思いつかなかった場合は解の公式を適用すればあとは地道に計算するだけである. 大人になってもその結果だけは覚えている人はそれなりの人数になると思われるが, これを導出できる人数となると少ないのではないだろうか.

xについての二次方程式 ax²+bx+c=0 の解は次のように求められる.

1. ax²+bx+c=0
2. a(x+b/(2a))²-b²/(4a)+c=0 (これが平方完成)
3. a(x+b/(2a))²=b²/(4a)-c (移項処理)
4. (x+b/(2a))²=(b²-4ac)/(4a²) (右辺を通分して両辺を a で割る)
5. x+b/(2a)=±√(b²-4ac)/(2a) (両辺の平方根をとる. 平方根には正のものと負のものが存在することに留意する.)
6. x=(-b±√(b²-4ac))/(2a) (左辺の定数部分を移項する)

これが二次方程式の解の公式の導出である. x の一次の項をうまく消去して処理する平方完成により解くことができたのである.

ここまでは 1 つの変数について考えてきたが, ここからは 2 つ以上の変数について考えてみよう. ただし、2 次式を考えるものとする.

例えば, 2 変数の場合は次のような二次曲線が得られることが知られている.

• 楕円 (例: x²–xy+y²=1)
• 双曲線 (例: x²–xy–y²=1)
• 放物線 (例: 2x²-4xy+2y²–x–y=1)

3 変数であれば二次曲面が得られる. 以下はその一例である.

• 錐面 (例: X^2 + Z^2 + 2 X Y + 2 Y Z – 2 Z X = 0)
• 楕円面 (例: 5 X^2 – 2 sqrt(2) Y (X + Z) + 2 X Z + 6 Y^2 + 5 Z^2 = 1)
• 一葉双曲面 (例: X^2 + 2 sqrt(2) X Y – 2 X Z + 2 sqrt(2) Y Z + Z^2 = 1)
• 二葉双曲面 (例: X^2 + 2 Y Z = -1)
• 楕円放物面 (例: (5 X^2 – 2 sqrt(2) Y (X + Z) + 2 X Z + 6 Y^2 + 5 Z^2) = 1)
• 双曲放物面 (例: (3 X^2 – 6 sqrt(2) Y (X + Z) – 2 X Z + 2 Y^2 + 3 Z^2) = 1)
• 楕円放物面 (例: (-X^2 – 6 sqrt(2) Y (X + Z) + 6 X Z + 2 Y^2 – Z^2) = 1)

4 変数以上ともなると図形的な表現は困難であるが、平方完成の話題からは離れてしまうので言及はやめておこう.

さて、このような多変数を扱うツールとしてベクトルや行列がある. ここまで扱った曲線はすべて次のような形で表わすことができる.

x^TAx+x^Tb+c=0 (x: 変数ベクトル, A: 行列, B: 定数ベクトル, c: 定数, ^T: 転置)

これを平行移動することにより x^Tb の項がない形に変形することを考える. その結果は (x-h)^TA(x-h)+k=0 となるはずである. 係数比較により h, k を A, b, c で表してみよう. ただし, A は正則行列 (すべての固有値が 0 とは異なる) であるものとする.

1. (x-h)^TA(x-h)+k=0
2. x^TAx–x^TAh-h^TAx+h^TAh+k=0 (第1項の展開)
3. x^TAx–x^TAh-((h^TAx)^T)^T+h^TAh+k=0 (第3項に転置を2回作用させる)
4. x^TAx–x^TAh-(x^TA^Th)^T+h^TAh+k=0 (第3項の内側の転置を作用. 積の順が逆転することに注意.)
5. x^TAx–x^TAh-x^TA^Th+h^TAh+k=0 (x^TA^Th はスカラー(1次正方行列)なので転置しても同じ値である)
6. x^TAx–x^T(A+A^T)h+h^TAh+k=0 (一部を結合法則でまとめる)
7. b=-(A+A^T)h, c=h^TAh+k (係数比較から得た結果)
8. h=-(A+A^T)^-1b, k=c-b^T(A+A^T)^-1A(A+A^T)^-1b (A+A^T は対称行列となるため (A+A^T)^T=A+A^T)
   • さらに, A^T=A (A そのものが対称行列) ならば h=-(1/2)A^-1b, k=c-(1/4)b^TA^-1b
9. 以上をまとめると, x^TAx+x^Tb+c=(x+(A+A^T)^-1b)^TA(x+(A+A^T)^-1b)+c-b^T(A+A^T)^-1A(A+A^T)^-1b となる.
   • さらに, A^T=A ならば x^TAx+x^Tb+c=(x+(1/2)A^-1b)^TA(x+(1/2)A^-1b)+c-(1/4)b^TA^-1b

以上が多変数の場合の平方完成となる.