Определение передаваемых частот Опции

[Kseniia]

Всем доброго времени суток!

Задача: есть DTMF метка, необходимо определить частоты которые передавались, с точностью до 3го знака после запятой. (DTMF - передается сумма двух гармоник, НЧ и ВЧ)

Что проделано:
выделяю из сигнала DTMF, fd=44100 Гц, для того, чтобы N была кратна двум, и равнялась 2048, немного обрезаю сигнал до 0.046 секунд (было 0.05 сек), таким образом T=1/fd=2.268*10^5.
Беру БПФ, строю график зависимости уровня(?)(20 log(sqwr(ReБПФ^2+ImБПФ^2)/fd-N), дБ от частоты, Гц (прикрепляю скрин, freqn=n*fd/N, n=0,...,N-1).

Пики приходятся на 947,46Гц и 1205,9Гц, хотя должны быть 941Гц и 1209Гц. такая разница по видимому из-за огромного ШАГА по частоте, в 21,55 Гц, то есть передаваемые частоты не кратны нашему "шагу".

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

вопрос: Как математически можно определить частоты?

[serjj]

я вот какой выход нашла, и без всякого БПФ... если ошибаюсь - подправьте пожауйста, если Вы именно это и советовали, то я пропустила - прошу прощения...
имеем сигнал s(x), далее находим составляющие коэф корреляции по синусу и по косинусу, и коэф корреляции, строим график, определяем максимумы, видно, что небольшая погрешность есть, но результат очень близок к истине.
...
шаг меняла, от шага зависит только график, а значения Xmax не изменяются.

Приведенные вами формулы не что иное как сэмплирование DTFT вашего сигнала. DFT и FFT тоже сэмплированный DTFT только с другими параметрами сэмплирования. Почитать можно например здесь. Если вы будете уменьшать частотный шаг, вы будете получать все больше точек DTFT функции, но разрешающая способность определяется только длиной исходной выборки вашего сигнала. Интерполяция DFT(FFT) в окрестностях максимума для уточнения результата - это восстановление значений функции DTFT в окрестностях этого самого максимума. Так что приведенные формулы корреляций не дадут результата принципиально иного в сравнении с Фурье преобразованием, т.к. математика там абсолютно одинаковая (внимательно сравните ваши формулы с формулировкой преобразования Фурье). Окно у вас действительно прямоугольное, т.к. данные перед преобразованием вы никак не взвешиваете.