Сжатие импульсов

редактировать

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

Содержание
  • 1 Простой импульс
    • 1.1 Описание сигнала
    • 1.2 Разрешение по дальности
    • 1.3 Энергия, необходимая для передачи этого сигнала
  • 2 Сжатие импульсов с помощью линейной частотной модуляции (или чирпирования)
    • 2.1 Основные принципы
    • 2.2 Взаимная корреляция между переданным и принятым сигналами
    • 2.3 Улучшение отношения сигнал / шум за счет сжатия импульсов
    • 2.4 Обработка растяжения
    • 2.5 Форма волны со ступенчатым изменением частоты
  • 3 Сжатие импульсов путем фазового кодирования
  • 4 Примечания
  • 5 Дополнительная литература
  • 6 См. Также
Простой импульс

Описание сигнала

Простейший сигнал, который может передать импульсный радар, - это импульс синусоидальной амплитуды, A {\ displaystyle A}A и несущая частота, f 0 {\ displaystyle f_ {0}}f_ {0} , усеченная прямоугольной функцией ширины, T {\ displaystyle \ scriptstyle T }\ scriptstyle T . Импульс передается периодически, но это не основная тема данной статьи; мы будем рассматривать только один импульс, s {\ displaystyle s}s . Если мы предположим, что импульс начинается в момент времени t = 0 {\ displaystyle t \, = \, 0}t \, = \, 0 , сигнал можно записать следующим образом, используя комплекс обозначение:

