Отправлено ВН 02 апреля 2003 г. 16:57 В ответ на: Расскажите как делать обнаружение сигнала. отправлено Newka 02 апреля 2003 г. 15:37

Посчитайте БПФ. Возьмите квадрат модуля. Найдите среднее этого квадрата модуля от 0 до fd/2. Возьмите значение квадрата модуля на интерес. частоте и отнесите к среднему. Отношение сравните с порогом,
~20-25. Превышает - значит попался, не превышает - ловим дальше.
Если длина синуса известна, то БПФ по той же длине. Окошко перед БПФ можно наложить.
А корреляция сигнала с Вашей синусоидой это есть ДПФ в одной точке.
Считается она: Ac(k)=summa(SIG(i+k)*cos(2*pi*1KHz*i/fd)); As=то же самое, только вместо косинуса - синус.
Параллельно: E(k)=summa(SIG(i+k)^2).
Далее либо (Ac^2+As^2)/E, либо корень из этого. И сравнение с порогом.
Или как SM рекомендовал.

• Ресурсов многовато кушает — AntZ   (02.04.2003 17:33, 616 байт)
  • Ответ: — ВН   (02.04.2003 17:54, 370 байт)
    • Автокорреляция требуется только в 2х точках, а в принципе можно обойтись и одной. Если автокорреляция в точке T (время полного периода) ~равна энергии сигнала (а это значение автокорреляции в точеке 0). Итого требуется N умножений и N-1 сложений, где N длина буфера. Фича в том, что гарантируется что есть периодическая составляющая, но нет гарантии что это синус — AntZ   (02.04.2003 18:34, пустое)
      • Ответ: — ВН   (02.04.2003 19:30, 793 байт)
  • Про второй метод - зачем корреляция-то? Герцель это то-же самое но менее ресурсоемко. — SM   (02.04.2003 17:38, пустое)
    • Возможно, я знаю что Герцель хороший :-), но сам его писать не пробовал :-) — AntZ   (02.04.2003 18:35, пустое)
• А зачем нам деление на Е. Можно ли просто оценить (Ac^2+As^2) — Newka   (02.04.2003 17:31, 1 байт)