数值分析原理课程实验——拉格朗日（Lagrange）插值

凉丶梦

已于 2022-04-26 21:57:46 修改

阅读量5.4k

点赞数 6

分类专栏：数值分析原理文章标签：算法 matlab c语言线性代数抽象代数

于 2021-04-29 22:47:34 首次发布

本文链接：https://blog.csdn.net/KissMoon_/article/details/116278449

版权

数值分析原理专栏收录该内容

4 篇文章 4 订阅

订阅专栏

本文介绍了拉格朗日插值法的原理和程序实现，探讨了插值次数、插值区间选择的影响，并通过Matlab和C语言代码展示了插值过程。同时，提到了牛顿、高斯列主元消去法和四阶龙格-库塔方法，讨论了内插与外推的可靠性。

摘要由CSDN通过智能技术生成

拉格朗日（Lagrange）插值

方法概要

在这里插入图片描述

待求问题

1.拉格朗日插值多项式的次数n越大越好吗？

在这里插入图片描述

2.插值区间越小越好吗？

在这里插入图片描述

3.内插比外推更可靠吗？

在这里插入图片描述

程序流程

在这里插入图片描述

程序代码

/*Matlab函数
function Result = Lagrange(n, x_in, a, b, f)
    h = (b-a)/n;
    M = zeros(n+1, 2);
    i = 0;
    while i <= n
        xi = a+i*h;
        yi = subs(f, symvar(f), xi);
        M(i+1, 1) = xi;
        M(i+1, 2) = yi;
        i = i+1;
    end
    y = 0.0;
    k = 1;
    while k <= n+1
       l = 1.0;
       for j = 1:k-1
           l = l*(x_in - M(j, 1))/(M(k, 1) - M(j, 1));
       end
       for j = k+1:n+1
           l = l*(x_in - M(j, 1))/(M(k, 1) - M(j, 1));
       end
       y = y + l*M(k, 2);
       k = k+1;
    end
    Result = [x_in, y];
end*/

/*C语言程序
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define N1 3   // n amount
#define N2 4   // x amount
#define N3 20  // n max

int Ns[N1] = {5, 10, 20};
double x[N2] = {-0.95, -0.05, 0.05, 0.95};
double l = -1.0;
double r = 1.0;

double X(int k, int n) {
    double h = (r - l) / n;
    return l + k * h;
}

double Y(double x) { return 1 / (1 + x * x); }

int main() {
    for (int i = 0; i < N2; i++) printf("\tx=%.2lf", x[i]);
    printf("\n");
    for (int i = 0; i < N1; i++) {
        double a[N3 + 1], b[N3 + 1];
        int n = Ns[i];
        for (int k = 0; k <= n; k++) {
            a[k] = X(k, n);  // x
            b[k] = Y(a[k]);  // y
        }
        printf("n=%d", n);
        for (int p = 0; p < N2; p++) {
            double y = 0.0;
            for (int k = 0; k <= n; k++) {
                double l = 1.0;
                for (int j = 0; j <= n; j++) {
                    if (j != k) l *= (x[p] - a[j]) / (a[k] - a[j]);
                }
                y += l * b[k];
            }
            printf("\t%.6lf", y);
        }
        printf("\n");
    }
    printf("Actual");
    for (int p = 0; p < N2; p++) printf("\t%.6lf", Y(x[p]));
    return 0;
}*/