数值计算方法编写lagrange插值的程序.docVIP

实验课程名称 数值计算方法 实验项目名称编写lagmnge插值的程序 年 级 专 业 学生 姓名 学 00学院 实验时间： 年 月 姓 名 学号 实验组 实验吋间 指导教师 成缋 实验项H名称 编写lagrange插值的程序 实验目的及要求： 掌握lagrange插值方法的思想，学会应用1 agrange插值方法对数值进行估算； 要求使用C语言或mat lab编程求解。 实验(或算法〉原理： 通过Ltzgm/职插值法的学习，掌握如何根据己知函数表 X X,…X” y y0 yt ??- x, ” n X— Y . 构造Ug醐ge插值多项式= 7 )V.用/(x) = Ux). J^xrXj 1 实验硬件及软件平台： MATLAB 7.10. 0 (R2010a) 实验步骤： 1、 分析问题； 2、 确定程序主线，进行模块化; 3、 编写程序; 4、 调试，优化。 实验内容(包括实验具体内容、算法分析、源代码等等): 脚本文件lan. m function zou=lan(j,x,n) syms z; s = l; p=l; for i=l:n if i==j continue else S=s*(x(j)-x (i)); p=p* (Z-X (i))； end end zou=p/s; 脚本文件chazhi. m clear,clc; answer=0; n=input (请输入点的个数：*); for i=l:n disp (请输入第 1) z disp (i) , x (i) =input (’ 个点的横坐标：’)； y (i) =input ( 1 纵坐标：1 ); end for j=l:n answer=answer+y(j)*lan(j, x, n); end answer a=x (1); b=x(1); for i=l:n if a〉x (i) a=x (i); end if bx (i) b=x(i); end end plot(x,y,▼r*1); hold on ezplot(answer,[a-l,b+l]); title (1 Lagrange插值。1); 实验结果与讨论： 1、在命令窗口中输入chazhi,回车; ▲ UITU4 I* bu 0^9 kW ?J J 龜今駟 M 3T r- W omumt ojbmt* 蠢壽 mr 好TTQMTua ? UID ?I 鱸叫? haw cl ??lit 产 SHIHH ?TM ■ ■u?*.rec? - rbM4 泰 iUK *— 5、根裾提示输入已知点的横纵坐标; V?bv iM ?? Eile Edit Debug BMalld Desktop Window Udp faai篇》龜今cm 抽3珀 似 Cuncnt Fokkc D:\360daU\BB教格衣的文 MATIA8 ， 3^5 MATJ13 710.0 (R2010a) a Current f old? X * 0 X * ? MAILAB ，fl令， L N4W Shortcuts ￡J How to Add ￡i WhM sNcw ibcd-mangle-m ch?hu?v3 ch*zhuninvestm larvm *J iunhM?nthuru-fnIjxoum么 zouzou^svO zouzou m Deuus 4 Start 3、根据提示输入数据，回车。这样就得到/插值多项式和图象 Coovnand Window ?*入?的个B: 11 濤 ??入麻 1 个点的?坐?: -5 ttfttFi 0.0016 ??入甬 6 2 个k5的購生林：- 纵坐抹：0.002S ?編入《 3 个点的櫬去株：-3 级轚标：0. 0044 ??入霣 4 个的?坐抹：-2 U至特：0.00的 入兩 5 个?的?生标：-1 纰全梓：0.0385 ??入* A 6 W * 1 ， ? Worksite ISUmc * Value M 君:MMmr ?5 ?5 Sb 5 5 B- U 11 ±b 11 11 Bn 11 11 1U1 double* -5 fly 1x11 double* 0A ? ?r Command Httfcxy -H □ X ■H □ X Sdect dM? to plot ▼ 0.0385 0 1 0.0385 2 0.0099 3 0.00“ 4 0.0025 5 0.0016 cs OVR 签名: 签名: 指导教师意见: