原创力文档- 35
- 0
- 约1.15万字
- 约 30页
- 2016-05-01 发布于湖北
- 举报
拉格朗日插值法
拉格朗日插值法基本原理：
通过平面上不同两点可以确定一条直线，这就是拉格朗日线性插值问题，对于不在同一条直线的三个点得到的插值多项式则为抛物线。
拉格朗日插值的基多项式（即基函数）为：
有了基函数以后就可以直接构造如下多项式：
该多项式就是拉格朗日插值法所求得的插值多项式。
拉格朗日插值法算法：
1、根据所给点的坐标依次写出其差值基函数（用循环可以轻易解决）
2、将差值基函数与其对应的点的函数值相乘得：
3、将2中各项累加即得差值多项式：
拉格朗日插值法程序：
function lagrange(A)%A为一个只有两行的矩阵，第一行为插值点，第二行为插值点对应的函数值
[m,n]=size(A);
f=1;
p=0;%两个用到的变量
syms x
for i=1:n
f=(x-A(1,i))*f;
end
for j=1:n
g(j)=f/(x-A(1,j));%求插值基函数的分母
h(j)=subs(g(j),x,A(1,j));%求插值基函数的分子
s(j)=g(j)/h(j)*A(2,j);%插值基函数
end
for k=1:n
s(j)=collect(s(j));%合并同类项
end
for i=1:n
p=p+s(i);
end
fprintf(拉格朗日插值法可得多项式:)
collect(p) % 可用lagr
文档评论（0）