вопрос по БПФ Опции

[Sneg_87]

"Частота дискретизации сигнала равна 44100 Гц. Размер БПФ равен 4096. Продолжительность сигнала из 4096 точек при данной частоте дискретиза-ции составляет 0.0929 секунды. Период первой синусоиды равен 4096 точкам, что по времени состав-ляет примерно 0.0929 секунды. Значит, частота первой синусоиды будет 10.77 Гц. Частота второй будет 21.53 Гц. И так далее."
а если частота сигнала, ну скажем, равна 15Гц, то как это отразится на преобразовании Фурье? разрешения по частоте не хватит для выявления этой частоты или же составляющая на этой частоте преобразованием будет определена. Помогите решить данный вопрос

[SPACUM]

Привет.
1.Очень советую написать самому програмку с графическим выводом и попробовать.
2.Или разобраться с Матлабом и попробовать.
3.Или попробовать на странице http://www.fourier-series.com/
Вообще БПФ = чистая математика и сильно отличается от бытовых представлений об частотах спектрах итд.

[fontp]

Вообще БПФ = чистая математика и сильно отличается от бытовых представлений об частотах спектрах итд.

Не особенно. В этом смысле и дискретный сигнал сильно отличается от непрерывного.
ДПФ (БПФ это реализация ДПФ) - это дискретизация непрерывного физического спектра взятого с функцией окна. Дискретные отсчеты свертки физического спектра с функцией окна в спектральной области. Если окно не применяется, то вставится прямоугольное окно (спектральные sin(x)/x) помимо нашего желания.

Поэтому каждый отсчет ДПФ для непрерывного спектра формируется суммой от всех частот, но с весом от функции окна. Функция окна быстро затухает по разности частот. А значит каждый отсчет ДПФ преобразования формируется непрерывным физическим спектром Фурье только от близких частот к данному бину . Ситуация ничем не отличается от той которая возникнет, если поставить банк узкополостных фильтров на соответствующих бинах ДПФ. Но ДПФ ещё и обладает полезными математическими свойствами. Спектральное разрешение в этих случаях непрерывного спектра сигнала всегда будет определяться шириной функции окна, которая в спектральной области в соответствии с критерием "неопределённости" Рэлея будет 1/N (или 1/T в размерных единицах)

Если мы знаем, что комплексная синусоида вообще одна на фоне шума - то в ДПФ мы в точности имеем отдискретизированую в спектральной области функцию окна, центрированую на частоте синусоиды. Поэтому всегда, независимо от того пападает ли частота кратно на бины ДПФ, мы можем увидев эту функцию окна в полученом ДПФ, провести интерполяцию и найти частоту, амплитуду и фазу этой синусоиды. Причем Райф и Бурстин доказали, что в случае одиночной синусоиды наибольшую точность даёт прямоугольное (т.е. никакое) окно. Они же предложили проводить интерполяцию посредством добавления нулей в данные и квадратичной интерполяцией в окрестности максимума.

Если спект линейчатый и гармоники находятся далеко друг от друга, эта же методика позволяет получать очень точные оценки этих синусоид, но с применением функций окон, изолирующих эти линии в спектре.Если на каждую линию спектра поставить функцию окна с соответствующей амплитудой и просумировать, то это то что мы получим в ДПФ и мы снова сможем проводить интерполяцию в том случае, если эти оконные отклики перекрываются слабо.

Походу интерполяция добавлением нулей и подгонки параболы фиттингом - не единственый способ интерполяции спектра вблизи максимума спектральной линии. Есть методы производящие "внутреннюю интерполяцию", без всякого добавления нулей. Вот как это всё выглядит в Матлабе

http://home.comcast.net/~kootsoop/EricJ2/index.htm
Лучший из них не очень давно предложен в работе МакЛеода (не путать с горцем Маклаудом и не размахивать здесь саблей!)

Прикрепленные файлы

 Macleod_Fast_N_ML_Estimate.pdf ( 196.95 килобайт ) Кол-во скачиваний: 312

 

[SPACUM]

Лучший из них не очень давно предложен в работе МакЛеода (не путать с горцем Маклаудом и не размахивать здесь саблей!)

Спасибо за МакЛеода. А нет ли у Вас [24] из него. Я делаю сейчас прибор с применением похожего принципа(чистая отсебятина, но работает) , мне иинтересно.

[fontp]

Спасибо за МакЛеода. А нет ли у Вас [24] из него. Я делаю сейчас прибор с применением похожего принципа(чистая отсебятина, но работает) , мне иинтересно.

Quinn II

Прикрепленные файлы

 Quinn_Estimation.pdf ( 174.48 килобайт ) Кол-во скачиваний: 207

 

[SPACUM]

Quinn II

Приятно читать матерого математика. Если бы прочел раньше - ничего бы не сделал, а теперь свои алгоритмы есть и к процессору адаптированы. Может выложите [1] из Quinn II.(Это последняя просьба). Еще мне интересно в каких приборах все эти алгоритмы применяются и зачем Кьюин на свой европатент получил. А можно ли по патенту определить где его применили?

[fontp]

Приятно читать матерого математика. Если бы прочел раньше - ничего бы не сделал, а теперь свои алгоритмы есть и к процессору адаптированы. Может выложите [1] из Quinn II.(Это последняя просьба). Еще мне интересно в каких приборах все эти алгоритмы применяются и зачем Кьюин на свой европатент получил. А можно ли по патенту определить где его применили?

Везде для синхронизации несущей в связи.
Голая несущая или модулированый сигнал после нелинейного преобразования - это комплексная экспонента в шуме.
Её, конечно, можно захватить ФАПЧ, но это во-первых долго, а во-вторых - не оптимально по отношению к шумам, поскольку нужно для захвата расширять шумовую полосу, а для штатной работы - сужать. Значительно эффективней быстро измерить частоту и фазу оптимальным feed forward методом, а потом уже запустить узкополосный ФАПЧ без частотной и фазовой расстройки. Для burst модемов TDMA, когда несущая меняется на каждом фрейме - это вообще единственный способ захвата несущей по преамбуле, никакой ФАПЧ не поможет
Только совершенно всё равно не понятно, как Quinn собирался узнать какой модем использует его вычислительный метод при синхронизации. Видимо wishfull thinking, вдруг кто сам заплатит
Quinn-I

Прикрепленные файлы

 Quinn1_EstimateFrequency.pdf ( 477.49 килобайт ) Кол-во скачиваний: 363

 

[SPACUM]

А если что-нибудь новое создать или много лучше старого. Можно своим именем назвать или надо регистрировать?

[AndrewN]

А если что-нибудь новое создать или много лучше старого. Можно своим именем назвать или надо регистрировать?

Не, по традиции имя посерёдке вставляют: Жан Батист Жозеф Кули Тьюки Виноград Темпертон Фурье. Добавьте ближе к концу...
Кстати, имя изобретателя китайской теоремы об остатках не вставили, забыли :)

[SPACUM]

Замечательно. Еще Шопенгауэр придумал (ЭРИСТИКА, ИЛИ ИСКУССТВО ПОБЕЖДАТЬ В СПОРАХ). Обязательно воспользуюсь. С китайцами действительно непонятки.