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

 
 
3 страниц V   1 2 3 >  
Reply to this topicStart new topic
> Вычислить разницу фазы.
DASM
сообщение Jan 7 2017, 03:36
Сообщение #1


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Дано. Сигнал и его копия, копия задержана 0.1 сек и пришла другим путем. Источник 455 кГц, монохроматичен и стабилен. Fd 10 msps, 8 bit. Пик фактор сигнала <5. Доплера пока нет в принятых. Нужно найти сдвиг фазы поточнее. Сейчас делаю что, просто тупо два Фурье от обоих, смотрю дельту angle на заданном "бине". Вроде получается но есть вопросы. Разрешение по времени значительно выше чем по Найквисту, получаю несколько бинов со значительной амлитудой. Что с ними делать? Может оконная ф-ция навроде хамминга нужна? Но для фазы. Второе - в половине случаев фурье дает угол меньше на 2пи, приходится прибавлять, чядн? unwrap как то совсем ни к месту, тупо if завел. Ну и есть способы лучше? Может xcorr? Сорри я нуб. .. Наверное надо сделать опорный синус и посчитать интеграл от угла с квадратурой опорной.. только не совсем понимаю как


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
Lmx2315
сообщение Jan 7 2017, 08:36
Сообщение #2


отэц
*****

Группа: Свой
Сообщений: 1 471
Регистрация: 18-09-05
Из: Москва
Пользователь №: 8 684



Говорю на основе теоретических фантазий:
1. Перемножить два сигнала раз они когерентны и отфильтровать - в остатке будет нечто пропорциональное разности фаз.
2. Перевести оба сигнала в комплексный вид и завести на два cordic-а с них снять мгновенные фазы и вычесть их друг с друга - получив разность фаз.


--------------------
"..не нравятся мои выборы? ..приходите в мой суд."
Узурпатор П.
Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 7 2017, 09:15
Сообщение #3


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Цитата(Lmx2315 @ Jan 7 2017, 11:36) *
Говорю на основе теоретических фантазий:
1. Перемножить два сигнала раз они когерентны и отфильтровать - в остатке будет нечто пропорциональное разности фаз.
2. Перевести оба сигнала в комплексный вид и завести на два cordic-а с них снять мгновенные фазы и вычесть их друг с друга - получив разность фаз.

Да зачем мне Кордик, ресурсов полно пока что. Смотрю пока на Фурье, он все равно в итоге тоже самое делает. А зачем в комплексный то переводить? Комплексный спектр и так выйдет после DFT


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
_pv
сообщение Jan 7 2017, 09:52
Сообщение #4


Гуру
******

Группа: Свой
Сообщений: 2 019
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



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


Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 7 2017, 10:42
Сообщение #5


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Цитата(_pv @ Jan 7 2017, 12:52) *
просто корреляция для нескольких значений сдвига между двумя сигналами может оказаться быстрее даже чем FFT.
а потом какую-нибудь параболу наименьшими квадратами на результат натянуть, чтобы положение максимума найти.

да это все замечательно, это будет следующее. Мне сейчас надо вообще как можно качественнее сдвиг этот найти, оптимизации все пока что побоку. Я просто не уверен что делаю все правильно через DFT. Тут эффекты от кратности длины выборки и частоты семплирования вылезают забавные. В общем я похоже просто пока путаюсь при переходе от FTк DFT. Попробую нулей что ли добить.

Цитата(_pv @ Jan 7 2017, 12:52) *
просто корреляция для нескольких значений сдвига между двумя сигналами может оказаться быстрее даже чем FFT.
а потом какую-нибудь параболу наименьшими квадратами на результат натянуть, чтобы положение максимума найти.

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


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
_pv
сообщение Jan 7 2017, 14:23
Сообщение #6


Гуру
******

Группа: Свой
Сообщений: 2 019
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



тестовые данные бы сюда прикрепили.

с Фурье есть вопрос на какой именно частоте разность фаз смотреть.
ну и плюс отсутсвие окна действительно положение максимума спектра будет двигать за счёт наличия "всяких забавных эффектов".
можно конечно и автокорреляцией сначала найти "основную частоту" и потом смотреть сколько там фазы у обоих сигналов просто посчитав интеграл Фурье для этой конкретной частоты, но тогда проще сразу считать взаимную корреляцию двух сигналов и искать на ней максимум.
Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 7 2017, 19:16
Сообщение #7


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



У меня не два разных сигнала, а один, пришедший разными путям и, причем его задержка мне известна, я же его и излучаю. Мне неизвестна разность фаз. А автокорреляция емнип через Фурье и делается


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
_pv
сообщение Jan 7 2017, 19:31
Сообщение #8


Гуру
******

Группа: Свой
Сообщений: 2 019
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



с таким объяснением про известную задерку но неизвестные фазы стало ещё непонятнее sm.gif
если бы данные, по которым ДВЕ спектральные картинки выше построены, сюда прикрепили, было бы гораздо понятнее.
Go to the top of the page
 
