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

 
 
2 страниц V  < 1 2  
Reply to this topicStart new topic
> Не могу понять ДПФ, вопрос чайника
kumle
сообщение Aug 25 2011, 08:09
Сообщение #16


Частый гость
**

Группа: Участник
Сообщений: 149
Регистрация: 15-12-09
Из: Москва
Пользователь №: 54 280



Цитата(SPACUM @ Aug 25 2011, 11:56) *
k - это нормированная частота. Если длина выборки = Т, то частота F = k / T.
Если Т = 1 секунда, то к - это частота в Герцах.
Разумеется к = 0, 1, 2, 3...(N/2 - 1) целое число.


Тогда в формуле аргумент косинуса у меня будет 2*pi*f*T*n/N
Go to the top of the page
 
+Quote Post
bahurin
сообщение Aug 25 2011, 08:40
Сообщение #17


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



Цитата(kumle @ Aug 25 2011, 11:34) *
И если я хочу вычислить амплитуду какой нибудь составляющей спектра я должен свои оцифрованные выборки
подставить в эту формулу.


Ничего подобного!!! Эта формула для реальной части спектра есть еще формула для мнимой части там вместо косинуса синус и чтобы получить амплитудный спектр вы должны взять корень квадратный из суммы квадратов реальной и мнимой частей спектра.
Go to the top of the page
 
+Quote Post
kumle
сообщение Aug 25 2011, 08:49
Сообщение #18


Частый гость
**

Группа: Участник
Сообщений: 149
Регистрация: 15-12-09
Из: Москва
Пользователь №: 54 280



Цитата(bahurin @ Aug 25 2011, 12:40) *
Ничего подобного!!! Эта формула для реальной части спектра есть еще формула для мнимой части там вместо косинуса синус и чтобы получить амплитудный спектр вы должны взять корень квадратный из суммы квадратов реальной и мнимой частей спектра.


Точно !
Это тоже было в книжке, я поторопил события и неуспел прочитать !
Спавсибо !
Теперь точно смогу вычислить.
Прикрепленное изображение
Go to the top of the page
 
+Quote Post
Дмитрий_Б
сообщение Aug 25 2011, 17:15
Сообщение #19


Местный
***

Группа: Участник
Сообщений: 211
Регистрация: 25-10-09
Пользователь №: 53 195



Цитата(Alexey Lukin @ Aug 25 2011, 00:27) *
Это классическая книга Рабинера-Гоулда. Никто не говорит, что X(k) в верхней формуле равно X(k) в нижней. Формула просто переписывается в другом виде, а значения коэффициентов Фурье при этом могут измениться.

Вот как? Означает ли это, что единственному сигналу соответствуют 2 различных спектра? Не нравится мне математика нетрадиционной ориентации. Сигналу соответствует единственный набор коэффициентов ряда Фурье. Они меняются в том и только в том случае, если меняется сигнал. Что до книги - мне, помнится, попадались в ней неточности, хотя книга хорошая.
Go to the top of the page
 
+Quote Post
bahurin
сообщение Aug 26 2011, 04:08
Сообщение #20


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



Цитата(Дмитрий_Б @ Aug 25 2011, 21:15) *
Вот как? Означает ли это, что единственному сигналу соответствуют 2 различных спектра? Не нравится мне математика нетрадиционной ориентации. Сигналу соответствует единственный набор коэффициентов ряда Фурье. Они меняются в том и только в том случае, если меняется сигнал. Что до книги - мне, помнится, попадались в ней неточности, хотя книга хорошая.


Да ну не цепляйтесь там по сути написано правильно, хотя и запутано.

Что касается
Цитата
Сигналу соответствует единственный набор коэффициентов ряда Фурье


ряд Фурье справедлив для периодического сигнала, тогда его спектр линейчатый. Но если после АЦП я получил N точек сигнала, то его спектр не является линейчатым, он непрерывный, хотя периодический, потому что сигнал цифровой. И чтобы сделать спектр линейчатым я должен свои N точек бесконечно повторить (сделать сигнал периодическим), но это уже совсем другой сигнал, поскольку о периодичности исходного сигнала, который я оцифровал ничего не известно. Вот и получается что на входе АЦП мы имеем один сигнал, а спектр рассчитываем для совершенно другого сигнала, со всеми эффектами возникающими при искусственной периодизации исходного оцифрованного сигнала.
Go to the top of the page
 
