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

 
 
2 страниц V  < 1 2  
Reply to this topicStart new topic
> FFT импульсного сигнала, как уменьшить искажения
_Anatoliy
сообщение Jun 10 2009, 18:08
Сообщение #16


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 10 2009, 18:47) *
Вы можете использовать perfect reconstruction critically-sampled filter banks, подполосы будут иметь вид не sin/x а могут иметь любое подавление за пределами своей полосы, перед вашим FFT добавится полифазный фильтр с определённой импульсной характеристикой.

Спасибо за подсказку!Буду искать информацию.
И тогда можно будет применить прямоугольное окно?
Go to the top of the page
 
+Quote Post
petrov
сообщение Jun 10 2009, 18:17
Сообщение #17


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(_Anatoliy @ Jun 10 2009, 22:08) *
Спасибо за подсказку!Буду искать информацию.
И тогда можно будет применить прямоугольное окно?


Окна с FFT это частный случай банка фильтров. Конструируя определённым образом фильтр прототип можно получить большие возможности, например большую избирательность при небольшом количестве подполос.
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 10 2009, 18:25
Сообщение #18


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 10 2009, 19:17) *
Окна с FFT это частный случай банка фильтров. Конструируя определённым образом фильтр прототип можно получить большие возможности, например большую избирательность при небольшом количестве подполос.

Спасибо!Я правильно понял что нужно отказаться от БПФ?
Go to the top of the page
 
+Quote Post
petrov
сообщение Jun 11 2009, 08:37
Сообщение #19


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(_Anatoliy @ Jun 10 2009, 22:25) *
Спасибо!Я правильно понял что нужно отказаться от БПФ?


FFT используется в банках фильтров. Я так понимаю что сейчас у вас количество подполос избыточно и выбирается исходя из требований по подавлению соседних каналов, с банком фильтров можно будет уменьшить в несколько раз количество подполос, сохранив избирательность по соседним каналам, при этом на выходе банка фильтров у вас будет блок отсчётов частоты на который будет попадать пакет целиком а не на край окна как сейчас.
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 11 2009, 10:17
Сообщение #20


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 11 2009, 10:37) *
FFT используется в банках фильтров. Я так понимаю что сейчас у вас количество подполос избыточно и выбирается исходя из требований по подавлению соседних каналов, с банком фильтров можно будет уменьшить в несколько раз количество подполос, сохранив избирательность по соседним каналам, при этом на выходе банка фильтров у вас будет блок отсчётов частоты на который будет попадать пакет целиком а не на край окна как сейчас.


Почему Вы считаете что у меня количество подполос избыточно?У меня есть ещё такой параметр как разрешение по частоте и он равен Fs/512,поэтому я выполняю БПФ размером 1024 точки и работаю с первой зоной Найквиста. Если я уменьшу количество подполос то потеряю разрешение.
Go to the top of the page
 
+Quote Post
petrov
сообщение Jun 11 2009, 10:51
Сообщение #21


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(_Anatoliy @ Jun 11 2009, 14:17) *
Почему Вы считаете что у меня количество подполос избыточно?У меня есть ещё такой параметр как разрешение по частоте и он равен Fs/512,поэтому я выполняю БПФ размером 1024 точки и работаю с первой зоной Найквиста. Если я уменьшу количество подполос то потеряю разрешение.


В общем FFT это всего лишь быстрый алгоритм, по сути у вас куча гетеродинов и ФНЧ фильтров спрямоугольной ИХ(или что-то типа колокольчика в зависимости от окна) и последующей децимацией, с банком фильтров вы не ограничиваете себя в виде ИХ анализирующего фильтра, она может быть подобной sinx/х, нужно вам это или нет решайте сами.

