2基于MATLAB的科学计算—函数逼近2.docVIP

数值分析—最佳逼近 ━基于ＭＡＴＬＡＢ的实现与分析 §２　函数最均方逼近 函数最佳均方逼近准则 　 　　　（９） 与下面的准则等价 　 　　　（１０） 为了讨论问题时方便，在下面的讨论中我们采用准则（１０）。 一般人们习惯于把一个次多项式写成的线性组合，即 　　　　　　　　（１１）次多项式做最佳均方逼近时，积分 　　 　　　（１２） 的结果依赖于次多项式系数，，即 　　 　　　（１３） 所以最佳平方逼近多项式必须满足如下条件： ，　　　　（１４） 即 ，　　　　　（１５） 由式（１５）可知，最佳平方逼近多项式的系数，，是阶线性方程组（１５）的解，但是线性方程组（１５）通常是病态方程组，当，时，对于，取，.由于 ， 则法方程的系数矩阵 . （6.12） 记 ，则 （6.13） 的解即为所求. 例5 在中求定义在[0，1]上的函数的最佳平方逼近多项式. 解 由于，而由上面结论，有 ，，，， 得法方程 ， 用列主元消去法解得.因此所求解为 . 平方误差为： . 【注】 用作基求最佳平方逼近多项式，得到法方程的系数矩阵为Hilbert矩阵，当较大时，是高度病态的.对于例5的三阶Hilbert矩阵，它的条件数为. 　　为避免解病态方程组，需要引入 函数内积的概念：对于，，定义内积如下： 　　　　　　　　（１６） 由内积诱导的范数（距离） 　　　　　　　（１７） 正交的概念： 　　　　　　（１８） 基于函数的内积与正交的概念，如果我们能找到一组正交多项式，： 　　　　　　（１９） 并且将次多项式表示成这组正交多项式的线性组合： 　　　　　　　　　　　（２０） 那么采用内积的记号，这时正规方程组（１５） 　　　　　（２１） 就具有下面的特殊形式： 　　　　　（２２） 正交多项式的构造： ，　　　　（２３） 其中 ，　 　（２４） 利用区间平移和伸缩变换：，． 　　　　　　　　　　　（２５） 　　　　　　　　　　　（２６） 可以原来定义在区间上的正交多项式族转化成区间上的正交多项式族。 Legendre多项式及其应用 在[-1，1]上，取权函数=1，由上述递推关系得勒让德多项式 ，，，，. 这是勒让德于1785年引进的.1814年罗德利克(Rodrigul)给出了简单的表达式 . 由于是次多项式，求阶导数后得 ， 于是得首项的系数.显然最高项系数为1的勒让德多项式为 . 勒让德多项式有下述几个重要性质. 【性质1】证明省略. 【性质2】. 【性质3】， 由，三项递推可得 . 证明省略.勒让德多项式的图形见图6-4. Plot_Legendre clear all clf x=sym(x); K=10; for n=1:K; if n==1 LN(n)=diff(((2^(n-1))*prod([1:n-1]))\(x^2-1)^(n-1),n-1)+1; else LN(n)=diff(((2^(n-1))*prod([1:n-1]))\(x^2-1)^(n-1),n-1); end end L(1)=0; P=; for i=1:K; if i==2 LN(i)=x^1; L(i+1)=L(i)+length(char(LN(i))); P=strcat(P,char(LN(i))); else L(i+1)=L(i)+length(char(LN(i))); P=strcat(P,char(LN(i))); end end for i=1:K; fplot(P(L(i)+1:L(i+1)),[-1,1]) hold on title(LEGENDRE POLYNORMIAL) plot([-1:0.05:1],0*[-1:0.05:1],r) pause hold off end 【性质4】在区间[-1，1]内有个不同的实零点. 图6-4 勒让德多项式图形 图6-5 切比雪夫多项式图形 例２　求函数在区间[-1，1]上的二次最佳平方逼近多项式. 解 利用勒让德多项式，，求二次最佳平方逼近多项式. ， ， ， 因此逼近多项式 . 逼近多项式的平方误差为 . 若所给的区间是，可以通过变换将其转化为[-1，1]. 例３　利用Legendre多项式求函数 在区间上的满足给定精度要求的最佳平方逼近多项式。 clear all syms x t; F=t*exp(t)*sin(t); N=input(Input the degree of