+Quote Post
Driver_GV
сообщение Aug 26 2011, 19:31
Сообщение #21


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-05-06
Из: Lugansk
Пользователь №: 17 011



Цитата(bahurin @ Aug 26 2011, 07:08) *
ряд Фурье справедлив для периодического сигнала, тогда его спектр линейчатый. Но если после АЦП я получил N точек сигнала, то его спектр не является линейчатым, он непрерывный, хотя периодический, потому что сигнал цифровой. И чтобы сделать спектр линейчатым я должен свои N точек бесконечно повторить (сделать сигнал периодическим), но это уже совсем другой сигнал, поскольку о периодичности исходного сигнала, который я оцифровал ничего не известно. Вот и получается что на входе АЦП мы имеем один сигнал, а спектр рассчитываем для совершенно другого сигнала, со всеми эффектами возникающими при искусственной периодизации исходного оцифрованного сигнала.


Используйте функции окна. Входные отсчеты по времени умножаются на соответствующую функцию окна, что влечет за собой обнуление сигнала на краях выборки. Это компромис между расширением центрального лепеска и подавлением боковых. я использую обычно Хамминга. Подавление боковых более 40 Дб.


--------------------
-------------------
picping.lg.ua
Go to the top of the page
 
+Quote Post
bahurin
сообщение Aug 27 2011, 06:54
Сообщение #22


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



Цитата(Driver_GV @ Aug 26 2011, 23:31) *
Используйте функции окна. Входные отсчеты по времени умножаются на соответствующую функцию окна, что влечет за собой обнуление сигнала на краях выборки. Это компромис между расширением центрального лепеска и подавлением боковых. я использую обычно Хамминга. Подавление боковых более 40 Дб.


Это понятно но мы не об этом. Умножение на функцию окна отличного от прямоугольного это уже третий сигнал со своим спектром. Выборка N точек сигнала это тоже умножение на функцию окна, только прямоугольного.
Go to the top of the page
 
+Quote Post
Дмитрий_Б
сообщение Aug 28 2011, 10:27
Сообщение #23


Местный
***

Группа: Участник
Сообщений: 211
Регистрация: 25-10-09
Пользователь №: 53 195



Цитата(bahurin @ Aug 26 2011, 08:08) *
Да ну не цепляйтесь там по сути написано правильно, хотя и запутано.

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


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

Сообщение отредактировал Дмитрий_Б - Aug 28 2011, 10:31
Go to the top of the page
 
+Quote Post
Zelepuk
сообщение Aug 28 2011, 14:59
Сообщение #24


Знающий
****

Группа: Участник
Сообщений: 634
Регистрация: 27-10-10
Пользователь №: 60 464



Цитата(Дмитрий_Б @ Aug 28 2011, 14:27) *
Не стал бы, поскольку и так это знаю. Но человек не может разобраться, поскольку математика хромает.
Рабинер и Голд не годится как учебник.


поэтому и рекомендую для ТС почитать лучше Сато или Кестера. rolleyes.gif
Go to the top of the page
 
+Quote Post
Serhiy_UA
сообщение Aug 29 2011, 06:01
Сообщение #25


Знающий
****

Группа: Свой
Сообщений: 721
Регистрация: 23-10-08
Из: next to Odessa
Пользователь №: 41 112



Цитата(Zelepuk @ Aug 28 2011, 17:59) *
поэтому и рекомендую для ТС почитать лучше Сато или Кестера. rolleyes.gif

Сообщите названия книг и год выпуска, пожалуйста. Эти авторы наверняка выпустили в свет не одну книгу...
Go to the top of the page
 
+Quote Post

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

 


RSS Текстовая версия Сейчас: 15th June 2025 - 17:41
Рейтинг@Mail.ru


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