Кстати а не слишком ли короткие пакеты получаются у вас, при которых у вас возникают проблемы, из скольки символов QPSK они состоят, несколько десятков, такие в жизни встречаются?
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 11 2009, 11:48
Сообщение #22


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 11 2009, 11:51) *
В общем FFT это всего лишь быстрый алгоритм, по сути у вас куча гетеродинов и ФНЧ фильтров спрямоугольной ИХ(или что-то типа колокольчика в зависимости от окна) и последующей децимацией, с банком фильтров вы не ограничиваете себя в виде ИХ анализирующего фильтра, она может быть подобной sinx/х, нужно вам это или нет решайте сами.

Кстати а не слишком ли короткие пакеты получаются у вас, при которых у вас возникают проблемы, из скольки символов QPSK они состоят, несколько десятков, такие в жизни встречаются?

1).Я думаю что простое изменение ИХ не решит проблему в комлексе.Вытащу амплитуду - загроблю избирательность и наоборот.
2).Встречаются.Подробности сказать не могу,но сигналы будут не только QPSK.
Go to the top of the page
 
+Quote Post
petrov
сообщение Jun 11 2009, 12:21
Сообщение #23


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(_Anatoliy @ Jun 11 2009, 15:48) *
1).Я думаю что простое изменение ИХ не решит проблему в комлексе.Вытащу амплитуду - загроблю избирательность и наоборот.


Возможно как раз наоборот банк фильтров позволит получит новое качество анализа.



Прикрепленное изображение




Цитата(_Anatoliy @ Jun 11 2009, 15:48) *
2).Встречаются.Подробности сказать не могу,но сигналы будут не только QPSK.


Какой-нибудь frequency-hopping для вояк анализируете? smile.gif

Просто тогда непонятно что вы хотите сделать и соответственно ничего и не посоветуешь...
Go to the top of the page
 
+Quote Post
samurad
сообщение Jun 11 2009, 16:05
Сообщение #24


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

Группа: Свой
Сообщений: 121
Регистрация: 9-05-08
Из: Япония
Пользователь №: 37 385



Цитата(petrov @ Jun 11 2009, 15:21) *
...Какой-нибудь frequency-hopping для вояк анализируете? smile.gif
Просто тогда непонятно что вы хотите сделать и соответственно ничего и не посоветуешь...

Задача похожа на UWB радар. Если я правильно понимаю,
1) точное время прихода нужного импульса неизвестно, а известно только окно (диапазон задержек), причем в таком окне другие импульсы не ожидаются (что характерно для радара),
2) амплитуда импульса слишком мала, даже после LNA, чтобы можно было его принять неоптимальным методом, (что характерно для UWB)
3) невозможно построить согласованный фильтр или коррелятор на каждый из возможных импульсов,
4) невозможно построить какой-л. фильтр с полосой примерно равной полосе принимаемого импульса, а можно только фильтр с много более узкой полосой (что странно для активного радара, но возможно для пассивного),
5) импульсы, приходящие с разной задержкой, в общем случае некоррелированы друг с другом ни по форме комплексной огибающей, ни по полосе (множественные цели?).

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

Можно посоветовать "загонять" импульсы в рециркулятор с активным усилителем в петли обратной связи и размыканием входной цепи рециркулятора после первого прохождения импульса. Выход рециркулятора модулировать известной последовательностью узкополосных импульсов с частотой немного отличающейся от частоты следования. Пропуская такую квантованную последовательность через скользящий коррелятор, можно различить пик с определенной переодичностью, зависящей от соотношения частот следования входного импульса и модулирующей последовательности. Полоса и амплитуда входного импульса будут определяться по полосе и амплитуде комплексной огибающей на выходе скользящего коррелятора.
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 16 2009, 10:10
Сообщение #25


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 11 2009, 14:21) *


С короткими импульсами разобрался,теперь у меня длительность импульса может быть 0,05...1*(время наблюдения) без существенных искажений амплитуды на выходе БПФ и боковые лепестки задавлены на 60дБ.
Теперь нужно побороть следствие окна - расширение основного лепестка.Это нужно сделать без увеличения времени наблюдения.Коллеги укажите,плз,правильное направление поиска по теме!Вообще это решаемо?
Go to the top of the page
 
