Блог

Я создал простой код для интерполяции Лагранжа векторов данных x и y и , но мне нужна помощь в построении графиков и отображении ошибок.

Я хочу построить равноудаленные точки вместе с интерполированной кривой Лагранжа на одном и том же графике и отобразить ошибку аппроксимации. Построение точек тривиально, но я борюсь с построением интерполированной полиномиальной кривой поверх графика (x,y).

 clear
clc

x = -1:0.25:1;
y = 1./(1 25.*x.^2);

N = 1:length(y);
M = zeros(length(y),length(x));
for n = 1:length(x)
    for i = 1:length(y)
        l = 1;
        index = find(N ~= i);
        for jj = 1:length(index)
            l = l.*(x(n)-y(index(jj)))./(y(i)-y(index(jj)));
        end
        M(i,n) = l;
    end
end

plot(x,y)
 

Мы ценим любую помощь