Есть опорный ЛЧМ сигнал, и есть исследуемый ЛЧМ сигнал.
Исследуемый сигнал имеет какую то задержку в семплах.

Нужно определелить эту задержку?

Ниже приведены несколько рисунков.
В графике ФЧХ отображается разница фаз исследуемого сигнала и опорного.

Сообщение отредактировал ivan219 - Jul 22 2011, 18:39

1. Можно посчитать взаимную корреляционную функцию опорного и сдвинутого сигналов и ее максимум будет смещен ровно на величину задержки
2. В ЛЧМ сигнале временная задержка соответствует частотному рассогласованию опорного и принятого сигналов. Вам надо раскрыть неоднозначность по фазе в ваших графиках и угол наклона даст временную задержку (при условии, что нет частотных сдвигов)

Можно по подробнее.

Да будут частотные сдвиги. Это для примера взят чистый сигнал.

Поэкспериментировал с корреляцией, и вот что получилось, смотрите рисунок.
Получается ошибка в 2 семпла.
При прохождении задержки от 0 до n-1 ошибка постепенно нарастает сначала 1 семпл, потом 2 бывает и больше потом также уходит до 0
Как с этим бороться???

Сообщение отредактировал ivan219 - Jul 22 2011, 21:31

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

R = abs( ifft( fft(s) * conj(fft(s0)) ) )




В этом случае частотный сдвиг принятого и опорного сигналов будет выглядеть как временная задержка. Почитайте книгу Кук Бернфельд Радиолокационные сигналы там про лчм очень подробно расписано.

PS ваши картинки взаимной корреляции какие-то неубидительные, возможно есть ошибка в модели.

Я делаю так, беру R = abs(IFFT(FFT(S) * FFT(S0))) (S0 повёрнут во времени во временной области 0 отсчёт идёт на место отсчёта N-1, а отсчёт N-1 идёт на место 0) размер блока FFT в 2 раза больше суммы выборок S + S0

Если не правильно объясните что означает функция conj()???

Вроде сейчас правильно определяется положение задержки.

Сообщение отредактировал ivan219 - Jul 23 2011, 12:05

Conj() это комплексно сопряженное число. Т.е. взяв бпф от опоры надо преобразовать ещё его в комплексно сопряженное число, чтобы получить отклик как функция неопределенности данного сигнала.

Что-то такое и предполагал.
А как это сделать?
Я предполагаю, что нужно умножить на j = Sqrt(-1) но как это сделать???
У меня есть вещественная часть Re и мнимая Im

Разобрался комплексно сопряженный спектр это спектр, умноженный на мнимую единицу j = Sqrt(-1) т.е. получается если брать вещественную и мнимую часть спектра то происходит:

происходит поворот фазы на 90 град. всех спектральных составляющих.
С корреляцией то же разобрался. Теперь всё точно вычисляет
Всем спасибо.

Да нет же, просто знак у мнимой части поменять надо пруф

Хоть один здравомыслящий человек нашелся... Похоже, еще не все потеряно.

Спасибо. Теперь точно разобрался

если вы развернули опорный сигнал, то операцию комплексного сопряжения conj делать не надо R = abs(IFFT(FFT(S) * FFT(S0))) в этом случае дает правильный результат, и картинки тоже стали более похожи на правду.

Значит, комплексное сопряжение, и поворот во времени это равноценные операции или всё-таки есть какая-то разница?
На картинке корреляция с комплексно сопряженным сигналом.

Выяснилась маленькая не приятность с фазовым сдвигом частотных компонентов если он больше длительности половины семпла Fi = 360 * F * Deley / N; 0.5 <= Deley < 1 то алгоритм определяет сдвиг на один семпл. И получается ошибочный результат.
Есть ли методы коррекции или другие алгоритмы???

Сообщение отредактировал ivan219 - Jul 27 2011, 08:09