+Quote Post
blackfin
сообщение Jan 7 2017, 19:43
Сообщение #9


Гуру
******

Группа: Свой
Сообщений: 2 590
Регистрация: 18-04-05
Пользователь №: 4 261



Цитата(DASM @ Jan 7 2017, 22:16) *
У меня не два разных сигнала, а один, пришедший разными путям и, причем его задержка мне известна, я же его и излучаю. Мне неизвестна разность фаз.

Вроде бы, уже решали эту задачу: Измерение фазового сдвига Sin..
Если сигналов два, вероятно, нужно будет вычислить фазу дважды, для обоих сигналов относительно опорных sin и cos.
Искомая фаза будет равна разности измеренных фаз.. ИМХО..
Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 7 2017, 21:10
Сообщение #10


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Цитата(blackfin @ Jan 7 2017, 22:43) *
Вроде бы, уже решали эту задачу: Измерение фазового сдвига Sin..
Если сигналов два, вероятно, нужно будет вычислить фазу дважды, для обоих сигналов относительно опорных sin и cos.
Искомая фаза будет равна разности измеренных фаз.. ИМХО..

Ну так Фурье оно и есть, только оно по всему спектру. Это и сделал и это уже работает. Удручен дрейфом фазы гетеродина, хоть и кварцеванного. Буду думать.
Задержка известна примерно, там один и тотже сигнал переключается с фиксированной частотой на две антенны. А на приемной имеем.. ну поняли ?
Это я прошелся уже скользящим Фурье по всей большой выборке с перекрывающимся окошком. Видно как плывет фаза sad.gif
У меня три периода, 1-ый работает одна антенна, второй - вторая, третий - молчание в эфире. Вот и видно их и шум между ними некоррелированный. Везде фаза сравнивается по одной частотной полоске из разложения Фурье и везде сравнение с самым первым "окном" в графике, то есть это long term drift. Код приводить не буду он ужасен. Но в принципе все работает
Эскизы прикрепленных изображений
Прикрепленное изображение
 


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
jorikdima
сообщение Jan 7 2017, 23:13
Сообщение #11


тут может быть ваша реклама
*****

Группа: Свой
Сообщений: 1 146
Регистрация: 15-03-06
Из: Санкт-Петербург/CA
Пользователь №: 15 280



Все таки не могу понять, чем не устраивает два измерения фазы (ПЛЛ) и разность между ними? Очень просто и быстро. Зачем фурье приплетать?
Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 8 2017, 02:59
Сообщение #12


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Ну зато сообразил откуда "дрейф", просто я частоту неточно смотрю. Да и дрейф все равно есть Главное что метод работает, и недалек от оптимального


--------------------
Любитель красивых вещей.
Go to the top of the page
 
+Quote Post
Lmx2315
сообщение Jan 8 2017, 08:59
Сообщение #13


отэц
*****

Группа: Свой
Сообщений: 1 471
Регистрация: 18-09-05
Из: Москва
Пользователь №: 8 684



Тему надо было назвать - "хочу измерять фазу через Фурье".


--------------------
"..не нравятся мои выборы? ..приходите в мой суд."
Узурпатор П.
Go to the top of the page
 
+Quote Post
Artem_Petrik
сообщение Jan 8 2017, 14:28
Сообщение #14


Местный
***

Группа: Свой
Сообщений: 443
Регистрация: 22-07-06
Из: Украина, г. Харьков
Пользователь №: 19 006



Посчитать свертку сигналов друг с другом (ака ВКФ), и найти ближайший максимум - его положение = задержка, если надо, из нее фазу уже вычислить
Go to the top of the page
 
+Quote Post
DASM
сообщение Jan 9 2017, 04:47
Сообщение #15


Школьник
******

Группа: Свой
Сообщений: 3 040
Регистрация: 28-05-05
Из: Санкт-Петербург
Пользователь №: 5 493



Цитата(Artem_Petrik @ Jan 8 2017, 17:28) *
Посчитать свертку сигналов друг с другом (ака ВКФ), и найти ближайший максимум - его положение = задержка, если надо, из нее фазу уже вычислить


Все немного не так оказалось. У меня по сути своей BPSK, очень низкоскоростной, порядка 50 Гц ключевание на 90 гр. Пока что считаю Фурье по "окнам", которые скользят по сигналу. Проблемы в неопределенности частоты реальной, точно ее можно определить только на бесконечном промежутке времени sad.gif Цель конечная на самом деле не получение информации от манипуляции, а вычисление сдвига фазы между битами с максимальной точностью. Как бы по Фурье все работает, ищу максимумы амплитуды и на этой частоте фазу. Но сигнал шум не совсем хорошо исходный, около 0 Дб. Конечно буду улучшать все параметры, но хочется выжать максимум в любом случае. О, вроде же любительская связь есть похожего принципа, надо будет посмотреть


--------------------
Любитель красивых вещей.
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 Текстовая версия Сейчас: 23rd May 2017 - 12:39
Рейтинг@Mail.ru


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