拉格朗日插值法_第1页
拉格朗日插值法_第2页
拉格朗日插值法_第3页
拉格朗日插值法_第4页
拉格朗日插值法_第5页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

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

文档简介

1、拉格朗日抛物线插值法1、定义若多项式lj（j=0,1,2.n）在n+1个节点x0x1.x=100 121 144 ; y=10 11 12 ; y2=lagrage(x,y,115) ; 输出 y2=10.72282)计算解：xy0.50.7070.866在Matlab窗口输入x= ; y=0.5 0.707 0.866 ; y2=lagrage(x,y,2PI/9) ; 输出 y2=0.6380均差与牛顿插值多项式1、 1）定义称 为函数关于的一阶均差， 称为的二阶均差。一般的，称为的k阶均差（均差也称为差商）。 2）牛顿插值公式推导根据均差定义，把x看成a,b上一点，则有把后一式带入前一式

2、可得最后一项中，均差部分含有x，为余项部分，记为而前面n+1项中，均差部分都不含有x，因而前面n+1项是关于x的n次多项式这就是牛顿插值公式。于是上式记。2、 Matlab文件 Newton int.mfunction n=length(x)D=Zeros(n.n);D(:,1)=for k=1:n-1;for i=1:n-kD(i+k,k+1)=(D(i+k,k)-D(i+k+1,k)/(x(i+k)-x(i);endendfor i=1:n;z=1for k=1:i-1; z=z*();endend3、 例题x1-12y0-34求二次插值解： x y 一次差商 二次差商 -1 -3 1 0

3、 2 4 4 N(x)=-3+3/2(x+1)+1/6(x+1)(x-1)在Matlab窗口输入x=-1,1,2;Y=-3 0 4;曲线拟合的最小二乘法1、 Matlab文件 polyfit.mfunction p=polyfit(x,y,n)A=Zeros(n+1,n+1)for t=0:nfor j=0:n A(i+1,j+1)=sum(x.(i+j)endb(i+1)=sum(I,*y)end a=Ab;p=fliplr(a)2、 例题设数据由表给出，表中第4行为，可以看出数学模型为，用最小二乘法确定a和b。i012341.005.101.6291.255.791.7561.506.531.8761.757.452.0082.008.462.135解：根据给定数据描图可确定拟合曲线方程为，它不是线性形式。两边取对数得，若令，则得。为确定A，b，先将转化为，数据见上表。根据最小二乘法，取，得 故有法方程5A+7.50b=9.404750A+11.875b=14.422解得A=1.122,b=0.505,a=于是得最小二乘拟合曲线为在Matlab文件输入x=1:0.25:2y=1

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

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

最新文档

评论

0/150

提交评论