matlab实现拉格朗日函数,拉格朗日插值多项式

Posted 邢逸

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了matlab实现拉格朗日函数,拉格朗日插值多项式相关的知识,希望对你有一定的参考价值。

%拉格朗日插值多项式  利用矩阵求解
x=1:0.2:3;%已知数据点x坐标向量:x
y=sin(x);%已知数据点x坐标向量:y
x1=1.1:0.2:3.1;%插值点的x坐标:x1
L=zeros(11,11);%另L矩阵为0        
        for i=1:11
            A=ones(10,1);%另A矩阵为10行1列的矩阵
            x2=x;
            x2(i)=[];
            x2‘;%10行一列 
            B=ones(1,11);%另B矩阵为1行11列的矩阵
            A*x1;%10行11列
            (x2‘)*B;%10行11列
             A*x1-(x2‘)*B;%11行11列
            ones(10,11);
            x(i);%提取x的第i个元素
            ones(10,11)*x(i);%10行11列的矩阵
            prod(A*x1-(x2‘)*B);%基函数的分子
            ones(10,11)*(x(i))-(x2‘)*B;%基函数的分母
            C=prod(A*x1-(x2‘)*B)./prod(ones(10,11)*(x(i))-(x2‘)*B);%对x2进行转置%C矩阵是一个1行11列的矩阵
            L(i,:)=C; %将A的第一行元素全部变为10 %将得到的矩阵赋值基函数的1,2,3。。。。11行
 
        end
    L;%11行11列
    y;%1行11列
    y1=y*L

  结果如下:

y1 =

  1 至 7 列

    0.8912    0.9636    0.9975    0.9917    0.9463    0.8632    0.7457

  8 至 11 列

    0.5985    0.4274    0.2392    0.0416

  

以上是关于matlab实现拉格朗日函数,拉格朗日插值多项式的主要内容,如果未能解决你的问题,请参考以下文章

12月学习进度4/31——计算机图形学拉格朗日插值的MATLAB实现

拉格朗日插值方法

拉格朗日插值Python代码实现

拉格朗日插值法

多项式函数插值:全域多项式插值单项式基插值拉格朗日插值牛顿插值 [MATLAB]

拉格朗日(Lagrange)插值算法