数值计算方法教案

1、3 最佳平方逼近31法方程设已知，且选择一函数类，其中且设在上线性无关（例如取或等）。研究最佳平方逼近问题：寻求（31）或写为 这里我们主要研究最佳平方逼近函数存在性，唯一性，计算等问题。设有，即使（31）式成立，来考查应满足什么条件。对于任一，即有，于是 (3.2) (3.2)式说明均方误差是多元函数（为二次函数），由设存在是极值问题（3.1）解，即说明存在使由多元函数取极值的必要条件，则有 （3.3）计算由（3。3）式，则有应满足方程组或或总结上述讨论有结论：（1）如果是最佳平方逼近函数，则系数满足方程组 （3.4）其中系数矩阵G是由基函数作内积构成，方程组称为法方程组。 （b）误差函数与

2、基函数正交，即事实上，由（3.4）式有 即 所以 （2）由线性无关得|G|0，则法方程组有唯一解在S中最佳平方逼近函数。事实上，即有 （3.5）如果能证明，对任何，则有那么，满足考查（记） 因 ） （因为），及（3.5）式有总结上述讨论有结论：定理6 （最佳平方逼近）（1）设；（2）选择函数类其中且于线性无关。则(a)在S中最佳平方逼近函数存在且唯一，即存在使 (b)可由解法方程组求得，于是的最佳平方逼近函数为 （因为）32 用多项式作最佳平方逼近已知。（1） 选取寻求使 显然，计算（2）求解法方程组： 即得：特别，设，则法方程组为： （或）求解，则可得上述矩阵称为Hilbert矩阵，是一个著

3、名病态矩阵（对解而言），且随增大，病态愈严重，求得比较准确的计算解就愈困难。因此，取中基，求是佳平方逼近多项式当较大时用一般计算方法求得的解是不可靠的，当增加时，这种方程组计算解精度由舍入误差影响而迅速恶化。一个回避病态矩阵的的办法是取中正交基来做最佳平方逼近。33用正交多项式作最佳平方逼近设。（1）选取中正交基权函数，寻求使 由设，。（2）求解法方程组 于是，得到在最佳平方逼近多项式 定理7 （用正交多项式作最佳平方逼近）（1） 设；（2） 选取中正交基即，为权函数，则在中最佳平方逼近多项式其中，均方误差= 由此，用正交多项式求得最佳平方逼近多项式，其中计算Gram矩阵时，只需计算对角元素，

4、这大大降低了法方程组系数矩阵所涉及的定积分计算量。例4 求在-1，1上3次最佳平方逼近多项式。解 取中正交基其中为Legendre多项式。于-1，1在中3次最佳逼近多项式为：其中且表3-1所以由表3-1： 下面画出图形和近似函数的图形，直观感受3次最佳平方逼近多项式对于的逼近效果。x=-1:0.01:1;y1=exp(x);y2=0.9963+0.9979*x+0.5367*x.2+0.1761*x.3;plot(x,y1,k-.,x,y2,r-)xlabel(x)ylabel(y)legend(f(x)=ex,p_3*(x)title(fontname new roman比较f(x)=ex与

5、其平方逼近函数p_3*(x)=0.9963+0.9979x+0.5367x2+0.1761x3图形);以下介绍如何应用Chebyshev多项式作最佳平方逼近。设。选取中正交基，其中为Chebyshev多项式。，权函数，寻求使由定理7，在中最佳平方逼近多项式为或均方误差 如果，要求在上最佳平方逼近多项式：作变换 于是， 且可用Legendre多项式求在-1，1的最佳平方逼近，其中最后，利用可得函数在上最佳平方逼近多项式例5 用Chebyshev多项式求在-1，1上3次最佳平方逼近多项式。解 3次最佳逼近多项式为其中，令于是，可用数值积分计算积分（见表3.2） 表3。21234562.1.0.0.0.0.0. 在上3次最佳逼近多项式为：且或 用Chebyshev多项式求得的3次最佳平方逼近的最大误差接近最佳一致逼近的误差（见图3-3），且误差函数的分布很相似（见本章6）。x=-1:0.01:1;y1=exp(x);y2=0.+0.*x+0.*x.2+0.*x.3;plot(x,y1,k-.,x,y2,r-)xlabel(x)ylabel(y)legend(f(x)