最佳平方逼近
函数逼近是使用一种简单易算的函数来近似表示一个复杂函数。
该问题可转化为求解线性方程组
G n C = F n G_{n}C=F_{n} GnC=Fn
其中,系数 C = ( c 0 , c 1 , ⋯ , c n ) T , F n = ( ( f , φ 0 ) , ( f , φ 1 ) , ⋯ , ( f , φ n ) ) T C=(c_{0},c_{1},\cdots,c_{n})^{\mathrm{T}},F_{n}=((f,\varphi_{0}),(f,\varphi_{1}),\cdots,(f,\varphi_{n}))^{\mathrm{T}} C=(c0,c1,⋯,cn)T,Fn=((f,φ0),(f,φ1),⋯,(f,φn))T
G n G_n Gn是格拉姆矩阵。称该线性方程组为法方程组或正规方程组。
最佳平方逼近的解函数为 φ ∗ = ∑ i = 0 n c i ∗ φ i \varphi^*=\sum_{i=0}^nc_i^*\varphi_i φ∗=∑i=0nci∗φi。
最佳平方逼近函数,继承内积,即 ( φ ∗ , φ ∗ ) = ( φ ∗ , f ) (\varphi^*,\varphi^*)=(\varphi^*,f) (φ∗,φ∗)=(φ∗,f)。
取逼近区间[a,b]为[0,1]时,其平方误差为:
∥
φ
∗
−
f
∥
2
2
=
(
f
,
f
)
−
F
n
T
C
∗
=
∫
0
1
f
2
(
x
)
d
x
−
F
n
T
C
∗
.
\parallel\varphi^*-f\parallel_2^2=(f,f)-F_n^\mathrm{T}C^*=\int_0^1f^2(x) \mathrm{d}x-F_n^\mathrm{T} C^* .
∥φ∗−f∥22=(f,f)−FnTC∗=∫01f2(x)dx−FnTC∗.
正交系
内积空间
V
V
V上的两个元素
f
f
f和
g
g
g,如果有内积
(
f
,
g
)
=
0
(f,g)=0
(f,g)=0,则称
f
f
f和
g
g
g关于内积
(
⋅
,
⋅
)
(\cdot,\cdot)
(⋅,⋅)正交。若内积空间上的元素系
{
f
i
}
\{f_{i}\}
{fi}满足两两正交
{
(
f
i
,
f
j
)
=
0
(
i
≠
j
)
,
(
f
i
,
f
i
)
=
γ
i
>
0
,
\begin{cases}(f_i,f_j)=0\quad(i\neq j) ,\\(f_i,f_i)=\gamma_i>0 ,\end{cases}
{(fi,fj)=0(i=j),(fi,fi)=γi>0,
则称
{
f
i
}
\{f_{i}\}
{fi}为正交系,若有
(
f
i
,
f
i
)
=
1
(
i
=
0
,
1
,
2
,
3...
)
(f_i,f_i)=1(i=0,1,2,3...)
(fi,fi)=1(i=0,1,2,3...),则称
{
f
i
}
\{f_{i}\}
{fi}为标准正交系。
给定一组正交基,法方程组系数矩阵
G
n
G_n
Gn为对角矩阵,其解向量为:
C
∗
=
(
(
φ
0
,
f
)
(
φ
0
,
φ
0
)
,
(
φ
1
,
f
)
(
φ
1
,
φ
1
)
,
⋯
,
(
φ
n
,
f
)
(
φ
n
,
φ
n
)
)
T
.
C^* = \left(\frac{(\varphi_0,f)}{(\varphi_0,\varphi_0)},\frac{(\varphi_1,f)}{(\varphi_1,\varphi_1)},\cdots,\frac{(\varphi_n,f)}{(\varphi_n,\varphi_n)}\right)^\mathrm{T}.
C∗=((φ0,φ0)(φ0,f),(φ1,φ1)(φ1,f),⋯,(φn,φn)(φn,f))T.
函数
f
(
x
)
f(x)
f(x)的最佳平方逼近函数为
φ
∗
=
(
φ
0
,
f
)
(
φ
0
,
φ
0
)
φ
0
+
(
φ
1
,
f
)
(
φ
1
,
φ
1
)
φ
1
+
⋯
+
(
φ
n
,
f
)
(
φ
n
,
φ
n
)
φ
n
.
\varphi^*=\frac{(\varphi_0,f)}{(\varphi_0,\varphi_0)}\varphi_0+\frac{(\varphi_1,f)}{(\varphi_1,\varphi_1)}\varphi_1+\cdots+\frac{(\varphi_n,f)}{(\varphi_n,\varphi_n)}\varphi_n.
φ∗=(φ0,φ0)(φ0,f)φ0+(φ1,φ1)(φ1,f)φ1+⋯+(φn,φn)(φn,f)φn.
平方误差为
∥
f
−
φ
∗
∥
2
2
=
(
f
,
f
)
−
∑
i
=
0
n
(
f
,
φ
i
)
2
(
φ
i
,
φ
i
)
.
\parallel f-\varphi^*\parallel_2^2=(f,f)-\sum_{i=0}^n\frac{(f,\varphi_i)^2}{(\varphi_i,\varphi_i)}.
∥f−φ∗∥22=(f,f)−i=0∑n(φi,φi)(f,φi)2.