Интерполяция полиномом Лагранжа и кусочная интерполяция

@MaryCantante · Регистрация: 14.11.2016

Author24 — интернет-сервис помощи студентам

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

Matlab M

    x=[0,0.2,0.4,0.6,0.8];
    y=2*cos(x);
function [p]=inter(a,x,y)
n=length(x);
p=0;
for i=1:n
    l=1;
    for k=1:n
        if k<>i then
            l=l .*((a-x(k)) ./(x(i)-x(k)));
        end;
     end;
     p=p+y(i) .*l;
end;
endfunction;
a=linspace(1,13,101);
plot(a,inter(a,x,y),x,y,'x')

вроде работает, но график странный выдает, график строит похожим больше на экспоненту, не могу понять, с чем связана эта ошибка...
+ нужно еще написать кусочную интерполяцию, где каждые 3 точки будут соединятся полиномом Лагранжа, а как эти точки выделять из общего числа понять не могу (примеры смотрела, находила только те, где всего по трем точкам интерполируют)
ХЕЛП!

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

Комментарий модератора

Правила форума: 5.5. Запрещено размещать тему в нескольких подразделах одного раздела одновременно (кросспостинг), а также дублировать тему в одном разделе.

@olbond · 17.11.2016, 03:59

Сообщение от MaryCantante

вроде работает, но график странный выдает

Так у Вас не интерполяция а экстрополяция - точки заданы на интервале [0;0.8] а Вы пытаетесь "интерполировать" на интервале [1;13]?

задайте точки для интерполяции правильно:

Matlab M

1	a=linspace(0,1.3,101);

Добавлено через 9 минут

Сообщение от MaryCantante

Не знаю, как брать по 3 точки, когда их гораздо больше

Нужно в цикле индексы увеличивать на 3, например:

Matlab M

for j = [1:3:103]
x_1 = x(j:j+2);
y_2 = y(j:j+2);
.........
end

@JRmsk · 11.04.2018, 00:35

Откуда взяли 1, 13, 101 ?

@olbond · 11.04.2018, 02:14

Сообщение от JRmsk

Откуда взяли 1, 13, 101 ?

Почитайте справку по linspace

@JRmsk 0 / 0 / 1 Регистрация: 17.03.2014 Сообщений: 10
	11.04.2018, 00:35	3
	Откуда взяли 1, 13, 101 ? 0

@olbond 143 / 133 / 42 Регистрация: 11.03.2013 Сообщений: 274
	11.04.2018, 02:14	4
	Сообщение от JRmsk Откуда взяли 1, 13, 101 ? Почитайте справку по linspace 0

Опции темы