измерение частоты аналогового сигнала, mega32 Опции

[Spym]

Доброго времени суток.

Требуется измерять частоту аналогового сигнала при помощи АЦП в mega32, точнее - искать сигнал с частотой 4850 - 4990Гц.

Какие существуют методы обнаружения известной частоты?
Чем определяется необходимая частота дискретизации входящего сигнала?

Мне это представляется как отсчет времени между всплесками на заданную амплитуду, и поиском всплесков с заданным периодом, с дальнейшей статистической обработкой полученных данных.
С ЦОС практически не знаком...

[Dopler]

Как быстро необходимо дать ответ об обнаруженной частоте?
Можно применить синхронный детектор.
В памяти контроллера нужно хранить два сигнала - sin и cos требуемой частоты (две квадратуры) с такой же частотой дискретизации, какая будет у входного сигнала. Опорные сигналы должны содержать целое число периодов. Входной сигнал перемножается на опорные. Если разложить по формуле sin(a)*sin(, то получится сумма синусов с частотами a-b и a+b, т.е. если частота a=b, то получается постоянная составляющая, плюс сигнал удвоенной частоты. Таким образом, если пропустить полученный сигнал (входной, умноженный на опорный) через ФНЧ, то сигнал двойной частоты отфильтруется, останется только постоянная составляющая. Ее амплитуда определяется как корень квадратный из суммы квадратов синусной и косинусной составляющей.

В данном случае можно сделать ФНЧ на 100 Гц (а можно и на 20, в зависимости от скорости, с которой надо давать ответ), корень квадратный можно попробовать вычислять таблицей, хотя если не нужны точные значения амплитуды сигнала искомой частоты, а нужно только детектировать ее наличие - то корень брать вообще не обязательно.

[_Pasha]

Как быстро необходимо дать ответ об обнаруженной частоте?
Можно применить синхронный детектор.
************************************************************
Ее амплитуда определяется как корень квадратный из суммы квадратов синусной и косинусной составляющей.

Держу в руках старинную книжку
"Микропроцессорный комплект БИС серии К1815 для цировой обработки сигналов" М.:"Радио и связь",1992
Я ее не выбросил из-за одной цитаты (с.232)

Yвых = max (|Y_Re|, 7/8 |Y_Re| + 1/2 |Y_Im|, 1/2|Y_Re|+7/8|Y_Im|, |Y_Im|)

Это касается приближенного вычисления амплитуды огибающей на выходе фильтра.
Заявленная погрешность 2.98%.
В формуле Y_Re, Y_Im - накопленные суммы в квадратурных каналах,
|x| - это я обозначил модуль значения (по-маткадовски)
Очень благородно! Пользовал лично, еще когда об аппаратном умножителе на аврках,пиках и мечтать не приходилось.

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