program for lagranges interpolation method using single for loop

3 次查看(过去 30 天)
how to program for lagranges interpolation method using single for loop
  4 个评论
Geoff Hayes
Geoff Hayes 2018-8-23
Please rename your local variable sum since it conflicts with the MATLAB built-in function of the same name.
Out of curiosity, why do you want to eliminate the inner for loop? For performance reasons?
PJS KUMAR
PJS KUMAR 2018-8-23
yes, for performance reasons can we write the program by using vector operations.

请先登录,再进行评论。

回答(1 个)

Geoff Hayes
Geoff Hayes 2018-8-23
PJS - if we assume that x is a vector/array then u is simply
u = (a - x(1)) * (a - x(2)) * ... * (a - x(i-1)) * (a - x(i+1)) * ... * (a - x(n));
where n is the number of elements in x. If this is true, then couldn't we do something like
for i = 1:length(x)
u = 1;
l = 1;
xt = x;
xt(i) = []; % remove the ith element
u = prod(a - xt);
% etc.
end
We use prod to calculate the product of all elements in the array. Similarly, for l
l = (x(i) - x(1)) * (x(i) - x(2)) * ... * (x(i) - x(i-1)) * (x(i) - x(i+1)) * ... * (ax(i) - x(n));
which means that the above can be replaced with
l = prod(x(i) - xt));
Try the above and see what happens!

类别

Help CenterFile Exchange 中查找有关 Interpolation 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by