Franky написал(а):ак вот, возможно нужно делать так если интересует лагранж 1-й степени: берём по две точки из выборки в N точек. И между каждыми двумя точками находим полином первой степени. Причём последняя точка из двух становиться первой для второй пары точек и так далее...
Правильно
Franky написал(а):если рассматривается вторая степень, то берём уже по три точки из выборки и строим полином между каждыми тремя точками полином второй степени, причем каждая последняя из трёх становиться первой бля следующей тройки... и так далее....
А вот тут ошибка. Склейка если так можно выразится стык в стык только при 1-ой степени. При второй идет как бы перекрытие, т.е. берём уже по три точки из выборки и строим между каждыми тремя точками полином второй степени, но делаем следующим образом. Есть скажем 5 отсчетов x0,x1,x2,x3,x4 и x5. Первая парабола назовем ее P0 по первым трем точкам x0,x1,x2, а вторая парабола P1 по трем точкам x1 x2 x3, т.е. интервал x1 x2 будет как в первой параболе P0 так и во второй P1 и на интервале x1 x2 будет сумма парабол P0+P1. Другими словами сдвиг только на один отсчет. третья парабола соответсвенно x2 x3 x4 и четвертая x3 x4 x5.
Вам до какой степени необходимо проанализировать? Если степень не очень высокая, то лучше всего построить фильтры фарроу всех порядков от 1 до K и проанализировать их. Есть в матлабе еще функция intfilt по ней написано в справке: b = intfilt(l,n,'Lagrange') designs an FIR filter that performs nth-order Lagrange polynomial interpolation on a sequence interleaved with l-1 consecutive zeros every r samples. b has length (n + 1)*l for n even, and length (n + 1)*l-1 for n odd. If both n and l are even, the filter designed is not linear phase.
. Ну денек повозится придется конечно, но мне кажется что лучше день потерять и потом за 5 минут долететь.