реклама на сайте
подробности

 
 
3 страниц V  < 1 2 3 >  
Reply to this topicStart new topic
> FFT в db, БПФ
Xenia
сообщение Jun 20 2014, 20:25
Сообщение #16


Гуру
******

Группа: Модератор FTP
Сообщений: 4 479
Регистрация: 20-02-08
Из: Москва
Пользователь №: 35 237



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

Откладывать на спектре следует не sin- и cos- составлющие (они дёргаются, переходя друг в друга, из-за того, что фаза плывёт), а их "гипотенузу" sqrt(sin^2+cos^2) - это будет амплитуда, от фазы независящая. Тогда эффект "шершавости", возможно, сам пропадет, без всяких окон.
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 20 2014, 20:31
Сообщение #17


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(ViKo @ Jun 20 2014, 23:16) *
Вижу - приблизились к нужной картинке. Вижу амплитуду в 2 раза меньше, чем без окна. Вижу не палку в спектре, а узкий колокол.

Не приблизились к нужной картинке, а сделали всё согласно правилам, за ширину "палки в спектре" и будет отвечать количество точек FFT. На любом спектроанализаторе Вы увидите такой же результат.
Go to the top of the page
 
+Quote Post
ViKo
сообщение Jun 20 2014, 20:40
Сообщение #18


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



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

На всех показанных картинках именно она, амплитуда, и есть.

Цитата(doom13 @ Jun 20 2014, 23:31) *
Не приблизились к нужной картинке, а сделали всё согласно правилам, за ширину "палки в спектре" и будет отвечать количество точек FFT. На любом спектроанализаторе Вы увидите такой же результат.

Задайте в вашем примере целое число периодов исследуемой частоты, и покажите результат БПФ, без окна и с окном.
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 20 2014, 21:05
Сообщение #19


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(Xenia @ Jun 20 2014, 23:25) *
Правую половину спектра надо отрезать - это зеркальная часть, побочный продукт алгоритма FFT. Тем более что по теории число коэффициентов не может превышать частоту Найквиста (половины числа точек в исходном массиве).

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

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

Вообще-то там и показана амплитуда:
Код
sqrt((Re{FFT})^2 + (Im{FFT})^2)

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

Цитата(ViKo @ Jun 20 2014, 23:40) *
Задайте в вашем примере целое число периодов исследуемой частоты, и покажите результат БПФ, без окна и с окном.

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

Эскизы прикрепленных изображений
Прикрепленное изображение
Прикрепленное изображение
Прикрепленное изображение
 
Go to the top of the page
 
+Quote Post
Xenia
сообщение Jun 20 2014, 21:14
Сообщение #20


Гуру
******

Группа: Модератор FTP
Сообщений: 4 479
Регистрация: 20-02-08
Из: Москва
Пользователь №: 35 237



Цитата(doom13 @ Jun 21 2014, 01:05) *
Вообще-то там и показана амплитуда:
Код
sqrt((Re{FFT})^2 + (Im{FFT})^2)


А чем тогда синенький от красненького отличается? sm.gif А то я поняла так, что один из них sin, а другой cos.
Go to the top of the page
 
+Quote Post
ViKo
сообщение Jun 20 2014, 21:16
Сообщение #21


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(doom13 @ Jun 21 2014, 00:05) *
Ну вот здесь получается, видимо, частный случай, не всегда в нужную выборку влезет целое число периодов. Так что на это, думаю, полагаться не стоит.

Возвращаемся к началу дискуссии... Полагаться, как вы говорите, не на что. В одном случае будут красивые "палки", в другом - не очень. Как повезет.
Теоретики придумали десятка два оконных функций, каждая имеет некую особенность, и еще нужно вспомнить, какую когда включать. И ни одна не гарантирует идеальную достоверность амплитуды и частоты. laughing.gif
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 20 2014, 21:17
Сообщение #22


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(Xenia @ Jun 21 2014, 00:14) *
А чем тогда синенький от красненького отличается? sm.gif А то я поняла так, что один из них sin, а другой cos.

Картинки слева на право:
Цитата(doom13 @ Jun 20 2014, 23:05) *
1) Сигнал и сигнал после применения оконной функции (окно Чебышева из Matlab - chebwin(1024,80)).
2) Спектр сигнала без использования оконной функции.
3) Спектр сигнала с использованием оконной функции.

Ну и цветом конечно же sm.gif
Go to the top of the page
 
+Quote Post
ViKo
сообщение Jun 20 2014, 21:17
Сообщение #23


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(Xenia @ Jun 21 2014, 00:14) *
А чем тогда синенький от красненького отличается? sm.gif А то я поняла так, что один из них sin, а другой cos.

В первом на сигнал не наложена оконная функция (плоское окно), в другом - ... наложена sm.gif
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 20 2014, 21:28
Сообщение #24


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(ViKo @ Jun 21 2014, 00:16) *
Возвращаемся к началу дискуссии... Полагаться, как вы говорите, не на что. В одном случае будут красивые "палки", в другом - не очень. Как повезет.
Теоретики придумали десятка два оконных функций, каждая имеет некую особенность, и еще нужно вспомнить, какую когда включать. И ни одна не гарантирует идеальную достоверность амплитуды и частоты. laughing.gif

Красивые, как Вы говорите, палки без использования окна будут в одном частном случае, если в рассматриваемую выборку помещается целое число периодов сигнала (из всего множества анализируемых частот - одна совпадёт и то, если повезёт).
С оконной функцией мы получаем красивый спектр во всех случаях.
Go to the top of the page
 
