Предельный по точности способ оценки частоты, фазы и амплитуды синусоиды (пусть даже и с гармониками) в присутствии белого шума - это квадратичная интерполяция спектра
https://ccrma.stanford.edu/STANM/stanms/sta...14/stanm114.pdfс продолжениями до 118
https://ccrma.stanford.edu/STANM/stanm/node3.htmlПредельный в смысле минимальности числа используемых отсчетов при заданной точности, поскольку он выходит на предел Крамера-Рао максимального правдоподобия. Ненужную полосу и ненужные гармоники предварительно нужно просто отфильтровать.
Ryfe&Boorstyn когда-то доказали, что метод выходит на предельные в статистическом смысле точности оценки
параметров для одиночной синусоиды, причем с использованием прямоугольного (никакого) окна. Они первыми и вывели предельную оценку Крамера-Рао для этой задачи. Потом они же показали, что для многих хорошо разрешимых (в смысле критерия Релея, мы говорим ещё - "принципа неопределённости" 1/T) спектральных линий вариация метода с хорошей функцией окна близка к оптимальной. Функция окна подбирается здесь, чтобы не просто разрешить,а изолировать спектральные линии.
Импульсный шум ничего не изменит. Цветной - может. Метод становится неоптимальным только когда шум сильно неоднороден по спектру (гармоническая помеха) или когда шум стремится к нулю (детерминированый сигнал).
Могут быть и другие методы. В соседней теме обсуждают возможность ресамплинга в исходной области на частоту кратную реальной измеренной, для точного измерения амплитуды основной частоты и гармоник посредством ДПФ.
http://electronix.ru/forum/index.php?showt...45&start=45Понятно, что если взять частоту дискретизации достаточно высокой (например 3*частота наивысшей гармоники),
отфильтровав предварительно все более высокие гармоники, можно провести ресамплинг на нужную частоту
интерполятором(например фарроу). В принципе, если интерполятор фарроу не обеспечит нужную точность (ему нужно будет минимум примерно 30*частоту наивысшей гармоники), всегда можно поднять частоту дискретизации хорошим полифазным фильтром на фиксированую более высокую частоту, после чего перейти на нужную (уже переменную) частоту интерполятором фарроу. Потом ДПФ (БПФ) Для любой точности задача чисто практическая