第3章函数逼近与计算(4551KB).ppt

于是得到关于c1,c2,…,cn的方程组 法方程组(或正规方程组) （8） 记 ? 多项式拟合问题 若取 称为多项式拟合。 则得到的拟合函数 ? 非线性拟合 （1）双曲拟合 （2）对数拟合 （3）指数拟合 （4）幂函数拟合 t=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]; y=[4.00,6.40,8.00,8.80,9.22,9.50,9.70,9.86, 10.00,10.20,10.32,10.42,10.50,10.55,10.58,10.60]; Plot(t,y,‘.’); y1=11.3253*10^(-3)*exp(-1.0567./t); Plot(t,y1); 作散点图 举例 给定 的一组数据 求拟合函数 。 解： 作图可知所有点都分布在一条直线的附近，即拟合函数 近似为一个线性函数，故可设 将数据带回，即可得拟合函数 4.用正交函数作最小二乘拟合 给定 与正交函数系 可利用正交化手续进行构造 则最小二乘拟合函数为： 其中， * 连续函数 在 上线性无关的 充分必要条件是它们的克莱姆(Gram)行列式 定理 其中， 是任意实数，则 并称 是生成集合 的一个基底。 的全体是 的一个子集，记为 设 是 上线性无关的连 续函数, 对任意的 为 的最佳平方逼近元。 2、最佳平方逼近元的定义 设 为线性内积空间， 为 上 个 线性无关元， 记由 张成的 的子空间为 ， 即 定义 在 的子空间 中， 求 的在2-范数 意义下的最佳逼近元 ， 即求 ，使不等式 对任意 成立. 若满足上式的 存在， 称 3.最佳平方逼近元的存在性 定理1 设 为线性内积空间， 由线性无关组 张成的线性空间 为 的子空间， 存在 为 的最佳平方逼近元. 则对任意的 Remark: 线性内积空间 的子空间 的线性无关组 选取不同， 在 中求得的对 的最佳 平方逼近元 也不同，求解 的难易程度也不同。 4. 最佳平方逼近元的特征 定理2 内积空间) 为 的最佳平方逼近元的充要条件是： (线性 与一切 正交。 其中， 为 的 个线性无关元。 REMARK： 定理2中所说的 与一切 正交， 与一切 的内积等于零， 是指 即 证： 必要性. 用反证法. 设 为 的最佳平方逼近元， 不与所有的 正交. 但 即存在 使得 则 令 所以 必须与一切 正交. 且 这说明 不是对 的最佳平方逼近元， 与假设条件矛盾， 充分性. 仍记 则对任意的 ，有 而 对任意 成立, 即 为 的最佳平方逼近元。 所以 进而有 5.最佳平方逼近元的惟一性 定理3 线性内积空间 的子空间 中若存在对 的 最佳平方逼近元，则惟一. 6. 最佳平方逼近元的求解 现假定线性内积空间 上的内积已定义， 并且 的 子空间的一组基底 也确定， 最佳平方逼近元. 那么，对具体的被逼近元 如何求 使其为 的 由最佳平方逼近元的充要条件， 若假定 则可以得出 其中， 为待定系数。 恒等变形为 用矩阵式表示这个方程组为 此方程组称为法方程组。 若所选取的一组基底 满足 则称其为正交基，此时 7. 最佳平方逼近元的误差估计 设 的最佳平方逼近元，则其平方误差为： 为 均方误差为： 五、连续函数的最佳平方逼近 1. 对于给定的函数 要求函数 使 若这样的 存在， 上的最佳平方逼近函数。 则称为 在区间 特别地，若 则称 为 在 上的 次最佳平方逼近多项式。 求最佳平方逼近函数 的问题 可归结为求它的系数 ,使多元函数 取得极小值。 由于 是关于 的二次函数， 故利用多元函数取得极值的必要条件，可得 (k = 0, 1, 2, …, n) 于是可得方程组 即 如采用函数内积记号，记 则方程组可以简写为 写成矩阵形式为 法方程组！ 由于?