+Quote Post
ViKo
сообщение Jun 20 2014, 21:39
Сообщение #25


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(doom13 @ Jun 21 2014, 00:28) *
Красивые, как Вы говорите, палки без использования окна будут в одном частном случае, если в рассматриваемую выборку помещается целое число периодов сигнала (из всего множества анализируемых частот - одна совпадёт и то, если повезёт).
С оконной функцией мы получаем красивый спектр во всех случаях.

Не одна, а целая куча. Их же много может быть, кратных.
Красивый спектр - именно синенький. Такой и хотел видеть топикстартер. А красненький - только приближение к нему.
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 20 2014, 21:54
Сообщение #26


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(ViKo @ Jun 21 2014, 00:39) *
Не одна, а целая куча. Их же много может быть, кратных.
Красивый спектр - именно синенький. Такой и хотел видеть топикстартер. А красненький - только приближение к нему.

То, что хотел увидеть топикстартер, он увидел выше и выразил свою благодарность за помощь в данном вопросе. Вам советую почитать хотя бы это. Думаю, спор далее не уместен и нужно закрывать тему.
Go to the top of the page
 
+Quote Post
ViKo
сообщение Jun 21 2014, 07:18
Сообщение #27


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(doom13 @ Jun 21 2014, 00:54) *
То, что хотел увидеть топикстартер, он увидел выше и выразил свою благодарность за помощь в данном вопросе. Вам советую почитать хотя бы это. Думаю, спор далее не уместен и нужно закрывать тему.

А мы спорим? sm.gif Мы говорим об одном и том же, только по-разному оцениваем степень "волшебства" оконных функций. Определяется упертостью оппонента. За совет спасибо, только я это читал и применял раньше вас (я так думаю). rolleyes.gif
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 21 2014, 08:28
Сообщение #28


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(ViKo @ Jun 21 2014, 10:18) *
А мы спорим? sm.gif Мы говорим об одном и том же, только по-разному оцениваем степень "волшебства" оконных функций. Определяется упертостью оппонента.

По всем признакам похоже на спор, и говорили мы, на сколько я помню, о разном: Вы о том, что нет необходимости в использовании оконной функции, я о том, что есть. На счёт степени упёртости оппонента позвольте с Вами согласиться, вообще нет ничего призрительного в признании своей неправоты и это не только моя оценка:
Цитата(Fat Robot @ Jun 20 2014, 19:37) *
Желание быть в оппозиции подчас "превозмогает и рассудок".


Цитата(ViKo @ Jun 21 2014, 10:18) *
За совет спасибо, только я это читал и применял раньше вас (я так думаю). rolleyes.gif

Читать то Вы читали, но не совсем то, о чём шла речь с самого начала темы. Ваши слова противоречат друг другу:
Цитата(ViKo @ Jun 20 2014, 22:11) *
Говорю, основываясь на собственном опыте. Точек, правда у меня 1024. Картинку, как на среднем рисунке у ТС, имел. А как на правом рисунке - нет. Перебираю 4 различных оконных функции + плоское. Если взять точек побольше, будет ближе к идеалу. А так, даже об амплитуде можно только гадать.

Цитата(ViKo @ Jun 20 2014, 13:34) *
Никогда не получит, ни при каких окнах. laughing.gif


То, о чём был начальный вопрос, и получается с использованием окна.
Go to the top of the page
 
+Quote Post
prig
сообщение Jun 23 2014, 07:36
Сообщение #29


Знающий
****

Группа: Свой
Сообщений: 869
Регистрация: 30-01-08
Из: СПб
Пользователь №: 34 595



Цитата(ViKo @ Jun 20 2014, 18:01) *
А как его узнать заранее, "правильное"...? Если частота сигнала известна, то и спектр не нужен. rolleyes.gif
А если в окно укладывается 20,5 периодов?

В реальных задачах, например OFDM, решают задачу синхронизации.
Соответственно, условия выборки определяются далеко не произвольно.
Зачем нужен спектр при известной частоте, за меня уже ответили.


Цитата(doom13 @ Jun 21 2014, 12:28) *
...Вы о том, что нет необходимости в использовании оконной функции, я о том, что есть.
...

Вообще-то, от конкретной задачи зависит. Если вспомнить тот же OFDM, то после применения оконной функции, хоть сколько-нибудь существенно отличной от прямоугольной, о собственно OFDM можно забыть.
Go to the top of the page
 
+Quote Post
doom13
сообщение Jun 23 2014, 09:08
Сообщение #30


Профессионал
*****

Группа: Свой
Сообщений: 1 404
Регистрация: 11-03-11
Из: Минск, Беларусь
Пользователь №: 63 539



Цитата(prig @ Jun 23 2014, 10:36) *
Вообще-то, от конкретной задачи зависит. Если вспомнить тот же OFDM, то после применения оконной функции, хоть сколько-нибудь существенно отличной от прямоугольной, о собственно OFDM можно забыть.

С OFDM не работал, видимо это так, в моих задачах применял, и для оценки ДД АЦП и при обработке. Тут, наверное, не совсем правильно было сказано и имелось ввиду использование окна применительно к поставленному вопросу. Вы верно дополнили.
Go to the top of the page
 
+Quote Post

3 страниц V  < 1 2 3 >
Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 19th July 2025 - 09:47
Рейтинг@Mail.ru


Страница сгенерированна за 0.01514 секунд с 7
ELECTRONIX ©2004-2016