哈工大计算方法实验拉格朗日.docx

1、哈工大计算方法实验拉格朗日实验题目1 Lagrange插值摘要给定平面上n+1个不同的数据点:则满足条件的n次拉格朗日插值多项式是存在唯一的。若，且充分光滑，则当时，有误差估计式前言利用拉格朗日插值多项式求的近似值程序设计流程问题1 （1）N = 5时，程序运行如下：TestLag(inline(1./(1+x.2), -5, 5, 5, :;将区间-5,5分为了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 10时，程序运行如下：TestLag(inline(1./(1+x.2), -5, 5, 10, :;将区间-5,5分为

2、了10段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 20时，程序运行如下：TestLag(inline(1./(1+x.2), -5, 5, 20, :;将区间-5,5分为了20段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = 问题1 （2）N = 5时，程序运行如下：TestLag(inline(exp(x), -1, 1, 5, );将区间-1,1分为了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N = 1

3、0时，程序运行如下：TestLag(inline(exp(x), -1, 1, 10, );将区间-1,1分为了10段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N = 20时，程序运行如下：TestLag(inline(exp(x), -1, 1, 20, );将区间-1,1分为了20段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * 0 问题2 （1）N = 5时，程序运行如下：TestLag(inline(1./(1+x.2), -1, 1, 5, );将区间-1,1分为

4、了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 10时，程序运行如下：TestLag(inline(1./(1+x.2), -1, 1, 10, );将区间-1,1分为了10段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 20时，程序运行如下：TestLag(inline(1./(1+x.2), -1, 1, 20, );将区间-1,1分为了20段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * 实验2 （2

5、）N = 5时，程序运行如下：TestLag(inline(exp(x), -5, 5, 5, );将区间-5,5分为了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 10时，程序运行如下：TestLag(inline(exp(x), -5, 5, 10, );将区间-5,5分为了10段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = N = 20时，程序运行如下：TestLag(inline(exp(x), -5, 5, 20, );将区间-5,5分为了20段计算插值的点 x

6、i = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * 问题3 （1）N =5时，程序运行如下：TestLag2(inline(1./(1+x.2), -1, 1, 5, );将区间-1,1分为了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N =10时，程序运行如下：TestLag2(inline(1./(1+x.2), -1, 1, 10, );将区间-1,1分为了10段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N =20时，程序

7、运行如下：TestLag2(inline(1./(1+x.2), -1, 1, 20, );将区间-1,1分为了20段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * 问题3 （2）N =5时，程序运行如下：TestLag2(inline(exp(x), -1, 1, 5, );将区间-1,1分为了5段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N =10时，程序运行如下：TestLag2(inline(exp(x), -1, 1, 10, );将区间-1,1分为了10段计算插

8、值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * N =20时，程序运行如下：TestLag2(inline(exp(x), -1, 1, 20, );将区间-1,1分为了20段计算插值的点 xi = 计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = * 问题4(1) 程序运行如下：TestLag3(1 4 9, 5 50 115 185)计算插值的点 xi = 5 50 115 185计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = (2) 程序运行如下：TestLag3(36 49 6

9、4, 5 50 115 185)计算插值的点 xi = 5 50 115 185计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = (3) 程序运行如下：TestLag3(100 121 144, 5 50 115 185)计算插值的点 xi = 5 50 115 185计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = (4) 程序运行如下：TestLag3(169 196 225, 5 50 115 185)计算插值的点 xi = 5 50 115 185计算出的插值 yi = 插值点处函数值 yFact = 计算误差 err = 实验所

10、用函数function yh = LagInterp(x, y, xh)% LagInterp 计算拉格朗日插值% Synopsis: yh = LagInterp(x, y, xh)% Input: x = 一维向量，将要做插值x的值% y = 一维向量，将要做插值y的值% xh = 数值或一维向量，计算插值的位置，支持计算一列xh的值% Output: yh = 数值或一维向量，通过计算插值的位置算出的插值 if min(size(x) 1 | min(size(y) 1 %判断x，y是否为一维向量 error(x,y must be vectors!); elseif length(x)

11、 = length(y) %判断x，y是否有同样多的元素 error(x and y must agree!); end yh = zeros(size(xh); L = zeros(length(x) - 1); for j = 1:length(xh) for i = 1:length(x) xCal = x; xCal(i) = ;%prod(xh(j) - xCal)/prod(x(i) - xCal)为拉格朗日基函数 L(i) = prod(xh(j) - xCal)/prod(x(i) - xCal); yh(j) = yh(j) + L(i) * y(i); %yh = sum(

12、L(i) * y(i) endendfunction TestLag(fx, a, b, n, xi)% TestLag 实验题目1 1,2% Synopsis: TestLag(fun, a, b, n, xi)% Input: fx = 用来验证插值计算准确率的函数% a,b = 节点选取上下限% n = 多项式次数，固定区间-a,b分段数% xi = 要计算插值的点 x = linspace(a, b, n); y = feval(fx, x); yi = LagInterp(x, y, xi); yFact = feval(fx, xi); err = yFact - yi; fpri

13、ntf(将区间%d,%d分为了%d段n, a, b, n); fprintf(计算插值的点 xi =n); disp(xi); fprintf(计算出的插值 yi =n); disp(yi); fprintf(插值点处函数值 yFact =n); disp(yFact); fprintf(计算误差 err =n); disp(err);function TestLag2(fx, a, b, n, xi)% TestLag2 实验题目1 3% Synopsis: TestLag2(fun, a, b, n, xi)% Input: fx = 用来验证插值计算准确率的函数% a,b = 节点选取上

14、下限% n = 多项式次数，固定区间-a,b分段数% xi = 要计算插值的点 x = zeros(1,n); for k = 1:n x(k) = cos( (2*k-1)*pi/(2*n) ); %构造非等距节点 end y = feval(fx, x); yi = LagInterp(x, y, xi); yFact = feval(fx, xi); err = yFact - yi; fprintf(将区间%d,%d分为了%d段n, a, b, n); fprintf(计算插值的点 xi =n); disp(xi); fprintf(计算出的插值 yi =n); disp(yi); f

15、printf(插值点处函数值 yFact =n); disp(yFact); fprintf(计算误差 err =n); disp(err);function TestLag3(x, xi)% TestLag3 实验题目1 4% Synopsis: TestLag3(fun, a, b, n, xi)% Input: x = 构造Lagrange插值的节点% xi = 要计算插值的点 fx = inline(sqrt(x); y = feval(fx, x); yi = LagInterp(x, y, xi); yFact = feval(fx, xi); err = yFact - yi;

16、fprintf(计算插值的点 xi =n); disp(xi); fprintf(计算出的插值 yi =n); disp(yi); fprintf(插值点处函数值 yFact =n); disp(yFact); fprintf(计算误差 err =n); disp(err);思考题1. 拉格朗日插值多项式的次数并不是越大越好，根据定义，插值式可以在节点处与实际函数匹配，但不能保证在节点之间很好的逼近实际函数。这个现象就是多项式摆动Runge现象，有时多项式摆动可以通过谨慎选择基础函数的取样点来减小。通常采用分段插值可以很好的消除多项式摆动现象。2. 在分段段数相同的情况下，插值区间越大，误差越大。原因是大部分情况下，相对于比较大的区间，函数在比较小的区间上的函数值变化较缓和，因此即使出现摆动也不会偏离原函数太大。3. 第一问中已经提到多项式摆动可以通过谨慎选择基础函数的取样点来减小。我们来看下的分布和的图像：可以看出，的分布是两端密集中间稀疏，的趋势是两边陡峭中间平缓，函数变化陡峭时节点增多正好可以增加插值的准确性。4. 一般来说，内插时插值收敛于实际函数，一旦超出内插的范围，插值函数会发散，且离插值区间越远外推误差越大。使用不用的插值方法在同一点外推的值也会相差很多，这说明外推本身就存在很大的不确定性。