А другие физические факторы. В жизни не бывает чистого синуса идеально синхронизированного с оцифровкой сигнала. Всегда есть небольшая амплитудная и частотная модуляция например из-за рыскания системы синхронизации. И мне очень хотелось бы иметь АЦП для оцифровки 62Гц имеющий хотя бы 22 честных бита при больших сигналах. Причем погрешность младших разрядов имеющихся в продаже не чистый белый шум, а величина сложным образом связанная с биениями сигнала и частоты оцифровки. А всякие расчеты основанные на том,что наши АЦП завтра будут выдавать результат в формате double лучше оставим потомкам и в запечатанной капсуле чтобы не видели нашего идиотизма. Я уже давно не моделирую ничего в Матлабе. При некотором опыте все идеи легче проверить на реальном процессоре, чем учесть все факторы при математическом моделировании.
А что касается поставленной задачи в жизни правильнее применить адаптивный метод постоянно подстраивающийся к мешающему фактору. А формально при идеальной синхронизации окна и идеальной синусоиде мешающего фактора метод предложенный автором несомненно наилучший, можно даже рассчитать эту гармонику один раз и вечно ее вычитать.

Устранять неопределенность можно увеличением общей длины выборки. Ничего не стоит записать входной и выходной сигналы в канале связи и проанализировать их многократно и неторопливо. Можно узнать насколько портится эта самая гармоника 62 Гц и сделать вывод как ее устранить. Наверняка она получит модуляцию (об этом писали), а значит и размытый спектр. То есть, в чистом виде ее (гармоники) на выходе канала связи уже нет. Если часть этой модуляции не случайная, ее можно компенсировать. Но следы все равно останутся. А режекторный фильтр - достаточно простое и устойчивое решение (без обратных связей).

Небольшое замечание. Если есть какая-то модуляция, то увеличение общей длины выборки приводит к размытию спектра. Чем больше бин занимает мешающий фактор тем хуже с ним бороться. Лучше уж адаптивно следить без всяких спектров.
Но вопрос был формальным: сам задал длину - довольно большую и сам задал помеху - абсолютно синхронную и не ограничил в ресурсах. Вопрос в правильной методике.

Здесь попросишь бумажку, а дадут томик война и мир, да ещё и вернуть скажут в целости и сохранности. А ты сидишь, читаешь, но в конце всё равно пальцем


Сообщение отредактировал ivan219 - Sep 20 2010, 18:24

Ну вот дошли руки и до эксперемента.
Сигнал синтезируется на ПК полезный сигнал и помеха. После чего они складываются и отправлястся на выход звуковой карты LinOut c LinOut сигнал поступает на LinIn и оцифровывается дальше делается ДПФ на частоту помехе она известна определяется амплитуда и фаза помехи синтезизуется сигнал помехи и из входного сигнала вычитается синтезируемый сигнал помехи (при вычитании среднего значения из сигнала результатов не заметил возможно они минимальны) после чего проделовалось FFT 1024 и получали чистый сигнал без помехи.
В приложении программа если есть навыки работы с паяльником можно посмотреть.

Сообщение отредактировал ivan219 - Sep 27 2010, 15:47