实验四用matlab实现拉格朗日插值、分段线性插值.pdf

实验四 用 MATLAB 实现拉格朗日插值、分段线性插值 一、实验目的： 1）学会使用MATLAB 软件； 2）会使用MATLAB 软件进行拉格朗日插值算法和分段线性差值算法； 二、实验内容： 1 用 MATLAB 实现 y = 1./(x.^2+1); （-1=x=1）的拉格朗日插值、分段线性 2.选择以下函数，在 n 个节点上分别用分段线性和三次样条插值的方法，计算 m 个插值点的 函数值，通过数值和图形的输出，将插值结果与精确值进行比较，适当增加n ，再作比较， 由此作初步分析： (1).y=sinx;( 0≤x≤2 π) (2).y=(1-x^2) （-1≤x≤1） 三、实验方法与步骤： 问题一 用拉格朗日插值法 1）定义函数：y = 1./(x.^2+1);将其保存在 f.m 文件中，程序如下： function y = f1(x) y = 1./(x.^2+1); 2）定义拉格朗日插值函数：将其保存在 lagrange.m 文件中，具体实现程序编程如下： function y = lagrange(x0,y0,x) m = length(x); /区间长度/ n = length(x0); for i = 1:n l(i) = 1; end for i = 1:m for j = 1:n for k = 1:n if j == k continue; end l(j) = ( x(i) -x0(k))/( x0(j) - x0(k) )*l(j); end end end y = 0; for i = 1:n y = y0(i) * l(i) + y; end 3）建立测试程序，保存在 text.m 文件中，实现画图： x=-1:0.001:1; y = 1./(x.^2+1); p=polyfit(x,y,n); py=vpa(poly2sym(p),10) plot_x=-5:0.001:5; f1=polyval(p,plot_x); figure plot(x,y,‘r,plot_x,f1) 二 分段线性插值： 建立 div_linear.m 文件。具体编程如下 /*分段线性插值函数：div_linear.m 文件*/ function y = div_linear(x0,y0,x,n) %for j = 1:length(x) for i = 1:n-1 if (x = x0(i)) (x = x0(i+1)) y = (x - x0(i+1))/(x0(i) - x0(i+1))*y0(i) + ( x - x0(i))/(x0(i+1) - x0(i))*y0(i+1); else continue; end end %end 测试程序（text2.m）: x0 = linspace( -5,5,50); y0= 1./(x0.^2+1); y=interp1(x0,y0,x0,linear) plot(x0,y0,x0,y,p); 2）运行测试程序 问题二：(1).分段线性插值 Matlab 命令如下： x=linspace(0,2*pi,100); y=sin(x); x1=linspace(0,2*pi,5); y1=sin(x1); plot(x,y,x1,y1,x1,y1,o,LineWidth,1.5), gtext(n=4) 图形如下： (2).三次样条插值选取7个基点计算插值Matlab命令如下 x0=linspace(-1,1,7); y0=(1-x0.^2); x=linspace(-1,1,100); y=interp1(x0,y0,x,spline); x1=linspace(-1,1,100); y1=(1-x1.^2); plot(x1,y1,k,x0,y0,+,x,y,r); 图形如下: 5、实验总结： 通过本次课程设计，我初步掌握了 MATLAB 运用，加深了对于各种线性插值的理解；培养了 独立工作能力和创造力；综合运用专业及基础知识，解决实际数学问题的能力；在本次课程 设计中，在老师的精心指导下，收益匪浅。同时对数学的研究有了更深入的认识。