数值分析 8.Approximation Theory

1.Discrete Least Squares Approximation

给定点$(x_1,y_1),...(x_m,y_m)$，求$P_n(x)=a_0+a_{1}x+...a_n{x}^n$,使得$E_2=\sum _{k=1}^m[P_n(x_k)-y_k{]}^2$最小。

需要E对所有a的导数都为0。

化简的结果为：$b_i=\sum _{i=1}^m{x}_i^k,c_k=\sum _{i=1}^m{y}_i{x}_i^k$

$\left[\begin{array}{ccc}{b}_{0+0}& ...& b_{0+n}\\ .& .& .\\ .& .& .\\ .& .& .\\ b_{n+0}& ...& b_{n+n}\end{array}\right]\left[\begin{array}{c}{a}_0\\ .\\ .\\ .\\ a_n\end{array}\right]=\left[\begin{array}{c}{c}_0\\ .\\ .\\ .\\ c_n\end{array}\right]$

可以通过取对数，取倒数再换元的方法将非多项式的函数转化为多项式函数。

如果是想用别的多项式基底表示也可采用：

$(f,g)=\sum _{i=1}^m{w}_{i}f(x_i)g(x_i)$

$(f,g)=\underset{a}{\overset{b}{\int }}w(x)f(x)g(x)dx$

$P(x)=a_0{\varphi }_0(x)+a_1{\varphi }_1(x)+...a_n{\varphi }_n(x)$

需要逼近的多项式为$f(x)$.

$\left[\begin{array}{ccc}({\varphi }_0,{\varphi }_0)& ...& ({\varphi }_0,{\varphi }_n)\\ .& .& .\\ .& .& .\\ .& .& .\\ ({\varphi }_n,{\varphi }_0)& ...& ({\varphi }_n,{\varphi }_n)\end{array}\right]\left[\begin{array}{c}{a}_0\\ .\\ .\\ .\\ a_n\end{array}\right]=\left[\begin{array}{c}({\varphi }_0,f)\\ .\\ .\\ .\\ ({\varphi }_n,f)\end{array}\right]$

2.Construction of the orthogonal polynomials

使用斯密特正交化法：

$$\begin{gathered} {\varphi }_0(x)=1, \\ {\varphi }_1(x)=x-B_1, \\ {\varphi }_k(x)=(x-B_k){\varphi }_{k-1}-C_k{\varphi }_{k-2} \end{gathered}$$

$$\begin{gathered} B_k=\frac{(x{\varphi }_{k-1},{\varphi }_{k-1})}{({\varphi }_{k-1},{\varphi }_{k-1})} \\ {C}_k==\frac{(x{\varphi }_{k-1},{\varphi }_{k-2})}{({\varphi }_{k-2},{\varphi }_{k-2})} \end{gathered}$$

3.Chebyshev Polynomials

第n个切比雪夫多项式：$T_n(x)=\cos (n\cdot \arccos x)$

递推公式：$T_0(x)=1,T_1(x)=x,T_{n+1}(x)=2x{T}_n(x)-T_{n-1}(x)$

1.f(x)在[a,b]上的n次插值多项式。

需要n个点构造拉格朗日多项式进行插值，n个点的选取如下：

切比雪夫多项式的零点：

$t_i=\cos (\frac{2k-1}{2n}\pi ),k=1,2,...,n$

从[-1,1]线性映射到[a,b]

$x_i=\frac{(b-a)t+(a+b)}{2}$

使用$x_i$进行线性插值。

2.要求切比雪夫多项式插值结果小于m

切比雪夫多项式的误差为$R_n=|\frac{f^{(n+1)}(ϵ)}{(n+1)!}\prod _{i=0}^n(x-x_i)|$

在[a,b]上的误差为$R_n=\frac{M}{(n+1)!(b-a{)}^{n+1}}$

M为导数的上确界，找一个较为合适的值。

然后代出n，按照1的方法做。

3.对泰勒公式降阶

$P_{n-1}(x)=P_n(x)+a_n\frac{T_{n-1}(x)}{2^{n-1}}$

其中$P_n(x)$是n阶泰勒展开，$T_n(x)$是第n个切比雪夫多项式。$a_n$是第n个泰勒展开的系数