а я откуда знаю... я с этими методами даже не знаком

Так попробуйте познакомиться со статейкой из IEEE. Кому нужно Вам или нам?
http://electronix.ru/forum/index.php?showt...&start=1215

Там как раз интерполяция по Лагранжу
Статейке в Springerlink из китайского журнала кажется защищенная от копирования, если я правильно понял, скопировать бесплатно её не получится

Вот эта?

В каком смысле - "Потом удалю"?

Вот та статья VU Jiekang из IEEE более соответствующая контексту оценки частоты силовой сети.
IEEE transactions on power delivery ISSN 0885-8977, 2005, vol. 20, no1, pp. 366-374 [9 page(s) (article)] (22 ref.)
Там нет однако оценок спектра
http://rapidshare.com/files/315859486/Powe...qEstimation.pdf

Большое спасибо за статью.
Я хотел бы попросить у TigerSHARK какую-нибудь реализацию с отклонением по частоте не менее 0.1 Гц. Я попробую определить нужные параметры и скажу как.
Сэмплы менее 24бит не интересны.
Эффект "растекания спектра" не относится ни к Вашим данным (согласно Котельникову Вы имеете все что нужно) ни к тому, что Вы желаете получить. Это эффект БПФ и окна,
Следовательно применив для реальных сигналов более прямой способ выделения известных гармоник - умножение на синус, на косинус, сглаживание... - Мы получим идеализированный невозможный для БПФ результат. (для ослабления эффекта Гиббса потребуются данные о предыстории и послеистории сигнала, всего на длину сглаживающего фильтра)
В случае реальных сигналов результат БПФ менее точный из-за неизбежной амплитудной и еще хуже незначительной частотной модуляции и из-за несинхронной части шумов. И чем больше число выборок тем хуже, уменьшается только шум.
Наверное существует оптимальный объем выборки для каждого сигнала, при котором шум достаточно уменьшен, а модуляции не сильно затронули соседние частоты.
Ресэмплинг мне не понравился. Самый лучший это подстановка sinx/x для некоторых значений выборки около требуемого значения времени - приводило к слишком длительным вычислениям, а применение Фэрроу слишком обогащало спектр Моего сигнала. А сигнал был принципиально несинхронным. В целом знание пред и послеистории значительно повышает точность ресэмплинга, делайте запас сразу.
Самое простое - попробовать. А вдруг у Вас ни шумов ни модуляций ни посторонних сигналов, а процессор Атом, тогда все должно быть шоколадно.
А еще мне встретилась задача где точность 10% = очень хорошо, а 6% = идеально. Удачи.

Готов предоставить всё, что в моих силах.
Немного не понял вопроса. Что вы имеете ввиду под "реализацией"...

Честно признаюсь, что Я - дилетант, и вокруг никто курс ЦОС не прослушал. Жаргон такой образовался: реализация = набор сэмплов полученных одновременно с нескольких датчиков при одном испытании, например "вчерашняя реапизация" или "реализация где датчик отвалился". Если у Вас есть правильный термин постараюсь выучить.
Постараюсь помочь, хотя времени мало - сегодня начальник выбрал корпус для прибора, с понедельника работа под давлением до макетного экземпляра прибора. Математика и фреймы написаны.
Мне нужно:
1. Какой АЦП.
2. Реализация(файл содержащий результаты АЦП без промежутков и заголовков)

