라그랑주 보간법(Lagrangian interpolation)이란 서로 다른
x 1 , ⋯ , x n + 1 x_{1},\cdots,x_{n+1} x 1 , ⋯ , x n + 1 에 대하여
n + 1 n+1 n + 1 개의 점
( x 1 , y 1 ) , ⋯ , ( x n + 1 , y n + 1 ) (x_{1},y_{1}),\cdots,(x_{n+1},y_{n+1}) ( x 1 , y 1 ) , ⋯ , ( x n + 1 , y n + 1 ) 이 주어져 있을때, 이 점을 모두 지나는
n n n 차 이하의
다항식 을 구하는 공식을 말한다. 이름대로
조제프루이 라그랑주 가 만들었다.
서로 다른
x 1 , ⋯ , x n + 1 x_{1},\cdots,x_{n+1} x 1 , ⋯ , x n + 1 에 대하여 아래의 다항식
p i ( x ) = ∏ j ≠ i x − x j x i − x j = ( x − x 1 ) ⋯ ( x − x i − 1 ) ( x − x i + 1 ) ⋯ ( x − x n + 1 ) ( x i − x 1 ) ⋯ ( x i − x i − 1 ) ( x i − x i + 1 ) ⋯ ( x i − x n + 1 ) p_{i}(x)=\displaystyle\prod_{j \neq i} \frac{x-x_{j}}{x_{i}-x_{j}}=\frac{(x-x_{1})\cdots(x-x_{i-1})(x-x_{i+1})\cdots(x-x_{n+1})}{(x_{i}-x_{1})\cdots(x_{i}-x_{i-1})(x_{i}-x_{i+1})\cdots(x_{i}-x_{n+1})} p i ( x ) = j = i ∏ x i − x j x − x j = ( x i − x 1 ) ⋯ ( x i − x i − 1 ) ( x i − x i + 1 ) ⋯ ( x i − x n + 1 ) ( x − x 1 ) ⋯ ( x − x i − 1 ) ( x − x i + 1 ) ⋯ ( x − x n + 1 ) 을 라그랑주 기저 다항식이라 한다. 또한
p ( x ) = y 1 p 1 ( x ) + ⋯ + y n + 1 p n + 1 ( x ) p(x)=y_{1}p_{1}(x)+\cdots+y_{n+1}p_{n+1}(x) p ( x ) = y 1 p 1 ( x ) + ⋯ + y n + 1 p n + 1 ( x ) 을 라그랑주 보간 다항식이라 하며, 이 다항식은 점
( x 1 , y 1 ) , ⋯ , ( x n + 1 , y n + 1 ) (x_{1},y_{1}),\cdots,(x_{n+1},y_{n+1}) ( x 1 , y 1 ) , ⋯ , ( x n + 1 , y n + 1 ) 을 모두 지나는 유일한
n n n 차 이하의 다항식이다.
집합
P n ( F ) \mathcal{P}_{n}(F) P n ( F ) 을
n n n 차 이하의 다항식 집합이라 하자.
P n ( F ) \mathcal{P}_{n}(F) P n ( F ) 가 벡터공간이므로,
쌍대공간 P n ∗ \mathcal{P}_{n}^{*} P n ∗ 가 존재한다. 임의의 자연수
i ≤ n + 1 i\leq n+1 i ≤ n + 1 에 대하여
선형범함수 (linear functional)
L i : P n ( F ) → F L_{i}:\mathcal{P}_{n}(F)\to F L i : P n ( F ) → F 을
L i ( p ) = p ( x i ) L_{i}(p)=p(x_{i}) L i ( p ) = p ( x i ) 라고 정의하자.
L 1 , ⋯ , L n + 1 L_{1},\cdots,L_{n+1} L 1 , ⋯ , L n + 1 이
P n ∗ \mathcal{P}_{n}^{*} P n ∗ 의 기저가 된다는 것은 다음 두 식
L j ( p i ) = p i ( x j ) = 0 L_{j}(p_{i})=p_{i}(x_{j})=0 L j ( p i ) = p i ( x j ) = 0 for
i ≠ j i \neq j i = j L i ( p i ) = p i ( x i ) = 1 L_{i}(p_{i})=p_{i}(x_{i})=1 L i ( p i ) = p i ( x i ) = 1 을 만족하는
p 1 , ⋯ , p n + 1 ∈ P n ( F ) p_{1},\cdots,p_{n+1}\in \mathcal{P}_{n}(F) p 1 , ⋯ , p n + 1 ∈ P n ( F ) 이 존재한다는 것을 보이면 되는데, 그 이유는, 그것이 존재한다면,
{ L 1 , ⋯ , L n + 1 } \{L_{1},\cdots,L_{n+1}\} { L 1 , ⋯ , L n + 1 } 이
{ p 1 , ⋯ , p n + 1 } \{p_{1},\cdots,p_{n+1}\} { p 1 , ⋯ , p n + 1 } 의 쌍대기저가 되기 때문이다. 물론
p i p_{i} p i 는
x 1 , ⋯ , x i − 1 , x i + 1 , ⋯ , x n + 1 x_{1},\cdots, x_{i-1}, x_{i+1},\cdots,x_{n+1} x 1 , ⋯ , x i − 1 , x i + 1 , ⋯ , x n + 1 을 근으로 가지며,
p i ( x i ) = 1 p_{i}(x_{i})=1 p i ( x i ) = 1 이므로,
p i ( x ) = ∏ j ≠ i x − x j x i − x j = ( x − x 1 ) ⋯ ( x − x i − 1 ) ( x − x i + 1 ) ⋯ ( x − x n ) ( x i − x 1 ) ⋯ ( x i − x i − 1 ) ( x i − x i + 1 ) ⋯ ( x i − x n ) ∈ P n ( F ) p_{i}(x)=\displaystyle\prod_{j \neq i} \frac{x-x_{j}}{x_{i}-x_{j}}=\frac{(x-x_{1})\cdots(x-x_{i-1})(x-x_{i+1})\cdots(x-x_{n})}{(x_{i}-x_{1})\cdots(x_{i}-x_{i-1})(x_{i}-x_{i+1})\cdots(x_{i}-x_{n})}\in \mathcal{P}_{n}(F) p i ( x ) = j = i ∏ x i − x j x − x j = ( x i − x 1 ) ⋯ ( x i − x i − 1 ) ( x i − x i + 1 ) ⋯ ( x i − x n ) ( x − x 1 ) ⋯ ( x − x i − 1 ) ( x − x i + 1 ) ⋯ ( x − x n ) ∈ P n ( F ) 임을 쉽게 알 수 있다. 따라서,
{ p 1 , ⋯ , p n + 1 } \{p_{1},\cdots,p_{n+1}\} { p 1 , ⋯ , p n + 1 } 은
P n ( F ) \mathcal{P}_{n}(F) P n ( F ) 기저이고, 임의의 다항식
p ∈ P n ( F ) p \in \mathcal{P}_{n}(F) p ∈ P n ( F ) 에 대하여,
p ( x ) = y 1 p 1 ( x ) + ⋯ + y n + 1 p n + 1 ( x ) p(x)=y_{1}p_{1}(x)+\cdots+y_{n+1}p_{n+1}(x) p ( x ) = y 1 p 1 ( x ) + ⋯ + y n + 1 p n + 1 ( x ) 인
y 1 , ⋯ , y n + 1 y_{1} ,\cdots,y_{n+1} y 1 , ⋯ , y n + 1 이 유일하게 존재하는데,
p ( x i ) = ∑ j ≠ i y j p j ( x i ) + y i p i ( x i ) = y i p(x_{i})=\displaystyle\sum_{j \neq i}y_{j}p_{j}(x_{i})+y_{i}p_{i}(x_{i})=y_{i} p ( x i ) = j = i ∑ y j p j ( x i ) + y i p i ( x i ) = y i 가 성립함을 확인할 수 있다.
서로 다른
x 1 , ⋯ , x n + 1 x_{1},\cdots,x_{n+1} x 1 , ⋯ , x n + 1 에 대하여, 다항식
x i x^{i} x i ( 0 ≤ i ≤ n ) (0\leq i \leq n) ( 0 ≤ i ≤ n ) 은 점
( x j , x j i ) (x_{j},x_{j}^{i}) ( x j , x j i ) 를 지나므로, 라그랑주 보간법에 의해
x i = ∑ j = 1 n + 1 x j i p j x^{i}=\displaystyle\sum_{j=1}^{n+1}x_{j}^{i}p_{j} x i = j = 1 ∑ n + 1 x j i p j 라고 쓸 수 있다. 그런데,
β = { 1 , x , x 2 , ⋯ , x n } \beta=\{1,x,x^{2},\cdots,x^{n}\} β = { 1 , x , x 2 , ⋯ , x n } 과
β ′ = { p 1 , ⋯ , p n + 1 } \beta^\prime=\{p_{1},\cdots,p_{n+1}\} β ′ = { p 1 , ⋯ , p n + 1 } 이 모두
P n ( F ) \mathcal{P}_{n}(F) P n ( F ) 의 기저이므로 방데르몽드 행렬(Vandermonde matrix)
( 1 x 1 x 1 2 ⋯ x 1 n 1 x 2 x 2 2 ⋯ x 2 n ⋮ ⋮ ⋮ ⋱ ⋮ 1 x n + 1 x n + 1 2 ⋯ x n + 1 n ) \begin{pmatrix}1&x_{1}&x_{1}^{2}&\cdots&x_{1}^{n}\\1&x_{2}&x_{2}^{2}&\cdots&x_{2}^{n}\\ \vdots&\vdots&\vdots&\ddots&\vdots\\1&x_{n+1}&x_{n+1}^{2}&\cdots&x_{n+1}^{n} \end{pmatrix} 1 1 ⋮ 1 x 1 x 2 ⋮ x n + 1 x 1 2 x 2 2 ⋮ x n + 1 2 ⋯ ⋯ ⋱ ⋯ x 1 n x 2 n ⋮ x n + 1 n 은
β \beta β 에서
β ′ \beta^{\prime} β ′ 으로의
기저변환행렬 이라고 할 수 있다.