RU2485526C2 - Способ оценки параметров и демодуляции случайных сигналов - Google Patents

Способ оценки параметров и демодуляции случайных сигналов Download PDF

Info

Publication number
RU2485526C2
RU2485526C2 RU2011133116/28A RU2011133116A RU2485526C2 RU 2485526 C2 RU2485526 C2 RU 2485526C2 RU 2011133116/28 A RU2011133116/28 A RU 2011133116/28A RU 2011133116 A RU2011133116 A RU 2011133116A RU 2485526 C2 RU2485526 C2 RU 2485526C2
Authority
RU
Russia
Prior art keywords
signal
amplitude
phase
parameters
demodulation
Prior art date
Application number
RU2011133116/28A
Other languages
English (en)
Other versions
RU2011133116A (ru
Inventor
Петр Васильевич Семкин
Александр Витальевич Кузовников
Виктор Григорьевич Сомов
Original Assignee
Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" filed Critical Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва"
Priority to RU2011133116/28A priority Critical patent/RU2485526C2/ru
Publication of RU2011133116A publication Critical patent/RU2011133116A/ru
Application granted granted Critical
Publication of RU2485526C2 publication Critical patent/RU2485526C2/ru

Links

Images

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

Изобретение относится к радиотехнике и может быть использовано для целей радиоконтроля, радиомониторинга, определения характеристик источников радиоизлучения. Способ основан на совмещении вейвлет-преобразования и статистической оценки параметров сигнала. Согласно изобретению оценка параметров и демодуляция сигналов проводится в два этапа с применением непрерывного вейвлет-преобразования с фиксированным значением масштаба в качестве первичной обработки, по результатам которой определяют моменты времени тактовых переходов амплитуды и фазы, зная которые, находят промежутки сигнала, на которых значения амплитуды и фазы постоянны. На втором этапе на этих промежутках используют статистический метод максимального правдоподобия для вычисления оценок амплитуды и фазы и демодуляции сигнала. Вместе с тем, подсчитав количество тактовых переходов, определяют техническую скорость передачи. При этом способ работоспособен при отношениях сигнал/шум более 5 дБ. Технический результат заключается в возможности определения технической скорости передачи и демодуляция сигналов с амплитудной, фазовой и амплитудно-фазовой модуляцией в отсутствии априорной информации о параметрах входного сигнала. 4 ил.

Description

Изобретение относится к радиотехнике и может быть использовано для целей радиоконтроля, радиомониторинга, определения параметров источников радиоизлучения. Способ основан на совмещении вейвлет-обработки и статистической оценки параметров сигналов.
В настоящее время известны различные способы оценки параметров и демодуляции сигналов с неизвестной структурой (см. патенты Российской Федерации №№2386165, 2341024, 2351005). Известные способы используют методы обработки сигналов, основанные на обработке в нейронных сетях, алгоритмах спектрального анализа (MUSIC, EV) и корреляционном анализе сигналов.
Существует способ обнаружения и определения частотных и пространственных параметров сигналов радиоэлектронных средств. Заявленное изобретение относится к области радиотехники и заключается в приеме сигналов плоской антенной решеткой, входящей в состав многоэлементной антенной системы, их усилении, перемножении и низкочастотной фильтрации выходного напряжения каждого антенного элемента с выходными напряжениями всех остальных антенных элементов и представлении результатов перемножения и фильтрации в виде корреляционной матрицы сигналов, задержке этих сигналов на время задержки τ3, значение которого обеспечивает выполнение следующего условия: вероятность изменения числа воздействующих на многоэлементную антенную систему за время τ3 сигналов более чем на один пренебрежимо мала, при поэлементном вычитании сигналов текущей и задержанной матриц сигналов и представлении результатов вычитания в виде разностной корреляционной матрицы сигнала, которая является ненулевой матрицей в случае изменения количества излучающих радиоэлектронных средств за период времени, равный τ3, и представляет собой корреляционную матрицу появившегося или исчезнувшего сигнала, и определении по виду разностной корреляционной матрицы сигнала значений рабочей частоты и направления прихода сигнала радиоэлектронного средства [Патент RU 2341024 C1, СПОСОБ ОБНАРУЖЕНИЯ РАДИОЭЛЕКТРОННЫХ СРЕДСТВ, опубликованный 10.12.2008].
Недостатком этого способа является невозможность определения типа модуляции обнаруженного сигнала и его демодуляции.
Также известен способ оценки параметров сигнала, основанный на спектральном анализе сигнала, выполняемом в два последовательных этапа, различающихся алгоритмами обработки сигнала и построением псевдоспектров сигнала на каждом этапе. На первом этапе выполняется вычисление и анализ собственных чисел и собственных векторов корреляционной матрицы сигнала (алгоритм MUSIC). На втором этапе выполняется вычисление и анализ весовых коэффициентов, обратно пропорциональных соответствующим собственным числам корреляционной матрицы сигнала (алгоритм EV). Устройство оценки параметров сигнала включает совокупность блоков спектрального анализа сигналов, включающий блок построения псевдоспектра, а также аналого-цифровой преобразователь и блок оценки частоты сигнала [Патент RU 2351005 C1, СПОСОБ ОЦЕНКИ ПАРАМЕТРОВ СИГНАЛА И УСТРОЙСТВО ДЛЯ РЕАЛИЗАЦИИ СПОСОБА (ВАРИАНТЫ), опубликованный 27.03.2009].
Недостатком данного способа является невозможность определения типа модуляции сигнала. Кроме того, в предложенном изобретении не решается задача оценки амплитуды и мгновенной фазы сигнала, то есть отсутствует возможность его демодуляции.
Наиболее близким к предлагаемому способу является способ определения структуры и демодуляции сигнала с неизвестной структурой, основанный на двухэтапной обработке сигнала с использованием первичной вейвлет-обработки для грубой оценки параметров сигнала (максимальная, минимальная амплитуда и частота, присутствие фазовых искажений) и вторичного анализа в нейронной сети для точного определения параметров сигнала, в которой проводится параллельная обработка сигнала, и автоматическая подстройка под каждый тип входного сигнала [Патент RU 2386165 C2, СПОСОБ ОПРЕДЕЛЕНИЯ СТРУКТУРЫ И ДЕМОДУЛЯЦИИ СИГНАЛА С НЕИЗВЕСТНОЙ СТРУКТУРОЙ, опубликованный 10.04.2010]. Данное изобретение выбрано в качестве прототипа.
Недостатком указанного способа является использование сложной обработки (многомасштабного вейвлет-преобразования, обработки в нейронной сети), которая требует значительных вычислительных и временных затрат.
Задачами настоящего изобретения является оценка параметров случайного сигнала и его демодуляция и устранение недостатков прототипа - упрощение обработки сигнала и уменьшение вычислительных затрат. Эту задачу в предлагаемом способе решают в два этапа, последовательно применяя непрерывное вейвлет-преобразование исследуемого сигнала с одним фиксированным значением масштаба вейвлет-функции, а затем проводя оценку амплитуды и фазы методом максимального правдоподобия, как это показано на блок-схеме на фиг.1.
Выбор непрерывного вейвлет-преобразования в качестве первичной обработки обусловлен тем, что оно не требует априорной информации о сигнале.
Непрерывное вейвлет-преобразование (НВП) сигнала определяют выражением:
Figure 00000001
где:
W(a, τ) - непрерывное вейвлет-преобразование сигнала,
а - масштаб вейвлет-функции,
τ - сдвиг вейвлет-функции по оси времени,
S(t) - входной сигнал,
t - текущее значение времени,
Ψ(t) - базисный вейвлет,
* - обозначает комплексное сопряжение.
Главным преимуществом использования вейвлет-преобразования является то, что оно позволяет исследовать во времени различные частотные составляющие принятого сигнала. Непрерывное вейвлет-преобразование дает набор коэффициентов, зависящих от масштаба вейвлет-функции а и параметра сдвига τ. За изменение рассматриваемой частоты отвечает параметр масштаба а: чем больше масштаб, тем ниже исследуемая частотная составляющая. Параметр сдвига τ отвечает за развертку по времени. Поскольку модулирующий сигнал всегда более низкочастотный, чем несущее колебание, в предлагаемом способе используют коэффициент вейвлет-преобразования с большим значением параметра а, соответствующим низкочастотным составляющим сигнала.
Экспериментальным путем выяснено, что для решаемой задачи наиболее оптимален вейвлет Добеши 8-го порядка. При вычислении вейвлет-преобразования используют только одно значение масштаба a, равное 560, что значительно уменьшает вычислительные затраты по сравнению с многомасштабным преобразованием, используемым в прототипе.
Достигаемым техническим результатом применения предварительной вейвлет-обработки является определение моментов времени ti, в которые происходят тактовые переходы фазы и амплитуды сигнала. Моменты тактовых переходов амплитуды и фазы по времени соответствуют максимумам коэффициента вейвлет-преобразования как это показано на фиг.2. Подсчитав количество максимумов коэффициента вейвлет-преобразования, определяют техническую скорость передачи (скорость манипуляции).
В промежутках между тактовыми переходами (t1, t2), (t2, t3), …, (ti-1, ti), амплитуда и фаза сигнала постоянны, поэтому на этих промежутках применяют метод максимального правдоподобия для оценки априорно неизвестных значений амплитуды и фазы. Для этого сигнал дискретизируют с частотой дискретизации Fs. В предлагаемом способе для достижения заявленного результата частоту дискретизации выбирают исходя из условия:
Figure 00000002
где fн - несущая частота исследуемого сигнала.
Оценки максимального правдоподобия амплитуды и фазы проводят по полученным дискретным выборкам, используя формулы:
Figure 00000003
Figure 00000004
где
Figure 00000005
Figure 00000006
Figure 00000007
N - количество отсчетов на участке (ti-1, ti),
S[n], n=1, 2, …, N - выборка отсчетов сигнала.
В результате для каждого интервала (ti-1, ti) получают две оценки
Figure 00000008
,
Figure 00000009
. По полученным зависимостям
Figure 00000010
,
Figure 00000011
определяют тип модуляции: находят количество значений, которые принимают амплитуда и фаза, после чего строят диаграмму рассеяния (см. фиг.3), по которой определяют тип модуляции. Например, одно значение амплитуды и два значения фазы соответствуют фазовой модуляции ФМ-2, одно значение амплитуды и четыре значения фазы соответствуют модуляции ФМ-4, четыре значения амплитуды и четыре значения фазы соответствуют модуляции КАМ-16 и так далее. В то же время, зависимости
Figure 00000012
,
Figure 00000013
в совокупности с диаграммой рассеяния представляют собой оценки демодулированной исходной битовой последовательности.
Достигаемым техническим результатом предлагаемого способа является определение технической скорости передачи и типа модуляции случайных сигналов с амплитудной, фазовой и амплитудно-фазовой манипуляцией в отсутствии априорной информации о сигнале, кроме того, данный способ позволяет демодулировать сигнал с априорно неизвестными параметрами. Данный результат получен путем последовательного применения вейвлет-обработки и метода максимального правдоподобия. Предлагаемый способ обеспечивает результат и при воздействии шумов на исследуемый сигнал (см. фиг.4). Пороговое значение отношения сигнал/шум, при котором выполняется заявленный результат, составляет 5 дБ. При этом обработка сигнала проводится в режиме максимально приближенном к реальному времени.

Claims (1)

  1. Способ оценки параметров и демодуляции случайных сигналов, основанный на двухэтапной обработке сигнала, при которой в качестве первичной используют вейвлет-обработку, отличающийся тем, что первичную вейвлет-обработку проводят с использованием только одного значения масштаба базисной вейвлет-функции, по результатам первичной вейвлет-обработки определяют моменты времени, в которые происходят тактовые переходы амплитуды и фазы, промежутки на которых амплитуда и фаза постоянны, а также техническую скорость передачи, после чего проводят вторичную обработку с использованием статистического метода максимального правдоподобия для оценки амплитуды и фазы, построения диаграммы рассеяния, определения типа модуляции и демодуляции сигнала.
RU2011133116/28A 2011-08-05 2011-08-05 Способ оценки параметров и демодуляции случайных сигналов RU2485526C2 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2011133116/28A RU2485526C2 (ru) 2011-08-05 2011-08-05 Способ оценки параметров и демодуляции случайных сигналов

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2011133116/28A RU2485526C2 (ru) 2011-08-05 2011-08-05 Способ оценки параметров и демодуляции случайных сигналов

Publications (2)

Publication Number Publication Date
RU2011133116A RU2011133116A (ru) 2013-02-10
RU2485526C2 true RU2485526C2 (ru) 2013-06-20

Family

ID=48786631

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2011133116/28A RU2485526C2 (ru) 2011-08-05 2011-08-05 Способ оценки параметров и демодуляции случайных сигналов

Country Status (1)

Country Link
RU (1) RU2485526C2 (ru)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020003843A1 (en) * 1999-05-28 2002-01-10 Wj Communications, Inc. Method and apparatus for high data rate wireless communications over wavefield spaces
RU2227306C2 (ru) * 2002-05-16 2004-04-20 Череповецкий военный инженерный институт радиоэлектроники Компенсационный способ измерения разности фаз случайных сигналов
US20080119716A1 (en) * 2006-05-17 2008-05-22 Olga Boric-Lubecke Determining presence and/or physiological motion of one or more subjects with quadrature doppler radar receiver systems
RU2386165C2 (ru) * 2008-06-16 2010-04-10 Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнева" Способ определения структуры и демодуляции сигнала с неизвестной структурой

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020003843A1 (en) * 1999-05-28 2002-01-10 Wj Communications, Inc. Method and apparatus for high data rate wireless communications over wavefield spaces
RU2227306C2 (ru) * 2002-05-16 2004-04-20 Череповецкий военный инженерный институт радиоэлектроники Компенсационный способ измерения разности фаз случайных сигналов
US20080119716A1 (en) * 2006-05-17 2008-05-22 Olga Boric-Lubecke Determining presence and/or physiological motion of one or more subjects with quadrature doppler radar receiver systems
RU2386165C2 (ru) * 2008-06-16 2010-04-10 Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнева" Способ определения структуры и демодуляции сигнала с неизвестной структурой

Also Published As

Publication number Publication date
RU2011133116A (ru) 2013-02-10

Similar Documents

Publication Publication Date Title
JP6026531B2 (ja) レーダー用デジタル受信機を用いるレーダーパルス検出
US9660772B2 (en) Detecting and processing weak signals using an array of antennas
CN102377452B (zh) 高速取样有限精度量化脉冲超宽带信号到达时间估计方法
RU2549207C2 (ru) Устройство обнаружения шумовых гидроакустических сигналов на основе квадратурного приемника
CN104218973A (zh) 基于Myriad滤波的跳频信号参数估计方法
US9244156B1 (en) Orthogonal polarization signal agnostic matched filter
CN105866748B (zh) 一种基于检测先验的固定窗长恒虚警检测方法
RU2485526C2 (ru) Способ оценки параметров и демодуляции случайных сигналов
RU2731130C1 (ru) Способ многоканального обнаружения источника шумоподобного радиосигнала
RU2282209C1 (ru) Способ и устройство обнаружения сложных широкополосных частотно-модулированных сигналов с фильтрацией в масштабно-временной области
US9157985B1 (en) Signal agnostic matched filter parameter estimator
CN108718223B (zh) 一种非合作信号的盲频谱感知方法
Kolchev et al. Use of hazard function for signal detection on ionograms
RU2608363C1 (ru) Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу
CN105764153A (zh) 一种上行随机接入的签名序列检测方法及装置
RU2292558C1 (ru) Способ определения энергетического спектра шумового электрического сигнала
RU2395158C1 (ru) Способ цифровой фильтрации сигналов
KR101633187B1 (ko) 레이더의 도플러 주파수 추정 장치 및 방법
RU2559734C1 (ru) Способ определения параметров модели замирания радиоканала по закону райса по информационному многочастотному сигналу
RU2548032C2 (ru) Способ оценивания отношения сигнал/шум при использовании сигналов с фазовой модуляцией
Pace et al. Nyquist folding analog-to-information receiver: Autonomous information recovery using quadrature mirror filtering
RU2555194C1 (ru) Способ обработки гидроакустического сигнала шумоизлучения объекта
RU2768238C1 (ru) Способ двухэтапной селекции спектральных компонент радиосигналов в многоканальной аппаратуре радиомониторинга
RU2285937C2 (ru) Способ обнаружения и определения координат источника радиоизлучения
RU2700798C2 (ru) Устройство обнаружения широкополосных полигармонических сигналов на фоне аддитивной помехи

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20180806