s (t) = {A e 2 i π f 0 t if 0 ≤ t < T 0 otherwise {\displaystyle s(t)=\left\{{\begin{array}{ll}Ae^{2i\pi f_{0}t}{\text{if}}\;0\leq ts (t) = \ left \ {{\ begin {array} {ll} Ae ^ {{2i \ pi f_ {0} t}} {\ text {if}} \; 0 \ leq t <T \\ 0 {\ text {иначе}} \ end {array}} \ right.

Разрешающая способность по дальности

Определим разрешающую способность по дальности, которая может быть получена с таким сигналом. Ответный сигнал, записанный r (t) {\ displaystyle \ scriptstyle r (t)}\ scriptstyle r (t) , представляет собой ослабленную и сдвинутую во времени копию исходного переданного сигнала (в действительности Доплеровский Эффект тоже может играть роль, но здесь это не важно.) Во входящем сигнале также присутствует шум, как в воображаемом, так и в реальном канале, который мы будем считать белым и Гауссовский (это обычно верно); мы пишем B (t) {\ displaystyle \ scriptstyle B (t)}\ scriptstyle B (t) , чтобы обозначить этот шум. Для обнаружения входящего сигнала обычно используется согласованная фильтрация. Этот метод является оптимальным, когда известный сигнал должен быть обнаружен среди аддитивного белого гауссовского шума.

Другими словами, вычисляется взаимная корреляция принятого сигнала с переданным сигналом. Это достигается посредством свертки входящего сигнала с сопряженной и обращенной во времени версией переданного сигнала. Эта операция может быть выполнена программно или аппаратно. Мы пишем < s, r>(t) {\ displaystyle \ scriptstyle (t)}\scriptstyle <s,\,r>(t) для этой взаимной корреляции. У нас есть:

⟨s, r⟩ (t) = ∫ t ′ = 0 + ∞ s ⋆ (t ′) р (t + t ′) dt ′ {\ displaystyle \ langle s, r \ rangle (t) = \ int _ {t '\, = \, 0} ^ {+ \ infty} s ^ {\ star} (t ') r (t + t') dt '}\langle s,r\rangle (t)=\int _{{t'\,=\,0}}^{{+\infty }}s^{\star }(t')r(t+t')dt'

Если отраженный сигнал возвращается к приемнику в момент tr {\ displaystyle \ scriptstyle t_ {r}}\ scriptstyle t_ {r} и ослабляется на коэффициент K {\ displaystyle \ scriptstyle K}\ scriptstyle K , это дает:

r (t) = {KA e 2 i π f 0 (t - tr) + B (t), если tr ≤ t < t r + T B ( t) otherwise {\displaystyle r(t)=\left\{{\begin{array}{ll}KAe^{2i\pi f_{0}(t\,-\,t_{r})}+B(t){\t_dv{if}}\;t_{r}\leq tr (t) = \ left \ {{\ begin {a rray} {ll} KAe ^ {{2i \ pi f_ {0} (t \, - \, t_ {r})}} + B (t) {\ t_dv {if}} \; t_ {r} \ leq t <t_ {r} + T \\ B (t) {\ t_dv {иначе}} \ end {array}} \ right.

Поскольку мы знаем передаваемый сигнал, получаем:

⟨s, r⟩ (t) = KA 2 Λ (t - tr T) e 2 i π f 0 ( t - тр) + В '(t) {\ displaystyle \ langle s, r \ rangle (t) = KA ^ {2} \ Lambda \ left ({\ frac {t-t_ {r}} {T}} \ справа) е ^ {2i \ pi f_ {0} (t \, - \, t_ {r})} + B '(t)}\langle s,r\rangle (t)=KA^{2}\Lambda \left({\frac {t-t_{r}}{T}}\right)e^{{2i\pi f_{0}(t\,-\,t_{r})}}+B'(t)

где B' (t) {\ displaystyle \ scriptsty le B '(t)}\scriptstyle B'(t), является результатом взаимной корреляции между шумом и передаваемым сигналом. Функция Λ {\ displaystyle \ Lambda}\ Lambda - функция треугольника, ее значение равно 0 на [- ∞, - 1 2] ∪ [1 2, + ∞] {\ displaystyle \ стиль сценария [- \ infty, \, - {\ frac {1} {2}}] \, \ cup \, [{\ frac {1} {2}}, \, + \ infty]}\ scriptstyle [- \ infty, \, - {\ frac {1} {2}}] \, \ cup \, [{\ frac {1} {2} }, \, + \ infty] , он увеличивается линейно на [- 1 2, 0] {\ displaystyle \ scriptstyle [- {\ frac {1} {2}}, \, 0]}\ scriptstyle [- {\ frac {1} {2}}, \, 0] , где достигает своего максимума 1, и оно линейно уменьшается на [0, 1 2] {\ displaystyle \ scriptstyle [0, \, {\ frac {1} {2}}]}\ scriptstyle [0, \, {\ frac {1} {2}}] , пока снова не достигнет 0. На рисунках в конце этого абзаца показана форма взаимной корреляции для выборочного сигнала (красным), в данном случае реального усеченного синуса, длительностью T = 1 {\ displaystyle \ scriptstyle T \, = \, 1 }\ scriptstyle T \, = \, 1 секунд единичной амплитуды и частоты f 0 = 10 {\ displaystyle \ scriptstyle f_ {0} \, = \, 10}\ scriptstyle f_ {0} \, = \, 10 герц. Два эхо-сигнала (отмечены синим цветом) возвращаются с задержкой 3 и 5 секунд и амплитудами, равными 0,5 и 0,3 амплитуды переданного импульса, соответственно; это просто случайные значения для примера. Поскольку сигнал реальный, взаимная корреляция взвешивается дополнительным коэффициентом ⁄ 2.

Если два импульса возвращаются (почти) одновременно, взаимная корреляция равна сумме взаимных корреляций двух элементарных сигналов. Чтобы отличить одну «треугольную» огибающую от огибающей другого импульса, четко видно, что время прихода двух импульсов должно быть разделено по крайней мере T {\ displaystyle \ scriptstyle T}\ scriptstyle T так что максимумы обоих импульсов могут быть разделены. Если это условие не выполняется, оба треугольника смешаются, и их невозможно будет разделить.

Поскольку расстояние, пройденное волной в течение T {\ displaystyle \ scriptstyle T}\ scriptstyle T , равно c T {\ displaystyle \ scriptstyle cT}\ scriptstyle cT (где c - скорость волны в среде), и поскольку это расстояние соответствует времени обхода, мы получаем:

Результат 1
Разрешение диапазона с синусоидальным импульсом составляет 1 2 c T {\ displaystyle \ scriptstyle {\ frac {1} {2}} cT}\ scriptstyle {\ frac {1} {2}} cT , где T {\ displaystyle \ scriptstyle T}\ scriptstyle T - длительность импульса, а c {\ displaystyle \ scriptstyle c}\ scriptstyle c , скорость волны.

Вывод: для увеличения разрешения необходимо уменьшить длину импульса.

Пример (простой импульс): переданный сигнал отображается красным цветом (несущая 10 Гц, амплитуда 1, длительность 1 секунда) и два эхо-сигнала (синий).
До согласованной фильтрацииПосле согласованной фильтрации
Если цели достаточно разнесены... ... можно различить эхо.
Если цели слишком близко... ... эхо смешивается.

Требуемая энергия для передачи этого сигнала

Мгновенная мощность переданного импульса составляет P (t) = | с | 2 (т) {\ displaystyle \ scriptstyle P (t) \, = \, | s | ^ {2} (t)}\ scriptstyle P (t) \, = \, | s | ^ {2} (t) . Энергия, вложенная в этот сигнал, равна:

E = ∫ 0 TP (t) dt = A 2 T {\ displaystyle E = \ int _ {0} ^ {T} P (t) dt = A ^ {2} T}E = \ int _ {0} ^ {T } P (t) dt = A ^ {2} T

Аналогично, энергия принятого импульса равна E r = K 2 A 2 T {\ displaystyle \ scriptstyle E_ {r} \, = \, K ^ {2} A ^ {2} T }\ scriptstyle E_ {r} \, = \, K ^ {2} A ^ {2} T . Если σ {\ displaystyle \ scriptstyle \ sigma}\ scriptstyle \ sigma - стандартное отклонение шума, отношение сигнал / шум (SNR) в приемнике будет:

SNR = E r σ 2 знак равно К 2 A 2 T σ 2 {\ Displaystyle SNR = {\ frac {E_ {r}} {\ sigma ^ {2}}} = {\ frac {K ^ {2} A ^ {2} T} { \ sigma ^ {2}}}}SNR = {\ frac {E_ {r}} {\ sigma ^ {{2}}}} = {\ frac {K ^ {2} A ^ {2} T} {\ sigma ^ {{2}}}}

SNR пропорционально длительности импульса T {\ displaystyle T}T , если другие параметры остаются постоянными. Это вводит компромисс: увеличение T {\ displaystyle T}T улучшает SNR, но снижает разрешение, и наоборот.

Сжатие импульсов с помощью линейной частотной модуляции (или чирпирования)

Основные принципы

Как можно получить достаточно большой импульс (чтобы по-прежнему иметь хорошее SNR в приемнике) без плохое разрешение? Здесь на сцену выходит сжатие импульсов. Основной принцип заключается в следующем:

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

В приложениях радар или сонар линейный щебетание - наиболее часто используемые сигналы для сжатия импульсов. Поскольку импульс имеет конечную длину, его амплитуда представляет собой функцию прямоугольника . Если переданный сигнал имеет продолжительность T {\ displaystyle \ scriptstyle T}\ scriptstyle T , начинается с t = 0 {\ displaystyle \ scriptstyle t \, = \, 0}\ scriptstyle t \, = \, 0 и линейно перемещает полосу частот Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f с центром на несущей f 0 {\ displaystyle \ scriptstyle f_ {0}}\ scriptstyle f_0 , это можно записать:

sc (t) = {A ei 2 π ((f 0 - Δ f 2) t + Δ f 2 T t 2) if 0 ≤ t < T 0 otherwise {\displaystyle s_{c}(t)=\left\{{\begin{array}{ll}Ae^{i2\pi \left(\left(f_{0}\,-\,{\frac {\Delta f}{2}}\right)t\,+\,{\frac {\Delta f}{2T}}t^{2}\,\right)}{\t_dv{if}}\;0\leq ts_ {c} (t) = \ left \ {{\ begin {array} {ll} Ae ^ {{i2 \ pi \ left (\ left (f_ {0} \, - \, {\ frac {\ Delta f} {2}} \ right) t \, + \, {\ frac {\ Delta f} {2T}} t ^ {2} \, \ right)}} {\ t_dv {if}} \; 0 \ leq t <T \\ 0 {\ t_dv {иначе}} \ end {array}} \ right.

Определение chirp, приведенное выше означает, что фаза чирпированного сигнала (то есть аргумент комплексной экспоненты) является квадратичной:

ϕ (t) = 2 π ((f 0 - Δ f 2) t + Δ f 2 T t 2) {\ displaystyle \ phi (t) = 2 \ pi \ left (\ left (f_ {0} \, - \, {\ frac {\ Delta f} {2}} \ right) t \, + \, {\ frac {\ Delta f} {2T}} t ^ {2} \, \ right)}\ phi (t) = 2 \ pi \ left (\ left (f_ {0} \, - \, {\ frac {\ Delta f} {2}} \ right) t \, + \, {\ frac {\ Delta f} {2T}} t ^ {2} \, \ right)

таким образом, мгновенная частота равна (по определению):

f (t) = 1 2 π [d ϕ dt] T знак равно f 0 - Δ е 2 + Δ е T T {\ displaystyle f (t) = {\ frac {1} {2 \ pi}} \ left [{\ frac {d \ phi} {dt}} \ right] _ {t} = f_ {0} - {\ frac {\ Delta f} {2}} + {\ frac {\ Delta f} {T}} t}f (t) = {\ frac {1} {2 \ pi}} \ left [{\ frac {d \ phi} {dt}} \ right] _ {t} = f_ {0} - {\ frac {\ Delta f} {2}} + {\ frac {\ Delta f} {T}} t

whi ch - предполагаемый линейный наклон, идущий от f 0 - Δ f 2 {\ displaystyle \ scriptstyle f_ {0} \, - \, {\ frac {\ Delta f} {2}}}\ scriptstyle f_ {0} \, - \, {\ frac {\ Delta f} {2}} от t = 0 {\ displaystyle \ scriptstyle t \, = \, 0}\ scriptstyle t \, = \, 0 до f 0 + Δ f 2 {\ displaystyle \ scriptstyle f_ {0} \, + \, {\ frac {\ Delta f} {2}}}\ scriptstyle f_ {0} \, + \, {\ frac {\ Delta f} {2}} при t = T {\ displaystyle \ scriptstyle t \, = \, T}\ scriptstyle t \, = \, T .

Отношение фазы к частоте часто используется в другом направлении, начиная с желаемого f (t) {\ displaystyle f (t)}f (t) и записывая фазу щебета через интегрирование частоты:

ϕ (t) = 2 π ∫ 0 tf (u) du {\ displaystyle \ phi (t) = 2 \ pi \ int _ {0} ^ {t} f (u) \, du}\ phi (t) = 2 \ pi \ int _ {0} ^ {t} f (u) \, du

Взаимная корреляция между переданным и принятый сигнал

Что касается «простого» импульса, давайте вычислим взаимную корреляцию между переданным и принятым сигналом. Для упрощения будем считать, что щебетание записывается не так, как указано выше, а в этой альтернативной форме (окончательный результат будет таким же):

sc ′ (t) = {A e 2 i π ( f 0 + Δ f 2 T t) t if - T 2 ≤ t < T 2 0 otherwise {\displaystyle s_{c'}(t)=\left\{{\begin{array}{ll}Ae^{2i\pi \left(f_{0}\,+\,{\frac {\Delta f}{2T}}t\right)t}{\t_dv{if}}\;-{\frac {T}{2}}\leq t<{\frac {T}{2}}\\0{\t_dv{otherwise}}\end{array}}\right.}s_{{c'}}(t)=\left\{{\begin{array}{ll}Ae^{{2i\pi \left(f_{0}\,+\,{\frac {\Delta f}{2T}}t\right)t}}{\t_dv{if}}\;-{\frac {T}{2}}\leq t<{\frac {T}{2}}\\0{\t_dv{otherwise}}\end{array}}\right.

Поскольку эта взаимная корреляция равна (за исключением коэффициента затухания K {\ displaystyle K}К ), чтобы автокорреляционная функция sc ′ {\ displaystyle \ scriptstyle s_ {c '}}\scriptstyle s_{{c'}}, это то, что мы рассматриваем:

⟨sc ′, sc ′⟩ (t) = ∫ - ∞ + ∞ sc ′ ⋆ (- t ′) sc ′ (t - t ′) dt ′ {\ displaystyle \ langle s_ {c '}, s_ {c'} \ rangle (t) = \ int _ {- \ infty} ^ {+ \ infty} s_ {c '} ^ {\ star} (- t') s_ {c '} (t-t') dt '}\langle s_{{c'}},s_{{c'}}\rangle (t)=\int _{{-\infty }}^{{+\infty }}s_{{c'}}^{\star }(-t')s_{{c'}}(t-t')dt'

Можно показать, что функция автокорреляции sc ′ {\ displaystyle s_ {c '}}s_{{c'}}равно:

⟨sc ′, sc ′⟩ (t) = A 2 T Λ (t T) sinc [Δ ft Λ (t T) ] е 2 я π е 0 T {\ Displaystyle \ langle s_ {c '}, s_ {c'} \ rangle (t) = A ^ {2} T \ Lambda \ left ({\ frac {t} {T} } \ right) \ mathrm {sinc} \ left [\ Delta ft \ Lambda \ left ({\ frac {t} {T}} \ right) \ ri ght] e ^ {2i \ pi f_ {0} t}}{\displaystyle \langle s_{c'},s_{c'}\rangle (t)=A^{2}T\Lambda \left({\frac {t}{T}}\right)\mathrm {sinc} \left[\Delta ft\Lambda \left({\frac {t}{T}}\right)\right]e^{2i\pi f_{0}t}}

Максимум функции автокорреляции sc ′ {\ displaystyle \ scriptstyle s_ {c '}}\scriptstyle s_{{c'}}достигается при 0. Около 0 эта функция ведет себя как член sinc (или кардинальный синус), определяемый здесь как sinc (x) = sin (π x) / (π x) {\ displaystyle sinc ( х) = грех (\ пи х) / (\ пи х)}{\ displaystyle sinc (x) = sin (\ pi x) / (\ pi x)} . Временная ширина этого кардинального синуса -3 дБ более или менее равна T '= 1 Δ f {\ displaystyle \ scriptstyle T' \, = \, {\ frac {1} {\ Delta f}}}\scriptstyle T'\,=\,{\frac {1}{\Delta f}}. Все происходит так, как если бы после согласованной фильтрации у нас было разрешение, которое было бы достигнуто с помощью простого импульса длительностью T ′ {\ displaystyle \ scriptstyle T '}\scriptstyle T'. Для общих значений Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f , T ′ {\ displaystyle \ scriptstyle T '}\scriptstyle T'меньше, чем T {\ displaystyle \ scriptstyle T}\ scriptstyle T , отсюда и название сжатия импульса.

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

Результат 2
Разрешение по расстоянию, достижимое с помощью линейной частотной модуляции импульса в полосе пропускания Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f , составляет: c 2 Δ f {\ displaystyle \ scriptstyle {\ frac {c} {2 \ Delta f}}}\ scriptstyle {\ frac {c} {2 \ Delta f}} , где c {\ displaystyle \ scriptstyle c}\ scriptstyle c - скорость волна.
Определение
Соотношение TT ′ = T Δ f {\ displaystyle \ scriptstyle {\ frac {T} {T ^ {\ prime}}} \, = \, T \ Delta f}\ scriptstyle {\ frac {T} {T ^ {\ prime}}} \, = \, T \ Delta f - степень сжатия импульса. Обычно он больше 1 (обычно от 20 до 30).
Пример (чирпированный импульс): переданный сигнал выделен красным цветом (несущая 10 Гц, модуляция 16 Гц, амплитуда 1, длительность 1 секунда) и два эхо-сигнала (синим цветом).
До согласованной фильтрации После согласованной фильтрации: эхо-сигналы короче по времени.

Улучшение отношения сигнал / шум посредством сжатия импульсов

Энергия сигнала не изменяется во время сжатия импульса. Однако теперь он находится в главном лепестке кардинального синуса, ширина которого составляет приблизительно T ′ ≈ 1 Δ f {\ displaystyle \ scriptstyle T '\, \ приблизительно \, {\ frac {1} {\ Delta f}}}\scriptstyle T'\,\approx \,{\frac {1}{\Delta f}}. Если P {\ displaystyle \ scriptstyle P}\ scriptstyle P - мощность сигнала до сжатия, а P ′ {\ displaystyle \ scriptstyle P '}\scriptstyle P'- мощность сигнал после сжатия, мы имеем:

P × T = P ′ × T ′ {\ displaystyle P \ times T = P '\ times T'}P\times T=P'\times T'

, что дает:

P ′ = P × TT ′ {\ Displaystyle P '= P \ times {\ frac {T} {T'}}}P'=P\times {\frac {T}{T'}}

Как следствие:

Результат 3
После сжатия импульса можно учитывать мощность принятого сигнала как усиленный T Δ f {\ displaystyle \ scriptstyle T \ Delta f}\ scriptstyle T \ Delta f . Это дополнительное усиление можно ввести в уравнение радара.
Пример: те же сигналы, что и выше, плюс аддитивный гауссовский белый шум (σ = 0,5 {\ displaystyle \ scriptstyle \ sigma \, = \, 0,5}\ scriptstyle \ sigma \, = \, 0,5 )
До согласованной фильтрации: сигнал скрыт в шуме После согласованной фильтрации: эхо-сигналы становятся видимыми.

Обработка растяжения

В то же время сжатие импульсов может обеспечить хорошее соотношение сигнал / шум и точное разрешение по диапазону одновременно цифровая обработка сигнала в такой системе может быть затруднена из-за высокой мгновенной ширины полосы сигнала (Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f может составлять сотни мегагерц или даже превышают 1 ГГц.) Обработка растяжения - это метод согласованной фильтрации широкополосной формы колебательного сигнала, который подходит для приложений, которым требуется очень точное разрешение диапазона в относительно коротких интервалах.

Обработка растяжения

На рисунке выше показан сценарий анализа обработки растяжения. Центральная контрольная точка (CRP) находится в середина интересующего окна диапазона в диапазоне R 0 {\ displaystyle \ scriptstyle R_ {0}}\ scriptstyle R_0 , что соответствует временной задержке t 0 {\ displaystyle \ scriptstyle t_ {0}}{\ displaystyle \ scriptstyle t_ {0}} .

Если переданный сигнал является формой сигнала ЛЧМ:

x (t) = exp ⁡ (j π Δ f T (t) 2) exp ⁡ (j 2 π f 0 (t)), 0 ≤ T ≤ T {\ displaystyle x (t) = \ exp \ left (j \ pi {\ frac {\ Delta f} {T}} (t) ^ {2} \ right) \ exp (j2 \ pi f_ {0} (t)), 0 \ leq t \ leq T}{\ displaystyle x (t) = \ exp \ left (j \ pi {\ frac {\ Delta f} {T}} (t) ^ { 2} \ right) \ exp (j2 \ pi f_ {0} (t)), 0 \ leq t \ leq T}

, то эхо-сигнал от цели на расстоянии R b {\ displaystyle \ scriptstyle R_ {b}}{\ displaystyle \ scriptstyle R_ {b}} может быть выражается как:

x ¯ (t) = ρ exp ⁡ (j π Δ f T (t - tb) 2) exp ⁡ (j 2 π f 0 (t - tb)), 0 ≤ t - tb ≤ T {\ displaystyle {\ bar {x}} (t) = \ rho \ exp \ left (j \ pi {\ frac {\ Delta f} {T}} (t-t_ {b}) ^ {2} \ right) \ ехр (j2 \ pi f_ {0} (t-t_ {b})), 0 \ leq t-t_ {b} \ leq T}{\ displaystyle {\ bar {x}} (t) = \ rho \ exp \ left (j \ pi {\ frac {\ Delta f}) {T}} (t-t_ {b}) ^ {2} \ right) \ exp (j2 \ pi f_ {0} (t-t_ {b})), 0 \ leq t-t_ {b} \ leq T}

где ρ {\ displaystyle \ scriptstyle \ rho}\ scriptstyle \ rho пропорционален отражательной способности рассеивателя. Затем мы умножаем эхо на exp ⁡ (- j 2 π f 0 t) exp ⁡ (- j π Δ f T (t - t 0) 2) {\ displaystyle \ scriptstyle \ exp (-j2 \ pi f_ {0} t) \ exp \ left (-j \ pi {\ frac {\ Delta f} {T}} (t-t_ {0}) ^ {2} \ right)}{\ displaystyle \ scriptstyle \ exp (-j2 \ pi f_ {0} t) \ exp \ left (-j \ pi {\ frac {\ Delta f} {T}} ( t-t_ {0}) ^ {2} \ right)} и эхо станет следующим:

y (t) = ρ exp ⁡ (- j 4 π R b λ) exp ⁡ (- j 2 π Δ f T δ tb (t - t 0)) exp ⁡ (j π Δ f T (δ tb) 2), t 0 ≤ t - δ tb ≤ t 0 + T {\ displaystyle y (t) = \ rho \ exp \ left (-j {\ frac {4 \ pi R_ {b}} { \ lambda}} \ right) \ exp \ left (-j2 \ pi {\ frac {\ Delta f} {T}} \ delta t_ {b} (t-t_ {0}) \ right) \ exp \ left ( j \ pi {\ frac {\ Delta f} {T}} (\ delta t_ {b}) ^ {2} \ right), t_ {0} \ leq t- \ delta t_ {b} \ leq t_ {0 } + T}{\ displaystyle y (t) = \ rho \ exp \ left (-j {\ frac {4 \ pi R_ {b}} {\ lambda}} \ right) \ exp \ left (-j2 \ pi {\ frac {\ Delta f} {T}} \ delta t_ {b} (t-t_ {0}) \ right) \ exp \ left (j \ pi {\ frac {\ Delta f} {T}} (\ delta t_ {b}) ^ {2} \ right), t_ {0} \ leq t- \ delta t_ {b} \ leq t_ {0} + T}

где λ {\ displaystyle \ scriptstyle \ lambda}\ scriptstyle \ lambda - длина волны электромагнитной волны в воздухе.

После проведения выборки и дискретного преобразования Фурье по y (t) частота синусоиды F b {\ displaystyle \ scriptstyle F_ {b}}{\ displaystyle \ scriptstyle F_ {b}} может быть решена:

F b = - δ tb Δ f T (H z) {\ displaystyle F_ {b} = - \ delta t_ {b} {\ frac {\ Delta f} {T}} (Гц)}{\ displaystyle F_ {b} = - \ delta t_ {b} {\ frac {\ Delta f} {T}} (Гц)}

и диапазон разности δ R b {\ displaystyle \ scriptstyle \ delta R_ {b}}{\ displaystyle \ scriptstyle \ delta R_ {b}} можно получить:

δ R b = - c TF b 2 Δ f {\ displaystyle \ delta R_ {b } = - {\ frac {cTF_ {b}} {2 \ Delta f}}}{\ displaystyle \ delta R_ {b} = - { \ frac {cTF_ {b}} {2 \ Delta f}}}

Чтобы показать, что полоса пропускания y (t) меньше исходной полосы пропускания сигнала Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f , мы предполагаем, что окно диапазона равно R w = c T w 2 {\ displaystyle \ scriptstyle R_ {w} = {\ frac {cT_ {w}} {2} }}{\ displaystyle \ scriptstyle R_ {w} = {\ frac {cT_ {w}} {2}}} длинный. Если цель находится на нижней границе окна диапазона, эхо будет достигнуто t 0 - T w / 2 {\ displaystyle \ scriptstyle t_ {0} -T_ {w} / 2}{\ displaystyle \ scriptstyle t_ {0} -T_ {w} / 2} секунды после передачи; аналогично, если цель находится на верхней границе окна диапазона, эхо будет достигнуто t 0 + T w / 2 {\ displaystyle \ scriptstyle t_ {0} + T_ {w} / 2}{\ displaystyle \ scriptstyle t_ {0} + T_ {w} / 2} секунд после передачи. Дифференциальное время прибытия δ tb {\ displaystyle \ scriptstyle \ delta t_ {b}}{\ displaystyle \ scriptstyle \ delta t_ {b}} для каждого случая составляет - T w / 2 {\ displaystyle \ scriptstyle -T_ {w} / 2}{\ displaystyle \ scriptstyle -T_ {w} / 2} и T w / 2 {\ displaystyle \ scriptstyle T_ {w} / 2}{\ displaystyle \ scriptstyle T_ {w} / 2} соответственно.

Затем мы можем получить полосу пропускания, учитывая разницу в частоте синусоиды для целей на нижней и верхней границе окна дальности:

Δ fs = F b, ближний - F b, дальний = - Δ е T (- T w / 2 - T w / 2) знак равно T w T Δ f {\ displaystyle \ Delta f_ {s} = F_ {b, near} -F_ {b, far} = - {\ frac {\ Дельта f} {T}} (- T_ {w} / 2-T_ {w} / 2) = {\ frac {T_ {w}} {T}} \ Delta f}{\ displaystyle \ Delta f_ {s} = F_ {b, near} -F_ {b, far} = - {\ frac {\ Delta f} {T}} (- T_ {w} / 2-T_ {w} / 2) = {\ frac {T_ {w}} {T}} \ Delta f} 

Как следствие:

Результат 4
Благодаря обработке растяжения полоса пропускания на выходе приемника меньше, чем исходная полоса пропускания сигнала, если T w < T {\displaystyle \scriptstyle T_{w}{\ displaystyle \ scriptstyle T_ {w} <T} , тем самым облегчая реализацию системы DSP в радиолокационной системе с линейной частотной модуляцией.

Чтобы продемонстрировать, что обработка растяжения сохраняет разрешение по дальности, нам нужно понимать, что y (t) на самом деле представляет собой последовательность импульсов с длительностью импульса T и периодом T trans {\ displaystyle \ scriptstyle T_ {trans}}{\ displaystyle \ scriptstyle T_ { транс}} , который равен периоду переданной последовательности импульсов. В результате преобразование Фурье y (t) на самом деле является функцией sinc с разрешением Рэлея 1 T {\ displaystyle \ scriptstyle {\ frac {1} {T}}}\ scriptstyle {\ frac {1} {T}} . То есть процессор сможет разрешать рассеиватели, F b {\ displaystyle \ scriptstyle F_ {b}}{\ displaystyle \ scriptstyle F_ {b}} как минимум Δ F b = 1 / T {\ displaystyle \ стиль сценария \ Delta F_ {b} = 1 / T}{\ displaystyle \ scriptstyle \ Delta F_ {b} = 1 / T} на расстоянии.

Следовательно,

1 T = | Δ f T Δ (δ t b) | ⇒ | Δ (δ t b) | Знак равно 1 Δ е {\ displaystyle {\ frac {1} {T}} = \ left \ vert {\ frac {\ Delta f} {T}} \ Delta (\ delta t_ {b}) \ right \ vert \ Rightarrow \ left \ vert \ Delta (\ delta t_ {b}) \ right \ vert = {\ frac {1} {\ Delta f}}}{\ displaystyle {\ frac {1} {T}} = \ left \ vert {\ frac {\ Delta f} {T }} \ Delta (\ delta t_ {b}) \ right \ vert \ Rightarrow \ left \ vert \ Delta (\ del ta t_ {b}) \ right \ vert = {\ frac {1} {\ Delta f}}}

и

Δ (δ R b) = c Δ ( δ tb) 2 знак равно с 2 Δ е {\ displaystyle \ Delta (\ delta R_ {b}) = {\ frac {c \ Delta (\ delta t_ {b})} {2}} = {\ frac {c} {2 \ Delta f}}}{\ displaystyle \ Delta (\ delta R_ {b}) = {\ frac {c \ Delta (\ delta t_ {b })} {2}} = {\ frac {c} {2 \ Delta f}}}

, что совпадает с разрешением исходного сигнала линейной частотной модуляции.

Форма волны ступенчатого изменения частоты

Хотя обработка растяжения может уменьшить полосу пропускания принимаемого сигнала основной полосы частот, все аналоговые компоненты во входной РЧ-схеме должны поддерживать мгновенную ширину полосы Δ е {\ Displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f . Кроме того, эффективная длина волны электромагнитной волны изменяется во время развертки частоты ЛЧМ-сигнала, и поэтому направление взгляда антенны неизбежно будет изменено в системе с фазированной решеткой.

Сигналы со ступенчатым изменением частоты - это альтернативный метод, позволяющий сохранить точное разрешение по диапазону и ОСШ принятого сигнала без большой мгновенной ширины полосы. В отличие от чирпирующего сигнала, который линейно перемещается по общей полосе пропускания Δ f {\ displaystyle \ scriptstyle \ Delta f}\ scriptstyle \ Delta f за один импульс, сигнал ступенчатой ​​частоты использует последовательность импульсов, где частота каждый импульс увеличивается на Δ F {\ displaystyle \ scriptstyle \ Delta F}{\ displaystyle \ scriptstyle \ Delta F} от предыдущего импульса. Сигнал основной полосы частот может быть выражен как:

x (t) = ∑ m = 0 M - 1 xp (t - m T) ej 2 π m Δ F (t - m T) {\ displaystyle x (t) = \ sum _ {m = 0} ^ {M-1} x_ {p} (t-mT) e ^ {j2 \ pi m \ Delta F (t-mT)}}{ \ Displaystyle х (т) = \ сумма _ {м = 0} ^ {M-1} x_ {p} (t-mT) e ^ {j2 \ pi m \ Delta F (t-mT)}}

где xp (t) {\ displaystyle \ scriptstyle x_ {p} (t)}{\ displaystyle \ scriptstyle x_ {p} (t)} - прямоугольный импульс длиной τ {\ displaystyle \ scriptstyle \ tau}\ scriptstyle \ tau , а M - количество импульсы в одной последовательности импульсов. Общая ширина полосы сигнала по-прежнему равна Δ f = M Δ F {\ displaystyle \ scriptstyle \ Delta f = M \ Delta F}{\ displaystyle \ scriptstyle \ Delta f = M \ Delta F} , но аналоговые компоненты можно сбросить для поддержки частота следующего импульса во время между импульсами. В результате можно избежать упомянутой выше проблемы.

Для вычисления расстояния до цели, соответствующего задержке tl + δ t {\ displaystyle \ scriptstyle t_ {l} + \ delta t}{\ displaystyle \ scriptstyle t_ {l} + \ delta t} , отдельные импульсы обрабатываются через простой импульсный согласованный фильтр:

hp (t) = xp ∗ (- t) {\ displaystyle h_ {p} (t) = x_ {p} ^ {*} (- t)}{\ displaystyle h_ {p } (t) = x_ {p} ^ {*} (- t)}

и выход согласованного фильтра:

ym (t) = sp ∗ (t - (tl + δ t) - m T) ej 2 π m Δ F (t - (tl + δ t) - m T) {\ displaystyle y_ {m} (t) = s_ {p} ^ {*} (t- (t_ {l} + \ delta t) -mT) e ^ {j2 \ pi m \ Delta F (t- (t_ {l } + \ delta t) -mT)}}{\ displaystyle y_ { m} (t) = s_ {p} ^ {*} (t- (t_ {l} + \ delta t) -mT) e ^ {j2 \ pi m \ Delta F (t- (t_ {l} + \ дельта t) -mT)}}

где

sp ∗ (t - (tl + δ t) - m T) = xp (t - (tl + δ t) - m T) ∗ hp (T) {\ Displaystyle s_ {p} ^ {*} (t- (t_ {l} + \ delta t) -mT) = ​​x_ {p} (t- (t_ {l} + \ delta t) -mT) * h_ {p} (t)}{\ displaystyle s_ {p} ^ {*} ( t- (t_ {l} + \ delta t) -mT) = ​​x_ {p} (t- (t_ {l} + \ delta t) -mT) * h_ {p} (t)}

Если взять образец ym (t) {\ displaystyle \ scriptstyle y_ {m} (t)}{\ displaystyle \ scriptstyle y_ {m} (t)} at t = tl + m T {\ displaystyle \ scriptstyle t = t_ {l} + mT}{\ displaystyle \ scriptstyle t = t_ {l} + mT} , мы можем получить:

y [l, m] = sp ∗ (δ t) ej 2 π m Δ F δ t {\ displaystyle y [l, m] = s_ {p} ^ {*} (\ delta t) e ^ {j2 \ pi m \ Delta F \ delta t}}{\ displaystyle y [l, m] = s_ {p } ^ {*} (\ delta t) e ^ {j2 \ pi m \ Delta F \ delta t}}

, где l означает диапазон l. Проведите DTFT (здесь m используется как время), и мы можем получить:

Y [l, ω] = ∑ m = 0 M - 1 y [l, m] e - j ω m = sp ∗ (δ t) ∑ м знак равно 0 M - 1 ej (ω - 2 π Δ F δ t) m {\ displaystyle Y [l, \ omega] = \ sum _ {m = 0} ^ {M-1} y [l, m] e ^ {- j \ omega m} = s_ {p} ^ {*} (\ delta t) \ sum _ {m = 0} ^ {M-1} e ^ {j (\ omega -2 \ pi \ Delta F \ delta t) m}}{\ displaystyle Y [ l, \ omega] = \ sum _ {m = 0} ^ {M-1} y [l, m] e ^ {- j \ omega m} = s_ {p} ^ {*} (\ delta t) \ сумма _ {м = 0} ^ {M-1} e ^ {j (\ omega -2 \ pi \ Delta F \ delta t) m}}

, и пик суммирования происходит, когда ω = 2 π Δ F δ t {\ displaystyle \ scriptstyle \ omega = 2 \ pi \ Delta F \ delta t}{\ displaystyle \ scriptstyle \ omega = 2 \ pi \ Delta F \ delta t} .

Следовательно, ДВПФ y [l, m] {\ displaystyle \ scriptstyle y [l, m]}{\ displaystyle \ scriptstyle y [l, m]} обеспечивает меру задержки цели относительно задержки ячейки диапазона tl {\ displaystyle \ scriptstyle t_ {l}}{\ displaystyle \ scriptstyle t_ {l}} :

δ t = ω p 2 π Δ F = fp Δ F {\ displaystyle \ delta t = {\ frac {\ omega _ {p}} {2 \ pi \ Delta F}} = {\ frac {f_ {p}} {\ Delta F}}}{\ displaystyle \ delta t = {\ frac { \ omega _ {p}} {2 \ pi \ Delta F}} = {\ frac {f_ {p}} {\ Delta F}}} 

и можно получить дифференциальный диапазон:

δ R = cfp 2 Δ F {\ displaystyle \ delta R = {\ frac {cf_ {p}} {2 \ Delta F}}}{\ displaystyle \ delta R = {\ frac {cf_ {p}} {2 \ Delta F}}}

где c - скорость света.

Чтобы продемонстрировать ступенчатую форму сигнала с сохранением разрешения по диапазону, следует отметить, что Y [l, ω] {\ displaystyle \ scriptstyle Y [l, \ omega]}{\ displaystyle \ scriptstyle Y [l, \ omega]} равно функция, подобная sinc, и поэтому она имеет разрешение Рэлея Δ fp = 1 / M {\ displaystyle \ scriptstyle \ Delta f_ {p} = 1 / M}{\ displaystyle \ scriptstyle \ Delta f_ {p} = 1 / M} . В результате:

Δ (δ t) = 1 M Δ F = 1 Δ f {\ displaystyle \ Delta (\ delta t) = {\ frac {1} {M \ Delta F}} = {\ frac { 1} {\ Delta f}}}{\ displaystyle \ Delta (\ delta t) = { \ frac {1} {M \ Delta F}} = {\ frac {1} {\ Delta f}}}

и, следовательно, разрешение дифференциального диапазона равно:

Δ (δ R) = c 2 Δ f {\ displaystyle \ Delta (\ delta R) = {\ frac {c} {2 \ Delta f}}}{\ displaystyle \ Delta (\ delta R) = {\ frac {c} {2 \ Delta f}}}

, что совпадает с разрешением исходного сигнала линейной частотной модуляции.

Сжатие импульсов с помощью фазового кодирования

Существуют и другие способы модуляции сигнала. Фазовая модуляция - широко используемый метод; в этом случае импульс делится на N {\ displaystyle \ scriptstyle N}\ scriptstyle N временных интервалов длительностью TN {\ displaystyle \ scriptstyle {\ frac {T} {N}}}\ scriptstyle {\ гидроразрыв {T} {N}} , для которого фаза в исходной точке выбирается в соответствии с заранее установленным соглашением. Например, можно не изменять фазу в течение некоторых временных интервалов (что сводится к тому, чтобы просто оставить сигнал таким, какой он есть в этих слотах) и сдвигать фазу сигнала в других слотах на π {\ displaystyle \ scriptstyle \ pi}\ scriptstyle \ pi (что эквивалентно изменению знака сигнала). Точный способ выбора последовательности фаз {0, π} {\ displaystyle \ scriptstyle \ {0, \, \ pi \}}\ scriptstyle \ {0, \, \ pi \} осуществляется в соответствии с техникой, известной как Коды Баркера. Возможно кодирование последовательности более чем на двух фазах (многофазное кодирование). Как и в случае с линейным чирпом, сжатие импульсов достигается за счет взаимной корреляции.

Преимущества кодов Баркера заключаются в их простоте (как указано выше, π {\ displaystyle \ scriptstyle \ pi}\ scriptstyle \ pi смещение фазы - это простая смена знака), но степень сжатия импульса ниже, чем в случае щебета, и сжатие очень чувствительно к изменениям частоты из-за эффекта Доплера, если это изменение больше, чем 1 T {\ displaystyle \ scriptstyle {\ frac {1} {T}}}\ scriptstyle {\ frac {1} {T}} .

Примечания
Дополнительная литература
См. Также
Последняя правка сделана 2021-06-02 10:39:58
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте