Обнаружение частоты, Нужно вычеслить наличие 1000Гц Опции

[Aleksandr Barano...]

А что, вычисление величины гармоники 1000 Гц - это не вычисление коэффициента корреляции сигнала с синусоидой частоты 1000 Гц?

[GetSmart]

А что, вычисление величины гармоники 1000 Гц - это не вычисление коэффициента корреляции сигнала с синусоидой частоты 1000 Гц?

Не пойму, почему все так считают. Корреляция - скалярная сумма, ДПФ - комплексная. Корреляция равна мощности сигнала, ДПФ равна вектору, содержащему амплитуду и фазу. Для определения уровня сигнала (синуса) методом корреляции нужно сделать множество свёрток с разными относительными сдвигами. Для ДПФ достаточно посчитать один раз и сразу выяснить относительный сдвиг между синусами и одновременно амплитуду. Кроме того, ДПФ обнаруживает только синусоиды, корреляция - вообще любые сигналы. И главная фича - корреляцию можно сделать сигнала самим с собой.

Сообщение отредактировал GetSmart - Nov 26 2008, 00:16

[Microwatt]

Да всё это звучит грозно, а на деле ... Гигобайты это на ПиСишках, а здесь их нет и сотнями байтов укладываемся.
Тут ещё одной информации не хватает - а именно превышение сигнала над фоном. Насколько сигнал чёткий будет?
Конечно фильтр перестраевыемый сделать значительно сложнее.

Так надо матчасть учить. Что умеет 16-разрядный процессор приблизительно понятно. И сколько времени программер будет уродоваться над всякими БПФ, пока 1000Гц от 2000Гц в идеальном спектре так-сяк программа отличать начнет.
А почему все-таки стандартный тональный декодер не годится? При уровне сигнала -6 дБ уверенно обнаруживает. Когда на экране самого сигнала за помехами не видно. За несколько периодов всего.
И перестраивается практически одним резистором. Это можно и программно делать.
Просто программеры не ищут легких путей в жизни?

[SasaVitebsk]

Так надо матчасть учить. Что умеет 16-разрядный процессор приблизительно понятно. И сколько времени программер будет уродоваться над всякими БПФ, пока 1000Гц от 2000Гц в идеальном спектре так-сяк программа отличать начнет.
А почему все-таки стандартный тональный декодер не годится? При уровне сигнала -6 дБ уверенно обнаруживает. Когда на экране самого сигнала за помехами не видно. За несколько периодов всего.
И перестраивается практически одним резистором. Это можно и программно делать.
Просто программеры не ищут легких путей в жизни?

Да нет, наоборот. Я просто говорю, что этот параметр надо учитывать при выборе алгоритма. А то наши советы бессмысленны.
Цифровой фильтр у меня также устойчиво обнаруживал сигнал DTMF на фоне разговора. При этом никогда не подрабатывал во время разговора. И, действительно, обнаруживал сигнал буквально по нескольким периодам. Мы добились качества работы неотличимого от ВЖ18. Точных параметров я вам не скажу, так как там всё несколько сложнее было (сигнал восстанавливался из цифры сначала).

У корелятора эти параметры раз в 10 хуже были.

[Alex_1811]

Вот до чего я докотился .
Пока на Delphi для простоты реализации.

Кто скажет, похоже на правду?

Прикрепленные файлы

 Project1.rar ( 165.34 килобайт ) Кол-во скачиваний: 26

 

[GetSmart]

Кто скажет, похоже на правду?

Очень похоже. Даже не отличишь
Только свёртку делают с нормированным (еденичным) синусом и косинусом, чтобы C (C:=sqrt(A*A+B*) содержало сразу же реальную амплитуду.

[Stanislav]

Возникла такая проблема:
Имеется Мега8 на вход АЦП которой через усилитель подключен микрофон.
Нужно детектировать наличее сигнала определенной частоты (в моем случаи 1000Гц).

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

ЗЫ. Для начала, можно попробовать полосовой фильтр конечной длины. Возьмите частоту выборки равной, скажем, 4 кГц, и реализуйте разностное отношение: Y(k)=Y(k-1)+X((k)-X(k-2)-X(k-40)+X(k-42). Это рекурсивный КИХ-фильтр, не содержащий операций умножения, и оптимальный для куска (ко)синусоиды частотой в 1 кГц и длиной в 10 периодов в белом шуме. Y(k) - текущий выходной отсчёт, Y(k-1) - предыдущий выходной отсчёт, X(k) - текущий входной отсчёт и т.д.
В реальности всё, конечно, сложнее.
В полосе фильтра может оказаться, например, форманта.

[Alex_1811]

Только свёртку делают с нормированным (еденичным) синусом и косинусом, чтобы C (C:=sqrt(A*A+B*) содержало сразу же реальную амплитуду.

Розжуйте пожалуйста белее популярно (доступно для меня).

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

Форму сигнала я вам не скажу так как сам не представляю.
Но получаем его приблизительно так:
1. Радиопередатчик передает синусоиду 1КГц в эфир.
2. Со стороны радиоприемника стоит мое устройство с подключ. микрофоном к АЦП Мега8.
3. Мне нужно определить есть ли на АЦП сигнал частотою 1КГц.

[МП41]

Это для синхронизации по 6-и точкам что-ли?

[GetSmart]

Розжуйте пожалуйста белее популярно (доступно для меня).

В конце обоих строк кода лишние умножения на 30.

A:=A+(Signal*sin(cnt*2*Pi*FREQ/512)*30);
B:=B+(Signal*cos(cnt*2*Pi*FREQ/512)*30);

[Alex_1811]

А от корня квадратного никак избавиться нельзя, для Меги8 напряг.

Сообщение отредактировал Alex_1811 - Nov 27 2008, 08:02

[GetSmart]

А от корня квадратного никак избавиться нельзя, для Меги8 напряг.

А зачем? Вам нужно задать порог срабатывания каким-нить числом. Задавайте порог сразу квадратом числа и сравнивайте. Выше - значит синусоида присутствует, ниже - отсутствует. Или нужно знать реальную амплитуду?

ЗЫ. Корень можно вычислить 3..5 итерациями методом Ньютона.

[Alex_1811]

Или нужно знать реальную амплитуду?

Для данного проэкта не нужно.

[Stanislav]

Форму сигнала я вам не скажу так как сам не представляю.
Но получаем его приблизительно так:
1. Радиопередатчик передает синусоиду 1КГц в эфир.
2. Со стороны радиоприемника стоит мое устройство с подключ. микрофоном к АЦП Мега8.

Непонятно что-то.
Сталыть, передатчик передаёт синусоиду в эфир.
А ловите его Вы микрофоном.
Если учесть, что плотность эфира =0, Ваша система работать не должна.

...3. Мне нужно определить есть ли на АЦП сигнал частотою 1КГц.

Я предлагаю Вам начать с простейших фильтров, на которых можно изучить элементарные свойства цифровой фильтрации, а затем уже двигаться дальше.