数值分析函数逼近公开课一等奖优质课大赛微课获奖课件.pptx

第六章函数迫近§1数据拟合最小二乘法§3函数最佳平方迫近§2正交多项式1第1页第1页

?Lagrange插值与最小二乘迫近图像描述2第2页第2页

?办法1:用3次Lagrange插值多项式近似x,y函数关系.?为何要用最小二乘迫近.xiyi24681.12.84.97.2例给定一组试验数据下列求x,y函数关系.?办法2:用直线来近似x,y函数关系.?3第3页第3页

?用直线y=a0+a1x来反应x,y之间函数关系.如何选取a0,a1?才干使直线最好地反应数据点基本趋势??残差向量?残差4第4页第4页

?衡量近似函数好坏原则：残差向量大小(1)使残差绝对值之和最小,即(2)使残差最大绝对值最小,即(3)使残差平方和最小,即最佳平方迫近或数据拟合最小二乘法最佳一致迫近5第5页第5页

问题:给定n个数据点(xi,yi)(i=1,2,…,n)求直线y=a0+a1x使得达到最小.?最小二乘一次多项式拟合§1数据拟合最小二乘法6第6页第6页

?令则原问题等价于求a0,a1使F(a0,a1)达到最小.利用多元函数取极值必要条件得正则方程组7第7页第7页

?由上式求得a0,a1,代入y=a0+a1x得到最小二乘拟合(直线)一次多项式.8第8页第8页

xiyi24681.12.84.97.2例给定一组试验数据下列求x,y函数关系.解正则方程组9第9页第9页

直线拟合误差很大抛物线拟合效果更加好10第10页第10页

问题:给定n个数据点(xi,yi)(i=1,2,…,n)求使得达到最小.?最小二乘二次多项式拟合11第11页第11页

?令则原问题等价于求a0,a1,a2,使F(a0,a1,a2)达到最小.利用多元函数取极值必要条件得12第12页第12页

?用Cholesky分解法求此对称正定阵?用MATLAB函数z=A\r?由上式求得a0,a1,a2,得到最小二乘拟合二次多项式正则方程组13第13页第13页

?最小二乘三次多项式拟合正则方程组14第14页第14页

?最小二乘m次多项式拟合(mn)正则方程组15第15页第15页

?指数拟合假如数据点(xi,yi)(i=1,2,…,n)分布近似指数曲线,则可考虑用指数函数去拟合数据.但是这是一个关于a,b非线性模型,故应通过适当变换,将其化为线性模型,然后利用最小二乘法求解.为此,对指数函数两端取对数,得16第16页第16页

则数据组(xi,yi)(i=1,2,…,n)最小二乘拟合指数曲线为这表明(xi,lnyi)(i=1,2,…,n)分布近似于直线,求出此数据组最小二乘拟合直线17第17页第17页

xiyi例给定一组试验数据下列求x,y函数关系.12346782367532?(1)作散点分布图点分布近似为抛物线18第18页第18页

?(2)拟定近似表示式设拟合曲线为二次多项式?(3)建立正则方程组19第19页第19页

故正则方程组为?(4)求解正则方程组得故所求拟合曲线为20第20页第20页

xiyi例给定一组试验数据下列求x,y函数关系.12346782367532?Matlab解法:polyfit([1,2,3,4,6,7,8],[2,3,6,7,5,3,2],2)ans=-0.38643.4318-1.318221第21页第21页

例测得一发射源发射强度I与时间t一组数据下列tiIi0.20.30.40.50.60.70.83.162.381.751.341.000.7