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

 
 
2 страниц V   1 2 >  
Reply to this topicStart new topic
> Найти разность фаз одинаковых сигналов
Sergey Krivonoso...
сообщение Feb 7 2018, 12:47
Сообщение #1


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

Группа: Участник
Сообщений: 89
Регистрация: 10-07-08
Из: Днепропетровск
Пользователь №: 38 859



Добрый день!
Есть 128 выборок сигнала и эталона. 128 выборок является периодом. Необходимо вычислить сдвиг по фазе 2 сигналов. Как это сделать без перебора?

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


На графике пример 2 сигналов с 36 выборками. Если 36 выборок я могу обработать методом перебора (сдвиг на 1 и сравнение сигналов), то 128 занимает в 10 раз больше времени.


--------------------
Go to the top of the page
 
+Quote Post
Rst7
сообщение Feb 7 2018, 14:09
Сообщение #2


Йа моск ;)
******

Группа: Модераторы
Сообщений: 4 322
Регистрация: 7-07-05
Из: Kharkiv-city
Пользователь №: 6 610



Ключевое слово для гугля будет "вычисление корреляционной функции через преобразование Фурье".


--------------------
"Практика выше (теоретического) познания, ибо она имеет не только достоинство всеобщности, но и непосредственной действительности." - В.И. Ленин
Go to the top of the page
 
+Quote Post
novikovfb
сообщение Feb 8 2018, 05:35
Сообщение #3


Местный
***

Группа: Участник
Сообщений: 477
Регистрация: 29-09-11
Пользователь №: 67 450



Цитата(Rst7 @ Feb 7 2018, 18:09) *
Ключевое слово для гугля будет "вычисление корреляционной функции через преобразование Фурье".

Правильно ли я понял, что предлагаете посчитать корреляцию эталона (один раз) и фактического сигнала с синусом и косинусом основной гармоники и по ним считать разность фаз?
Go to the top of the page
 
+Quote Post
ViKo
сообщение Feb 8 2018, 06:46
Сообщение #4


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 251
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Преобразование фурье ведь вычисляет фазовый спектр. Вот и найти фазу первой гармоники.
Go to the top of the page
 
+Quote Post
_pv
сообщение Feb 8 2018, 07:15
Сообщение #5


Гуру
******

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



Цитата(novikovfb @ Feb 8 2018, 11:35) *
Правильно ли я понял, что предлагаете посчитать корреляцию эталона (один раз) и фактического сигнала с синусом и косинусом основной гармоники и по ним считать разность фаз?

нет,
произведению сигналов в частотной области соответствует свертка сигналов во временной области
но даже для 128 отсчётов сделать три (два на самом деле, эталон можно посчитать один раз заранее) преобразования Фурье, хоть и быстрых, в зависимости от того на чём считается, может оказаться не быстрее прямого вычисления корреляции.
через Фурье выигрыш будет для более длинных выборок.
Go to the top of the page
 
+Quote Post
novikovfb
сообщение Feb 8 2018, 07:53
Сообщение #6


Местный
***

Группа: Участник
Сообщений: 477
Регистрация: 29-09-11
Пользователь №: 67 450



Цитата(_pv @ Feb 8 2018, 11:15) *
нет,
произведению сигналов в частотной области соответствует свертка сигналов во временной области
но даже для 128 отсчётов сделать три (два на самом деле, эталон можно посчитать один раз заранее) преобразования Фурье, хоть и быстрых, в зависимости от того на чём считается, может оказаться не быстрее прямого вычисления корреляции.
через Фурье выигрыш будет для более длинных выборок.

зачем полное Фурье? Достаточно посчитать корреляцию два раза (по синусу и по косинусу).
Go to the top of the page
 
+Quote Post
_pv
сообщение Feb 8 2018, 08:38
Сообщение #7


Гуру
******

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



Цитата(novikovfb @ Feb 8 2018, 14:53) *
зачем полное Фурье? Достаточно посчитать корреляцию два раза (по синусу и по косинусу).

какой частоты?
Go to the top of the page
 
+Quote Post
blackfin
сообщение Feb 8 2018, 08:47
Сообщение #8


Гуру
******

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



Цитата(_pv @ Feb 8 2018, 11:38) *
какой частоты?

Частоты основной гармоники сигнала и основной гармоники эталона. Они ведь равны, эти частоты, если верить графику из первого поста?..

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

PPS. И соглашусь, выигрыша по количеству вычислений этот метод не даст..
Go to the top of the page
 
+Quote Post
Sergey Krivonoso...
сообщение Feb 8 2018, 08:50
Сообщение #9


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

