《最小二乘法第八部分课件》由会员分享，可在线阅读，更多相关《最小二乘法第八部分课件（11页珍藏版）》请在装配图网上搜索。
1、最小二乘法第八部分课件 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望oyx 偏差)(iiixfyr有正有负,值都较小且便于计算,可由偏差平方和最小 min)(20iinixfy为使所有偏差的绝对来确定近似函数 f(x).最小二乘法原理最小二乘法原理:设有一列实验数据分布在某条曲线上,通过偏差平方和最小求该曲线的方法称为最小二乘法最小二乘法,找出的函数关系称为经验公式经验公式.),1,0(),(nkyxkk,它们大体 机动 目录 上页 下页 返回 结束 特别,当数
2、据点分布近似一条直线时,问题为确定 a,b 令min)(20bxayknkk),(baMaM0)(20kknkkxbxaybM0)(20bxayknkkbxay满足:使oyx得axnkk02bxnkk0nkkkyx0axnkk0bn)1(nkky0解此线性方程组即得 a,b称为法方程组机动 目录 上页 下页 返回 结束 例例1.为了测定刀具的磨损速度,每隔 1 小时测一次刀具的厚度,得实验数据如下:找出一个能使上述数据大体适合的经验公式.解解:通过在坐标纸上描点可看出它们大致在一条直线上,列表计算:故可设经验公式为(P67 例1)bxayoyt27.0 26.8 26.5 26.3 26.1
3、25.7 25.3 24.8 0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7(mm)iy(h)iti机动 目录 上页 下页 返回 结束 得法方程组a140b28a285.2088 b717解得,125.27,3036.0ba故所求经验公式为125.273036.0)(ttfyit0i2itiyiity70 0 27.0 07 49 24.8 137.628 140 208.5 717.0机动 目录 上页 下页 返回 结束 为衡量上述经验公式的优劣,计算各点偏差如下:oyt称为均方误差,Mn1对本题均方误差124.071M它在一定程度上反映了经验函数的好坏.机动 目录 上页 下页
4、 返回 结束 108165.0)(270iiitfyM偏差平方和为27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 0 1 2 3 4 5 6 7/mm)(itf算得的/mmiy实测的i27.125 26.518 25.911 25.30326.821 26.214 25.607 25.000 0.125 0.018 0.189 0.0030.021 0.086 0.093 0.200)(iitfy oyt称为均方误差,Mn1对本题均方误差124.071M它在一定程度上反映了经验函数的好坏.机动 目录 上页 下页 返回 结束 108165.0)(270iiitfy
5、M偏差平方和为i27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 0 1 2 3 4 5 6 7/mm)(itf算得的/mmiy实测的27.125 26.518 25.911 25.30326.821 26.214 25.607 25.000 0.125 0.018 0.189 0.0030.021 0.086 0.093 0.200)(iitfy 例例2.在研究某单分子化学反应速度时,得到下列数据:57.6 41.9 31.0 22.7 16.6 12.2 8.9 6.5 3 6 9 12 15 18 21 241 2 3 4 5 6 7 8iyii其中 表示从
6、实验开始算起的时间,y 表示时刻 反应 物的量.试根据上述数据定出经验公式).(fy(P70例2)解解:由化学反应速度的理论知,经验公式应取meky 其中k,m 为待定常数.对其取对数得kmylnlnkbmaXyYln,ln令bXaY(线性函数)(书中取的是常用对数)机动 目录 上页 下页 返回 结束 oy因此 a,b 应满足法方程组:akk812bkk8181lnkkkyakk81b881lnkky经计算得 994.2801081836ba714.238108ba解得:364.4,104.0ba,57.78,104.0bekm所求经验公式为104.057.78ey其均方误差为135.071M机动 目录 上页 下页 返回 结束 观测数据:用最小二乘法确定a,b 通过计算确定某些经验公式类型的方法通过计算确定某些经验公式类型的方法:),1,0(),(niyxii),2,1(,11niyyyxxxiiiiii令,)1(定值若iixybxay则考虑,lnln)2(定值若iixybxay 则考虑,ln)3(定值若iixyxbeay 则考虑axbylnlnln转化为axbylnln转化为机动 目录 上页 下页 返回 结束 作业 (习题8-10)P72 1,2习题课 目录 上页 下页 返回 结束
最小二乘法第八部分课件
