フーリエ級数・フーリエ変換・フーリエ解析

補間多項式

相異なる点 x0,x1, において、関数値 f(x0),f(x1), が与えられている。

(1) 2点 (x0,f(x0)),(x1,f(x1)) を通る高々1次の関数 y=a0+a1x の係数 a0,a1 を決定せよ。

(2) 3点 (x0,f(x0)),(x1,f(x1)),(x2,f(x2)) を通る高々2次の関数 y=a0+a1x+a2x2 の係数 a0,a1a2 を決定せよ。

(3) n+1(x0,f(x0)),(x1,f(x1)),,(xn,f(xn)) を通る高々 n 次の関数 y=a0+a1x++anxn が一意的に決まることを示せ。

(1) 単に答えを求めるなら
yf(x0)=f(x1)f(x0)x1x0(xx0)
となるが、(2), (3) を見越して、より一般的な方法で導くことにする。

今、求めたい a0,a1 は次の連立方程式の解である。
a0+a1x0=f(x0)a0+a1x1=f(x1)
これを行列で書けば
(1x01x1)(a0a1)=(f(x0)f(x1))
このような連立1次方程式の解を求める表式はクラメルの公式として知られている。
すなわち
A=(1x01x1), A0=(f(x0)x0f(x1)x1), A1=(1f(x0)1f(x1))
とするとき
a0=detA0detA, a1=detA1detA
と求められる。実際に a0,a1 をこの表式から求めると
detA=x1x0a0=x1f(x0)x0f(x1)a1=f(x1)f(x0)
となり、これは最初に示した一次式と一致する。

(2) (1) の方法を3行3列の場合に適用することにより
a0=1detA|f(x0)x0x02f(x1)x1x12f(x2)x2x22|, a1=1detA|1f(x0)x021f(x1)x121f(x2)x22|, a2=1detA|1x0f(x0)1x1f(x1)1x2f(x2)|
と求まる。ここで
A=(1x0x021x1x121x2x22)
である。

(3) 全く同様の式が (n+1)×(n+1)行列について成り立つが、ここで
A=(1x0x02x0n1x1x12x1n1xnxn2xnn)
が逆行列を持つかどうかが、一意的に a0,a1,,an が決まるかを決定する。
ここで、x0,x1,,xn は全て異なるので、detA0 であるので、A は逆行列を持つ。
従って、a0,a1,,an はただ1つに決まることが分かる。