Способ измерения времени прихода сигнала и устройство для его реализации

Изобретение относится к радиотехнике и может использоваться в приемниках для измерения времени прихода сигналов с двухпозиционной угловой манипуляцией. Технический результат - повышение вычислительной эффективности за счет исключения неопределенности частоты приема, позволяющее обойтись поиском аргумента максимума одномерной дискретной кросс-корреляционной функции. Устройство для измерения времени прихода сигнала с двухпозиционной угловой манипуляцией содержит средство приема сигналов, аналого-цифровой преобразователь (АЦП), первый процессор БПФ, вычислитель кросс-спектров, второй процессор БПФ, устройство определения аргументов максимума ДККФ сигнала, третий процессор БПФ, а также последовательно включенные между выходом АЦП и входом первого процессора БПФ устройство квадратурного разложения сигнала, первый фильтр нижних частот (ФНЧ), постоянное запоминающее устройство (ПЗУ) и устройство вычитания по модулю 2π, между вторым выходом устройства квадратурного разложения сигнала и вторым входом ПЗУ включен второй ФНЧ, вычитающий вход устройства вычитания соединен с выходом ПЗУ через элемент задержки на длительность символа сигнала, а вход третьего процессора БПФ подключен к выходу формирователя двухуровневого модулирующего сигнала через последовательно соединенные дифференциальный декодер и экспандер частоты дискретизации. 1 з.п. ф-лы, 3 ил.

 

Изобретение относится к радиотехнике и может использоваться для измерения времени прихода сигналов с двухпозиционной угловой манипуляцией на приемной позиции.

Измерение времени прихода (ВП) сигналов на приемной позиции с известными координатами имеет большое значение в системах дальней космической связи с космическими аппаратами (КА) типа Mars Polarlander, Mars Pathfinder [1] или с КА типа Voyager и Galileo [2], в том числе для определения параметров движения таких КА.

Известен ряд способов измерения ВП сигналов, в том числе с разнесенных приемных позиций [3-6], основанных на нахождении аргументов максимума двумерной кросс-корреляционной функции, потенциально дающих статистически оптимальные максимально правдоподобные оценки. Однако потенциальная точность способов [3-6] нереализуема на практике ввиду отсутствия плавно регулируемых эталонов времени и частоты [7, 8].

Известен ряд цифровых способов измерения ВП сигналов с разнесенных приемных позиций [8], основанных на нахождении аргумента максимума кросс-корреляционной функции или аргумента минимума разностных дискретных кросс-корреляционных функций (ДККФ), позволяющих реализовать потенциальную точность способов [3-6] за счет исключения погрешности дискретности путем параболической интерполяции окрестностей максимума или минимума соответствующих функций. Однако цифровые способы измерения ВП сигналов [8-9] неприменимы при наличии неопределенности по частоте приема, что имеет место в системах дальней космической связи с космическими аппаратами. Другим недостатком цифровых способов измерения, представленных в работе [8], является низкая вычислительная эффективность на больших размерах выборки данных, поскольку они основаны на прямом вычислении кросс-корреляционных функций (при таком методе нахождения кросс-корреляционных функций количество операций умножения пропорционально квадрату длины (n) выборки данных, подобная пропорциональность обычно обозначается как О(n2)).

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

В патентах [10-12] представлен ряд цифровых способов совместного измерения разности времени прихода (РВП) и РЧП сигналов с разнесенных приемных позиций, основанных на нахождении аргументов максимума двумерной ДККФ (ДДККФ), которые позволяют производить измерения при наличии движения источника сигнала или приемной позиции. Основная цель, поставленная автором патентов [10-12], посвящена снижению потоков информации передаваемых между разнесенными приемными позициями и компенсации систематической погрешности, возникающей при корреляции компрессированного и опорного сигналов. Однако в способах совместного измерения РВП и РЧП [10-12] не решены вопросы исключения погрешности дискретности (корреляционная зависимость между РВП и РЧП [13] не позволяет очевидным способом распространить на совместные измерения подходы, изложенные в [8]) и вопросы повышения вычислительной эффективности, так как кросс-корреляционные функции в способах [10-12] порождаются прямым методом через свертку двух последовательностей данных, как и в [8], а за счет n-кратного повторения (для обеспечения поиска аргумента максимума в частотной области) количество операций умножения увеличивается и пропорционально кубу длины (n) выборки данных, т.е. O(n3).

Наиболее близким к предлагаемому способу измерения ВП сигналов, по совокупности используемых действий над сигналом, является способ [14], основанный на нахождении аргументов максимума ДДККФ (вычислительно существенно более эффективный, чем аналоги), принятый за прототип.

Согласно этому способу:

1. Принимают сигнал на двух разнесенных приемных позициях.

2. Преобразуют сигнал, полученный от одной приемной позиции, в первый цифровой поток данных, который в цифровой форме представляет сигнал как ряд значений функции времени.

3. Преобразуют сигнал, полученный от другой приемной позиции, во второй цифровой поток данных, который в цифровой форме представляет сигнал как ряд значений функции времени.

4. Преобразуют с использованием быстрого преобразования Фурье (БПФ) упомянутый первый цифровой поток данных в значения первого спектра, который представляет сигнал как ряд значений функции частоты S1(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных.

5. Преобразуют с использованием БПФ упомянутый второй цифровой поток данных в значения второго спектра, который представляет сигнал как ряд значений функции частоты S2(k).

6. Взаимно перемножают значения спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала для порождения кросс-спектра.

7. Спектр одного из сигналов Si(k), где i=1, 2 номер одного из двух спектров, преобразуют по частоте на величины m, выбираемые согласно требованиям по разрешению РЧП, из целочисленного ряда значений от единицы до (N-1), где N - длина обрабатываемого цифрового потока данных, с созданием последовательности преобразованных по частоте спектров Si,m(k), полученной из исходных спектральных составляющих сигнала по следующему правилу

,

где i - номер одного из двух принимаемых сигналов.

8. Значения спектров Si,m(k) преобразованных по частоте сигналов взаимно перемножают с комплексно сопряженными значениями спектра второго из двух исходных сигналов для порождения множества кросс-спектров.

9. Вычисляют ДДККФ сигнала с использованием обратного БПФ множества кросс-спектров.

10. Определяют разности времени прихода и частоты приема сигналов как аргументы максимума ДДККФ сигнала.

По существу способ, описанный выше [14], реализует вычислительно эффективное по сравнению с аналогами [10-12] нахождение аргументов максимума ДДККФ с использованием при вычислении ДДККФ быстрых преобразований Фурье на основании теоремы Винера-Хинчина [15], определяющей взаимосвязь между спектром и корреляционной функцией сигнала.

В описании прототипа [14] представлена оценка объема вычислительных затрат на нахождение двумерной кросс-корреляционной функции в способе-прототипе и в аналогах, показывающая, что для выборки данных длиной n=1024 при k=n в способе-прототипе на нахождение ДДККФ требуется порядка О(k·n·log2n+k·n) или 10 миллионов операций (страница 10 описания способа-прототипа), а в аналогах требуется от O(3·k·n·log2n+k·n) до O(n3) или от 30 миллионов до миллиарда операций (страница 5 описания способа-прототипа).

Несмотря на существенное снижение вычислительных затрат по сравнению с аналогами, одним из недостатков способа-прототипа является недостаточная вычислительная эффективность. Другим недостатком способа-прототипа является, как и у аналогов [10-12], длительное время нахождения аргумента максимума, то есть длительное время измерения.

Устройство-прототип [14] содержит первое средство приема сигналов, подключенное к устройству определения аргументов максимума ДДККФ, через последовательно включенные первый аналого-цифровой преобразователь (АЦП), первый процессор БПФ, вычислитель кросс-спектров и второй процессор БПФ. Выход устройства определения аргументов максимума ДДККФ является выходом устройства измерения. Между вторым средством приема сигналов и вторым входом вычислителя кросс-спектров последовательно включены второй АЦП и третий процессор БПФ. Между выходом первого процессора БПФ и третьим входом вычислителя кросс-спектров включено арифметическое устройство.

Недостатком устройства-прототипа является недостаточное быстродействие определения РЧП и РВП сигналов. Другим недостатком устройства-прототипа является длительное время нахождения аргумента максимума, то есть длительное время измерения.

Техническим результатом изобретения является повышение вычислительной эффективности за счет исключения неопределенности частоты приема, позволяющее исключить необходимость двумерного поиска аргументов максимума ДДККФ для сигнала с двухпозиционной угловой манипуляцией, и обойтись, таким образом, поиском аргумента максимума одномерной ДККФ.

Технический результат достигнут тем, что в способе измерения времени прихода сигнала с двухпозиционной угловой манипуляцией, включающем прием сигнала, аналого-цифровое преобразование принятого сигнала, в первый цифровой поток данных, представляющий сигнал как ряд значений функции времени, преобразованных в цифровую форму, использование быстрого преобразования Фурье (БПФ) для двух сигналов, представляющего оба сигнала в виде ряда дискретных значений спектра Si(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных, а i=1, 2 - индексная переменная, означающая порядковый номер преобразуемого с использованием БПФ цифрового потока данных в значения спектра, взаимное перемножение значений спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала, порождающее кросс-спектр, вычисление дискретной кросс-корреляционной функции (ДККФ) сигнала с использованием обратного БПФ, определение времени прихода (ВП) сигнала как аргумента максимума ДККФ сигнала согласно изобретению осуществляют квадратурное разложение первого цифрового потока данных относительно номинальной центральной частоты модулированного сигнала в соответствующее указанному потоку множество синфазных и квадратурных отсчетов, полученные синфазные и квадратурные отсчеты независимо низкочастотно фильтруют с частотой среза, соответствующей скорости манипуляции модулирующего сигнала, получают множество текущих фаз сигнала как аргументов комплексных чисел, в качестве действительной части которых используют соответствующие отфильтрованные синфазные отсчеты, а в качестве мнимой - соответствующие отфильтрованные квадратурные отсчеты сигнала, задерживают множество полученных текущих фаз сигнала на длительность символа модулирующего сигнала, вычитают по модулю 2π из каждой полученной текущей фазы соответствующее значение задержанной текущей фазы сигнала, а полученный разностный цифровой поток данных преобразуют с использованием БПФ в значения первой функции частоты S1(k), формируют двухуровневый модулирующий сигнал и дифференциально декодируют, отображают каждый символ дифференциально декодированного сигнала на соответствующее количество отсчетов преобразованного в цифровую форму принятого сигнала относительно шкалы времени приемной позиции и формируют, таким образом, второй цифровой поток данных, который преобразуют с использованием БПФ в значения второй функции частоты S2(k).

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

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

Другим техническим результатом изобретения является исключение погрешности дискретности за счет использования того факта, что окрестность главного максимума дискретной кросс-корреляционной функции сигнала имеет форму параболы [13, стр.8-9], и за счет разработки оптимального в среднеквадратическом смысле аналитического метода оценки параметров параболы и определения аргумента ее максимума.

Технический результат достигнут тем, что в способе измерения времени прихода сигнала согласно изобретению слева и справа от аргумента максимума ДККФ сигнала выбирают по L точек указанной функции, значения ДККФ для выбранных точек и точки максимума объединяют в порядке временного следования в вектор-столбец с, из единиц формируют вектор-столбец v0 такой же размерности, как вектор с, соответствующие выбранным точкам ДККФ аргументы, выражают через индексную переменную l=0, ±1, …, ±L, где индексы l с отрицательным знаком соответствуют временам до точки максимума ДККФ, нулевой индекс l соответствует точке максимума ДККФ, а индексы l с положительным знаком соответствуют временам после точки максимума ДККФ, аргументы выбранных точек ДККФ упорядочивают так же, как и компоненты вектора с, и объединяют в вектор-столбец v1, квадраты компонент вектора v1 аналогично объединяют в вектор-столбец v2, а уточненное значение ВП определяют по следующему правилу

,

где N=2L+1 - размерность сформированных векторов, надстрочный индекс T обозначает операцию транспонирования вектора, Ts - период дискретизации сигнала, - скалярное произведение векторов v1 и с, коэффициент S определяется выражением

,

коэффициент Q определяется другим выражением

,

a (NSv2-S2v0)Tc - скалярное произведение разностного вектора (NSv2-S2v0) с вектором с.

Способ реализуется устройством измерения прихода сигнала с двухпозиционной угловой манипуляцией по пп.1 и 2, содержащим средство приема сигналов, вход которого является входом устройства измерения, соединенного с аналого-цифровым преобразователем (АЦП), последовательно включенные первый процессор БПФ, вычислитель кросс-спектров, второй процессор БПФ и устройство определения аргументов максимума дискретной кросс-корреляционной функции (ДККФ) сигнала, выход которого является выходом устройства измерения, при этом к второму входу вычислителя кросс-спектров подключен третий процессор БПФ, согласно изобретению между выходом АЦП и входом первого процессора БПФ последовательно включены устройство квадратурного разложения сигнала, первый фильтр нижних частот (ФНЧ), постоянное запоминающее устройство (ПЗУ) и устройство вычитания по модулю 2π, между вторым выходом устройства квадратурного разложения сигнала и вторым входом ПЗУ включен второй ФНЧ, вычитающий вход устройства вычитания соединен с выходом ПЗУ через элемент задержки на длительность символа сигнала, а вход третьего процессора БПФ подключен к выходу формирователя двухуровневого модулирующего сигнала через последовательно соединенные дифференциальный декодер и экспандер частоты дискретизации.

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

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

На фиг.1 представлен пример полной фазы сигнала с двухпозиционной фазовой манипуляцией при наличии частотной расстройки.

На фиг.2 представлен пример разности фаз сигналов со сдвигом на длительность символа Т.

На фиг.3 приведена структурная схема устройства, в котором реализуется предложенный способ.

Согласно предлагаемому способу:

1. Принимают сигнал.

2. Преобразуют принятый сигнал в первый цифровой поток данных, представляющий сигнал как ряд значений функции времени, преобразованных в цифровую форму.

3. Осуществляют квадратурное разложение первого цифрового потока данных относительно номинальной центральной частоты модулированного сигнала в соответствующее указанному потоку множество синфазных и квадратурных отсчетов.

4. Полученные синфазные и квадратурные отсчеты независимо низкочастотно фильтруют с частотой среза, соответствующей скорости манипуляции модулирующего сигнала.

5. Получают множество текущих фаз сигнала как аргументов комплексных чисел, в качестве действительной части которых используют соответствующие отфильтрованные синфазные отсчеты, а в качестве мнимой - соответствующие отфильтрованные квадратурные отсчеты сигнала.

6. Задерживают множество полученных текущих фаз сигнала на длительность символа модулирующего сигнала.

7. Вычитают по модулю 2π из каждой полученной текущей фазы соответствующее значение задержанной текущей фазы сигнала.

8. Преобразуют с использованием быстрого преобразования Фурье (БПФ) полученный разностный цифровой поток данных как первый цифровой поток данных в дискретные значения первого спектра S1(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных.

9. Формируют двухуровневый модулирующий сигнал.

10. Дифференциально декодируют сформированный сигнал.

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

12. Преобразуют с использованием БПФ упомянутый второй цифровой поток данных в значения второго спектра S2(k).

13. Взаимно перемножают значения спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала для порождения кросс-спектра.

14. Вычисляют дискретную кросс-корреляционную функцию (ДККФ) сигнала с использованием обратного БПФ кросс-спектра.

15. Определяют время прихода сигнала как аргумент максимума ДККФ сигнала.

Покажем, что в предложенном способе измерения времени прихода сигнала с двухпозиционной угловой манипуляцией возможно исключение необходимости двумерного поиска при нахождении аргумента максимума ДККФ.

Сигнал с двухпозиционной фазовой манипуляцией описывается следующим образом [24]:

где bk∈{-1,+1}, p(t) - прямоугольный импульс с единичной амплитудой длительностью Т, определенный на интервале [-T/2, T/2], А - амплитуда сигнала, а ωс - частота несущей. Полная фаза такого сигнала

где θk∈{0,π} и связана с bk соотношением

При квадратурном разложении сигнала (1) несущая частота, которая априорно неизвестна, в полной фазе сигнала практически всегда будет частотой ω0, относительно которой осуществляется квадратурное разложение, то есть в полной фазе сигнала будет присутствовать компонента, линейно зависящая от времени, которая соответствует второму слагаемому выражения (3), что не позволяет напрямую использовать выражение для полной фазы для формирования ДККФ.

На фиг.1 представлен пример полной фазы сигнала с двухпозиционной фазовой манипуляцией при наличии частотной расстройки, демонстрирующий влияние компоненты, линейно зависящей от времени.

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

Разность фаз со сдвигом на длительность символа Т позволяет исключить частотную расстройку как компоненту, линейно зависящую от времени

Как видно из выражения (5), разность фаз Ψ(t) за исключением постоянного смещения ω0Т, определяемого величиной частотной расстройки, однозначно соответствует дифференциально декодированной модулирующей функции сигнала (смотри фиг.2).

Поскольку постоянное смещение не влияет на аргумент кросс-корреляционной функции [19], то одномерная корреляционная обработка дифференциально декодированной модулирующей функции сигнала с функцией разности фаз сигнала Ψ(t) позволит определить время прихода сигнала как аргумент максимума кросс-корреляционной функции указанных сигналов. Аналогичный результат можно показать и для сигнала с манипуляцией минимальным сдвигом [24].

Согласно второму варианту предлагаемого способа, исключающему погрешность измерений, вызываемую некратностью длительности двухуровневого модулирующего символа и частоты выборки аналогоцифрового преобразования принятого сигнала:

1. Принимают сигнал.

2. Преобразуют принятый сигнал в первый цифровой поток данных, представляющий сигнал как ряд значений функции времени, преобразованных в цифровую форму.

3. Осуществляют квадратурное разложение первого цифрового потока данных относительно номинальной центральной частоты модулированного сигнала в соответствующее указанному потоку множество синфазных и квадратурных отсчетов.

4. Полученные синфазные и квадратурные отсчеты независимо низкочастотно фильтруют с частотой среза, соответствующей скорости манипуляции модулирующего сигнала.

5. Отфильтрованные синфазные и квадратурные отсчеты дополнительно децимируют и интерполируют с обеспечением кратности частоты отсчетов первого цифрового потока данных длительности символа сигнала.

6. Получают множество текущих фаз сигнала как аргументов комплексных чисел, в качестве действительной части которых используют соответствующие синфазные отсчеты, а в качестве мнимой - соответствующие квадратурные отсчеты сигнала.

7. Задерживают множество полученных текущих фаз сигнала на длительность символа модулирующего сигнала.

8. Вычитают по модулю 2π из каждой полученной текущей фазы соответствующее значение задержанной текущей фазы сигнала.

9. Преобразуют с использованием быстрого преобразования Фурье (БПФ) полученный разностный цифровой поток данных как первый цифровой поток данных в дискретные значения первого спектра S1(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных.

10. Формируют двухуровневый модулирующий сигнал.

11. Дифференциально декодируют сформированный сигнал.

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

13. Преобразуют с использованием БПФ упомянутый второй цифровой поток данных в значения второго спектра S2(k).

14. Взаимно перемножают значения спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала для порождения кросс-спектра.

15. Вычисляют дискретную кросс-корреляционную функцию (ДККФ) сигнала с использованием обратного БПФ кросс-спектра.

16. Определяют время прихода сигнала как аргумент максимума ДККФ сигнала.

Некратность длительности двухуровневого модулирующего символа и частоты выборки аналого-цифрового преобразования принятого сигнала можно устранить, к примеру, полифазной фильтрацией [18], позволяющей за счет интерполяции и децимации обеспечить целочисленное количество отсчетов сигнала на длительности символа на конечном интервале времени.

Согласно другому варианту предлагаемого способа, исключающему погрешность дискретности:

1. Принимают сигнал.

2. Преобразуют принятый сигнал в первый цифровой поток данных, представляющий сигнал как ряд значений функции времени, преобразованных в цифровую форму.

3. Осуществляют квадратурное разложение первого цифрового потока данных относительно номинальной центральной частоты модулированного сигнала в соответствующее указанному потоку множество синфазных и квадратурных отсчетов.

4. Полученные синфазные и квадратурные отсчеты независимо низкочастотно фильтруют с частотой среза, соответствующей скорости манипуляции модулирующего сигнала.

5. Получают множество текущих фаз сигнала как аргументов комплексных чисел, в качестве действительной части которых используют соответствующие отфильтрованные синфазные отсчеты, а в качестве мнимой - соответствующие отфильтрованные квадратурные отсчеты сигнала.

6. Задерживают множество полученных текущих фаз сигнала на длительность символа модулирующего сигнала.

7. Вычитают по модулю 2π из каждой полученной текущей фазы соответствующее значение задержанной текущей фазы сигнала.

8. Преобразуют с использованием быстрого преобразования Фурье (БПФ) полученный разностный цифровой поток данных как первый цифровой поток данных в дискретные значения первой функции частоты S1(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных.

9. Формируют двухуровневый модулирующий сигнал.

10. Дифференциально декодируют сформированный сигнал.

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

12. Преобразуют с использованием БПФ упомянутый второй цифровой поток данных в значения второго спектра S2(k).

13. Взаимно перемножают значения спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала для порождения кросс-спектра.

14. Вычисляют дискретную кросс-корреляционную функцию (ДККФ) сигнала с использованием обратного БПФ кросс-спектра.

15. Определяют время прихода сигнала как аргумент максимума ДККФ сигнала.

16. Слева и справа от аргумента максимума ДККФ сигнала выбирают по L точек указанной функции.

17. Значения ДККФ для выбранных точек и точки максимума объединяют в порядке временного следования в вектор-столбец с.

18. Из единиц формируют вектор-столбец v0 такой же размерности, как вектор с.

19. Соответствующие выбранным точкам ДККФ аргументы выражают через индексную переменную l=0, ±1, …, ±L, где индексы l с отрицательным знаком соответствуют временам до точки максимума ДККФ, нулевой индекс l соответствует точке максимума ДККФ, а индексы l с положительным знаком соответствуют временам после точки максимума ДККФ.

20. Аргументы выбранных точек ДККФ упорядочивают так же, как и компоненты вектора с, и объединяют в вектор-столбец v1.

21. Квадраты компонент вектора v1 аналогично объединяют в вектор-столбец v2.

22. Определяют по следующему правилу уточненное значение ВП

,

где N=2L+1 - размерность сформированных векторов, надстрочный индекс T обозначает операцию транспонирования вектора, Ts - период дискретизации сигнала, - скалярное произведение векторов v1 и c, коэффициент S определяется выражением

,

коэффициент Q определяется другим выражением

,

a (NSv2-S2v0)Tc - скалярное произведение разностного вектора (NSv2-S2v0) с вектором с.

При цифровом преобразовании непрерывного сигнала неизбежно появляются погрешности дискретности [17-18], их влияние может быть существенно снижено интерполяцией ДККФ в окрестности экстремума, что не влечет значительных вычислительных затрат, возникающих при интерполяции сигнала в целом.

Вариант предлагаемого способа, исключающий погрешность дискретности при оценке параметров ДККФ, основан на том, что в окрестности максимума указанная функция имеет форму параболы [13, стр.8-9], оценив параметры которой в результате решения системы линейных уравнений, можно уточнить значение времени прихода как аргумента максимума ДККФ, поскольку ДККФ является выборочной от кросс-корреляционной функции в непрерывном времени.

Уравнение, описывающее вершину кросс-корреляционной функции CXY(τ) как параболу, имеет вид

Для заданных отсчетов cn кросс-корреляционной функции CXY(τ)

где ξn - погрешность задания отсчетов КФ, Ts - период дискретизации.

Для параболы (6) оценка точки экстремума (в данном случае максимума исходя из смысла задачи и предположения о единственности экстремума в рассматриваемой окрестности) определяется соотношением

Оценки для коэффициентов a, b в (8) могут быть получены из линейной системы вида

где ξ - вектор погрешностей задания отсчетов кросс-корреляционной функции (7). Кроме того, в (9) предполагается нечетное число исходных данных (7) и принята симметричная индексация n=-K, -(K-1), …, -1, 0, 1, 2, …, K, и хотя данные предположения не существенны для общей оценки по методу наименьших квадратов (МНК), они не умаляют общность, а их смысл будет пояснен позднее.

Согласно [22] МНК оценка вектора коэффициентов для (9) и соответствующая (8) оценка задержки имеют вид

где V - псевдообратная матрица, которая при существовании W=(VTV)-1 имеет вид V=WVT [22-23].

Приняв оценку , определяемую в (8) за требуемое решение, заметим, что для ее получения в указанном виде должна быть использована операция обращения матрицы, что не всегда вычислительно эффективно, особенно при реализации слежения за задержкой в реальном времени. Матрица V не зависит от исходных измерений (7) и может быть рассчитана априорно, но она зависит от объема данных N. При таком подходе необходимо хранить набор матриц для различных N либо требовать работы с всегда фиксированным объемом, что тоже может оказаться не оптимальным. Кроме того, при вычислении оценки согласно первому выражению в (10) осуществляются затраты на оценку не используемую далее. С целью повышения вычислительной эффективности получим для величины тождественные соотношения вида

где αN, βN - весовые N-вектора, определяемые явными арифметическими зависимостями от N.

Рассматривая в (9) матрицу V=[v0 v1 v2] как совокупность ее столбцов, определения которых очевидно из (9), заметим, что в силу симметрии имеют место следующие свойства

Тогда используемая в оценке (10) и требующая обращения матрица W принимает вид

где D=det(VTV)=S(NQ-S2) - определитель, значение которого для дальнейшего можно и не вычислять.

Из определений (10) и (11) видно, что искомые весовые функции , являются второй и третьей строками матрицы V, соответственно получим

Тогда с учетом (14) из определения (11) можно найти уточненную искомую оценку времени прихода в виде

Таким образом, при соответствующем выборе точек в окрестности максимума ДККФ, являющейся, в дискретном времени, выборочной функцией для кросс-корреляционной функции CXY(τ), с соответствующим формированием из выбранных отсчетов ДККФ вектора с, с формированием из аргументов выбранных отсчетов и их квадратов векторов v1 и v2, а также с формированием вектора v0 из единиц можно найти уточненную оценку времени прихода при использовании правила (15).

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

Устройство для реализации способа измерения времени прихода сигнала с двухпозиционной угловой манипуляцией (смотри Фигуру 3) содержит средство приема сигналов 1, соединенное с аналого-цифровым преобразователем 2 (АЦП), последовательно включенные первый процессор БПФ 3, вычислитель кросс-спектров 4, второй процессор БПФ 5 и устройство определения аргумента максимума дискретной кросс-корреляционной функции (ДККФ) сигнала 6. Вход средства приема сигналов 1 является входом устройства измерения. Выход устройства 6 является выходом устройства измерения. Ко второму входу вычислителя кросс-спектров 4 подключен третий процессор БПФ 7. Между выходом АЦП 2 и входом первого процессора БПФ 3 последовательно включены устройство квадратурного разложения сигнала 8, первый фильтр нижних частот (ФНЧ) 9, постоянное запоминающее устройство (ПЗУ) 10 и устройство вычитания по модулю 2π 11. Между вторым выходом устройства квадратурного разложения сигнала 8 и вторым входом ПЗУ 10 включен второй ФНЧ 12. Вычитающий вход устройства вычитания 11 соединен с выходом ПЗУ 10 через элемент задержки на длительность символа сигнала 13. Вход третьего процессора БПФ 7 подключен к выходу формирователя двухуровневого модулирующего сигнала 14 через последовательно соединенные дифференциальный декодер 15 и экспандер частоты дискретизации 16.

В другом варианте устройства измерения времени прихода сигнала к выходу устройства определения аргумента максимума ДККФ сигнала 6 дополнительно подключено арифметическое устройство 17.

Предлагаемое устройство для реализации способа измерения времени прихода сигнала с двухпозиционной угловой манипуляцией работает следующим образом.

Средство приема сигналов 1 принимает сигнал. С выхода первого средства приема сигналов 1 принятый сигнал подают на вход АЦП 2 и преобразуют сигнал в первый цифровой поток данных, который в цифровой форме представляет сигнал как ряд значений функции времени. Устройство квадратурного разложения сигнала 8 осуществляет квадратурное разложение первого цифрового потока данных относительно номинальной центральной частоты модулированного сигнала в соответствующее указанному потоку множество синфазных и квадратурных отсчетов. Полученные синфазные и квадратурные отсчеты независимо низкочастотно фильтруют ФНЧ 9 и 12 с частотой среза, соответствующей скорости манипуляции модулирующего сигнала. В ПЗУ 10 получают множество текущих фаз сигнала как аргументов комплексных чисел, в качестве действительной части которых используют соответствующие отфильтрованные синфазные отсчеты, а в качестве мнимой - соответствующие отфильтрованные квадратурные отсчеты сигнала. Элемент задержки 13 задерживает множество полученных текущих фаз сигнала на длительность символа модулирующего сигнала. В устройстве вычитания 11 вычитают по модулю 2π из каждой полученной текущей фазы соответствующее значение задержанной текущей фазы сигнала. В первом процессоре БПФ 3 преобразуют полученный разностный цифровой поток данных как первый цифровой поток данных в дискретные значения первой функции частоты S1(k), где k - целочисленная индексная переменная, изменяющаяся в пределах длины данных. Формирователь 14 формирует двухуровневый модулирующий сигнал. Сформированный сигнал дифференциально декодируют декодером 15. Каждый символ дифференциально декодированного сигнала отображают на соответствующее количество отсчетов преобразованного в цифровую форму принятого сигнала относительно шкалы времени приемной позиции путем увеличения частоты дискретизации с использованием экспандера частоты дискретизации 16 и формируют, таким образом, второй цифровой поток данных. В третьем процессоре БПФ 7 преобразуют упомянутый второй цифровой поток данных в значения второго спектра, который представляет сигнал как ряд значений функции частоты S2(k). В вычислителе кросс-спектров 4 взаимно перемножают значения спектра одного из сигналов с комплексно сопряженными значениями спектра другого сигнала, поступившие из первого 3 и третьего 7 процессоров БПФ для порождения кросс-спектра. Выходной сигнал из вычислителя 4 подается на вход второго процессора БПФ 5, в котором с использованием обратного БПФ кросс-спектра вычисляют ДККФ сигнала. Полученная ДККФ сигнала с выхода второго процессора БПФ 5 поступает на вход устройства 2, где определяют время прихода сигнала как аргумент максимума ДККФ сигнала.

В другом варианте устройства измерения времени прихода сигнала выходной сигнал устройства определения аргумента максимума 6 дополнительно подается на арифметическое устройство 17, которое работает в соответствии с правилом (15) и определяет уточненное значение времени прихода сигнала.

ИСТОЧНИКИ ИНФОРМАЦИИ

1. Harcke, L., and G.Wood, Laboratory and Flight Performance of the Mars Pathfinder (15,1/6) Convolutionally Encoded Telemetry Link, TDA PR 42-129, January-March 1997, pp.1-11, May 15, 1997.

2. Rebold, T.A., M.Tinto, S.W.Asmar, and E.R.Kursinski, Neptune Revisited: Synthesizing Coherent Doppler From Voyager′s Noncoherent Downlink, TDA PR 42-131, July-September 1997, pp.1-19, November 15, 1997.

3. Knapp С.H. and Carter G.С., Estimation of time delay in the presence of source or receiver motion, J. Acoust. Soc. Amer., vol.61, no.6, pp.1545-1549, June 1977.

4. Adams W.В., Kuhn J.P., and Whyland W.P., Correlator compensation requirements for passive time-delay estimation with moving source or receivers, IEEE Trans. Acoust, Speech, Signal Processing, vol.ASSP-28, pp.158-168, Apr. 1980.

5. Stuller J.A., Maximum-likelihood estimation of time-varying delay - Part I, IEEE Trans. Acoust, Speech, Signal Processing, vol.ASSP-35, pp.300-313, Mar. 1987.

6. Jin Q., Wong K.M., and Luo Z.Q.Т., The estimation of time delay and doppler stretch of wideband signals, IEEE Trans. Signal Processing, vol.43, pp.904-916, Apr. 1995.

7. Одуан К., Гино Б. Измерение времени. Основы GPS, M.: Техносфера, 2002, 400 с.

8. Jacovitti G. and Scarano G., Discrete time techniques for time delay estimation, IEEE Trans. Signal Processing, vol.41, pp.525-533, Feb. 1993 (аналоги).

9. US Patent No. 6201499, Hawkes K. D., Koehler J. L., Time difference of arrival measurement system, МПК 7: G01S 001/24, March 13, 2001 (аналог).

10. US Patent No. 5570099, DesJardins G.A., TDOA/FDOA technique for locating a transmitter, МПК 7: G01S 003/16; G01S 005/04; G01S 015/00; G01S 003/80, October 29, 1996 (аналог).

11. US Patent No. 5874916, DesJardins G.A., Frequency selective TDOA/FDOA cross-correlation, February 23, 1999 (аналог).

12. European Patent No. 0649033 A2, DesJardins G.A., Improved TDOA/FDOA technique for locating a transmitter, МПК 7: G01S 5/06, April 19, 1995 (аналог).

13. Moura Jose M.F., Baggeroer A.B., Passive Systems Theory with Narrow-band and linear Constraints: Part 1 - Spatial Diversity, IEEE Journal on Ocean Engineering, vol.OE-3, pp.5-13, No.1 Jan. 1978.

14. Подчиненко Н.Е., Скрипкин А.А., Щербачев В.А. Патент РФ №2256192. Способ измерения разностей времени прихода и частоты приема сигналов и устройство для его реализации (Способ и устройство-прототип).

15. Радиотехнические цепи и сигналы. - Под ред. К.А.Самойло. - М.: Радио и связь, 1982. - 528 с.

16. Блейхут Р. Быстрые алгоритмы цифровой обработки сигналов. - М.: Мир, 1989, 448 с.

17. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. - М.: Мир, 1978, 848 с.

18. Оппенгейм А.В., Шафер Р.В. Цифровая обработка сигналов. - М.: Связь, 1979, 416 с.

19. Бендат Дж., Пирсол А. Прикладной анализ случайных данных. - М.: Мир, 1989, 540 с.

20. Френкс Л. Теория сигналов. - М.: Сов. радио, 1974, 344 с.

21. Беклемишев Д.В. Курс аналитической геометрии и линейной алгебры. - М.: Наука, 1974, 320 с.

22. Альберт А. Регрессия, псевдоинверсия, рекуррентное оценивание. - М.: Наука, 1977, 223 с.

23. Гантмахер Ф.Р. Теория матриц. - М.: Наука, 1988, 552 с.

24. F.Xiong, Digital Modulation Techniques, Artech House, Second Edition. 2006.

1. Устройство для измерения времени прихода сигнала с двухпозиционной угловой манипуляцией, содержащее средство приема сигналов, вход которого является входом устройства измерения, соединенное с аналого-цифровым преобразователем (АЦП), последовательно включенные первый процессор БПФ, вычислитель кросс-спектров, второй процессор БПФ и устройство определения аргументов максимума дискретной кросс-корреляционной функции (ДККФ) сигнала, выход которого является выходом устройства измерения, при этом ко второму входу вычислителя кросс-спектров подключен третий процессор БПФ, отличающееся тем, что между выходом АЦП и входом первого процессора БПФ последовательно включены устройство квадратурного разложения сигнала, первый фильтр нижних частот (ФНЧ), постоянное запоминающее устройство (ПЗУ) и устройство вычитания по модулю 2π, между вторым выходом устройства квадратурного разложения сигнала и вторым входом ПЗУ включен второй ФНЧ, вычитающий вход устройства вычитания соединен с выходом ПЗУ через элемент задержки на длительность символа сигнала, а вход третьего процессора БПФ подключен к выходу формирователя двухуровневого модулирующего сигнала через последовательно соединенные дифференциальный декодер и экспандер частоты дискретизации.

2. Устройство для измерения времени прихода сигнала с двухпозиционной угловой манипуляцией по п.1, отличающееся тем, что к выходу устройства определения аргумента максимума дискретной кросс-корреляционной функции сигнала дополнительно подключено арифметическое устройство.



 

Похожие патенты:

Изобретение относится к радиотехнике и может использоваться для измерения времени прихода /4-QPSK сигнала на приемной позиции. .

Изобретение относится к радиотехнике и может использоваться для измерения времени прихода сигналов с М-позиционной квадратурной амплитудной манипуляцией на приемной позиции.

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

Изобретение относится к контрольно-измерительной технике. .

Изобретение относится к навигации, а именно к системам определения положения объекта без использования отражения или вторичного излучения, и может быть использовано для коррекции инерциальных навигационных систем летательных аппаратов, систем прицеливания и предупреждения столкновений.
Изобретение относится к области измерительной техники и может быть использовано в обзорно-поисковых системах, в частности в пассивных РЛС. .

Изобретение относится к области автомобилестроения. .

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

Изобретение относится к радиотехнике и может использоваться для измерения времени прихода сигнала с четырехпозиционной (квадратурной) фазовой манипуляцией со сдвигом (OQPSK) на приемной позиции

Изобретение относится к области информационно-измерительных систем и предназначается для решения задач измерения дальности и линейных размеров объектов по их цифровым фотографическим изображениям

Устройство относится к радиотехнике, а именно к антенно-фидерным устройствам СВЧ бортового радиооборудования самолетов. Техническим результатом является обеспечение кругового обзора пространства приемопередатчиком и тремя радиоприемными устройствами с трехантенной системой и улучшение энергетических характеристик коммутационно-разделительного устройства. Коммутационно-разделительное устройство содержит три коммутатора, три четвертьволновых трансформатора сопротивлений, управляемый фазовращатель, тройник, гибридные устройства, которые выполнены в виде трех направленных ответвителей, первый коммутатор выполнен на основе моста из p-i-n-диодов, второй и третий коммутаторы выполнены на основе полумостов из p-i-n-диодов. 2 ил.

Изобретение относится к локационной технике. Достигаемый технический результат - уменьшение габаритов без уменьшения точности определения угловых координат и дальности до удаленных объектов. Указанный результат достигается благодаря использованию антенны, обеспечивающей моноимпульсную обработку моноимпульсного измерителя координат, элемента ИЛИ и блока автосопровождения по направлению, при этом выход антенны, обеспечивающей моноимпульсную обработку, соединен с входом моноимпульсного измерителя координат, имеющего группу выходов, соединенную с группой входов элемента ИЛИ и с четвертой группой входов блока вторичной обработки, имеющего вторую и третью группы входов, соответственно соединенные через блок автосопровождения по направлению с первой и второй группами выходов телевизионного координатора, и имеющим третью группу выходов, соединенную с третьей группой входов индикатора, при этом выход элемента ИЛИ соединен со вторым входом преобразователя дальности. 1 ил.

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

Изобретение относится к определению местоположения с использованием нескольких разнесенных источников излучения. Достигаемый технический результат - автоматизация процесса, повышение точности измерения. Указанный результат достигается за счет того, что по сигналам трех лазерных маяков, установленных вдоль профиля строящейся автотрассы, посредством двух разнесенных цифровых фотокамер, установленных на раме дорожно-строительной машины и выполненных в виде фоточувствительных матриц, размещенных в фокальной плоскости фотообъектива, осуществляют обработку оцифрованного изображения, снимаемого с фоточувствительных матриц для определения координат изображений лазерных маяков, вычисляют координаты трех лазерных маяков относительно дорожно-строительной машины, а затем определяют положение дорожно-строительной машины относительно автотрассы. 1 ил.

Изобретение относится к информационно-измерительной технике. Достигаемый технический результат изобретения - повышение точности определения расстояния от наблюдателя до объектов с работающими на излучение источниками радиоволн и определения координат данных источников. Указанный результат достигается за счет того, что способ определения расстояния от наблюдателя, являющегося техническим средством, до работающего на излучение источника радиоволн заключается в определении высоты пересечения линии радиогоризонта над уровнем моря hрг данного источника излучения наблюдателем, которая устанавливается по факту входа или выхода наблюдателя в область пространства, в которой энергия излучения передается прямой волной, по установленным в процессе наблюдения радиотехническим характеристикам источника излучения радиоволн определяют данный источник радиоволн, отождествляют данный источник радиоволн с его носителем, на котором он используется, определяют из известных геометрических характеристик определенного носителя его высоту, которая соответствует высоте источника излучения над уровнем моря hн, после чего расстояние от наблюдателя до работающего на излучение источника радиоволн определяют согласно формуле , где a - радиус земного шара. 2 з.п. ф-лы, 1 ил.

Изобретение относится к средствам для измерения времени прихода сигналов с двухпозиционной угловой манипуляцией на приемной позиции. Техническим результатом изобретения является повышение вычислительной эффективности и повышение точности измерения. Данный технический результат достигается за счет исключения неопределенности частоты приема, позволяющего исключить необходимость двумерного поиска аргументов максимума двумерной дискретной кросс-корреляционной функции (ДККФ) для сигнала с двухпозиционной угловой манипуляцией и обойтись, таким образом, поиском аргумента максимума одномерной ДККФ. Формируют двухуровневый модулирующий сигнал и дифференциально декодируют. Отображают каждый символ дифференциально декодированного сигнала на соответствующее количество отсчетов преобразованного в цифровую форму принятого сигнала относительно шкалы времени приемной позиции путем повтора значений дифференциально декодированного сигнала с повышением его частоты выборки до частоты выборки преобразованного в цифровую форму принятого сигнала. Формируют второй цифровой поток данных, который преобразуют с использованием БПФ в значения второго спектра. 2 н.з. и 3 з.п. ф-лы, 3 ил.

Предлагаемая группа изобретений относится к области вооружения и военной техники, в частности к стрельбе комплекса вооружения боевой машины (БМ) по цели. Предлагаемый способ стрельбы вооружения БМ по цели включает обнаружение и распознавание цели, взятие на сопровождение и сопровождение цели с одновременным дальнометрированием, определение угловых поправок стрельбы из математических выражений с использованием в качестве входных параметров, в частности, значений угловых скоростей, поступающих с органов управления наводчика или командира. Способ также включает постоянное отклонение с учетом угловых поправок стволов пушечной или пулеметной установки (ПУ) относительно линии визирования и стрельбу по цели. При определении угловых поправок стрельбы используют значения угловых скоростей, скорректированных с учетом предварительно полученной до стрельбы экспериментальной зависимости угловой скорости линии визирования. Согласно изобретению система снабжена последовательно соединенными блоком управления (БУ) и устройством корректирования угловой скорости линии визирования по горизонтальному и вертикальному каналам. Для определения экспериментальной угловой скорости линии визирования последовательно с использованием специально организуемого стенда для измерения угловой скорости для каждой заданной угловой скорости поочередно перемещают башню или блок оружия соответственно по горизонтальному или вертикальному каналам. При каждом перемещении по истечении заданного времени замеряют их углы поворота, определяют искомую угловую скорость линии визирования. По полученным значениям воспроизводят зависимость угловой скорости линии визирования от угловой скорости, поступающей с органов управления ПН или ПК или автомата сопровождения, и запоминают эту зависимость. Достигается повышение точности и соответственно эффективности стрельбы вооружения БМ по подвижным целям, в особенности по скоростным целям, а также при больших полетных временах снаряда, обусловленных, в частности, низкоскоростной баллистикой и большими дальностями стрельбы. 3 н. и 3 з.п. ф-лы, 5 ил.

Изобретение относится к радиотехнике и может быть использовано в системах контроля воздушного пространства с использованием прямых и рассеянных воздушными объектами радиосигналов, излучаемых множеством неконтролируемых и контролируемых передатчиков радиоэлектронных систем различного назначения. Достигаемый технический результат изобретения - расширение зоны действия системы обнаружения воздушных объектов, а также повышение точности оценки доплеровского сдвига частоты сигналов, отраженных от воздушного объекта, при наличии мешающих сигналов, переотраженных от объектов индустрии. Указанные результаты достигаются за счет применения антенн с широкой диаграммой направленности, а также за счет компенсации прямого сигнала передатчика и сигналов, отраженных от объектов индустрии, за счет усреднением оценки доплеровского сдвига частоты по оценкам, полученным в каждом из четырех приемных пунктов. 2 ил.

Изобретение относится к радиотехнике и может использоваться в приемниках для измерения времени прихода сигналов с двухпозиционной угловой манипуляцией

Наверх