实验报告七-一
计 计 算 机 课 程 实 验 报 告 2009—2010 学年度 第二学期 院(系):数学与计算机科学学院 实验课程 数值分析 班 级 10 计本 学 号 07 姓 名 赵霞 指导教师 陈耀庚 实验题目 拉格朗日多项式逼近日 期 2012 年 11 月 2 日 实验目的及要求 使用 MATLAB 求出拉格朗日多项式 一、实验 程序:
:
function [C,L,y]=lagran(X,Y)%拉格朗日多项式逼近%x,y 为差值节点 w=length(X);L=zeros(w,w);for k=1:w V=1;for j=1:w if k~=j V=conv(V,poly(X(j)))/(X(k)-X(j))end end L(k,:)=V;end C=Y*L;y=poly2sym(C,"x");二、实验结果:
:
x=[1 2 3] x = 1 2 3 >> y=[5 6 8] y =
5 6 8 >> lagran(x,y)V =-1 2 V = 0.5000-2.5000 3.0000 V = 1-1 V =-1 4-3 V = 0.5000-0.5000 V = 0.5000-1.5000 1.0000 ans = 0.5000-0.5000 5.0000 三、出现的问题:
四、实验总结:
注:页面不够可以附页
