|
|
Нужна еще тактовая частота Fs, на которой поступают отсчеты сигнала. Я буду считать, что она у Вас Fs=20 Гц.
Если Вам нужно разрешение dF Герц, то необходимо взять число отсчетов равное
N = Fs/dF (с округлением до старшего целого).
При Fs=20Гц и dF=0.5 Гц получаем N=40. Если тактовая частота намного превышает 20 Гц, то возможно имеет смысл понизить частоту дискретизации. Для этого в общем случае необходимо вначале выполнить фильтрацию, оставив только сигнал в полосе 0..10Гц. Однако в Вашем случае этого делать не обязательно, поскольку, как я понял, сигнал у Вас и так находится в данном диапазоне (если нет, то предварительная фильтрация необходима).
Положим, что Fs=20*k, где k - некоторое число большее 2. Тогда частоту дискретизации можно понизить в [k] раз, где операция [.] означает выделение целой части. Таким образом, минимальная длина БПФ, обеспечивающая разрешение dF будет равна
N'=Fs/([k]*dF) (с округлением до старшего целого).
Здесь еще замечу, что как правило, алгоритм БПФ реализован для длин кратных 2^k, например, 16, 32, 64, ... Поэтому на практике округлить результат в приведенной выше формуле придется до ближайшей старшей степени двойки.
Какой алгоритм Radix2 или Radix4 Вы будете использовать - роли никакой не играет. На результате это никак не отразится. Выбор алгоритма БПФ может лишь ускорить быстродействие и/или уменьшить затраты памяти. Во многом выбор алгоритма может быть обусловлен Вашим процессором (от процессора зависит, насколько эффективной окажется реализация, положим Radix4; может оказаться, что на Вашем DSP лучше будет применить Radix2). Но вообще, могу посоветовать заглянуть в библиотеку функций, которая поставляется вместе с Visual DSP. В ней Вы найдете реализацию, которая наверняка является одной из лучших для данного процессора (хотя по опыту могу сказать, что это бывает не всегда так).
Разрядность БПФ (16, 24 или 32 битная точность) выбирайте, исходя из задачи (допустимый уровень шумов).
|
|
