数值积分,最小二乘法详解.ppt

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
华长生制作 * 例3. 在某化学反应里,测得生成物浓度y%与时间t的 数据如下,试建立y关于t的经验公式 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 解: 具有图示的图形的曲线很多,本题特提供两种形式 华长生制作 * 两边取对数,得 得 即为拟合函数 基函数为 解法方程组得 平方误差为 华长生制作 * 用最小二乘法得 即 无论从图形还是从平方误差考虑 在本例中指数函数拟合比双曲线拟合要好 平方误差为 华长生制作 * 三、加权最小二乘法 各点的重要性可能是不一样的 重度: 即权重或者密度,统称为权系数 定义加权 平方误差为 ----- 9 华长生制作 * 使得 华长生制作 * 由多元函数取极值的必要条件 得 即 华长生制作 * 同样引入记号 定义加权内积 ----- 10 华长生制作 * 矩阵形式 法方程组 为 方程组 10 式化为 ----- 11 --- 12 华长生制作 * 平方误差为 作为特殊情形,用多项式作拟合函数的法方程组为 ----- 13 华长生制作 * 四、用正交多项式作最小二乘拟合 即 正交多项式如何选取呢 ----- 14 华长生制作 * 使得 华长生制作 * 由 可知 因此 华长生制作 * 而 因此 华长生制作 * 可知 最后可得正交多项式选取的方法: ----- 15 由 华长生制作 * 使得 由正交多项式的性质,法方程组 华长生制作 * ----- 16 ----- 17 可化为 即 得 即 为利用正交多项式的最小二乘解 华长生制作 * 平方误差为 华长生制作 * 例4. 是用最小二乘法求拟合这组数据的多项式 解: 从散点图可知 数据和二次多项式拟合较好 因此选用二次多项式作 这组数据的拟合函数 华长生制作 * 设拟合函数 取 华长生制作 * 华长生制作 * 因此拟合多项式为 平方误差为 * * * * 3.3最佳平方逼近 华长生制作 * 3.3.1最佳平方逼近问题的求解 连续函数空间C[a,b]上定义了内积就形成了一个内积空间, 对内积空间任一元素f x ∈ C[a,b],也可用线形无关的基表示。 则称 是f x 在 中的最佳平方逼近函数。 下面我们先讨论在区间[a,b]上 一般的最佳平方逼近问题。 设 是C[a,b]中的线性无关函数,记 对于f x ∈C[a,b],若存在 ,使得 求 等价于求多元函数 的极小值。利用多元函数求极小值的必要条件有 华长生制作 * 按内积的定义,上式可写为 这是关于 的线性方程组,称为法方程或正规方程组。   由于 线性无关,故法方程的系数矩阵非奇异,于是方程组有唯一解 。从而得到 即对任意 ,有 事实上,有 因此,对任意 ,有 ,从而也有 华长生制作 * 于是 这就证明了结论,从而也证明了f在 中的最佳平方逼近的存在唯一性。 称最佳逼近的误差 为平方误差。   考虑特殊情形,设[a,b] [0,1], 。对于f∈C[a,b], 在 中最佳平方逼近多项式可以表示为 华长生制作 * 相应于法方程中的系数矩阵为 称之为Hilbert矩阵 例1 设 ,求[0,1]上的一次最佳平方逼近多项式。 解 由于 得方程组 华长生制作 * 解得c0 0.394,c1 0.246。从而最佳平方逼近为 平方误差 由于Hilbert矩阵是病态的(见第5章),用 作基时,求法方程的解,舍入误差很大。实用的办法是采用正交多项式作基。 华长生制作 * 若 是正交多项式组,则法方程系数矩阵为 基于正交基函数的最佳平方逼近 华长生制作 * 于是f x 的最佳平方逼近多项式为 例2 设 f x ex ,在[-1,1]上用勒让德多项式作f的三次多次最佳平方逼近多项式。 解 用勒让德多项式 P x k 0,1,2,3 , 可得 华长生制作 * 于是最佳平方逼近为 。 平方误差 华长生制作 * 3.4曲线拟合的最小二乘方法 华长生制作 * 实例:考察某种纤维的强度与其拉伸倍数的关系,下表 是实际测定的24个纤维样品的强度与相应的拉伸倍数 是记录: 曲线拟合模型及其求解 华长生制作 * 纤维强度随拉伸 倍数增加而增加 并且24个点大致分 布在一条直线附近

文档评论(0)

1192212 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档