实验报告3(插值与拟合)..docVIP

淮海工学院实验报告 班级 精算111 姓名 张琪 学号 2011122373 指导教师 张滦云 课程名称 数学建模与实验 成绩 实验室 数学建模实验室 实验项目 3、插值与拟合 同组实验者 实验日期 2013-3-22 一、实验目的 二、实验所用软件及版本 MATLAB 7.1 三、实验内容1、 function y=lagr(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0; for k=1:n p=1; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 2. x0=[0 3 5 7 9 11 12 13 14 15]; y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; x=0:0.1:15; y1=lagr(x0,y0,x); y2=interp1(x0,y0,x); y3=interp1(x0,y0,x,spline); [x,y1,y2,y3] subplot(3,1,1) plot(x0,y0,k+,x,y1,r) grid,title(lagrange) subplot(3,1,2) plot(x0,y0,k+,x,y2,r) grid,title(piecewise linear) subplot(3,1,3) plot(x0,y0,k+,x,y3,r) grid,title(spline) ans = 0 0 0 0 0.1000 -4.9437 0.0400 0.0499 0.2000 -8.8200 0.0800 0.0990 0.3000 -11.7726 0.1200 0.1474 0.4000 -13.9310 0.1600 0.1951 0.5000 -15.4117 0.2000 0.2421 0.6000 -16.3190 0.2400 0.2884 0.7000 -16.7457 0.2800 0.3340 0.8000 -16.7747 0.3200 0.3788 0.9000 -16.4792 0.3600 0.4230 1.0000 -15.9238 0.4000 0.4665 1.1000 -15.1651 0.4400 0.5094 1.2000 -14.2524 0.4800 0.5515 1.3000 -13.2287 0.5200 0.5930 1.4000 -12.1307 0.5600 0.6338 1.5000 -10.9898 0.6000 0.6739 1.6000 -9.8326 0.6400 0.7134 1.7000 -8.6812 0.6800 0.7523 1.8000 -7.5539 0.7200 0.7904 1.9000 -6.4654 0.7600 0.8280 2.0000 -5.4272 0.8000 0.8649 2.1000 -4.4484 0.8400 0.9012 2.2000 -3.5353 0.8800 0.9368 2.3000 -2.6923 0.9200 0.9719 2.4000 -1.9220 0.9600 1.0063 2.5000 -1.2253 1.0000 1.0401 2.6000 -0.6019 1.0400 1.0732 2.7000 -0.0504 1.0800