FFT в db, БПФ Опции

[djhall]

Мне необходимо построить график БПФ в том виде как он обычно представляется в даташитах на АЦП, т.е. в децибелах.
У меня есть исходный файл с отсчётами АЦП на который подан синусоидальный сигнал. Нахожу БПФ по 2048 точкам и далее нахожу 10*log10(БПФ) и строю график. Далее показан сигнал и 10*log10(БПФ) (первые две картинки)- явно не совсем хороший спектр. В чём может быть ошибка? Подскажите как правильно построить такой график.

Ожидал увидеть нечто такое (третья картика)

Эскизы прикрепленных изображений

 

[Xenia]

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

Откладывать на спектре следует не sin- и cos- составлющие (они дёргаются, переходя друг в друга, из-за того, что фаза плывёт), а их "гипотенузу" sqrt(sin^2+cos^2) - это будет амплитуда, от фазы независящая. Тогда эффект "шершавости", возможно, сам пропадет, без всяких окон.

[doom13]

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

Вопрос был не в этом, тут сам смысл применения окна показан, а правую половину, если смущает, можно и выбросить.

Откладывать на спектре следует не sin- и cos- составлющие (они дёргаются, переходя друг в друга, из-за того, что фаза плывёт), а их "гипотенузу" sqrt(sin^2+cos^2) - это будет амплитуда, от фазы независящая. Тогда эффект "шершавости", возможно, сам пропадет, без всяких окон.

Вообще-то там и показана амплитуда:

Код

sqrt((Re{FFT})^2 + (Im{FFT})^2)

, так что если правую часть выбросить, то остальное всё верно.

Задайте в вашем примере целое число периодов исследуемой частоты, и покажите результат БПФ, без окна и с окном.

Ну вот здесь получается, видимо, частный случай, не всегда в нужную выборку влезет целое число периодов. Так что на это, думаю, полагаться не стоит.

Эскизы прикрепленных изображений

 

[Xenia]

Вообще-то там и показана амплитуда:

Код

sqrt((Re{FFT})^2 + (Im{FFT})^2)

А чем тогда синенький от красненького отличается? А то я поняла так, что один из них sin, а другой cos.