... у меня никакой аппаратной реализации((
АЦП по расчётам 12 бит минимум.
Моделирую всё в MATLAB. Подаю на вход алгоритма полигармонический синус содержащий 40 гармоник кратных 50 Гц (50, 100, ... до 2000) с одинаковыми амплитудами и нулевыми фазами (это наихудший теоретический случай, позволяющий оценить погрешность сверху) Ставлю кратную частоту дискретизации (минимум 6400Гц для БПФ). Получается эталонный спектр.
Затем меняю частоту основной гармоники (скажем 47.35 ставлю), не меняя частоты дискретизации. Смотрю на спектр - искажён.... вот тут, получается, и надо как-то частоту дискретизации поменять, и, соответсвенно, временное окно анализа.

В результате понимаю, что нужно как то взять с запасом частоту дискретизации, с учётом погрешности и порядка фильтра-интерполятора. Порядок расчёта длинны выборки БПФ у меня есть...

Так вот главный вопрос всё ещё актуален: фильтр фарроу меняет длину временного окна? как я понял он только выборки сдвигает во времени и всё((

Так он может сдвигать на разные величины сдвига. Поэтому можно получить любое количество промежуточных отсчетов в любых промежуточных точках в каких пожелаете. В частности, можно провести ресамплинг на любую частоту дискретизации c*Fd.

К-во отсчетов меняет. А продолжительность в секундах - нет. Это уже сами добавляете или отбрасываете.

И полифазные фильтры позволяют перейти на любую "рациональную" частоту (m/n)*Fd. Полифазные фильтры, как говорилось, практически не вносят искажений в АЧХ, в отличие от. Но Вам нужно было бы поднимать частоту дискретизации на очень высокий уровень, чтобы иметь возможность отслеживать основной тон с достаточной точностью.

Фильтр фарроу количество отсчетов не меняет. Он только сдвигает сигнал, а частота дискретизации остается та же. Меняет количество отсчетов ресэмплер, если он сделан на основе фильтра фарроу.

Фильтр фарроу выполняет в чистом виде интеполяцию лагранжа.

Вот есть отсчеты x(0) x(1) x(2) x(3) ...
Мы хотим увеличить частоту дискретизации в 1.1 раза
Получим отсчеты x(0) x(1.1) x(2.2) x(3.3) ...

Отсчет x(1.1) находится между x(1) и x(2)
Интерполятор 3-й степени вычисляет его так:
h(0)*x(0)+h(1)*x(1)+h(2)*x(2)+h(3)*x(3)

Для отсчета x(2.2)
h(0)*x(1)+h(1)*x(2)+h(2)*x(3)+h(3)*x(4)
итд

К-ты h(0...3) вычисляются для каждого нового отсчета исходя из дробной части его индекса.

Сам по себе фильтр фильтр фарроу это просто элемент дробной задержки, у которого есть параметр mu, который определяет, сколько эта задержка составляет относительно интвервала дискретизации входного сигнала. А на основе фильтра фарроу можно построить переменную задержку или ресэмплер. Чтобы сделатьресэмплер недостаточно просто фильтра, нужна еще хитрая схема синхронизации.

Это уже философия. Так Вы человека совсем запутаете.
Фильтр - не меняет. Постоянное mu.
А интерполятор - меняет. Хоть на каждом такте mu можно менять.
Самый обычный интерполятор B-сплайном.
Хоть для сдвига (фильтра), хоть ресамплинга.
Короче, если человек хочет поменять, то он может поменять с помощью интерполятора фарроу.
Трудно даже задачу представить, где mu не меняется. Скорее фильтр - это когда mu меняется медленно, а
интерполятор, когда mu меняется быстро

Хотя может если в FPGA то и хитрая синхронизация...

Прошу прощения Я в матлабе не работаю. Всегда на живом процессоре. Сейчас LPC2478.
Я думаю разработка или публикация по практической тематике должны производиться в следующем порядке:
1. При помощи своего или чужого прибора получить запись интересующего процесса.
2. На живом процессоре или в матлабе(по желанию) получить желаемое.
3. Проверить на других объектах.
4. Приступить к созданию прибора или публикации.
В Вашем случае п. 1 не выполнен.
А кто Вам сказал, что эти синусы там имеются и какой величины и что 12бит хватит?
У WU Jiekang 16бит. У меня 24 и меньше не получится.
А что касается Фарроу, то это способ пересчитать сэмплы в нужные моменты времени.(Ресэмплинг)
Просто берете выборку с запасом и пересчитываете.
Еще не ясно как реальный сигнал на это издевательство прореагирует.
Если выполните п.1 помогу с п.2.