Участок интерполяции - PullRequest
Новая
       25

Участок интерполяции

0 голосов
/ 01 апреля 2020

Проблема: у меня есть работающий скрипт интерполяции Лагранжа 

x0 = [-3 -2 -1 0 1 2 3];
y0 = [-15 1 5 3 1 5 21]; 
x=-1.65;

function y = lagrange(x, x0, y0)

n = length(x0); 
y = 0;

  for i=1:n
      p = 1;
      for j=1:n
          if j == i
              continue;
          endif;
          p *= (x-x0(j)) / (x0(i)-x0(j));
      endfor;
      y += y0(i) * p;   
  endfor;

endfunction

y = lagrange(x, x0, y0);

, но мне нужно нарисовать график этих значений на гладкой кривой (вероятно, в linspace). Предварительное условие: я не могу использовать встроенную функцию «полифит». Нужно ли указывать вектор со 100 значениями для построения аналогичной диаграммы?

...