Группа: Участник
Сообщений: 89
Регистрация: 10-07-08
Из: Днепропетровск
Пользователь №: 38 859



Цитата(novikovfb @ Feb 8 2018, 07:35) *
Правильно ли я понял, что предлагаете посчитать корреляцию эталона (один раз) и фактического сигнала с синусом и косинусом основной гармоники и по ним считать разность фаз?

Если у меня период равен 128 выборкам, то период основной гармоники тоже равен 128? Предложение понравилось, это ускорит вычисления, только будет ли это работать?


--------------------
Go to the top of the page
 
+Quote Post
blackfin
сообщение Feb 8 2018, 09:17
Сообщение #10


Гуру
******

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



Цитата(Sergey Krivonosov @ Feb 8 2018, 11:50) *
Предложение понравилось, это ускорит вычисления, только будет ли это работать?

Работать будет, если научитесь быстро вычислять арктангенс и делить комплексные числа. biggrin.gif
Go to the top of the page
 
+Quote Post
novikovfb
сообщение Feb 8 2018, 09:36
Сообщение #11


Местный
***

Группа: Участник
Сообщений: 477
Регистрация: 29-09-11
Пользователь №: 67 450



Цитата(blackfin @ Feb 8 2018, 12:47) *
PPS. И соглашусь, выигрыша по количеству вычислений этот метод не даст..

По сравнению с идеей 128 раз сравнивать эталон с реальным сигналом (см. первое сообщение)?
Go to the top of the page
 
+Quote Post
blackfin
сообщение Feb 8 2018, 09:52
Сообщение #12


Гуру
******

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



Цитата(novikovfb @ Feb 8 2018, 12:36) *
По сравнению с идеей 128 раз сравнивать эталон с реальным сигналом (см. первое сообщение)?

Судя по графику из первого поста, нет необходимости делать сравнение все 128 раз.

Во-первых, можно взять лишь один период эталонного сигнала и искать совпадение с этим укороченным эталоном и на укороченном интервале.
Это, ессно, в том случае, если на интервале в 128 точек укладывается несколько периодов эталонного сигнала.. /См. рис./

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

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

Зачем искать сигнал там, где заранее известно, что его нет? biggrin.gif
Go to the top of the page
 
+Quote Post
Sergey Krivonoso...
сообщение Feb 8 2018, 12:32
Сообщение #13


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

Группа: Участник
Сообщений: 89
Регистрация: 10-07-08
Из: Днепропетровск
Пользователь №: 38 859



Цитата(blackfin @ Feb 8 2018, 11:52) *
Во-первых, можно взять лишь один период эталонного сигнала и искать совпадение с этим укороченным эталоном и на укороченном интервале.
Это, ессно, в том случае, если на интервале в 128 точек укладывается несколько периодов эталонного сигнала.. /См. рис./

Вы наверное не поняли задачу. Все получаемые сигналы и все эталоны имеют один и тот же период 128 выборок. Если сигнал и эталон совпадают, они могут отличаются только сдвигом на n количество выборок (0-127).


--------------------
Go to the top of the page
 
+Quote Post
blackfin
сообщение Feb 8 2018, 12:50
Сообщение #14


Гуру
******

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



Цитата(Sergey Krivonosov @ Feb 8 2018, 15:32) *
Вы наверное не поняли задачу. Все получаемые сигналы и все эталоны имеют один и тот же период 128 выборок. Если сигнал и эталон совпадают, они могут отличаются только сдвигом на n количество выборок (0-127).

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

Все, конечно, зависит от уровня шума и спектров обоих сигналов.
Go to the top of the page
 
+Quote Post
Sergey Krivonoso...
сообщение Feb 8 2018, 12:56
Сообщение #15


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

Группа: Участник
Сообщений: 89
Регистрация: 10-07-08
Из: Днепропетровск
Пользователь №: 38 859



Цитата(blackfin @ Feb 8 2018, 14:50) *
Даже если и так, никто не мешает сделать сначала грубый поиск по 16-ти точкам отстоящим друг от друга на 8 выборок, а потом уточнить сдвиг между сигналами вычисляя корреляцию уже для всех 128-ми точек и для сдвигов на максимум 8 точек влево и вправо относительно найденной грубой оценки искомой точки.

Все, конечно, зависит от уровня шума и спектров обоих сигналов.

Это хорошее предложение, но проблема в том что в сигнале может быть шум и по такому методу можно пропустить нужный сигнал.


--------------------
Go to the top of the page
 
+Quote Post

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

 


RSS Текстовая версия Сейчас: 22nd May 2018 - 05:03
Рейтинг@Mail.ru


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