函数的数值逼近正规版资料_第1页
函数的数值逼近正规版资料_第2页
函数的数值逼近正规版资料_第3页
函数的数值逼近正规版资料_第4页
函数的数值逼近正规版资料_第5页
已阅读5页,还剩19页未读, 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第四章函数的数值(shùzí)逼近第一页,共24页。什么(shénme)是插值与拟合?有何区别?0.00.20.40.60.811.21.41.61.80.511.53420-1144.1代数(dàishù)插值多项式第二页,共24页。求一个(yīɡè)经验函数y=g(x),使g(xi)=f(xi),i=1,…n.oxy●●●●●

y0

x1

x2

xn

y1

y2

yn

x0y=f(x)g(x)第三页,共24页。Lagrange插值线性插值(n=1)求次数(cìshù)≤1的多项式L1(x). 满足条件L1(x0)=y0,L1(x1)=y1,y=f(x)y=L1(x)x0x1xy第四页,共24页。第五页,共24页。类似的可以(kěyǐ)得到l1(x),l2(x)l0(x),l1(x),l2(x)称为(chēnɡwéi)以x0,x1,x2为节点的插值基函数。二次插值(n=2)求次数(cìshù)≤2的多项式L2(x),使其满足条件L2(x0)=y0,L2(x1)=y1,L2(x2)=y2第六页,共24页。令Ln(x)=l0(x)y0+l1(x)y1+…+ln(x)yn

求n次多项式lj(x),(j=0,1,…,n)使其满足条件n次插值多项式:求次数(cìshù)≤n的多项式Ln(x),使其满足Ln(x0)=y0,Ln(x1)=y1,......,Ln(xn)=yn..(7)容易(róngyì)求得lj(x)(j=0,1,…,n)称为(chēnɡwéi)以x0,x1,...,xn为节点的Lagrange插值基函数。第七页,共24页。 ....(9)公式(gōngshì)(9)就是n次Lagrange插值多项式.特点(tèdiǎn):构造容易,L-型插值基函数理论上有意义,但增加节点要重新计算,不适合编程计算。实际应用:只用低次插值。第八页,共24页。第九页,共24页。第十页,共24页。如何(rúhé)估计插值的截断误差?4.2多项式插值的进一步分析(fēnxī)第十一页,共24页。第十二页,共24页。生成(shēnɡchénɡ)Langrange插值多项式:x=1:6;y=[161821171512];disp([x;y])u=.75:.05:6.25;v=polyinterp(x,y,u);plot(x,y,'o',u,v,'-');第十三页,共24页。4.5曲线拟合的最小二乘方法(fāngfǎ)4.3分段(fēnduàn)插值与保形插值(略)4.4样条函数(hánshù)插值(略)第十四页,共24页。最小二乘原理(yuánlǐ)第十五页,共24页。例4.3直线(zhíxiàn)拟合第十六页,共24页。第十七页,共24页。线性相关与线性无关(wúguān)第十八页,共24页。最小二乘原理(yuánlǐ)的实质第十九页,共24页。求一个(yīɡè)经验函数y=g(x),使g(xi)=f(xi),i=1,…n.类似的可以(kěyǐ)得到l1(x),l2(x)线性插值(n=1)求次数(cìshù)≤1的多项式L1(x).第二十三页,共24页。求n次多项式lj(x),(j=0,1,…,n)使其满足条件公式(gōngshì)(9)就是n次Lagrange插值多项式.满足条件L1(x0)=y0,L1(x1)=y1,第二十二页,共24页。4样条函数(hánshù)插值(略)第四章函数的数值(shùzí)逼近,Ln(xn)=yn.disp([x;y])lj(x)(j=0,1,…,n)称为(chēnɡwéi)以x0,x1,.如何(rúhé)估计插值的截断误差?3直线(zhíxiàn)拟合生成(shēnɡchénɡ)Langrange插值多项式:第二十页,共24页。第二十一页,共24页。第二十二页,共24页。4样条函数(hánshù)插值(略)1代数(dàishù)插值多项式3直线(zhíxiàn)拟合公式(gōngshì)(9)就是n次Lagrange插值多项式.3直线(zhíxiàn)拟合l0(x),l1(x),l2(x)称为(chēnɡwéi)以x0,x1,x2为节点的插值基函数。第二十三页,共24页。第二十二页,共24页。公式(gōngshì)(9)就是n次Lagrange插值多项式.第二十三页,共24页。y=[161821171512];,Ln(xn)=yn

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论