Lagrange多项式(附代码)
拉格朗日插值在解决数据的拟合过程中比较常用,以下主要分享关于拉格朗日插值方法的matlab编程,具体如下:
Lagrange函数编写
function p = Lagrange(X,Y)
n = length(X);
syms x;
p = 0;
for k = 1:n
q=1;
for j = 1:n
if j~= k
q = q*(x-X(j))/(X(k)-X(j));
end
end
p = p+q*Y(k);
end
输入:X = [0.5 0.7 0.9 1.1 1.3 1.5 1.7 1.9];
Y = [0.4794 0.6442 0.7833 0.8912 0.9636 0.9975 0.9917 0.9463];
p = Lagrange(X,Y)
运行结果如下:
(83125*(x - 1/2)(x - 7/10)(x - 9/10)(x - 11/10)(x - 13/10)(x - 17/10)(x - 19/10))/256 + (16105*(5x - 5/2)(x - 3/2)(x - 9/10)(x - 11/10)(x - 13/10)(x - 17/10)(x - 19/10))/1152 - (3995(5x - 7/2)(x - 3/2)(x - 9/10)(x - 11/10)(x - 13/10)(x - 17/10)(x - 19/10))/2688 - (13055((5x)/2 - 5/4)(x - 3/2)(x - 7/10)(x - 11/10)(x - 13/10)(x - 17/10)(x - 19/10))/128 + (13925((5x)/3 - 5/6)(x - 3/2)(x - 7/10)(x - 9/10)(x - 13/10)(x - 17/10)(x - 19/10))/48 - (20075((5x)/4 - 5/8)(x - 3/2)(x - 7/10)(x - 9/10)(x - 11/10)(x - 17/10)(x - 19/10))/48 - (49585((5x)/6 - 5/12)(x - 3/2)(x - 7/10)(x - 9/10)(x - 11/10)(x - 13/10)(x - 19/10))/384 + (47315((5x)/7 - 5/14)(x - 3/2)(x - 7/10)(x - 9/10)(x - 11/10)(x - 13/10)*(x - 17/10))/2304
感受
第一次发表博客,希望大家多多支持,加油,加油!