+Quote Post
petrov
сообщение Jun 16 2009, 11:57
Сообщение #26


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(_Anatoliy @ Jun 16 2009, 14:10) *
С короткими импульсами разобрался,теперь у меня длительность импульса может быть 0,05...1*(время наблюдения) без существенных искажений амплитуды на выходе БПФ и боковые лепестки задавлены на 60дБ.
Теперь нужно побороть следствие окна - расширение основного лепестка.Это нужно сделать без увеличения времени наблюдения.Коллеги укажите,плз,правильное направление поиска по теме!Вообще это решаемо?


Например можно сделать банк фильтров на основе анализирующего фильтра типа приподнятый косинус, коэффициент ската 1, ширина основного лепестка во времени 2T, ессно во времени будут ещё боковые лепестки, ширина основного лепестка по частоте по уровню 3 дБ 1/T, ширина основного лепестка по частоте 2/T, уровень боковых лепестков по частоте за пределами основного любой, ограничивается только сложностью полифазного фильтра перед FFT, анализируемые полосы находятся на расстоянии 1/T друг от друга. Но это так примерно, возможно не будет там ортогональности совершенной, подробнее надо исследовать, возможно есть лучшие импульсные характеристики для анализирующего фильтра.
Go to the top of the page
 
+Quote Post
litv
сообщение Jun 17 2009, 06:11
Сообщение #27


Местный
***

Группа: Свой
Сообщений: 401
Регистрация: 6-10-04
Из: Воронеж
Пользователь №: 806



Надо накладывать окна с перекрытием 50% и будет вам счастье.
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 17 2009, 06:27
Сообщение #28


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(petrov @ Jun 16 2009, 13:57) *
Например можно сделать банк фильтров на основе анализирующего фильтра типа приподнятый косинус, коэффициент ската 1, ширина основного лепестка во времени 2T, ессно во времени будут ещё боковые лепестки, ширина основного лепестка по частоте по уровню 3 дБ 1/T, ширина основного лепестка по частоте 2/T, уровень боковых лепестков по частоте за пределами основного любой, ограничивается только сложностью полифазного фильтра перед FFT, анализируемые полосы находятся на расстоянии 1/T друг от друга. Но это так примерно, возможно не будет там ортогональности совершенной, подробнее надо исследовать, возможно есть лучшие импульсные характеристики для анализирующего фильтра.

Ага,спасибо,попробую.

Цитата(litv @ Jun 17 2009, 08:11) *
Надо накладывать окна с перекрытием 50% и будет вам счастье.

Очень интересно,спасибо.А можно чуть подробней?
Go to the top of the page
 
+Quote Post
litv
сообщение Jun 17 2009, 08:55
Сообщение #29


Местный
***

Группа: Свой
Сообщений: 401
Регистрация: 6-10-04
Из: Воронеж
Пользователь №: 806



Обычно все анализаторы спектра работают с перекрытием.
Например tektronix.

http://www.afc-ingenieros.com/NotasTecn/RT...iendo%20FFT.pdf

http://www2.tek.com/cmswpt/tidetails.lotr?...14448&lc=EN

http://www.dliengineering.com/vibman/overlapprocessing.htm
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Jun 17 2009, 09:42
Сообщение #30


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832



Цитата(litv @ Jun 17 2009, 09:55) *
Обычно все анализаторы спектра работают с перекрытием.
Например tektronix.

http://www.afc-ingenieros.com/NotasTecn/RT...iendo%20FFT.pdf

http://www2.tek.com/cmswpt/tidetails.lotr?...14448&lc=EN

http://www.dliengineering.com/vibman/overlapprocessing.htm


Спасибо за информацию!
Сразу не готов сказать подходит к моей задаче или нет.Перевариваю.
Go to the top of the page
 
+Quote Post

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

 


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


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