Определение задержки Опции

[ivan219]

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

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

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

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

Эскизы прикрепленных изображений

 

[ivan219]

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

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

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

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

[bahurin]

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

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

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

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

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

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

[ivan219]

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

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

Я делаю так, беру 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

Эскизы прикрепленных изображений

 

[bahurin]

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

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

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

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

[ivan219]

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

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

Эскизы прикрепленных изображений