数值分析拉格朗日插值法上机实验报告_第1页
数值分析拉格朗日插值法上机实验报告_第2页
数值分析拉格朗日插值法上机实验报告_第3页
全文预览已结束

下载本文档

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

文档简介

1、课题一:拉格朗日插值法1. 实验目的1.学习和掌握拉格朗日插值多项式2运用拉格朗日插值多项式进行计算2. 实验过程作出插值点(1.00,0.00),( -1.00,-3.00),( 2.00,4.00)二、算法步骤已知:某些点的坐标以及点数。输入:条件点数以及这些点的坐标。输出:根据给定的点求出其对应的拉格朗日插值多项式的值。3. 程序流程:(1)输入已知点的个数;(2)分别输入已知点的X坐标;(3)分别输入已知点的Y坐标;程序如下:#in elude <iostream>#in elude vconi o.h>#i nclude <malloc.h>/*拉格朗日

2、插值算法float lagrange(float *x,float *y,float xx,int n)*/ int i,j;float *a,yy=0.0;/*a作为临时变量,记录拉格朗日插值多项*/a=(float*)malloc (n *sizeof(float);for(i=0;i<=n-1;i+)ai=yi;for(j=0;j< 二n-1;j+)if(j!=i)ai*=(xx-xj)/(xi-xj);yy+=ai;free(a);return yy;int mai n()int i;int n;float x20,y20,xx,yy;printf("Input

3、n:");sca nf("%d",&n);if(n <=0)printf("Error! The value of n must in (0,20)."); getch();return 1;for(i=0;i< 二n-1;i+)pri ntf("x%d:",i);sca nf("%f", &xi);prin tf("n");for(i=0;i< 二n-1;i+)prin tf("y%d:",i);sca nf("%f&q

4、uot;, &yi);prin tf("n");prin tf("I nput xx:");sca nf("%f", &xx);yy=lagra nge(x,y,xx ,n);prin tf("x=%f,y=%fn",xx,yy);getch();举例如下: 已知当x=1 , -1,2时f(x)=0 , -3,4 ,求f(1.5)的值。运行结果如下:In put n:3x0:1x1:-1x2:2y0:0y1:-3y2:4In put xx:1.5x=1.500000,y=1.7916673、实验总结拉格朗日插值模型简单,结构紧凑,是经典的插值法。但是由 于拉格朗日的插值多项式和每个节点都有关 ,当改变节点个数时, 需要重新计算。且当增大插值阶数时容易出现龙格现象。在物理化学,资产价值鉴定工作和计算某一时刻的卫星坐标和 钟差等这些方面可以应用Lagrange插值。采用拉格朗日插值法计算设 备等功能重置成本,计算精度较高,方法快捷。但是这方法只能针对 可比性较强的标准设备,方法本身也只考虑了单一功能参 数,它的应 用范围因此受到了一定的限制。作

温馨提示

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

评论

0/150

提交评论