信息与计算张椿木.docVIP

数值分析实验报告 学 院：数学与计算科学 专 业：信息与计算科学 班 级：信计一班 学 号： 2012750217 姓 名： 张椿木 指导老师：曹学年 2014年6月14日 湘潭大学数学与计算科学学院 学生实验报告 专业 数学与计算科学 班级 1 姓名 张椿木 学号 2012750217 实验日期 2014 年 3 月 20 日 至 2014 年 6 月 20 日 实验项目 多项式插值的震荡现象 实 验 内 容 与 要 求 实验内容： 考虑区间[-1，1]的一个等距划分，分点为： x(i)=-1+2i/n,i=0,1,2…,n 泽拉格朗日插值多项式为： L(x)=∑l(i)(x)/(1+25x(j)^2 ) i=0,1,…n 其中l(i)(x), i=0,1,…n,n是n次拉格朗日插值基函数。 实验要求： ⑴ 选择不断增大的分点数目n=2,3…，画出f(x)及插值多项式函数 L(x)在[-1，1]上的图象，比较分析实验结果。 （2）选择其它的函数，例如定义在区间[-5，5]上的函数 h(x)=x/(1+x^4) , g(x)=arctanx 重复上述的实验看其结果如何。 （3）区间[a,b]上切比雪夫点的定义为： xk=（b+a）/2+((b-a)/2)cos((2k-1)π/(2(n+1))),k=1,2,^,n+1， 以x1,x2^x(n+1)为插值节点构造上述各函数的拉格朗日插值多项 式，比较其结果。 注：实验结果、参考文献及相应程序附后 指导教师意见 成绩： 签名： 年 月 日 备注 注：实验结果、参考文献及相应程序附后 一 实验题目与设计 插值法 实验过程： 程序： 多项式插值的震荡现象（实验2.1） for m=1:6 subplot(2,3,m) %把窗口分割成2*3大小的窗口 largrang(6*m) %对largrang函数进行运行 if m==1 title(longn=6) elseif m==2 title(longn=12) elseif m==3 title(longn=18) elseif m==4 title(longn=24) elseif m==5 title(longn=30) elseif m==6 title(longn=36) end %对每个窗口分别写上标题为插值点的个数 end 保存为：chazhi.m function largrang(longn) mm=input(please input mm(运行第几个函数就输入mm为几):mm=) if mm==1 %d表示定义域的边界值 d=1; elseif mm==2||mm==3 d=5; end x0=linspace(-d,d,longn); %x的节点 if mm==1 y0=1./(1.+25.*x0.^2); elseif mm==2 y0=x0./(1.+x0.^4); elseif mm==3 y0=atan(x0); end x=sym(x);