[Alexey_N 0]

Господа, я пытаюсь сделать обнаружитель гармонического сигнала в его смеси с шумом. По-моему это - классическая постановка, и все давно это уже умеют.

А вот у меня не получается :(.

Беру входную смесь с шумом (массив чисел) и умножаю её на синус и на косинус той же частоты, что и сигнал (тоже массивы чисел). Полученные два массива почленно возвожу в квадрат, суммирую и беру корень из суммы. Полученный результат накапливаю, ну типа сложением...

Ожидаю, что в процессе обработки массивов на выходе значение будет нарастать. Так и есть, нарастает, но при этом не зависит от частоты входного сигнала. Ну то есть совпадает частота коррелятора со входным сигналом или нет, на выходе нарастающая пила одинакова.

Вопрос - я что-то не так делаю?

Или просто дикая ошибка где-то?

Спасибо.

[petrov 6]

После умножения на комплексную синусоиду накопить и только потом модуль вычислять.

 

[Herz 4]

А просто FFT нельзя сделать?

[Alexey_N 0]

После умножения на комплексную синусоиду накопить и только потом модуль вычислять.

Я делаю:

Вход = Сигнал + Шум

(Вход*Синус)^2 + (Вход*Косинус)^2 = Сумма

Результат = Результат + (Сумма)^0.5

 

В каком месте нужен модуль?

 

А просто FFT нельзя сделать?

FFT работает, без проблем.

Но мне нужен именно коррелятор! Вот тот самый, - из старинных советских учебников, именно он почему-то не выходит...

 

[Serg76 0]

В каком месте нужен модуль?

после свертки Вашего Входа с комплексной экспонентой

[Onkel 1]

Господа, я пытаюсь сделать обнаружитель гармонического сигнала в его смеси с шумом. По-моему это - классическая постановка, и все давно это уже умеют.

А вот у меня не получается :(.

Беру входную смесь с шумом (массив чисел) и умножаю её на синус и на косинус той же частоты, что и сигнал (тоже массивы чисел). Полученные два массива почленно возвожу в квадрат, суммирую и беру корень из суммы.

Спасибо.

 

синус квадрат плюс косинус квадрат будет 1, если я правильно понял вас, то и суммы почленные будут 1*входной сигнал.

Я для проверки самописанных всяких функци использую фурье- анализ экзела, там есть такое, кстати операции все комплексные.

Если у вас не просто таблица, а рельный эксперимент - не проще применить синхродетектирование?

 

 

(Вход*Синус)^2 + (Вход*Косинус)^2 = Сумма

сумма будет просто вход ^2

 

[Alexey_N 0]

сумма будет просто вход ^2

Да, согласен, в суете накосячил, нужно сначала интегрировать. Исправил, - всё заработало.

СПАСИБО!

:)