1 重み付き内積と最良近似
定義 1.1. a < b a<b a < b を実数とする。関数w : ( a , b ) → R w\colon(a,b)\to\R w : ( a , b ) → R が、任意のx ∈ ( a , b ) x\in(a,b) x ∈ ( a , b ) に対してw ( x ) > 0 w(x)>0 w ( x ) > 0 を満たし、a < c < d < b a<c<d<b a < c < d < b を満たす任意の実数c , d c,d c , d について[ c , d ] [c,d] [ c , d ] で Riemann 可積分であり、広義積分
μ 0 : = lim c → a + , d → b − ∫ c d w ( x ) d x \mu_0:=\lim_{c\to a+,\ d\to b-}\int_c^dw(x)\,dx μ 0 := c → a + , d → b − lim ∫ c d w ( x ) d x が有限であるとき、w w w を( a , b ) (a,b) ( a , b ) 上の 重み関数 (weight function ) という。
補題 1.2. c < d c<d c < d を実数とし、w : [ c , d ] → R w\colon[c,d]\to\R w : [ c , d ] → R を Riemann 可積分な関数とする。
h : [ c , d ] → R h\colon[c,d]\to\R h : [ c , d ] → R が連続ならば、h w hw h w は[ c , d ] [c,d] [ c , d ] で Riemann 可積分である。
任意のx ∈ [ c , d ] x\in[c,d] x ∈ [ c , d ] に対してw ( x ) > 0 w(x)>0 w ( x ) > 0 ならば、∫ c d w ( x ) d x > 0 \int_c^dw(x)\,dx>0 ∫ c d w ( x ) d x > 0 である。
証明. (1) を示す。§D1.13 定理 2.1 によりH : = max x ∈ [ c , d ] ∣ h ( x ) ∣ H:=\max_{x\in[c,d]}\lvert h(x)\rvert H := max x ∈ [ c , d ] ∣ h ( x )∣ は有限であり、Riemann 可積分なw w w は有界であるからK : = sup x ∈ [ c , d ] ∣ w ( x ) ∣ K:=\sup_{x\in[c,d]}\lvert w(x)\rvert K := sup x ∈ [ c , d ] ∣ w ( x )∣ は有限である。したがってh w hw h w は有界である。[ c , d ] [c,d] [ c , d ] の部分区間I I I と関数g g g に対してosc I g : = sup I g − inf I g \operatorname{osc}_Ig:=\sup_Ig-\inf_Ig osc I g := sup I g − inf I g と置く。x , y ∈ I x,y\in I x , y ∈ I ならば
∣ h ( x ) w ( x ) − h ( y ) w ( y ) ∣ ≤ ∣ h ( x ) ∣ ∣ w ( x ) − w ( y ) ∣ + ∣ w ( y ) ∣ ∣ h ( x ) − h ( y ) ∣ ≤ H osc I w + K osc I h \lvert h(x)w(x)-h(y)w(y)\rvert\le\lvert h(x)\rvert\,\lvert w(x)-w(y)\rvert+\lvert w(y)\rvert\,\lvert h(x)-h(y)\rvert\le H\operatorname{osc}_Iw+K\operatorname{osc}_Ih ∣ h ( x ) w ( x ) − h ( y ) w ( y )∣ ≤ ∣ h ( x )∣ ∣ w ( x ) − w ( y )∣ + ∣ w ( y )∣ ∣ h ( x ) − h ( y )∣ ≤ H osc I w + K osc I h であるから、osc I ( h w ) ≤ H osc I w + K osc I h \operatorname{osc}_I(hw)\le H\operatorname{osc}_Iw+K\operatorname{osc}_Ih osc I ( h w ) ≤ H osc I w + K osc I h である。ε > 0 \varepsilon>0 ε > 0 をとる。§D1.15 定理 3.1 によりh h h は[ c , d ] [c,d] [ c , d ] で一様連続であるから、δ > 0 \delta>0 δ > 0 であって、x , y ∈ [ c , d ] x,y\in[c,d] x , y ∈ [ c , d ] 、∣ x − y ∣ < δ \lvert x-y\rvert<\delta ∣ x − y ∣ < δ ならば∣ h ( x ) − h ( y ) ∣ < ε \lvert h(x)-h(y)\rvert<\varepsilon ∣ h ( x ) − h ( y )∣ < ε となるものが存在する。§D1.17 定理 2.4 により[ c , d ] [c,d] [ c , d ] の分割P P P であってU ( w , P ) − L ( w , P ) < ε U(w,P)-L(w,P)<\varepsilon U ( w , P ) − L ( w , P ) < ε を満たすものが存在する。P P P に分割点を加えて、各小区間の長さがδ \delta δ 未満である分割Q Q Q をとると、§D1.17 補題 1.3 によりU ( w , Q ) − L ( w , Q ) < ε U(w,Q)-L(w,Q)<\varepsilon U ( w , Q ) − L ( w , Q ) < ε であり、Q Q Q の各小区間I I I でosc I h ≤ ε \operatorname{osc}_Ih\le\varepsilon osc I h ≤ ε である。Q Q Q の小区間について和をとると
U ( h w , Q ) − L ( h w , Q ) ≤ H ( U ( w , Q ) − L ( w , Q ) ) + K ε ( d − c ) ≤ ( H + K ( d − c ) ) ε U(hw,Q)-L(hw,Q)\le H\bigl(U(w,Q)-L(w,Q)\bigr)+K\varepsilon(d-c)\le\bigl(H+K(d-c)\bigr)\varepsilon U ( h w , Q ) − L ( h w , Q ) ≤ H ( U ( w , Q ) − L ( w , Q ) ) + K ε ( d − c ) ≤ ( H + K ( d − c ) ) ε である。ε > 0 \varepsilon>0 ε > 0 は任意であるから、§D1.17 定理 2.4 によりh w hw h w は[ c , d ] [c,d] [ c , d ] で Riemann 可積分である。
(2) を示す。w > 0 w>0 w > 0 であるから、§D1.17 命題 3.4 により∫ c d w d x ≥ 0 \int_c^dw\,dx\ge0 ∫ c d w d x ≥ 0 である。∫ c d w d x = 0 \int_c^dw\,dx=0 ∫ c d w d x = 0 と仮定する。閉区間[ c k , d k ] [c_k,d_k] [ c k , d k ] (k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 )を、c k < d k c_k<d_k c k < d k 、[ c k + 1 , d k + 1 ] ⊆ [ c k , d k ] [c_{k+1},d_{k+1}]\subseteq[c_k,d_k] [ c k + 1 , d k + 1 ] ⊆ [ c k , d k ] 、w w w の[ c k , d k ] [c_k,d_k] [ c k , d k ] への制限は Riemann 可積分でその積分は0 0 0 、k ≥ 1 k\ge1 k ≥ 1 ならばsup [ c k , d k ] w < 1 / k \sup_{[c_k,d_k]}w<1/k sup [ c k , d k ] w < 1/ k を満たすように帰納的に定める。[ c 0 , d 0 ] : = [ c , d ] [c_0,d_0]:=[c,d] [ c 0 , d 0 ] := [ c , d ] とする。[ c k , d k ] [c_k,d_k] [ c k , d k ] が定まったとき、[ c k , d k ] [c_k,d_k] [ c k , d k ] 上のw w w の上積分は積分0 0 0 に等しいから、[ c k , d k ] [c_k,d_k] [ c k , d k ] の分割P P P であってU ( w , P ) < ( d k − c k ) / ( k + 1 ) U(w,P)<(d_k-c_k)/(k+1) U ( w , P ) < ( d k − c k ) / ( k + 1 ) を満たすものが存在する。P P P の小区間の長さの和はd k − c k d_k-c_k d k − c k であるから、P P P の小区間のうちw w w の上限が1 / ( k + 1 ) 1/(k+1) 1/ ( k + 1 ) 未満であるものが存在し、それを[ c k + 1 , d k + 1 ] [c_{k+1},d_{k+1}] [ c k + 1 , d k + 1 ] とする。§D1.17 命題 3.1 と§D1.17 定理 3.6 により、w w w のP P P の各小区間への制限は Riemann 可積分であり、それらの積分の和は0 0 0 である。§D1.17 命題 3.4 によりそれらの積分はすべて0 0 0 以上であるから、すべて0 0 0 であり、[ c k + 1 , d k + 1 ] [c_{k+1},d_{k+1}] [ c k + 1 , d k + 1 ] 上の積分も0 0 0 である。
区間が入れ子であるから、任意のk , j ∈ N ≥ 0 k,j\in\N k , j ∈ N ≥ 0 に対してc k ≤ d j c_k\le d_j c k ≤ d j である。x ∗ : = sup k c k x^*:=\sup_kc_k x ∗ := sup k c k と置くと、任意のk ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 に対してc k ≤ x ∗ ≤ d k c_k\le x^*\le d_k c k ≤ x ∗ ≤ d k である。したがって任意のk ∈ N ≥ 1 k\in\NN k ∈ N ≥ 1 に対してw ( x ∗ ) ≤ sup [ c k , d k ] w < 1 / k w(x^*)\le\sup_{[c_k,d_k]}w<1/k w ( x ∗ ) ≤ sup [ c k , d k ] w < 1/ k であり、w ( x ∗ ) ≤ 0 w(x^*)\le0 w ( x ∗ ) ≤ 0 である。これはw ( x ∗ ) > 0 w(x^*)>0 w ( x ∗ ) > 0 と両立しないから、∫ c d w d x > 0 \int_c^dw\,dx>0 ∫ c d w d x > 0 である。▨
補題 1.3. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、μ 0 : = ∫ a b w ( x ) d x \mu_0:=\int_a^bw(x)\,dx μ 0 := ∫ a b w ( x ) d x とする。h : [ a , b ] → R h\colon[a,b]\to\R h : [ a , b ] → R を連続関数とする。
広義積分∫ a b h ( x ) w ( x ) d x \int_a^bh(x)w(x)\,dx ∫ a b h ( x ) w ( x ) d x は絶対収束し、∣ ∫ a b h ( x ) w ( x ) d x ∣ ≤ μ 0 ∥ h ∥ ∞ \bigl|\int_a^bh(x)w(x)\,dx\bigr|\le\mu_0\lVert h\rVert_\infty ∫ a b h ( x ) w ( x ) d x ≤ μ 0 ∥ h ∥ ∞ である。
任意のx ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] に対してh ( x ) ≥ 0 h(x)\ge0 h ( x ) ≥ 0 であり、h ( x 0 ) > 0 h(x_0)>0 h ( x 0 ) > 0 を満たすx 0 ∈ [ a , b ] x_0\in[a,b] x 0 ∈ [ a , b ] が存在するならば、∫ a b h ( x ) w ( x ) d x > 0 \int_a^bh(x)w(x)\,dx>0 ∫ a b h ( x ) w ( x ) d x > 0 である。
証明. (1) を示す。a < c < d < b a<c<d<b a < c < d < b とする。w > 0 w>0 w > 0 であるから、§D1.17 定理 3.6 と§D1.17 命題 3.4 により∫ c d w d x \int_c^dw\,dx ∫ c d w d x はc c c について減少しd d d について増加し、その極限μ 0 \mu_0 μ 0 以下である。g : [ a , b ] → R g\colon[a,b]\to\R g : [ a , b ] → R を任意のx ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] に対して0 ≤ g ( x ) ≤ ∥ h ∥ ∞ 0\le g(x)\le\lVert h\rVert_\infty 0 ≤ g ( x ) ≤ ∥ h ∥ ∞ を満たす連続関数とし、a < e < b a<e<b a < e < b を一つ固定する。補題 1.2 (1) によりg w gw g w は( a , b ) (a,b) ( a , b ) に含まれる各閉区間で Riemann 可積分であるから、a < c < e a<c<e a < c < e に対してG − ( c ) : = ∫ c e g ( x ) w ( x ) d x G_-(c):=\int_c^eg(x)w(x)\,dx G − ( c ) := ∫ c e g ( x ) w ( x ) d x 、e < d < b e<d<b e < d < b に対してG + ( d ) : = ∫ e d g ( x ) w ( x ) d x G_+(d):=\int_e^dg(x)w(x)\,dx G + ( d ) := ∫ e d g ( x ) w ( x ) d x と置くことができる。g w ≥ 0 gw\ge0 g w ≥ 0 であるから、§D1.17 定理 3.6 と§D1.17 命題 3.4 によりG − G_- G − は減少しG + G_+ G + は増加する。g w ≤ ∥ h ∥ ∞ w gw\le\lVert h\rVert_\infty w g w ≤ ∥ h ∥ ∞ w であるから、§D1.17 命題 3.4 によりG − ( c ) ≤ ∥ h ∥ ∞ ∫ c e w d x ≤ ∥ h ∥ ∞ μ 0 G_-(c)\le\lVert h\rVert_\infty\int_c^ew\,dx\le\lVert h\rVert_\infty\mu_0 G − ( c ) ≤ ∥ h ∥ ∞ ∫ c e w d x ≤ ∥ h ∥ ∞ μ 0 、G + ( d ) ≤ ∥ h ∥ ∞ ∫ e d w d x ≤ ∥ h ∥ ∞ μ 0 G_+(d)\le\lVert h\rVert_\infty\int_e^dw\,dx\le\lVert h\rVert_\infty\mu_0 G + ( d ) ≤ ∥ h ∥ ∞ ∫ e d w d x ≤ ∥ h ∥ ∞ μ 0 である。単調で有界な関数の片側極限は存在するから、lim c → a + G − ( c ) \lim_{c\to a+}G_-(c) lim c → a + G − ( c ) とlim d → b − G + ( d ) \lim_{d\to b-}G_+(d) lim d → b − G + ( d ) は有限な値として存在する。a < c < e < d < b a<c<e<d<b a < c < e < d < b ならば§D1.17 定理 3.6 により∫ c d g w d x = G − ( c ) + G + ( d ) \int_c^dgw\,dx=G_-(c)+G_+(d) ∫ c d g w d x = G − ( c ) + G + ( d ) であるから、g w gw g w の広義積分は存在する。
h ± : = max { ± h , 0 } h^\pm:=\max\{\pm h,0\} h ± := max { ± h , 0 } と置くと、h + h^+ h + 、h − h^- h − 、∣ h ∣ = h + + h − \lvert h\rvert=h^++h^- ∣ h ∣ = h + + h − は[ a , b ] [a,b] [ a , b ] 上の連続関数であり、その値は0 0 0 以上∥ h ∥ ∞ \lVert h\rVert_\infty ∥ h ∥ ∞ 以下である。上の議論をg = ∣ h ∣ g=\lvert h\rvert g = ∣ h ∣ に適用すると∣ h w ∣ = ∣ h ∣ w \lvert hw\rvert=\lvert h\rvert w ∣ h w ∣ = ∣ h ∣ w の広義積分が存在し、g = h + g=h^+ g = h + 、g = h − g=h^- g = h − に適用するとh + w h^+w h + w とh − w h^-w h − w の広義積分が存在する。h w = h + w − h − w hw=h^+w-h^-w h w = h + w − h − w であるから、§D1.17 定理 3.2 によりa < c < d < b a<c<d<b a < c < d < b に対して∫ c d h w d x = ∫ c d h + w d x − ∫ c d h − w d x \int_c^dhw\,dx=\int_c^dh^+w\,dx-\int_c^dh^-w\,dx ∫ c d h w d x = ∫ c d h + w d x − ∫ c d h − w d x であり、h w hw h w の広義積分も存在する。したがってh w hw h w の広義積分は絶対収束する。§D1.17 定理 3.5 と§D1.17 命題 3.4 により∣ ∫ c d h w d x ∣ ≤ ∫ c d ∣ h ∣ w d x ≤ ∥ h ∥ ∞ ∫ c d w d x ≤ ∥ h ∥ ∞ μ 0 \bigl|\int_c^dhw\,dx\bigr|\le\int_c^d\lvert h\rvert w\,dx\le\lVert h\rVert_\infty\int_c^dw\,dx\le\lVert h\rVert_\infty\mu_0 ∫ c d h w d x ≤ ∫ c d ∣ h ∣ w d x ≤ ∥ h ∥ ∞ ∫ c d w d x ≤ ∥ h ∥ ∞ μ 0 であり、c → a + c\to a+ c → a + 、d → b − d\to b- d → b − として評価を得る。
(2) を示す。h h h はx 0 x_0 x 0 で連続であるから、a < c < d < b a<c<d<b a < c < d < b を満たすc , d c,d c , d であって、任意のx ∈ [ c , d ] x\in[c,d] x ∈ [ c , d ] に対してh ( x ) ≥ h ( x 0 ) / 2 h(x)\ge h(x_0)/2 h ( x ) ≥ h ( x 0 ) /2 となるものが存在する。補題 1.2 (1) によりh w hw h w は[ c , d ] [c,d] [ c , d ] で Riemann 可積分であり、[ c , d ] [c,d] [ c , d ] 上でh w ≥ ( h ( x 0 ) / 2 ) w hw\ge(h(x_0)/2)\,w h w ≥ ( h ( x 0 ) /2 ) w であるから、§D1.17 命題 3.4 と補題 1.2 (2) により∫ c d h w d x ≥ ( h ( x 0 ) / 2 ) ∫ c d w d x > 0 \int_c^dhw\,dx\ge(h(x_0)/2)\int_c^dw\,dx>0 ∫ c d h w d x ≥ ( h ( x 0 ) /2 ) ∫ c d w d x > 0 である。h ≥ 0 h\ge0 h ≥ 0 であるから、a < c ′ ≤ c < d ≤ d ′ < b a<c'\le c<d\le d'<b a < c ′ ≤ c < d ≤ d ′ < b ならば∫ c ′ d ′ h w d x ≥ ∫ c d h w d x \int_{c'}^{d'}hw\,dx\ge\int_c^dhw\,dx ∫ c ′ d ′ h w d x ≥ ∫ c d h w d x である。c ′ → a + c'\to a+ c ′ → a + 、d ′ → b − d'\to b- d ′ → b − として∫ a b h w d x ≥ ( h ( x 0 ) / 2 ) ∫ c d w d x > 0 \int_a^bhw\,dx\ge(h(x_0)/2)\int_c^dw\,dx>0 ∫ a b h w d x ≥ ( h ( x 0 ) /2 ) ∫ c d w d x > 0 を得る。▨
定義 1.4. a < b a<b a < b を実数とし、C ( [ a , b ] ) C([a,b]) C ([ a , b ]) を[ a , b ] [a,b] [ a , b ] 上の連続な実数値関数の全体、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数とする。f , g ∈ C ( [ a , b ] ) f,g\in C([a,b]) f , g ∈ C ([ a , b ]) に対して、補題 1.3 (1) により有限な値
⟨ f , g ⟩ w : = ∫ a b f ( x ) g ( x ) w ( x ) d x , ∥ f ∥ w : = ⟨ f , f ⟩ w 1 / 2 \langle f,g\rangle_w:=\int_a^bf(x)g(x)w(x)\,dx,\qquad\lVert f\rVert_w:=\langle f,f\rangle_w^{1/2} ⟨ f , g ⟩ w := ∫ a b f ( x ) g ( x ) w ( x ) d x , ∥ f ∥ w := ⟨ f , f ⟩ w 1/2 を置く。⟨ ⋅ , ⋅ ⟩ w \langle\cdot,\cdot\rangle_w ⟨ ⋅ , ⋅ ⟩ w をw w w に関する 重み付き内積 (weighted inner product ) という。
命題 1.5. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、μ 0 : = ∫ a b w ( x ) d x \mu_0:=\int_a^bw(x)\,dx μ 0 := ∫ a b w ( x ) d x とする。⟨ ⋅ , ⋅ ⟩ w \langle\cdot,\cdot\rangle_w ⟨ ⋅ , ⋅ ⟩ w はC ( [ a , b ] ) C([a,b]) C ([ a , b ]) 上の内積であり、任意のf ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) に対して∥ f ∥ w ≤ μ 0 1 / 2 ∥ f ∥ ∞ \lVert f\rVert_w\le\mu_0^{1/2}\lVert f\rVert_\infty ∥ f ∥ w ≤ μ 0 1/2 ∥ f ∥ ∞ である。さらに、n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して、[ a , b ] [a,b] [ a , b ] への制限P n → C ( [ a , b ] ) \mathcal P_n\to C([a,b]) P n → C ([ a , b ]) は単射であり、1 , x , … , x n 1,x,\dots,x^n 1 , x , … , x n の制限はC ( [ a , b ] ) C([a,b]) C ([ a , b ]) の一次独立な元である。
証明. 広義積分は被積分関数について線形であるから、⟨ ⋅ , ⋅ ⟩ w \langle\cdot,\cdot\rangle_w ⟨ ⋅ , ⋅ ⟩ w は対称な双線形形式である。f ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) に対してf 2 ≥ 0 f^2\ge0 f 2 ≥ 0 であるから、補題 1.3 (1) をh = f 2 h=f^2 h = f 2 に適用して0 ≤ ⟨ f , f ⟩ w ≤ μ 0 ∥ f 2 ∥ ∞ = μ 0 ∥ f ∥ ∞ 2 0\le\langle f,f\rangle_w\le\mu_0\lVert f^2\rVert_\infty=\mu_0\lVert f\rVert_\infty^2 0 ≤ ⟨ f , f ⟩ w ≤ μ 0 ∥ f 2 ∥ ∞ = μ 0 ∥ f ∥ ∞ 2 を得る。f ≠ 0 f\ne0 f = 0 ならばf ( x 0 ) 2 > 0 f(x_0)^2>0 f ( x 0 ) 2 > 0 を満たすx 0 ∈ [ a , b ] x_0\in[a,b] x 0 ∈ [ a , b ] が存在し、補題 1.3 (2) により⟨ f , f ⟩ w > 0 \langle f,f\rangle_w>0 ⟨ f , f ⟩ w > 0 である。したがって⟨ ⋅ , ⋅ ⟩ w \langle\cdot,\cdot\rangle_w ⟨ ⋅ , ⋅ ⟩ w は正定値である。
p ∈ P n p\in\mathcal P_n p ∈ P n の[ a , b ] [a,b] [ a , b ] への制限が0 0 0 ならば、p p p は[ a , b ] [a,b] [ a , b ] の相異なるn + 1 n+1 n + 1 点で0 0 0 になるから、§E20.12 補題 1.2 (2) をk i = 1 k_i=1 k i = 1 、K = n + 1 K=n+1 K = n + 1 、m = n m=n m = n として適用してp = 0 p=0 p = 0 を得る。1 , x , … , x n 1,x,\dots,x^n 1 , x , … , x n はP n \mathcal P_n P n の基底であるから、その制限は一次独立である。▨
定理 1.6. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 、f ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) とする。
p ∈ P n p\in\mathcal P_n p ∈ P n であって、任意のq ∈ P n q\in\mathcal P_n q ∈ P n に対して∥ f − p ∥ w ≤ ∥ f − q ∥ w \lVert f-p\rVert_w\le\lVert f-q\rVert_w ∥ f − p ∥ w ≤ ∥ f − q ∥ w を満たすものがただ一つ存在する。
p ∈ P n p\in\mathcal P_n p ∈ P n が(1) の元であることと、任意のq ∈ P n q\in\mathcal P_n q ∈ P n に対して⟨ f − p , q ⟩ w = 0 \langle f-p,q\rangle_w=0 ⟨ f − p , q ⟩ w = 0 であることは同値である。
p p p を(1) の元とすると、任意のq ∈ P n q\in\mathcal P_n q ∈ P n に対して∥ f − q ∥ w 2 = ∥ f − p ∥ w 2 + ∥ p − q ∥ w 2 \lVert f-q\rVert_w^2=\lVert f-p\rVert_w^2+\lVert p-q\rVert_w^2 ∥ f − q ∥ w 2 = ∥ f − p ∥ w 2 + ∥ p − q ∥ w 2 である。
φ 0 , … , φ n \varphi_0,\dots,\varphi_n φ 0 , … , φ n をP n \mathcal P_n P n の基底とする。c 0 , … , c n ∈ R c_0,\dots,c_n\in\R c 0 , … , c n ∈ R について、∑ j = 0 n c j φ j \sum_{j=0}^nc_j\varphi_j ∑ j = 0 n c j φ j が(1) の元であることと、
∑ j = 0 n ⟨ φ j , φ i ⟩ w c j = ⟨ f , φ i ⟩ w ( 0 ≤ i ≤ n ) \sum_{j=0}^n\langle\varphi_j,\varphi_i\rangle_wc_j=\langle f,\varphi_i\rangle_w\qquad(0\le i\le n) j = 0 ∑ n ⟨ φ j , φ i ⟩ w c j = ⟨ f , φ i ⟩ w ( 0 ≤ i ≤ n )
が成り立つことは同値である。
(4) の基底がi ≠ j i\ne j i = j ならば⟨ φ i , φ j ⟩ w = 0 \langle\varphi_i,\varphi_j\rangle_w=0 ⟨ φ i , φ j ⟩ w = 0 を満たすとき、(1) の元はp = ∑ i = 0 n ( ⟨ f , φ i ⟩ w / ∥ φ i ∥ w 2 ) φ i p=\sum_{i=0}^n\bigl(\langle f,\varphi_i\rangle_w/\lVert\varphi_i\rVert_w^2\bigr)\varphi_i p = ∑ i = 0 n ( ⟨ f , φ i ⟩ w / ∥ φ i ∥ w 2 ) φ i であり、
∥ f − p ∥ w 2 = ∥ f ∥ w 2 − ∑ i = 0 n ⟨ f , φ i ⟩ w 2 ∥ φ i ∥ w 2 \lVert f-p\rVert_w^2=\lVert f\rVert_w^2-\sum_{i=0}^n\frac{\langle f,\varphi_i\rangle_w^2}{\lVert\varphi_i\rVert_w^2} ∥ f − p ∥ w 2 = ∥ f ∥ w 2 − i = 0 ∑ n ∥ φ i ∥ w 2 ⟨ f , φ i ⟩ w 2
である。
証明. 命題 1.5 により1 , x , … , x n 1,x,\dots,x^n 1 , x , … , x n は内積空間( C ( [ a , b ] ) , ⟨ ⋅ , ⋅ ⟩ w ) (C([a,b]),\langle\cdot,\cdot\rangle_w) ( C ([ a , b ]) , ⟨ ⋅ , ⋅ ⟩ w ) の一次独立な元であり、これに§D3.14 定理 2.1 を適用すると、P n \mathcal P_n P n を張る正規直交系e 0 , … , e n e_0,\dots,e_n e 0 , … , e n を得る。p : = ∑ k = 0 n ⟨ f , e k ⟩ w e k p:=\sum_{k=0}^n\langle f,e_k\rangle_we_k p := ∑ k = 0 n ⟨ f , e k ⟩ w e k と置く。§D3.14 命題 3.2 により、f − p f-p f − p はP n \mathcal P_n P n のすべての元に直交し、任意のq ∈ P n q\in\mathcal P_n q ∈ P n に対して∥ f − p ∥ w ≤ ∥ f − q ∥ w \lVert f-p\rVert_w\le\lVert f-q\rVert_w ∥ f − p ∥ w ≤ ∥ f − q ∥ w であって、等号はq = p q=p q = p のときに限る。これで(1) は示された。
(2) を示す。(1) の元p p p についてf − p f-p f − p がP n \mathcal P_n P n に直交することは上で示した。逆にp ′ ∈ P n p'\in\mathcal P_n p ′ ∈ P n が任意のq ∈ P n q\in\mathcal P_n q ∈ P n に対して⟨ f − p ′ , q ⟩ w = 0 \langle f-p',q\rangle_w=0 ⟨ f − p ′ , q ⟩ w = 0 を満たすとする。p − p ′ ∈ P n p-p'\in\mathcal P_n p − p ′ ∈ P n であるから、∥ p − p ′ ∥ w 2 = ⟨ f − p ′ , p − p ′ ⟩ w − ⟨ f − p , p − p ′ ⟩ w = 0 \lVert p-p'\rVert_w^2=\langle f-p',p-p'\rangle_w-\langle f-p,p-p'\rangle_w=0 ∥ p − p ′ ∥ w 2 = ⟨ f − p ′ , p − p ′ ⟩ w − ⟨ f − p , p − p ′ ⟩ w = 0 であり、p ′ = p p'=p p ′ = p である。
(3) を示す。f − q = ( f − p ) + ( p − q ) f-q=(f-p)+(p-q) f − q = ( f − p ) + ( p − q ) であり、p − q ∈ P n p-q\in\mathcal P_n p − q ∈ P n であるから(2) により⟨ f − p , p − q ⟩ w = 0 \langle f-p,p-q\rangle_w=0 ⟨ f − p , p − q ⟩ w = 0 である。∥ f − q ∥ w 2 \lVert f-q\rVert_w^2 ∥ f − q ∥ w 2 を展開して等式を得る。
(4) を示す。p = ∑ j c j φ j p=\sum_jc_j\varphi_j p = ∑ j c j φ j と置く。P n \mathcal P_n P n はφ 0 , … , φ n \varphi_0,\dots,\varphi_n φ 0 , … , φ n で張られるから、(2) の条件は⟨ f − p , φ i ⟩ w = 0 \langle f-p,\varphi_i\rangle_w=0 ⟨ f − p , φ i ⟩ w = 0 (0 ≤ i ≤ n 0\le i\le n 0 ≤ i ≤ n )と同値であり、これは⟨ f , φ i ⟩ w = ∑ j c j ⟨ φ j , φ i ⟩ w \langle f,\varphi_i\rangle_w=\sum_jc_j\langle\varphi_j,\varphi_i\rangle_w ⟨ f , φ i ⟩ w = ∑ j c j ⟨ φ j , φ i ⟩ w (0 ≤ i ≤ n 0\le i\le n 0 ≤ i ≤ n )と同値である。
(5) を示す。基底の元は0 0 0 でないから∥ φ i ∥ w 2 > 0 \lVert\varphi_i\rVert_w^2>0 ∥ φ i ∥ w 2 > 0 である。(4) の方程式は∥ φ i ∥ w 2 c i = ⟨ f , φ i ⟩ w \lVert\varphi_i\rVert_w^2c_i=\langle f,\varphi_i\rangle_w ∥ φ i ∥ w 2 c i = ⟨ f , φ i ⟩ w となり、係数の式を得る。(3) をq = 0 q=0 q = 0 として適用すると∥ f ∥ w 2 = ∥ f − p ∥ w 2 + ∥ p ∥ w 2 \lVert f\rVert_w^2=\lVert f-p\rVert_w^2+\lVert p\rVert_w^2 ∥ f ∥ w 2 = ∥ f − p ∥ w 2 + ∥ p ∥ w 2 であり、φ i \varphi_i φ i が互いに直交するから∥ p ∥ w 2 = ∑ i c i 2 ∥ φ i ∥ w 2 = ∑ i ⟨ f , φ i ⟩ w 2 / ∥ φ i ∥ w 2 \lVert p\rVert_w^2=\sum_ic_i^2\lVert\varphi_i\rVert_w^2=\sum_i\langle f,\varphi_i\rangle_w^2/\lVert\varphi_i\rVert_w^2 ∥ p ∥ w 2 = ∑ i c i 2 ∥ φ i ∥ w 2 = ∑ i ⟨ f , φ i ⟩ w 2 / ∥ φ i ∥ w 2 である。▨
定義 1.7. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 とする。f ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) に対して、定理 1.6 (1) のただ一つの元をΠ n w f \Pi_n^wf Π n w f と書き、f f f のP n \mathcal P_n P n によるw w w に関する 二乗平均最良近似 (best approximation in the weighted mean square ) という。
2 直交多項式と三項漸化式
定義 2.1. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数とする。
k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 とする。x k x^k x k の係数が1 1 1 であるk k k 次の実係数多項式をk k k 次の モニック多項式 (monic polynomial ) という。
実係数多項式の列( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 が、各k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 についてπ k \pi_k π k はk k k 次のモニック多項式であり、k ≥ 1 k\ge1 k ≥ 1 ならば任意のq ∈ P k − 1 q\in\mathcal P_{k-1} q ∈ P k − 1 に対して⟨ π k , q ⟩ w = 0 \langle\pi_k,q\rangle_w=0 ⟨ π k , q ⟩ w = 0 であるという条件を満たすとき、( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 をw w w に関する モニック直交多項式列 (monic orthogonal polynomials ) といい、π k \pi_k π k をw w w に関するk k k 次の 直交多項式 (orthogonal polynomial ) という。
命題 2.2. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数とする。
w w w に関するモニック直交多項式列( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 はただ一つ存在し、π 0 = 1 \pi_0=1 π 0 = 1 、π k = x k − Π k − 1 w ( x k ) \pi_k=x^k-\Pi_{k-1}^w(x^k) π k = x k − Π k − 1 w ( x k ) (k ≥ 1 k\ge1 k ≥ 1 )である。
n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 とする。π 0 , … , π n \pi_0,\dots,\pi_n π 0 , … , π n はP n \mathcal P_n P n の基底であり、i ≠ j i\ne j i = j ならば⟨ π i , π j ⟩ w = 0 \langle\pi_i,\pi_j\rangle_w=0 ⟨ π i , π j ⟩ w = 0 である。任意のf ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) に対して
Π n w f = ∑ k = 0 n ⟨ f , π k ⟩ w ∥ π k ∥ w 2 π k \Pi_n^wf=\sum_{k=0}^n\frac{\langle f,\pi_k\rangle_w}{\lVert\pi_k\rVert_w^2}\,\pi_k Π n w f = k = 0 ∑ n ∥ π k ∥ w 2 ⟨ f , π k ⟩ w π k
であり、n ≥ 1 n\ge1 n ≥ 1 ならばΠ n w f = Π n − 1 w f + ( ⟨ f , π n ⟩ w / ∥ π n ∥ w 2 ) π n \Pi_n^wf=\Pi_{n-1}^wf+\bigl(\langle f,\pi_n\rangle_w/\lVert\pi_n\rVert_w^2\bigr)\pi_n Π n w f = Π n − 1 w f + ( ⟨ f , π n ⟩ w / ∥ π n ∥ w 2 ) π n である。
n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 とし、p ∈ P n p\in\mathcal P_n p ∈ P n が任意のq ∈ P n − 1 q\in\mathcal P_{n-1} q ∈ P n − 1 に対して⟨ p , q ⟩ w = 0 \langle p,q\rangle_w=0 ⟨ p , q ⟩ w = 0 を満たすとする。p p p のx n x^n x n の係数をc c c とするとp = c π n p=c\,\pi_n p = c π n である。
証明. (1) を示す。0 0 0 次のモニック多項式は1 1 1 だけであり、k = 0 k=0 k = 0 では直交の条件は課されない。k ≥ 1 k\ge1 k ≥ 1 とする。k k k 次のモニック多項式は、r ∈ P k − 1 r\in\mathcal P_{k-1} r ∈ P k − 1 によってx k − r x^k-r x k − r とただ一通りに書かれる。定理 1.6 (2) をf = x k f=x^k f = x k 、n = k − 1 n=k-1 n = k − 1 として適用すると、任意のq ∈ P k − 1 q\in\mathcal P_{k-1} q ∈ P k − 1 に対して⟨ x k − r , q ⟩ w = 0 \langle x^k-r,q\rangle_w=0 ⟨ x k − r , q ⟩ w = 0 であることとr = Π k − 1 w ( x k ) r=\Pi_{k-1}^w(x^k) r = Π k − 1 w ( x k ) であることは同値である。したがってπ k \pi_k π k はただ一つ存在し、π k = x k − Π k − 1 w ( x k ) \pi_k=x^k-\Pi_{k-1}^w(x^k) π k = x k − Π k − 1 w ( x k ) である。
(2) を示す。i < j i<j i < j ならばπ i ∈ P j − 1 \pi_i\in\mathcal P_{j-1} π i ∈ P j − 1 であるから⟨ π j , π i ⟩ w = 0 \langle\pi_j,\pi_i\rangle_w=0 ⟨ π j , π i ⟩ w = 0 である。π i ≠ 0 \pi_i\ne0 π i = 0 であるから、∑ i c i π i = 0 \sum_ic_i\pi_i=0 ∑ i c i π i = 0 とπ j \pi_j π j の内積をとってc j ∥ π j ∥ w 2 = 0 c_j\lVert\pi_j\rVert_w^2=0 c j ∥ π j ∥ w 2 = 0 、すなわちc j = 0 c_j=0 c j = 0 を得る。したがってπ 0 , … , π n \pi_0,\dots,\pi_n π 0 , … , π n はn + 1 n+1 n + 1 次元のP n \mathcal P_n P n の一次独立なn + 1 n+1 n + 1 個の元であり、基底である。展開式は定理 1.6 (5) であり、n ≥ 1 n\ge1 n ≥ 1 のときの等式はn n n とn − 1 n-1 n − 1 の展開式の差である。
(3) を示す。p − c π n p-c\pi_n p − c π n のx n x^n x n の係数は0 0 0 であるからp − c π n ∈ P n − 1 p-c\pi_n\in\mathcal P_{n-1} p − c π n ∈ P n − 1 である。p p p とπ n \pi_n π n はともにP n − 1 \mathcal P_{n-1} P n − 1 のすべての元に直交するから、∥ p − c π n ∥ w 2 = ⟨ p − c π n , p − c π n ⟩ w = 0 \lVert p-c\pi_n\rVert_w^2=\langle p-c\pi_n,p-c\pi_n\rangle_w=0 ∥ p − c π n ∥ w 2 = ⟨ p − c π n , p − c π n ⟩ w = 0 であり、p = c π n p=c\pi_n p = c π n である。▨
定理 2.3. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 をw w w に関するモニック直交多項式列とする。π − 1 : = 0 \pi_{-1}:=0 π − 1 := 0 と置き、k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 に対して多項式x π k x\pi_k x π k をx ⋅ π k ( x ) x\cdot\pi_k(x) x ⋅ π k ( x ) で定め、
α k : = ⟨ x π k , π k ⟩ w ∥ π k ∥ w 2 ( k ∈ N ≥ 0 ) , β 0 : = 0 , β k : = ∥ π k ∥ w 2 ∥ π k − 1 ∥ w 2 ( k ∈ N ≥ 1 ) \alpha_k:=\frac{\langle x\pi_k,\pi_k\rangle_w}{\lVert\pi_k\rVert_w^2}\quad(k\in\N),\qquad\beta_0:=0,\qquad\beta_k:=\frac{\lVert\pi_k\rVert_w^2}{\lVert\pi_{k-1}\rVert_w^2}\quad(k\in\NN) α k := ∥ π k ∥ w 2 ⟨ x π k , π k ⟩ w ( k ∈ N ≥ 0 ) , β 0 := 0 , β k := ∥ π k − 1 ∥ w 2 ∥ π k ∥ w 2 ( k ∈ N ≥ 1 ) と置く。
任意のk ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 に対してπ k + 1 = ( x − α k ) π k − β k π k − 1 \pi_{k+1}=(x-\alpha_k)\pi_k-\beta_k\pi_{k-1} π k + 1 = ( x − α k ) π k − β k π k − 1 であり、k ≥ 1 k\ge1 k ≥ 1 ならばβ k > 0 \beta_k>0 β k > 0 である。
k ∈ N ≥ 1 k\in\NN k ∈ N ≥ 1 とする。実数α , β \alpha,\beta α , β がπ k + 1 = ( x − α ) π k − β π k − 1 \pi_{k+1}=(x-\alpha)\pi_k-\beta\pi_{k-1} π k + 1 = ( x − α ) π k − β π k − 1 を満たすならば、α = α k \alpha=\alpha_k α = α k 、β = β k \beta=\beta_k β = β k である。
証明. (1) を示す。k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 とし、r : = x π k − π k + 1 r:=x\pi_k-\pi_{k+1} r := x π k − π k + 1 と置く。x π k x\pi_k x π k とπ k + 1 \pi_{k+1} π k + 1 はともにk + 1 k+1 k + 1 次のモニック多項式であるからr ∈ P k r\in\mathcal P_k r ∈ P k であり、命題 2.2 (2) によりr = ∑ j = 0 k d j π j r=\sum_{j=0}^kd_j\pi_j r = ∑ j = 0 k d j π j と書かれる。π i \pi_i π i (0 ≤ i ≤ k 0\le i\le k 0 ≤ i ≤ k )との内積をとるとd i ∥ π i ∥ w 2 = ⟨ r , π i ⟩ w d_i\lVert\pi_i\rVert_w^2=\langle r,\pi_i\rangle_w d i ∥ π i ∥ w 2 = ⟨ r , π i ⟩ w である。π k + 1 \pi_{k+1} π k + 1 はP k \mathcal P_k P k に直交し、被積分関数の等式x π k ( x ) π i ( x ) w ( x ) = π k ( x ) ⋅ x π i ( x ) w ( x ) x\pi_k(x)\pi_i(x)w(x)=\pi_k(x)\cdot x\pi_i(x)w(x) x π k ( x ) π i ( x ) w ( x ) = π k ( x ) ⋅ x π i ( x ) w ( x ) により
⟨ r , π i ⟩ w = ⟨ x π k , π i ⟩ w = ⟨ π k , x π i ⟩ w \langle r,\pi_i\rangle_w=\langle x\pi_k,\pi_i\rangle_w=\langle\pi_k,x\pi_i\rangle_w ⟨ r , π i ⟩ w = ⟨ x π k , π i ⟩ w = ⟨ π k , x π i ⟩ w である。i ≤ k − 2 i\le k-2 i ≤ k − 2 ならばx π i ∈ P k − 1 x\pi_i\in\mathcal P_{k-1} x π i ∈ P k − 1 であるからd i = 0 d_i=0 d i = 0 である。i = k i=k i = k ならばd k = α k d_k=\alpha_k d k = α k である。k ≥ 1 k\ge1 k ≥ 1 かつi = k − 1 i=k-1 i = k − 1 ならば、x π k − 1 − π k ∈ P k − 1 x\pi_{k-1}-\pi_k\in\mathcal P_{k-1} x π k − 1 − π k ∈ P k − 1 であるから⟨ π k , x π k − 1 ⟩ w = ∥ π k ∥ w 2 \langle\pi_k,x\pi_{k-1}\rangle_w=\lVert\pi_k\rVert_w^2 ⟨ π k , x π k − 1 ⟩ w = ∥ π k ∥ w 2 であり、d k − 1 = β k d_{k-1}=\beta_k d k − 1 = β k である。以上によりr = α k π k + β k π k − 1 r=\alpha_k\pi_k+\beta_k\pi_{k-1} r = α k π k + β k π k − 1 であり(k = 0 k=0 k = 0 ではπ − 1 = 0 \pi_{-1}=0 π − 1 = 0 によりr = α 0 π 0 r=\alpha_0\pi_0 r = α 0 π 0 )、漸化式を得る。k ≥ 1 k\ge1 k ≥ 1 ならば∥ π k ∥ w 2 > 0 \lVert\pi_k\rVert_w^2>0 ∥ π k ∥ w 2 > 0 、∥ π k − 1 ∥ w 2 > 0 \lVert\pi_{k-1}\rVert_w^2>0 ∥ π k − 1 ∥ w 2 > 0 であるからβ k > 0 \beta_k>0 β k > 0 である。
(2) を示す。(1) の等式との差をとると( α k − α ) π k = ( β − β k ) π k − 1 (\alpha_k-\alpha)\pi_k=(\beta-\beta_k)\pi_{k-1} ( α k − α ) π k = ( β − β k ) π k − 1 である。命題 2.2 (2) によりπ k − 1 \pi_{k-1} π k − 1 とπ k \pi_k π k は一次独立であるから、α = α k \alpha=\alpha_k α = α k 、β = β k \beta=\beta_k β = β k である。▨
3 直交多項式の零点
定理 3.1. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 とし、π n \pi_n π n をw w w に関するn n n 次の直交多項式とする。このとき、a < z 1 < ⋯ < z n < b a<z_1<\dots<z_n<b a < z 1 < ⋯ < z n < b を満たす実数z 1 , … , z n z_1,\dots,z_n z 1 , … , z n が存在してπ n = ∏ i = 1 n ( x − z i ) \pi_n=\prod_{i=1}^n(x-z_i) π n = ∏ i = 1 n ( x − z i ) である。特に、π n \pi_n π n の複素数の零点はz 1 , … , z n z_1,\dots,z_n z 1 , … , z n だけであり、π n ′ ( z i ) = ∏ j ≠ i ( z i − z j ) ≠ 0 \pi_n'(z_i)=\prod_{j\ne i}(z_i-z_j)\ne0 π n ′ ( z i ) = ∏ j = i ( z i − z j ) = 0 (1 ≤ i ≤ n 1\le i\le n 1 ≤ i ≤ n )である。
証明. π n ≠ 0 \pi_n\ne0 π n = 0 であるから、§E20.12 補題 1.2 (2) によりπ n \pi_n π n の( a , b ) (a,b) ( a , b ) における零点は高々n n n 個である。それらをy 1 < ⋯ < y r y_1<\dots<y_r y 1 < ⋯ < y r (0 ≤ r ≤ n 0\le r\le n 0 ≤ r ≤ n )とし、y 0 : = a y_0:=a y 0 := a 、y r + 1 : = b y_{r+1}:=b y r + 1 := b と置く。0 ≤ j ≤ r 0\le j\le r 0 ≤ j ≤ r について、π n \pi_n π n は開区間( y j , y j + 1 ) (y_j,y_{j+1}) ( y j , y j + 1 ) に零点をもたないから、中間値の定理によりこの区間で一定の符号ε j ∈ { 1 , − 1 } \varepsilon_j\in\{1,-1\} ε j ∈ { 1 , − 1 } をとる。1 ≤ j ≤ r 1\le j\le r 1 ≤ j ≤ r のうちε j − 1 ≠ ε j \varepsilon_{j-1}\ne\varepsilon_j ε j − 1 = ε j を満たすj j j に対応するy j y_j y j を小さい順にz 1 < ⋯ < z m z_1<\dots<z_m z 1 < ⋯ < z m とし、q : = ∏ i = 1 m ( x − z i ) q:=\prod_{i=1}^m(x-z_i) q := ∏ i = 1 m ( x − z i ) (m = 0 m=0 m = 0 のときq : = 1 q:=1 q := 1 )と置く。
1 ≤ j ≤ r 1\le j\le r 1 ≤ j ≤ r について、y j y_j y j をπ n \pi_n π n のk j k_j k j 重の零点とし、π n = ( x − y j ) k j g j \pi_n=(x-y_j)^{k_j}g_j π n = ( x − y j ) k j g j 、g j ( y j ) ≠ 0 g_j(y_j)\ne0 g j ( y j ) = 0 と書くと、g j g_j g j はy j y_j y j を含む開区間で符号が一定であるから、ε j − 1 ≠ ε j \varepsilon_{j-1}\ne\varepsilon_j ε j − 1 = ε j であることはk j k_j k j が奇数であることと同値である。したがってπ n \pi_n π n は( a , b ) (a,b) ( a , b ) で奇数重複の零点z 1 , … , z m z_1,\dots,z_m z 1 , … , z m でだけ符号を変える。q q q は各z i z_i z i で一位の零点をもち、q / ( x − z i ) = ∏ l ≠ i ( x − z l ) q/(x-z_i)=\prod_{l\ne i}(x-z_l) q / ( x − z i ) = ∏ l = i ( x − z l ) はz i z_i z i で0 0 0 でないから、q q q は( a , b ) (a,b) ( a , b ) でz 1 , … , z m z_1,\dots,z_m z 1 , … , z m でだけ0 0 0 になり、各z i z_i z i で符号を変える。π n \pi_n π n とq q q は( a , b ) (a,b) ( a , b ) の同じ点z 1 , … , z m z_1,\dots,z_m z 1 , … , z m でだけ符号を変えるので、π n q \pi_nq π n q は( a , b ) (a,b) ( a , b ) からy 1 , … , y r y_1,\dots,y_r y 1 , … , y r を除いた集合で一定の符号ε ∈ { 1 , − 1 } \varepsilon\in\{1,-1\} ε ∈ { 1 , − 1 } をとり、h : = ε π n q h:=\varepsilon\,\pi_nq h := ε π n q はこの集合で正である。h h h は[ a , b ] [a,b] [ a , b ] で連続であるから、任意のx ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] に対してh ( x ) ≥ 0 h(x)\ge0 h ( x ) ≥ 0 であり、h h h は( y 0 , y 1 ) (y_0,y_1) ( y 0 , y 1 ) で正である。補題 1.3 (2) により⟨ π n , q ⟩ w = ε ∫ a b h ( x ) w ( x ) d x ≠ 0 \langle\pi_n,q\rangle_w=\varepsilon\int_a^bh(x)w(x)\,dx\ne0 ⟨ π n , q ⟩ w = ε ∫ a b h ( x ) w ( x ) d x = 0 である。
m ≤ n − 1 m\le n-1 m ≤ n − 1 と仮定するとq ∈ P n − 1 q\in\mathcal P_{n-1} q ∈ P n − 1 であり、直交多項式の定義により⟨ π n , q ⟩ w = 0 \langle\pi_n,q\rangle_w=0 ⟨ π n , q ⟩ w = 0 である。これは⟨ π n , q ⟩ w ≠ 0 \langle\pi_n,q\rangle_w\ne0 ⟨ π n , q ⟩ w = 0 と両立しないから、m ≥ n m\ge n m ≥ n である。z 1 , … , z n z_1,\dots,z_n z 1 , … , z n は( a , b ) (a,b) ( a , b ) の相異なる零点であるから、§E20.12 補題 1.2 (3) をk i = 1 k_i=1 k i = 1 、K = n K=n K = n として適用すると、π n \pi_n π n のx n x^n x n の係数が1 1 1 であることからπ n = ∏ i = 1 n ( x − z i ) \pi_n=\prod_{i=1}^n(x-z_i) π n = ∏ i = 1 n ( x − z i ) を得る。この積表示から、π n \pi_n π n の複素数の零点はz 1 , … , z n z_1,\dots,z_n z 1 , … , z n だけであり、積の微分によりπ n ′ ( z i ) = ∏ j ≠ i ( z i − z j ) \pi_n'(z_i)=\prod_{j\ne i}(z_i-z_j) π n ′ ( z i ) = ∏ j = i ( z i − z j ) である。▨
例 3.2. ( 0 , 1 ) (0,1) ( 0 , 1 ) 上のw ( x ) : = x w(x):=x w ( x ) := x は重み関数であり、∫ 0 1 x k w ( x ) d x = 1 / ( k + 2 ) \int_0^1x^kw(x)\,dx=1/(k+2) ∫ 0 1 x k w ( x ) d x = 1/ ( k + 2 ) (k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 )、μ 0 = 1 / 2 \mu_0=1/2 μ 0 = 1/2 である。定理 2.3 (1) により
α 0 = ⟨ x , 1 ⟩ w ∥ 1 ∥ w 2 = 1 / 3 1 / 2 = 2 3 , π 1 = x − 2 3 , ∥ π 1 ∥ w 2 = 1 4 − 4 3 ⋅ 1 3 + 4 9 ⋅ 1 2 = 1 36 \alpha_0=\frac{\langle x,1\rangle_w}{\lVert1\rVert_w^2}=\frac{1/3}{1/2}=\frac23,\qquad\pi_1=x-\frac23,\qquad\lVert\pi_1\rVert_w^2=\frac14-\frac43\cdot\frac13+\frac49\cdot\frac12=\frac1{36} α 0 = ∥ 1 ∥ w 2 ⟨ x , 1 ⟩ w = 1/2 1/3 = 3 2 , π 1 = x − 3 2 , ∥ π 1 ∥ w 2 = 4 1 − 3 4 ⋅ 3 1 + 9 4 ⋅ 2 1 = 36 1 であり、β 1 = ( 1 / 36 ) / ( 1 / 2 ) = 1 / 18 \beta_1=(1/36)/(1/2)=1/18 β 1 = ( 1/36 ) / ( 1/2 ) = 1/18 、⟨ x π 1 , π 1 ⟩ w = 1 / 5 − ( 4 / 3 ) ( 1 / 4 ) + ( 4 / 9 ) ( 1 / 3 ) = 2 / 135 \langle x\pi_1,\pi_1\rangle_w=1/5-(4/3)(1/4)+(4/9)(1/3)=2/135 ⟨ x π 1 , π 1 ⟩ w = 1/5 − ( 4/3 ) ( 1/4 ) + ( 4/9 ) ( 1/3 ) = 2/135 、α 1 = ( 2 / 135 ) ⋅ 36 = 8 / 15 \alpha_1=(2/135)\cdot36=8/15 α 1 = ( 2/135 ) ⋅ 36 = 8/15 である。したがって
π 2 = ( x − 8 15 ) ( x − 2 3 ) − 1 18 = x 2 − 6 5 x + 3 10 \pi_2=\Bigl(x-\frac8{15}\Bigr)\Bigl(x-\frac23\Bigr)-\frac1{18}=x^2-\frac65x+\frac3{10} π 2 = ( x − 15 8 ) ( x − 3 2 ) − 18 1 = x 2 − 5 6 x + 10 3 である。π 2 \pi_2 π 2 の零点は3 / 5 ± 6 / 10 3/5\pm\sqrt6/10 3/5 ± 6 /10 、すなわち約0.3551 0.3551 0.3551 と約0.8449 0.8449 0.8449 であり、定理 3.1 のとおり( 0 , 1 ) (0,1) ( 0 , 1 ) の相異なる二点である。
4 Legendre 多項式と Chebyshev の重み
命題 4.1. ( − 1 , 1 ) (-1,1) ( − 1 , 1 ) 上の重み関数w P ( x ) : = 1 w_P(x):=1 w P ( x ) := 1 に関するモニック直交多項式列を( π n ) n ∈ N ≥ 0 (\pi_n)_{n\in\N} ( π n ) n ∈ N ≥ 0 、三項漸化式の係数をα n , β n \alpha_n,\beta_n α n , β n とする。n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対してP n P_n P n をn n n 次の Legendre 多項式とし、c n c_n c n をP n P_n P n のx n x^n x n の係数とする。このとき、任意のn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して
c n = ( 2 n ) ! 2 n ( n ! ) 2 , π n = P n c n , ∥ π n ∥ w P 2 = 2 ( 2 n + 1 ) c n 2 , α n = 0 c_n=\frac{(2n)!}{2^n(n!)^2},\qquad\pi_n=\frac{P_n}{c_n},\qquad\lVert\pi_n\rVert_{w_P}^2=\frac2{(2n+1)c_n^2},\qquad\alpha_n=0 c n = 2 n ( n ! ) 2 ( 2 n )! , π n = c n P n , ∥ π n ∥ w P 2 = ( 2 n + 1 ) c n 2 2 , α n = 0 であり、n ≥ 1 n\ge1 n ≥ 1 ならばβ n = n 2 / ( 4 n 2 − 1 ) \beta_n=n^2/(4n^2-1) β n = n 2 / ( 4 n 2 − 1 ) である。さらに、任意のf ∈ C ( [ − 1 , 1 ] ) f\in C([-1,1]) f ∈ C ([ − 1 , 1 ]) とn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して
Π n w P f = ∑ k = 0 n 2 k + 1 2 ( ∫ − 1 1 f ( x ) P k ( x ) d x ) P k \Pi_n^{w_P}f=\sum_{k=0}^n\frac{2k+1}2\Bigl(\int_{-1}^1f(x)P_k(x)\,dx\Bigr)P_k Π n w P f = k = 0 ∑ n 2 2 k + 1 ( ∫ − 1 1 f ( x ) P k ( x ) d x ) P k である。
証明. §E10.15 定理 2.1 によりP n = 1 2 n n ! d n d x n ( x 2 − 1 ) n P_n=\frac1{2^nn!}\frac{d^n}{dx^n}(x^2-1)^n P n = 2 n n ! 1 d x n d n ( x 2 − 1 ) n であり、n = 0 , 1 n=0,1 n = 0 , 1 としてP 0 = 1 P_0=1 P 0 = 1 、P 1 = x P_1=x P 1 = x を得る。二項展開により( x 2 − 1 ) n − x 2 n (x^2-1)^n-x^{2n} ( x 2 − 1 ) n − x 2 n の次数は2 n − 2 2n-2 2 n − 2 以下であるから、そのn n n 回微分の次数はn − 2 n-2 n − 2 以下であり、d n d x n x 2 n = ( 2 n ) ! n ! x n \frac{d^n}{dx^n}x^{2n}=\frac{(2n)!}{n!}x^n d x n d n x 2 n = n ! ( 2 n )! x n である。したがってP n P_n P n はn n n 次の多項式であり、
c n = 1 2 n n ! ⋅ ( 2 n ) ! n ! = ( 2 n ) ! 2 n ( n ! ) 2 c_n=\frac1{2^nn!}\cdot\frac{(2n)!}{n!}=\frac{(2n)!}{2^n(n!)^2} c n = 2 n n ! 1 ⋅ n ! ( 2 n )! = 2 n ( n ! ) 2 ( 2 n )! である。特にc 0 = c 1 = 1 c_0=c_1=1 c 0 = c 1 = 1 であり、任意のn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対してc n + 1 / c n = ( 2 n + 2 ) ( 2 n + 1 ) / ( 2 ( n + 1 ) 2 ) = ( 2 n + 1 ) / ( n + 1 ) c_{n+1}/c_n=(2n+2)(2n+1)/(2(n+1)^2)=(2n+1)/(n+1) c n + 1 / c n = ( 2 n + 2 ) ( 2 n + 1 ) / ( 2 ( n + 1 ) 2 ) = ( 2 n + 1 ) / ( n + 1 ) 、すなわち( n + 1 ) c n + 1 = ( 2 n + 1 ) c n (n+1)c_{n+1}=(2n+1)c_n ( n + 1 ) c n + 1 = ( 2 n + 1 ) c n である。
n ≥ 1 n\ge1 n ≥ 1 とする。§E10.15 定理 3.1 により、P n P_n P n は次数がn n n 未満の任意の多項式と⟨ ⋅ , ⋅ ⟩ w P \langle\cdot,\cdot\rangle_{w_P} ⟨ ⋅ , ⋅ ⟩ w P について直交する。命題 2.2 (3) によりP n = c n π n P_n=c_n\pi_n P n = c n π n である。n = 0 n=0 n = 0 ではP 0 = 1 = π 0 P_0=1=\pi_0 P 0 = 1 = π 0 、c 0 = 1 c_0=1 c 0 = 1 である。§E10.15 命題 4.2 により∥ π n ∥ w P 2 = ∥ P n ∥ w P 2 / c n 2 = 2 / ( ( 2 n + 1 ) c n 2 ) \lVert\pi_n\rVert_{w_P}^2=\lVert P_n\rVert_{w_P}^2/c_n^2=2/((2n+1)c_n^2) ∥ π n ∥ w P 2 = ∥ P n ∥ w P 2 / c n 2 = 2/ (( 2 n + 1 ) c n 2 ) である。
π 1 = x = ( x − 0 ) π 0 \pi_1=x=(x-0)\pi_0 π 1 = x = ( x − 0 ) π 0 であるから、定理 2.3 (1) のk = 0 k=0 k = 0 の等式と比べてα 0 = 0 \alpha_0=0 α 0 = 0 である。n ≥ 1 n\ge1 n ≥ 1 とする。§E10.15 定理 4.1 を( n + 1 ) c n + 1 = ( 2 n + 1 ) c n (n+1)c_{n+1}=(2n+1)c_n ( n + 1 ) c n + 1 = ( 2 n + 1 ) c n で割り、P k = c k π k P_k=c_k\pi_k P k = c k π k を代入すると
π n + 1 = x π n − n c n − 1 ( 2 n + 1 ) c n π n − 1 \pi_{n+1}=x\pi_n-\frac{n\,c_{n-1}}{(2n+1)c_n}\,\pi_{n-1} π n + 1 = x π n − ( 2 n + 1 ) c n n c n − 1 π n − 1 である。n c n = ( 2 n − 1 ) c n − 1 nc_n=(2n-1)c_{n-1} n c n = ( 2 n − 1 ) c n − 1 (n = 1 n=1 n = 1 ではc 0 = c 1 = 1 c_0=c_1=1 c 0 = c 1 = 1 による)であるからπ n − 1 \pi_{n-1} π n − 1 の係数はn 2 / ( ( 2 n + 1 ) ( 2 n − 1 ) ) n^2/((2n+1)(2n-1)) n 2 / (( 2 n + 1 ) ( 2 n − 1 )) であり、定理 2.3 (2) によりα n = 0 \alpha_n=0 α n = 0 、β n = n 2 / ( 4 n 2 − 1 ) \beta_n=n^2/(4n^2-1) β n = n 2 / ( 4 n 2 − 1 ) である。
最後の式は、命題 2.2 (2) の展開式にπ k = P k / c k \pi_k=P_k/c_k π k = P k / c k を代入し、⟨ f , π k ⟩ w P / ∥ π k ∥ w P 2 ⋅ π k = ⟨ f , P k ⟩ w P / ∥ P k ∥ w P 2 ⋅ P k \langle f,\pi_k\rangle_{w_P}/\lVert\pi_k\rVert_{w_P}^2\cdot\pi_k=\langle f,P_k\rangle_{w_P}/\lVert P_k\rVert_{w_P}^2\cdot P_k ⟨ f , π k ⟩ w P / ∥ π k ∥ w P 2 ⋅ π k = ⟨ f , P k ⟩ w P / ∥ P k ∥ w P 2 ⋅ P k と§E10.15 命題 4.2 を用いて得る。▨
補題 4.2. w T ( x ) : = ( 1 − x 2 ) − 1 / 2 w_T(x):=(1-x^2)^{-1/2} w T ( x ) := ( 1 − x 2 ) − 1/2 は( − 1 , 1 ) (-1,1) ( − 1 , 1 ) 上の重み関数であり、∫ − 1 1 w T ( x ) d x = π \int_{-1}^1w_T(x)\,dx=\pi ∫ − 1 1 w T ( x ) d x = π である。任意のg ∈ C ( [ − 1 , 1 ] ) g\in C([-1,1]) g ∈ C ([ − 1 , 1 ]) に対して
∫ − 1 1 g ( x ) w T ( x ) d x = ∫ 0 π g ( cos θ ) d θ \int_{-1}^1g(x)w_T(x)\,dx=\int_0^\pi g(\cos\theta)\,d\theta ∫ − 1 1 g ( x ) w T ( x ) d x = ∫ 0 π g ( cos θ ) d θ である。
証明. w T w_T w T は( − 1 , 1 ) (-1,1) ( − 1 , 1 ) で正であり、連続であるから、§D1.17 定理 2.5 により− 1 < c < d < 1 -1<c<d<1 − 1 < c < d < 1 を満たす各[ c , d ] [c,d] [ c , d ] で Riemann 可積分である。0 < δ < π / 2 0<\delta<\pi/2 0 < δ < π /2 とする。cos \cos cos は[ δ , π − δ ] [\delta,\pi-\delta] [ δ , π − δ ] でC 1 C^1 C 1 級かつ狭義単調減少であり、この区間を[ − cos δ , cos δ ] [-\cos\delta,\cos\delta] [ − cos δ , cos δ ] の上へ写し、θ ∈ ( 0 , π ) \theta\in(0,\pi) θ ∈ ( 0 , π ) に対してsin θ = ( 1 − cos 2 θ ) 1 / 2 > 0 \sin\theta=(1-\cos^2\theta)^{1/2}>0 sin θ = ( 1 − cos 2 θ ) 1/2 > 0 である。置換積分により、任意のg ∈ C ( [ − 1 , 1 ] ) g\in C([-1,1]) g ∈ C ([ − 1 , 1 ]) に対して
∫ − cos δ cos δ g ( x ) w T ( x ) d x = ∫ δ π − δ g ( cos θ ) sin θ ( 1 − cos 2 θ ) 1 / 2 d θ = ∫ δ π − δ g ( cos θ ) d θ \int_{-\cos\delta}^{\cos\delta}g(x)w_T(x)\,dx=\int_\delta^{\pi-\delta}g(\cos\theta)\frac{\sin\theta}{(1-\cos^2\theta)^{1/2}}\,d\theta=\int_\delta^{\pi-\delta}g(\cos\theta)\,d\theta ∫ − c o s δ c o s δ g ( x ) w T ( x ) d x = ∫ δ π − δ g ( cos θ ) ( 1 − cos 2 θ ) 1/2 sin θ d θ = ∫ δ π − δ g ( cos θ ) d θ である。g = 1 g=1 g = 1 とすると、− 1 < c < d < 1 -1<c<d<1 − 1 < c < d < 1 に対してe : = max ( ∣ c ∣ , ∣ d ∣ ) e:=\max(|c|,|d|) e := max ( ∣ c ∣ , ∣ d ∣ ) 、δ : = arccos e \delta:=\arccos e δ := arccos e と置けば∫ c d w T d x ≤ ∫ − e e w T d x = π − 2 δ < π \int_c^dw_T\,dx\le\int_{-e}^ew_T\,dx=\pi-2\delta<\pi ∫ c d w T d x ≤ ∫ − e e w T d x = π − 2 δ < π である。w T > 0 w_T>0 w T > 0 であるから∫ c d w T d x \int_c^dw_T\,dx ∫ c d w T d x はc c c について減少しd d d について増加し、上の評価によりπ \pi π 以下であるから、w T w_T w T の広義積分は有限であり、w T w_T w T は重み関数である。補題 1.3 (1) によりg w T gw_T g w T の広義積分は存在するから、その値はc = − cos δ c=-\cos\delta c = − cos δ 、d = cos δ d=\cos\delta d = cos δ 、δ → 0 + \delta\to0+ δ → 0 + とした極限に等しく、右辺のδ → 0 + \delta\to0+ δ → 0 + の極限∫ 0 π g ( cos θ ) d θ \int_0^\pi g(\cos\theta)\,d\theta ∫ 0 π g ( cos θ ) d θ に一致する。g = 1 g=1 g = 1 として∫ − 1 1 w T d x = π \int_{-1}^1w_T\,dx=\pi ∫ − 1 1 w T d x = π を得る。▨
命題 4.3. w T ( x ) : = ( 1 − x 2 ) − 1 / 2 w_T(x):=(1-x^2)^{-1/2} w T ( x ) := ( 1 − x 2 ) − 1/2 を( − 1 , 1 ) (-1,1) ( − 1 , 1 ) 上の重み関数とし、w T w_T w T に関するモニック直交多項式列を( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 、三項漸化式の係数をα k , β k \alpha_k,\beta_k α k , β k とする。k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 に対してT k T_k T k をk k k 次の Chebyshev 多項式とする。
j , k ∈ N ≥ 0 j,k\in\N j , k ∈ N ≥ 0 に対して、j ≠ k j\ne k j = k ならば⟨ T j , T k ⟩ w T = 0 \langle T_j,T_k\rangle_{w_T}=0 ⟨ T j , T k ⟩ w T = 0 であり、∥ T 0 ∥ w T 2 = π \lVert T_0\rVert_{w_T}^2=\pi ∥ T 0 ∥ w T 2 = π 、k ≥ 1 k\ge1 k ≥ 1 ならば∥ T k ∥ w T 2 = π / 2 \lVert T_k\rVert_{w_T}^2=\pi/2 ∥ T k ∥ w T 2 = π /2 である。
π 0 = 1 \pi_0=1 π 0 = 1 であり、k ≥ 1 k\ge1 k ≥ 1 ならばπ k = 2 1 − k T k \pi_k=2^{1-k}T_k π k = 2 1 − k T k である。
任意のk ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 に対してα k = 0 \alpha_k=0 α k = 0 であり、β 1 = 1 / 2 \beta_1=1/2 β 1 = 1/2 、k ≥ 2 k\ge2 k ≥ 2 ならばβ k = 1 / 4 \beta_k=1/4 β k = 1/4 である。
任意のf ∈ C ( [ − 1 , 1 ] ) f\in C([-1,1]) f ∈ C ([ − 1 , 1 ]) とn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して、a k : = ( 2 / π ) ∫ 0 π f ( cos θ ) cos k θ d θ a_k:=(2/\pi)\int_0^\pi f(\cos\theta)\cos k\theta\,d\theta a k := ( 2/ π ) ∫ 0 π f ( cos θ ) cos k θ d θ (k ∈ N ≥ 0 k\in\N k ∈ N ≥ 0 )と置くと、Π n w T f = a 0 / 2 + ∑ k = 1 n a k T k \Pi_n^{w_T}f=a_0/2+\sum_{k=1}^na_kT_k Π n w T f = a 0 /2 + ∑ k = 1 n a k T k である。
証明. (1) を示す。補題 4.2 と§E20.11 補題 3.2 (2) により
⟨ T j , T k ⟩ w T = ∫ 0 π cos j θ cos k θ d θ = 1 2 ∫ 0 π ( cos ( j − k ) θ + cos ( j + k ) θ ) d θ \langle T_j,T_k\rangle_{w_T}=\int_0^\pi\cos j\theta\cos k\theta\,d\theta=\frac12\int_0^\pi\bigl(\cos(j-k)\theta+\cos(j+k)\theta\bigr)\,d\theta ⟨ T j , T k ⟩ w T = ∫ 0 π cos j θ cos k θ d θ = 2 1 ∫ 0 π ( cos ( j − k ) θ + cos ( j + k ) θ ) d θ である。整数l l l について∫ 0 π cos l θ d θ \int_0^\pi\cos l\theta\,d\theta ∫ 0 π cos l θ d θ はl = 0 l=0 l = 0 のときπ \pi π 、l ≠ 0 l\ne0 l = 0 のときsin ( l π ) / l = 0 \sin(l\pi)/l=0 sin ( l π ) / l = 0 であるから、主張の値を得る。
(2) を示す。k ≥ 1 k\ge1 k ≥ 1 とする。§E20.11 補題 3.2 (1) によりT 0 , … , T k − 1 ∈ P k − 1 T_0,\dots,T_{k-1}\in\mathcal P_{k-1} T 0 , … , T k − 1 ∈ P k − 1 であり、(1) によりこれらは互いに直交する0 0 0 でないk k k 個の元であるから、k k k 次元のP k − 1 \mathcal P_{k-1} P k − 1 の基底である。したがってT k T_k T k はP k − 1 \mathcal P_{k-1} P k − 1 のすべての元に直交する。§E20.12 命題 5.2 (1) によりT k ∈ P k T_k\in\mathcal P_k T k ∈ P k のx k x^k x k の係数は2 k − 1 2^{k-1} 2 k − 1 であるから、命題 2.2 (3) によりT k = 2 k − 1 π k T_k=2^{k-1}\pi_k T k = 2 k − 1 π k である。
(3) を示す。π 1 = T 1 = x \pi_1=T_1=x π 1 = T 1 = x であるから、定理 2.3 (1) のk = 0 k=0 k = 0 の等式π 1 = ( x − α 0 ) π 0 \pi_1=(x-\alpha_0)\pi_0 π 1 = ( x − α 0 ) π 0 と比べてα 0 = 0 \alpha_0=0 α 0 = 0 である。π 2 = T 2 / 2 = ( 2 x 2 − 1 ) / 2 = x π 1 − π 0 / 2 \pi_2=T_2/2=(2x^2-1)/2=x\pi_1-\pi_0/2 π 2 = T 2 /2 = ( 2 x 2 − 1 ) /2 = x π 1 − π 0 /2 である。k ≥ 2 k\ge2 k ≥ 2 ならば、§E20.11 定義 3.1 の漸化式T k + 1 = 2 x T k − T k − 1 T_{k+1}=2xT_k-T_{k-1} T k + 1 = 2 x T k − T k − 1 に2 − k 2^{-k} 2 − k を掛け、π k − 1 = 2 2 − k T k − 1 \pi_{k-1}=2^{2-k}T_{k-1} π k − 1 = 2 2 − k T k − 1 を用いてπ k + 1 = x π k − π k − 1 / 4 \pi_{k+1}=x\pi_k-\pi_{k-1}/4 π k + 1 = x π k − π k − 1 /4 を得る。定理 2.3 (2) により主張の値を得る。
(4) を示す。命題 2.2 (2) の展開式の各項は、(2) により⟨ f , T k ⟩ w T / ∥ T k ∥ w T 2 ⋅ T k \langle f,T_k\rangle_{w_T}/\lVert T_k\rVert_{w_T}^2\cdot T_k ⟨ f , T k ⟩ w T / ∥ T k ∥ w T 2 ⋅ T k に等しい。補題 4.2 と§E20.11 補題 3.2 (2) により⟨ f , T k ⟩ w T = ∫ 0 π f ( cos θ ) cos k θ d θ \langle f,T_k\rangle_{w_T}=\int_0^\pi f(\cos\theta)\cos k\theta\,d\theta ⟨ f , T k ⟩ w T = ∫ 0 π f ( cos θ ) cos k θ d θ であり、(1) のノルムで割って係数a 0 / 2 a_0/2 a 0 /2 とa k a_k a k (k ≥ 1 k\ge1 k ≥ 1 )を得る。▨
例 4.4. f ( x ) : = ∣ x ∣ f(x):=|x| f ( x ) := ∣ x ∣ をC ( [ − 1 , 1 ] ) C([-1,1]) C ([ − 1 , 1 ]) の元とし、n = 2 n=2 n = 2 とする。
w P : = 1 w_P:=1 w P := 1 のとき、∫ − 1 1 ∣ x ∣ P 0 d x = 1 \int_{-1}^1|x|P_0\,dx=1 ∫ − 1 1 ∣ x ∣ P 0 d x = 1 、∫ − 1 1 ∣ x ∣ P 1 d x = 0 \int_{-1}^1|x|P_1\,dx=0 ∫ − 1 1 ∣ x ∣ P 1 d x = 0 、∫ − 1 1 ∣ x ∣ P 2 d x = 2 ∫ 0 1 x ( 3 x 2 − 1 ) / 2 d x = 1 / 4 \int_{-1}^1|x|P_2\,dx=2\int_0^1x(3x^2-1)/2\,dx=1/4 ∫ − 1 1 ∣ x ∣ P 2 d x = 2 ∫ 0 1 x ( 3 x 2 − 1 ) /2 d x = 1/4 である。命題 4.1 の式により
Π 2 w P f = 1 2 ⋅ 1 + 5 2 ⋅ 1 4 ⋅ 3 x 2 − 1 2 = 3 16 + 15 16 x 2 \Pi_2^{w_P}f=\frac12\cdot1+\frac52\cdot\frac14\cdot\frac{3x^2-1}2=\frac3{16}+\frac{15}{16}x^2 Π 2 w P f = 2 1 ⋅ 1 + 2 5 ⋅ 4 1 ⋅ 2 3 x 2 − 1 = 16 3 + 16 15 x 2
であり、定理 1.6 (5) により∥ f − Π 2 w P f ∥ w P 2 = 2 / 3 − ( 1 / 2 ) 2 ⋅ 2 − ( 5 / 8 ) 2 ⋅ ( 2 / 5 ) = 1 / 96 \lVert f-\Pi_2^{w_P}f\rVert_{w_P}^2=2/3-(1/2)^2\cdot2-(5/8)^2\cdot(2/5)=1/96 ∥ f − Π 2 w P f ∥ w P 2 = 2/3 − ( 1/2 ) 2 ⋅ 2 − ( 5/8 ) 2 ⋅ ( 2/5 ) = 1/96 である。
w = w T w=w_T w = w T のとき、命題 4.3 (4) の係数は、θ ↦ π − θ \theta\mapsto\pi-\theta θ ↦ π − θ の置換によりa 0 = ( 2 / π ) ⋅ 2 ∫ 0 π / 2 cos θ d θ = 4 / π a_0=(2/\pi)\cdot2\int_0^{\pi/2}\cos\theta\,d\theta=4/\pi a 0 = ( 2/ π ) ⋅ 2 ∫ 0 π /2 cos θ d θ = 4/ π 、a 1 = 0 a_1=0 a 1 = 0 、a 2 = ( 2 / π ) ∫ 0 π / 2 ( cos 3 θ + cos θ ) d θ = 4 / ( 3 π ) a_2=(2/\pi)\int_0^{\pi/2}(\cos3\theta+\cos\theta)\,d\theta=4/(3\pi) a 2 = ( 2/ π ) ∫ 0 π /2 ( cos 3 θ + cos θ ) d θ = 4/ ( 3 π ) である。したがって
Π 2 w T f = 2 π + 4 3 π ( 2 x 2 − 1 ) = 2 3 π + 8 3 π x 2 \Pi_2^{w_T}f=\frac2\pi+\frac4{3\pi}(2x^2-1)=\frac2{3\pi}+\frac8{3\pi}x^2 Π 2 w T f = π 2 + 3 π 4 ( 2 x 2 − 1 ) = 3 π 2 + 3 π 8 x 2
であり、∥ f ∥ w T 2 = ∫ 0 π cos 2 θ d θ = π / 2 \lVert f\rVert_{w_T}^2=\int_0^\pi\cos^2\theta\,d\theta=\pi/2 ∥ f ∥ w T 2 = ∫ 0 π cos 2 θ d θ = π /2 と定理 1.6 (5) により∥ f − Π 2 w T f ∥ w T 2 = π / 2 − ( 2 2 / π + ( 2 / 3 ) 2 / ( π / 2 ) ) = π / 2 − 44 / ( 9 π ) ≈ 0.01461 \lVert f-\Pi_2^{w_T}f\rVert_{w_T}^2=\pi/2-(2^2/\pi+(2/3)^2/(\pi/2))=\pi/2-44/(9\pi)\approx0.01461 ∥ f − Π 2 w T f ∥ w T 2 = π /2 − ( 2 2 / π + ( 2/3 ) 2 / ( π /2 )) = π /2 − 44/ ( 9 π ) ≈ 0.01461 である。
Π 2 w P f \Pi_2^{w_P}f Π 2 w P f の係数( 3 / 16 , 15 / 16 ) = ( 0.1875 , 0.9375 ) (3/16,15/16)=(0.1875,0.9375) ( 3/16 , 15/16 ) = ( 0.1875 , 0.9375 ) とΠ 2 w T f \Pi_2^{w_T}f Π 2 w T f の係数( 2 / ( 3 π ) , 8 / ( 3 π ) ) ≈ ( 0.2122 , 0.8488 ) (2/(3\pi),8/(3\pi))\approx(0.2122,0.8488) ( 2/ ( 3 π ) , 8/ ( 3 π )) ≈ ( 0.2122 , 0.8488 ) は異なる。誤差f − p f-p f − p は、p = Π 2 w P f p=\Pi_2^{w_P}f p = Π 2 w P f ではx = 0 x=0 x = 0 で− 3 / 16 -3/16 − 3/16 、x = ± 1 x=\pm1 x = ± 1 で− 1 / 8 -1/8 − 1/8 であり、p = Π 2 w T f p=\Pi_2^{w_T}f p = Π 2 w T f ではx = 0 x=0 x = 0 で− 2 / ( 3 π ) ≈ − 0.2122 -2/(3\pi)\approx-0.2122 − 2/ ( 3 π ) ≈ − 0.2122 、x = ± 1 x=\pm1 x = ± 1 で1 − 10 / ( 3 π ) ≈ − 0.0610 1-10/(3\pi)\approx-0.0610 1 − 10/ ( 3 π ) ≈ − 0.0610 である。x → ± 1 x\to\pm1 x → ± 1 で非有界な重みw T w_T w T に関する最良近似は、x = ± 1 x=\pm1 x = ± 1 での誤差の絶対値がw P w_P w P の場合より小さく、x = 0 x=0 x = 0 での誤差の絶対値が大きい。
5 二乗平均近似の収束
定理 5.1. a < b a<b a < b を実数とし、w w w を( a , b ) (a,b) ( a , b ) 上の重み関数、μ 0 : = ∫ a b w ( x ) d x \mu_0:=\int_a^bw(x)\,dx μ 0 := ∫ a b w ( x ) d x 、f ∈ C ( [ a , b ] ) f\in C([a,b]) f ∈ C ([ a , b ]) とする。n ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対してE n ( f ) : = inf q ∈ P n ∥ f − q ∥ ∞ E_n(f):=\inf_{q\in\mathcal P_n}\lVert f-q\rVert_\infty E n ( f ) := inf q ∈ P n ∥ f − q ∥ ∞ をf f f のP n \mathcal P_n P n による最良一様近似誤差とする。
任意のn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して∥ f − Π n w f ∥ w ≤ μ 0 1 / 2 E n ( f ) \lVert f-\Pi_n^wf\rVert_w\le\mu_0^{1/2}E_n(f) ∥ f − Π n w f ∥ w ≤ μ 0 1/2 E n ( f ) である。
任意のn ∈ N ≥ 0 n\in\N n ∈ N ≥ 0 に対して∥ f − Π n + 1 w f ∥ w ≤ ∥ f − Π n w f ∥ w \lVert f-\Pi_{n+1}^wf\rVert_w\le\lVert f-\Pi_n^wf\rVert_w ∥ f − Π n + 1 w f ∥ w ≤ ∥ f − Π n w f ∥ w である。
lim n → ∞ ∥ f − Π n w f ∥ w = 0 \lim_{n\to\infty}\lVert f-\Pi_n^wf\rVert_w=0 lim n → ∞ ∥ f − Π n w f ∥ w = 0 である。
証明. (1) を示す。q ∈ P n q\in\mathcal P_n q ∈ P n を任意にとる。定理 1.6 (1) と命題 1.5 により∥ f − Π n w f ∥ w ≤ ∥ f − q ∥ w ≤ μ 0 1 / 2 ∥ f − q ∥ ∞ \lVert f-\Pi_n^wf\rVert_w\le\lVert f-q\rVert_w\le\mu_0^{1/2}\lVert f-q\rVert_\infty ∥ f − Π n w f ∥ w ≤ ∥ f − q ∥ w ≤ μ 0 1/2 ∥ f − q ∥ ∞ である。q q q について下限をとって主張を得る。
(2) を示す。Π n w f ∈ P n ⊆ P n + 1 \Pi_n^wf\in\mathcal P_n\subseteq\mathcal P_{n+1} Π n w f ∈ P n ⊆ P n + 1 であるから、定理 1.6 (1) をn + 1 n+1 n + 1 について適用して主張を得る。
(3) を示す。ε > 0 \varepsilon>0 ε > 0 をとる。§E2.32 問題 5.2 により、実係数多項式p p p であって∥ f − p ∥ ∞ < ε \lVert f-p\rVert_\infty<\varepsilon ∥ f − p ∥ ∞ < ε を満たすものが存在する。N N N をp p p の次数(p = 0 p=0 p = 0 のときN : = 0 N:=0 N := 0 )とすると、n ≥ N n\ge N n ≥ N ならばp ∈ P n p\in\mathcal P_n p ∈ P n でありE n ( f ) ≤ ∥ f − p ∥ ∞ < ε E_n(f)\le\lVert f-p\rVert_\infty<\varepsilon E n ( f ) ≤ ∥ f − p ∥ ∞ < ε である。(1) により、n ≥ N n\ge N n ≥ N ならば∥ f − Π n w f ∥ w < μ 0 1 / 2 ε \lVert f-\Pi_n^wf\rVert_w<\mu_0^{1/2}\varepsilon ∥ f − Π n w f ∥ w < μ 0 1/2 ε である。▨
6 離散と連続の最小二乗の比較
例 6.1. f ( x ) : = ∣ x ∣ f(x):=|x| f ( x ) := ∣ x ∣ 、次数2 2 2 、区間[ − 1 , 1 ] [-1,1] [ − 1 , 1 ] 、( − 1 , 1 ) (-1,1) ( − 1 , 1 ) 上の重み関数w P : = 1 w_P:=1 w P := 1 とする。整数m ≥ 3 m\ge3 m ≥ 3 に対して等間隔の節点x i : = − 1 + 2 i / ( m − 1 ) x_i:=-1+2i/(m-1) x i := − 1 + 2 i / ( m − 1 ) (0 ≤ i ≤ m − 1 0\le i\le m-1 0 ≤ i ≤ m − 1 )をとり、A : = ( x i j ) ∈ R m × 3 A:=(x_i^{\,j})\in\R^{m\times3} A := ( x i j ) ∈ R m × 3 (0 ≤ i ≤ m − 1 0\le i\le m-1 0 ≤ i ≤ m − 1 、0 ≤ j ≤ 2 0\le j\le2 0 ≤ j ≤ 2 )、b : = ( ∣ x i ∣ ) i ∈ R m b:=(|x_i|)_i\in\R^m b := ( ∣ x i ∣ ) i ∈ R m と置く。x 0 , x 1 , x 2 x_0,x_1,x_2 x 0 , x 1 , x 2 は相異なるから、§E20.12 定理 1.3 (3) によりA A A の第0 0 0 行から第2 2 2 行までからなる3 × 3 3\times3 3 × 3 行列は正則であり、A A A の列は一次独立である。§E20.6 命題 1.2 (2) によりA x = b Ax=b A x = b の最小二乗解c ( m ) c^{(m)} c ( m ) はただ一つであり、p m : = ∑ j c j ( m ) x j p_m:=\sum_jc^{(m)}_jx^j p m := ∑ j c j ( m ) x j と置く。m = 3 m=3 m = 3 では§E20.12 定理 1.3 (2) の補間多項式x 2 x^2 x 2 が残差0 0 0 を与えるから、p 3 = x 2 p_3=x^2 p 3 = x 2 は節点− 1 , 0 , 1 -1,0,1 − 1 , 0 , 1 におけるf f f の補間多項式である。正規方程式を有理数の演算で解くと次の値を得る。
多項式p p p
p p p の係数( 1 , x , x 2 ) (1,x,x^2) ( 1 , x , x 2 )
∥ f − p ∥ w P 2 \lVert f-p\rVert_{w_P}^2 ∥ f − p ∥ w P 2
max [ − 1 , 1 ] ∣ f − p ∣ \max_{[-1,1]}\lvert f-p\rvert max [ − 1 , 1 ] ∣ f − p ∣
p 3 p_3 p 3 (補間)
( 0 , 0 , 1 ) (0,0,1) ( 0 , 0 , 1 )
1 / 15 ≈ 0.06667 1/15\approx0.06667 1/15 ≈ 0.06667
1 / 4 1/4 1/4
p 5 p_5 p 5
( 6 / 35 , 0 , 6 / 7 ) (6/35,0,6/7) ( 6/35 , 0 , 6/7 )
8 / 525 ≈ 0.01524 8/525\approx0.01524 8/525 ≈ 0.01524
6 / 35 ≈ 0.1714 6/35\approx0.1714 6/35 ≈ 0.1714
p 21 p_{21} p 21
( 594 / 3059 , 0 , 2750 / 3059 ) (594/3059,0,2750/3059) ( 594/3059 , 0 , 2750/3059 )
≈ 0.01076 \approx0.01076 ≈ 0.01076
594 / 3059 ≈ 0.1942 594/3059\approx0.1942 594/3059 ≈ 0.1942
Π 2 w P f \Pi_2^{w_P}f Π 2 w P f
( 3 / 16 , 0 , 15 / 16 ) (3/16,0,15/16) ( 3/16 , 0 , 15/16 )
1 / 96 ≈ 0.01042 1/96\approx0.01042 1/96 ≈ 0.01042
3 / 16 = 0.1875 3/16=0.1875 3/16 = 0.1875
定理 1.6 (1) により∥ f − p ∥ w P \lVert f-p\rVert_{w_P} ∥ f − p ∥ w P はP 2 \mathcal P_2 P 2 の中でp = Π 2 w P f p=\Pi_2^{w_P}f p = Π 2 w P f のときに最小であり、表の第三列の最小値はΠ 2 w P f \Pi_2^{w_P}f Π 2 w P f の行にある。第四列ではp 5 p_5 p 5 の値がΠ 2 w P f \Pi_2^{w_P}f Π 2 w P f の値より小さく、Π 2 w P f \Pi_2^{w_P}f Π 2 w P f はP 2 \mathcal P_2 P 2 の中で最大値ノルムの誤差を最小にしない。
m → ∞ m\to\infty m → ∞ でc ( m ) c^{(m)} c ( m ) はΠ 2 w P f \Pi_2^{w_P}f Π 2 w P f の係数( 3 / 16 , 0 , 15 / 16 ) (3/16,0,15/16) ( 3/16 , 0 , 15/16 ) に収束する。実際、h : = 2 / ( m − 1 ) h:=2/(m-1) h := 2/ ( m − 1 ) と置くと、正規方程式A T A c ( m ) = A T b A^{\mathsf T}Ac^{(m)}=A^{\mathsf T}b A T A c ( m ) = A T b にh h h を掛けた式の係数h ∑ i x i j + k h\sum_ix_i^{\,j+k} h ∑ i x i j + k と右辺h ∑ i ∣ x i ∣ x i j h\sum_i|x_i|x_i^{\,j} h ∑ i ∣ x i ∣ x i j は、連続関数x j + k x^{j+k} x j + k と∣ x ∣ x j |x|x^j ∣ x ∣ x j の刻みh h h の Riemann 和に項h ⋅ 1 j + k h\cdot1^{j+k} h ⋅ 1 j + k 、h ⋅ 1 h\cdot1 h ⋅ 1 を加えたものであるから、それぞれ⟨ x j , x k ⟩ w P \langle x^j,x^k\rangle_{w_P} ⟨ x j , x k ⟩ w P と⟨ f , x j ⟩ w P \langle f,x^j\rangle_{w_P} ⟨ f , x j ⟩ w P に収束する。極限の行列G : = ( ⟨ x j , x k ⟩ w P ) j , k G:=(\langle x^j,x^k\rangle_{w_P})_{j,k} G := (⟨ x j , x k ⟩ w P ) j , k は、G c = 0 Gc=0 G c = 0 ならば∥ ∑ j c j x j ∥ w P 2 = c T G c = 0 \lVert\sum_jc_jx^j\rVert_{w_P}^2=c^{\mathsf T}Gc=0 ∥ ∑ j c j x j ∥ w P 2 = c T G c = 0 となってc = 0 c=0 c = 0 であるから正則である。逆行列は正則行列の上で連続であるから、c ( m ) c^{(m)} c ( m ) はG − 1 ( ⟨ f , x j ⟩ w P ) j G^{-1}(\langle f,x^j\rangle_{w_P})_j G − 1 (⟨ f , x j ⟩ w P ) j に収束し、定理 1.6 (4) によりこの極限はΠ 2 w P f \Pi_2^{w_P}f Π 2 w P f の係数である。
7 演習
問題 7.1. w U ( x ) : = ( 1 − x 2 ) 1 / 2 w_U(x):=(1-x^2)^{1/2} w U ( x ) := ( 1 − x 2 ) 1/2 を( − 1 , 1 ) (-1,1) ( − 1 , 1 ) 上の重み関数とする。実係数多項式の列( U k ) k ∈ N ≥ 0 (U_k)_{k\in\N} ( U k ) k ∈ N ≥ 0 をU 0 : = 1 U_0:=1 U 0 := 1 、U 1 : = 2 x U_1:=2x U 1 := 2 x 、U k + 1 : = 2 x U k − U k − 1 U_{k+1}:=2xU_k-U_{k-1} U k + 1 := 2 x U k − U k − 1 (k ≥ 1 k\ge1 k ≥ 1 )で定める。w U w_U w U に関するモニック直交多項式列( π k ) k ∈ N ≥ 0 (\pi_k)_{k\in\N} ( π k ) k ∈ N ≥ 0 がπ k = 2 − k U k \pi_k=2^{-k}U_k π k = 2 − k U k であることを示し、三項漸化式の係数α k , β k \alpha_k,\beta_k α k , β k を求めよ。
解答. U k U_k U k は次数k k k 、x k x^k x k の係数2 k 2^k 2 k の多項式である。実際、k = 0 , 1 k=0,1 k = 0 , 1 では定義から成り立ち、k ≥ 1 k\ge1 k ≥ 1 についてU k U_k U k が次数k k k 、係数2 k 2^k 2 k でありU k − 1 U_{k-1} U k − 1 の次数がk − 1 k-1 k − 1 であれば、2 x U k − U k − 1 2xU_k-U_{k-1} 2 x U k − U k − 1 は次数k + 1 k+1 k + 1 、x k + 1 x^{k+1} x k + 1 の係数2 k + 1 2^{k+1} 2 k + 1 である。
θ ∈ ( 0 , π ) \theta\in(0,\pi) θ ∈ ( 0 , π ) に対してU k ( cos θ ) sin θ = sin ( k + 1 ) θ U_k(\cos\theta)\sin\theta=\sin(k+1)\theta U k ( cos θ ) sin θ = sin ( k + 1 ) θ である。実際、k = 0 k=0 k = 0 では両辺はsin θ \sin\theta sin θ 、k = 1 k=1 k = 1 では2 cos θ sin θ = sin 2 θ 2\cos\theta\sin\theta=\sin2\theta 2 cos θ sin θ = sin 2 θ であり、加法定理sin ( k + 2 ) θ + sin k θ = 2 cos θ sin ( k + 1 ) θ \sin(k+2)\theta+\sin k\theta=2\cos\theta\sin(k+1)\theta sin ( k + 2 ) θ + sin k θ = 2 cos θ sin ( k + 1 ) θ により、k − 1 k-1 k − 1 とk k k で等式が成り立てばU k + 1 ( cos θ ) sin θ = 2 cos θ sin ( k + 1 ) θ − sin k θ = sin ( k + 2 ) θ U_{k+1}(\cos\theta)\sin\theta=2\cos\theta\sin(k+1)\theta-\sin k\theta=\sin(k+2)\theta U k + 1 ( cos θ ) sin θ = 2 cos θ sin ( k + 1 ) θ − sin k θ = sin ( k + 2 ) θ である。
w U w_U w U は( − 1 , 1 ) (-1,1) ( − 1 , 1 ) で正であり、[ − 1 , 1 ] [-1,1] [ − 1 , 1 ] へ連続に延びるから、§D1.17 定理 2.5 により− 1 < c < d < 1 -1<c<d<1 − 1 < c < d < 1 を満たす各[ c , d ] [c,d] [ c , d ] で Riemann 可積分であって、広義積分は延長の[ − 1 , 1 ] [-1,1] [ − 1 , 1 ] 上の積分に等しく有限である。したがってw U w_U w U は重み関数であり、⟨ f , g ⟩ w U \langle f,g\rangle_{w_U} ⟨ f , g ⟩ w U は通常の積分∫ − 1 1 f g w U d x \int_{-1}^1fg\,w_U\,dx ∫ − 1 1 f g w U d x に等しい。x = cos θ x=\cos\theta x = cos θ (θ ∈ [ 0 , π ] \theta\in[0,\pi] θ ∈ [ 0 , π ] )と置換すると、d x = − sin θ d θ dx=-\sin\theta\,d\theta d x = − sin θ d θ 、w U ( cos θ ) = sin θ w_U(\cos\theta)=\sin\theta w U ( cos θ ) = sin θ であるから、
⟨ U j , U k ⟩ w U = ∫ 0 π U j ( cos θ ) U k ( cos θ ) sin 2 θ d θ = ∫ 0 π sin ( j + 1 ) θ sin ( k + 1 ) θ d θ = 1 2 ∫ 0 π ( cos ( j − k ) θ − cos ( j + k + 2 ) θ ) d θ \langle U_j,U_k\rangle_{w_U}=\int_0^\pi U_j(\cos\theta)U_k(\cos\theta)\sin^2\theta\,d\theta=\int_0^\pi\sin(j+1)\theta\,\sin(k+1)\theta\,d\theta=\frac12\int_0^\pi\bigl(\cos(j-k)\theta-\cos(j+k+2)\theta\bigr)\,d\theta ⟨ U j , U k ⟩ w U = ∫ 0 π U j ( cos θ ) U k ( cos θ ) sin 2 θ d θ = ∫ 0 π sin ( j + 1 ) θ sin ( k + 1 ) θ d θ = 2 1 ∫ 0 π ( cos ( j − k ) θ − cos ( j + k + 2 ) θ ) d θ である。整数l ≠ 0 l\ne0 l = 0 について∫ 0 π cos l θ d θ = 0 \int_0^\pi\cos l\theta\,d\theta=0 ∫ 0 π cos l θ d θ = 0 であるから、j ≠ k j\ne k j = k ならば⟨ U j , U k ⟩ w U = 0 \langle U_j,U_k\rangle_{w_U}=0 ⟨ U j , U k ⟩ w U = 0 、∥ U k ∥ w U 2 = π / 2 \lVert U_k\rVert_{w_U}^2=\pi/2 ∥ U k ∥ w U 2 = π /2 である。
k ≥ 1 k\ge1 k ≥ 1 とする。U 0 , … , U k − 1 U_0,\dots,U_{k-1} U 0 , … , U k − 1 はP k − 1 \mathcal P_{k-1} P k − 1 の互いに直交する0 0 0 でないk k k 個の元であるから、P k − 1 \mathcal P_{k-1} P k − 1 の基底であり、U k U_k U k はP k − 1 \mathcal P_{k-1} P k − 1 のすべての元に直交する。命題 2.2 (3) によりU k = 2 k π k U_k=2^k\pi_k U k = 2 k π k である。k = 0 k=0 k = 0 ではU 0 = 1 = π 0 U_0=1=\pi_0 U 0 = 1 = π 0 である。
π 1 = U 1 / 2 = x \pi_1=U_1/2=x π 1 = U 1 /2 = x であるから、定理 2.3 (1) のk = 0 k=0 k = 0 の等式と比べてα 0 = 0 \alpha_0=0 α 0 = 0 である。k ≥ 1 k\ge1 k ≥ 1 ならば、U k + 1 = 2 x U k − U k − 1 U_{k+1}=2xU_k-U_{k-1} U k + 1 = 2 x U k − U k − 1 に2 − k − 1 2^{-k-1} 2 − k − 1 を掛け、π k − 1 = 2 1 − k U k − 1 \pi_{k-1}=2^{1-k}U_{k-1} π k − 1 = 2 1 − k U k − 1 を用いてπ k + 1 = x π k − π k − 1 / 4 \pi_{k+1}=x\pi_k-\pi_{k-1}/4 π k + 1 = x π k − π k − 1 /4 を得る。定理 2.3 (2) により、k ∈ N ≥ 1 k\in\NN k ∈ N ≥ 1 に対してα k = 0 \alpha_k=0 α k = 0 、β k = 1 / 4 \beta_k=1/4 β k = 1/4 である。▨