数值计算方法实验报告1.doc_第1页
数值计算方法实验报告1.doc_第2页
数值计算方法实验报告1.doc_第3页
数值计算方法实验报告1.doc_第4页
全文预览已结束

下载本文档

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

文档简介

长春理工大学学生实验报告学院名称计算机科学与技术专业班级学号学生姓名实验日期2014.09.12成绩课程名称数值计算方法实验题目实验1 代数插值一、实验目的和要求使用拉格朗日插值法或牛顿插值法求解：已知f(x)在6个点的函数值如下表所示，运用插值方法，求f(0.54)的近似值。X0.400.500.600.700.80f(x)0.410750.578150.696750.888111.02652二、主要设备 PC，Windows操作系统，VC+6.0编程平台；三、实验内容和原理设函数在区间a,b上n+1互异节点x0,x1,xn上的函数值分别为y0,y1,yn,求n次插值多项式Pn(x),满足条件Pn(xj)=yj, j=0,1,n令Ln(x)=y0l0(x)+y1l1(x)+ynln(x)= yili(x)其中l0(x),l1(x), ln(x) 为以x0,x1,xn为节点的n次插值基函数，则Ln(x)是一次数不超过n的多项式，且满足Ln(xj)=yj, L=0,1,n再由插值多项式的唯一性，得Pn(x)Ln(x) 流程图如下：四、操作方法与实验步骤#include #include #include #include void difference(float *x,float *y,int n) float *f; int k,i; f=(float *)malloc(n*sizeof(float); for(k=1;k=n;k+) f0=yk; for(i=0;ik;i+)fi+1=(fi-yi)/(xk-xi); yk=fk; return; int main() int i,n; float x20,y20,xx,yy; printf(请输入数据个数n：); scanf(%d,&n);printf(n); for(i=0;i=0;i-)yy=yy*(xx-xi)+yi; printf(n近似值为：F(%f)=%fn,xx,yy); 五、实验结果与分析 分析： 拉格朗日插值的优点是插值多项式特别容易建立，缺点是增加节点是原有多项式不能利用，必须重新建立，即所有基函数都要重新计算，这就造成计算量的增加。牛顿插值法则很好地避免了上述问题。五、 讨论、心得 本实验有两种插值方法可以选用，由于时间关系，最终选用牛顿插值法。若是下去有时间的话，可以再用拉格朗日插值
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

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

最新文档

评论

0/150

提交评论