Ch3函数逼近与快速傅里叶变换精选.pptVIP

对方程y aebx 两边取对数得lny lna+bx，如果令z lny,A lna，则z A+bx,? 1,x . 为确定A和b，先将数据 xi, yi 可转化为 xi, zi ，见数据表. 根据最小二乘法，取 得 故有法方程 解得A 1.122, b 0.505, a eA 3.071. 于是得最小二乘拟合曲线为 现在很多计算机配有自动选择数学模型的程序，其方法与本例相同. 程序中因变量与自变量变换的函数类型较多，通过计算比较误差找到拟合得较好的曲线，最后输出曲线图形及数学表达式. 解 作散点图如右， 从右图可以看出这些点接近一条抛物线，因此设所求公式为 x 0 1 2 3 4 5 y 5 2 1 1 2 3 例 已知一组观测数据如表所示，试用最小二乘法求一个多项式拟合这组数据。 代入法方程 4.6 式有 解之可得 故所求拟合多项式为 由于 故cj 0对j 0,1, …,n成立.由此得?0,?1,…,?n线性无关，于是可直接得到正交多项式的以下性质. 　　 1 对任何多项式P x ∈Hn均可表示为函数族?0 x ,?1 x ,…,?n x 的线性组合，即 2 ?n x 与任何次数小于n的多项式P x ∈Hn-1都正交，即 3.2.2 勒让德 Legendre 多项式 当区间为 ，权函数 时， 并用 表示. 罗德利克 Rodrigul）给出了简单的表达式 （2.5） 正交化得到的多项式就称为勒让德 Legendre 多项式， 由 勒让德多项式重要性质： （2.7） 性质1 正交性 （2.8） 由于 是偶次多项式，经过偶次求导仍为 偶次多项式，经过奇次求导则为奇次多项式. 性质2 奇偶性 性质3 递推关系  在区间 内有 个不同的实零点. 性质4 性质2 （2.8） 由于 是偶次多项式，经过偶次求导仍为 偶次多项式，经过奇次求导则为奇次多项式，故 为偶数时 为偶函数， 为奇数时 为奇函数，于是（2.8）成 立. 奇偶性 性质3 递推关系  在区间 内有 个不同的实零点. 性质4 3.2.3 切比雪夫 Chebyshev 多项式 当权函数 ，区间为 时，由序 列 正交化得到的正交多项式就是切比雪夫 Chebyshev 多项式. 它可表示为 （2.10） 若令 ， 则 切比雪夫多项式有很多重要性质： 性质5 递推关系 （2.11） 性质6 切比雪夫多项式 在区间 上带权 （2.12） 正交，且 性质7 只含 的偶次幂， 只含 的奇次幂. 这个性质由递推关系直接得到. 可以用 的线性组合表示 ， 性质8 在区间 上有 个零点 其公式为 （2.13） 解 利用T5 x 的零点和区间变换可知节点 例4 求f x ex在[0, 1]上的4次拉格朗日插值多项式L4 x ，插值节点用T5 x 的零点，并估计误差 即 对应的拉格朗日插值多项式为 利用 2.15 式可得误差估计 而 于是有 　　在第2章中已经知道，由于高次插值出现龙格现象，一般Ln x 不收敛于f x ，因此它并不适用. 但若用切比雪夫多项式零点插值却可避免龙格现象，可保证整个区间上收敛. 见书p65例5. 3.2.5 其他常用的正交多项式 一般说，如果区间[a, b]及权函数? x 不同, 则由 1,x,…,xn,… 正交化得到的多项式也不同. 除上述两种最重要的正交多项式外，下面再给出三种较常用的正交多项式. 1. 第二类切比雪夫多项式 区间为[-1, 1]时, 取权函数 由序列 1,x,…,xn,… 正交化得到的多项式就称为第二类切比雪夫多项式，其表达式为 令x cosθ，可得 即 Un x 为[-1,1]上正交多项式族，且有递推关系式为 2. 拉盖尔多项式 区间为[0, +∞ 时, 取权函数 由序列 1,x,…,xn,… 正交化得到的多项式就称为拉盖尔 Laguerre 多项式，其表达式为 它也具有正交性质 和递推关系式 3. 埃尔米特多项式 区间为 -∞,+∞ 时, 取权函数 由序列 1,x,…,xn,… 正交化得到的多项式就称为埃尔米特 Hermite 多项式，其表达式为 它也具有正交性质 和递推关系式 3.3 最佳平方逼近 法方程组 /*normal equations */ Hilbert阵！ 解 由公式有 得法方程组 例 求 在 [0,1]上的一次最佳平方逼近多项式 注意与例4的区别 . 解出 故得所求的一次最佳平方逼近多项式为 平方误差为 最大误差为 若用 1,x,…,xn 做基，求最佳平方多项式，计算简单，但当n较大时, 系数矩阵H是高度病态的 见第5章 ,因此直接求解法方程是相当困难的，故通常是采用正交多项式做基底构造最佳平方多项式. 3.3.2 用正