|
|
ru.algorithms- RU.ALGORITHMS ---------------------------------------------------------------- From : George Shepelev 2:461/124 25 Oct 2001 12:55:31 To : Alegz Lee Subject : ДПФ, БПФ --------------------------------------------------------------------------------
Alegz, ты ещё здесь сидишь?
Среда Октябрь 24 2001 15:47, Alegz Lee wrote to Evgenij Masherov:
AL>>> Для дипломной работы понадобились вычисления ДПФ, БПФ. Какие
AL>>> есть алгоритмы на этот счёт и их реализации? (Ожидаемое
AL>>> количество точек 6000 - 10 000).
EM>> БПФ - неоднократно постились в эху. Могу выслать на С мылом.
AL> Да, пожалуйста. Можно и на baliset[at]magnetic.org.ru
Важно только не забывать о недостатках БПФ. Это преобразование
корректно применяется _ИСКЛЮЧИТЕЛЬHО_ к _периодическим_ сигналам
с периодом, равным периоду обработки. Т.е. на самом деле это
преобразование определяется значениями _всех_ выборок по времени
от нуля до бесконечности.
Иногда это преобразование (для простоты) применяют к "окну"
непериодической функции, что может приводить к _сколь угодно
большим_ ошибкам (см. теорему Агеева) :-\
Предлагавшиеся "хитрые приёмы обхода проблемы", вроде
"обнуления" значений всех остальных точек при расчёте
спектра или "зеркальных отражений окна" данную проблему
не решают, поскольку исходят из молчаливого предположения
о взаимной независимости "окон". Этот тезис не соответствует
истине, раз уж мы взялись работать со _спектральным_
представлением сигнала, нельзя игнорировать факт
"бесконечности" каждой из спектральных составляющих
во времени, её "не-локализуемости".
AL>>> Также интересуют алгоритмы построения огибающей заданной
AL>>> гладкости
EM>> Преобразование Гильберта, дающее мнимую часть аналитического
EM>> сигнала. Модуль комплексного аналитического сигнала дает
EM>> огибающую.
AL> А можно об этом подробнее?
Если _уже_ получен спектр сигнала, особых проблем не должно быть.
Считается (для удовлетворения нескольким разумным критериям),
что исходный сигнал задаётся комплексной функцией от времени
.
u(t) = u(t) + jv(t), где j - мнимая единица
так, что сопряжённый с u(t) сигнал v(t) строится из u(t) поворотом фаз
всех спектральных составляющих на -pi/2 в области положительных частот
и на pi/2 в области отрицательных частот.
После чего огибающая получается так:
______________
. / 2 2
A(t) = | u(t) | = |/ u (t) + v (t)
Ежли кто-то лучше разбирается в математике - пусть поправит/дополнит...
Георгий
[ignore-list "TWO SKUNKS" Team]
--- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- ---
* Origin: Должен же быть кто-то умнее? (2:461/124)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
Архивное /ru.algorithms/18153bd815f4.html, оценка из 5, голосов 10
|