版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上插值法和曲线拟合电子科技大学 摘要:理解拉格朗日多项式插值、分段线性插值、牛顿前插,曲线拟合,用matlab编程求解函数,用插值法和分段线性插值求解同一函数,比较插值余项;用牛顿前插公式计算函数,计算函数值;对于曲线拟合,用不同曲线拟合数据。关键字:拉格朗日插值多项式;分段线性插值;牛顿前插;曲线拟合引言:在数学物理方程中,当给定数据是不同散点时,无法确定函数表达式,求解函数就需要很大的计算量,我们有多种方法对给定的表格函数进行求解,我们这里,利用插值法和曲线拟合对函数进行求解,进一步了解函数性质,两种方法各有利弊,适合我们进行不同的散点函数求解。正文:一、插值法和分
2、段线性插值1拉格朗日多项式原理对某个多项式,已知有给定的k + 1个取值点:其中对应着的位置,而对应着函数在这个位置的取值。假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:其中每个为拉格朗日基本多项式(或称插值基函数),其表达式为:拉格朗日基本多项式的特点是在 上取值为1,在其它的点 上取值为0。2分段线性插值原理给定区间a,b, 将其分割成a=x0 <x1 <<xn =b, 已知函数y= f(x) 在这些插值结点的函数值为yk =f(xk)(k=0,1,n)求一个分段函数Ih(x), 使其满足:(1) Ih(xk )=yk ,(k=0
3、,1,n) ;(2) 在每个区间xk ,xk+1 上,Ih (x)是个一次函数。易知,Ih(x)是个折线函数, 在每个区间xk ,xk+1 上,(k=0,1,n) , 于是, Ih (x)在a,b上是连续的,但其一阶导数是不连续的。3拉格朗日插值多项式算法输入,令。对,计算4分段线性插值算法输入(x,y),k=0,1,n;计算5插值法和分段线性插值程序按下列数据分别作五次插值和分段线性插值,画出两条插值曲线以及给定数据点。求x1=0.32, x2=0.55, x3=0.68 时的函数近似值,并比较两种方法的插值余项。0.300.420.500.580.660.721.044031.084621
4、.118031.156031.198171,23223拉格朗日插值程序:function lagrintxi=0.32,0.55,0.68;%xi=0.2:0.001:0.8;x=0.3,0.42,0.50,0.58,0.66,0.72;y=1.04403,1.08462,1.11803,1.15603,1.19817,1.23223;L=zeros(size(y);m=length(xi);for i=1:m dxi=xi(i)-x; L(1)=prod(dxi(2:6)/prod(x(1)-x(2:6); L(6)=prod(dxi(1:6-1)/prod(x(6)-x(1:6-1); f
5、or j=2:6-1 num=prod(dxi(1:j-1)*prod(dxi(j+1:6); den=prod(x(j)-x(1:j-1)*prod(x(j)-x(j+1:6); L(j)=num/den; end yi(i)=sum(y.*L); fprintf('x=%f,y=%fn',xi(i),yi(i); endplot(xi,yi,'r');axis(0.2 0.8 1.03 1.24);hold on plot(x,y,'b.','markersize',20)grid on分段线性插值算法程序:function
6、y=div %xi=0.3:0.001:0.72;x0=0.3,0.42,0.50,0.58,0.66,0.72;y0=1.04403,1.08462,1.11803,1.15603,1.19817,1.23223;k=1;xi=0.32,0.55,0.68;for j=1:3 for i=1:5 if xi(j)>=x0(i) && xi(j)<=x0(i+1) && k<=3 lx(1)=(xi(j)-x0(i+1)/(x0(i)-x0(i+1); lx(2)=(xi(j)-x0(i)/(x0(i+1)-x0(i); y(k)=lx(1)*
7、y0(i)+lx(2)*y0(i+1); k=k+1; end endendplot(xi,y,'r');axis(0.2 0.8 1.03 1.24);hold on plot(x0,y0,'b.','markersize',20)grid on6运算结果拉格朗日插值结果x=0.,y=1.x=0.,y=1.x=0.,y=1.拉格朗日插值余项:分段插值结果ans =1.0508 1.1418 1.2095分段线性插值余项:由于拉格朗日插值的余项比分段线性插值的余项要求更为严格,点少、区间小的时候,拉格朗日插值要更好。但在区间较大、节点较多的时候,
8、分段线性插值要更好。二、牛顿前插1牛顿前插原理次牛顿前插公式:插值余项:,阶差分记作。 阶差商是差分和差商之间的关系是2牛顿前插算法输入。对,计算各阶差分计算函数值3牛顿前插程序:编写一个用牛顿前插公式计算函数值的程序,要求先输出差分表,再计算x点的函数值0.1250.2500.3750.5000.6250.7500.7960.7730.7440.7040.6560.602分别求x=0.158和x=0.636的三次插值的值,并比较二者的插值余项。这里以x=0.636为例function P=newtonchax0=0.636;X=0.125 0.250 0.375 0.500 0.625 0.
9、750;Y=0.796 0.773 0.744 0.704 0.656 0.602;h=abs(X(2)-X(1);n=find(abs(x0-X)<3*h);X=X(n(1):n(end);Y=Y(n(1):n(end);w=length(X);R=zeros(w,w);R(:,1)=Y(:);for k=2:w for j=k:w R(j,k)=R(j,k-1)-R(j-1,k-1); endendt=(x0-X(1)/h;T=1; for m=1:w-1 T=T*(t-m+1); N(m)=R(m+1,m+1)*T/factorial(m); endP=R(1,1)+sum(N);
10、4运行结果:差分表0.0000000000.0000-0.000000000.0000-0.0000-0.00000000.0000-0.0000-0.0000-0.0000000.0000-0.0000-0.00000.00000.000000.0000-0.0000-0.00000.0000-0.0000-0.0000X=0.636时ans =0.4000 x=0,158时ans =0.0000三、曲线拟合1曲线拟合原理:给定数据。记拟合函数的形式为(1.1),其中为已知的线性无关函数。求系数使得(1.2)取最小值。称(1.3)为拟合函数或经验公式。如果,则(1.3)为次最小二乘拟合多项式
11、2曲线拟合算法:已知数据对,求多项式,使得为最小。注意到此时,多项式系数满足下面的线性方程组:其中,然后只要调用线性方程组的函数程序即可3曲线拟合程序:试分别用抛物线y=a+bx2和指数曲线y=aebx拟合下列数据12.53.543.81.5026.033.0画出数据点和两条拟合曲线,并通过计算2个拟合函数残差向量的2范数来比较拟合优劣。用抛物线y=a+bx拟合程序:function ZXEx=1 2.52 3.52 42;y=3.8 1.50 26.0 33.0;m=1;S=zeros(1,2*m+1);T=zeros(m+1,1);for k=1:2*m+1 S(k)=sum(x.(k-1
12、);endfor k=1:m+1 T(k)=sum(x.(k-1).*y);endA=zeros(m+1,m+1);a=zeros(m+1,1);for i=1:m+1 for j=1:m+1 A(i,j)=S(i+j-1); endenda=AT;for k=1:m+1 fprintf('a%d=%fn',k,a(k);endp=polyfit(x,y,1);u=polyval(p,x);plot(sqrt(x),u,'b')hold onplot(sqrt(x),y,'b.')grid on指数曲线y=aebx拟合程序:function ZX
13、E2x=1 2.5 3.5 4;y=3.8 1.50 26.0 33.0;y=log(y);m=1;S=zeros(1,2*m+1);T=zeros(m+1,1);for k=1:2*m+1 S(k)=sum(x.(k-1);endfor k=1:m+1 T(k)=sum(x.(k-1).*y);endA=zeros(m+1,m+1);a=zeros(m+1,1);for i=1:m+1 for j=1:m+1 A(i,j)=S(i+j-1); endenda=AT;for k=1:m+1 fprintf('a%d=%fn',k,a(k);endp=polyfit(x,y,1);u=polyval(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 餐饮应急方案
- 老年人电脑培训内容
- 关于春节的调查报告15篇
- 地板铺设脚手架施工合同范本
- 通信线路以此合同为准
- 住宅小区造价管理聘用合同
- 污水处理防水防腐施工合同
- 家电公司办公区门窗定制合同
- 产业园办公用房托管合同
- 商业综合体简易工程施工合同
- 五年级数学上册8解方程课件
- 教学工作中存在问题及整改措施
- 内部项目跟投协议书(正)
- 钢管静压桩质量监理细则
- 5000头奶牛养殖场新建项目环境评估报告书
- 16飞机颠簸教学课件
- IATF16949-过程审核检查表-(含审核记录)-
- 散文阅读精练(100篇优秀散文每篇均附能力训练)-现代文精练
- 《成本会计》考试复习题库(浓缩300题)
- 工作成功案例分享模板
- 安全管理的几点做法1000字
评论
0/150
提交评论