Delphi版本插值Lagrange

最新推荐文章于 2023-03-29 16:41:52 发布

转载最新推荐文章于 2023-03-29 16:41:52 发布 · 289 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

原文链接：http://www.cnblogs.com/cinemaparadiso/p/11586255.html

收录于

本文介绍了一种在数值分析中常用的插值方法——拉格朗日插值法的实现过程。通过一个具体的算法函数，展示了如何利用输入的离散数据点，计算出任意点处的插值结果。

function Lagrange(x, y:ArrayDouble; N: integer; x1: Double): double;
var
I, J: Integer;
a: Double;
begin
Result := 0;
for i := 0 to N - 1 do
begin
a := 1;
for J := 0 to N - 1 do
begin
if (I <> J) then
a := a * (x1 - x[J]) / (x[I] - x[J]);
end;
Result := Result + a * y[I];
end;
end;