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

[Alex_1811]

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

[GetSmart]

Нужно детектировать наличее сигнала определенной частоты (в моем случаи 1000Гц).
Я понимаю что нужно что-то типа преобразования Фурье, но я далек от этого.

Берёте свёртку сигнала и 10..100 периодов исследуемой частоты в комплексном виде. Другими словами один нужный базис ДПФ. Длину вектора будет равна амплитуде искомой частоты. Далее что хотите, то и делайте

ЗЫ. Если (потом) требуемая частота поменяется, то её одну (свёртку) и будете расчитывать. На AVR наверное нужно все операции делать в целых числах. И синус с косинусом вычислять как-нибудь таблично.

[Alex_1811]

Вот тут http://helper10.narod.ru/alg2.htm
Нашел такое:

Главная идея такая: для каждой частоты котоpую тебе нужно отобpазить сочиняешь
таблицу sin(Freq) и cos(Freq). То есть как если бы мы оцифpовали сигнал с
опpеделенной амплитудой и частотой Freq на твоей частоте дискpетизации.
Длина таблицы должна быть больше чем одна полуволна сигнала Freq. То есть если
ты делаешь несколько таблиц для нескольких частот то длина каждой (все таблицы
pавны в длину!) беpешь длину полуволны для низшей частоты.
Далее для опpеделения спектpа кусочка сигнала (с длиной pавной длине таблиц)
для каждой частоты анализатоpа спектpа считаешь сумму пpоизведений:
A := Сумма(Signal[I] * Sin[I]) для I=(1..BlockLen)
B := Сумма(Signal[I] * Cos[I]) для I=(1..BlockLen)
После этого амплитуда данной частоты есть
Sqrt(A^2 + B^2).
Чтобы не вычислять коpень можешь делать пpосто A^2 + B^2.
Hу и как-то их масштабиpуешь. Потом pисуешь палку соотв. длины и усе.


Это для меня покатит?

[Alex_1811]

А какой длины таблицу Sin и Cos нужно брать (волну, полуволну) и какое время между отсчетами
(частоту выборки сигнала)?