RU2760744C1 - Способ измерения частоты модуляции - Google Patents

Способ измерения частоты модуляции Download PDF

Info

Publication number
RU2760744C1
RU2760744C1 RU2020139965A RU2020139965A RU2760744C1 RU 2760744 C1 RU2760744 C1 RU 2760744C1 RU 2020139965 A RU2020139965 A RU 2020139965A RU 2020139965 A RU2020139965 A RU 2020139965A RU 2760744 C1 RU2760744 C1 RU 2760744C1
Authority
RU
Russia
Prior art keywords
frequency
spectrum
radio
power spectrum
frequencies
Prior art date
Application number
RU2020139965A
Other languages
English (en)
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 RU2020139965A priority Critical patent/RU2760744C1/ru
Application granted granted Critical
Publication of RU2760744C1 publication Critical patent/RU2760744C1/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • G01R23/165Spectrum analysis; Fourier analysis using filters

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

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

Description

Изобретение относится к области измерения параметров радиосигналов и может быть использовано в системах радиоконтроля за использованием радиочастотного спектра.
Известен способ измерения частоты модуляции методом детектирования, в соответствии с которым частотно-модулированное колебание преобразуют с несущей на промежуточную частоту, детектируют и измеряют частоту продетектированного сигнала [1. Измерения в технике связи / под ред. М.А. Ракк. – М.: ГОУ «Учебно-методический центр по образованию на железнодорожном транспорте». 2008, с. 492-495].
Измерения данным способом сопровождаются погрешностями, обусловленными нелинейностью частотного детектирования, при этом необходима информации о несущей частоте и девиации частоты, измерение этих параметров увеличивает трудоёмкость процесса радиоконтроля.
Из известных наиболее близким к предлагаемому является способ измерения частоты модуляции по спектру радиосигнала. Способ основан на свойствах энергетического спектра частотно модулированного тоном колебания. Спектр «состоит из бесконечного числа боковых частот, расположенных попарно симметрично относительно несущей частоты и отличающихся от последней» на величину кратную частоте модуляции [2. Гоноровский И.С. Радиотехнические спектры и сигналы. – М.: Радио и связь, 1986, с 86].
Учёт закономерностей энергетического спектра, как теоретического объекта, который определяется при бесконечном времени наблюдения на непрерывном и бесконечном множестве частот принципиально позволяет определять частоту модуляции. Например, как частотное расстояние от несущей частоты до ближайшей гармоники. Недостаток в том, что в этом случае необходима информация о несущей частоте или её необходимо дополнительно измерить с соответствующими затратами.
Однако главный недостаток способа-прототипа состоит в том, что в нём не определены конкретные операции над сигналом, которые необходимо выполнять в процессе измерений при практических ограничениях, в частности по частоте и времени.
Технической задачей настоящего изобретения является определение частоты модуляции по спектру радиосигналов.
Поставленная техническая задача решается за счёт того, что в известном способе измерения частоты модуляции, включающем приём частотно-модулированного радиосигнала и преобразование его в спектр мощности, новым является то, что спектр мощности фильтруют, исключая дискретные составляющие по уровню менее заданного порога, определяют, по мере увеличения, частоту каждого локального максимума отфильтрованного спектра и уточняют её путём усреднения частот локального максимума и соседних с ним дискретных составляющих с весом равным значениям спектра мощности этих составляющих, затем по соседним уточнённым частотам определяют попарные разности и их среднее медианное значение.
В настоящем изобретении выполнен переход от теоретического энергетического спектра к его оценке на ограниченном временном интервале наблюдения и в полосе частот приёма: спектру мощности.
Спектр мощности определяется на дискретном ограниченном множестве частот, например, в виде квадратов модулей коэффициентов дискретного преобразования Фурье. Ему присущи основные свойства энергетического спектра, однако есть и особенности.
Дискретный характер спектра мощности позволяет определять локальные максимумы гармоник с шагом, равным дискрете спектрального анализа, обратно пропорциональной времени наблюдения (длительности реализации) сигнала. При ограничении времени наблюдения происходит размывание гармоник, часть энергии просачивается в соседние частотные позиции, что снижает разрешающую способность по частоте, но позволяет уточнить положение локального максимума на континуальном множестве частот весовым усреднением.
Другая особенность состоит в том, что воздействие шумов сопровождается образованием значительного количества дополнительных локальных максимумов с нерегулярным положением. В обеспечение помехозащищённости выполняется пороговая фильтрация спектра мощности.
Из цитированного ранее свойства энергетического спектра следует, что расстояния между соседними гармониками равны частоте модуляции. Однако часть гармоник в энергетическом спектре отсутствует, например, на несущей частоте при индексе модуляции равном 2,4; 5,52; 8,65; 11,79 [3. Измерения в технике связи/ Под ред. М.А. Ракк. – М.: ГОУ «Учебно-методический центр по образованию на железнодорожном транспорте». 2008, с. 490-491]. Поэтому применительно к фильтрованному спектру мощности можно утверждать, что расстояния между соседними, по мере увеличения частоты, гармониками кратны частоте модуляции.
Учёт свойства спектра в такой формулировке не предполагает знания несущей частоты сигнала, что упрощает процесс измерения, но возникает неопределённость кратности. Более детальный анализ показывает, что кратность имеет преимущественно единичное значение, что и позволяет преодолеть неопределённость путём медианной оценки попарных разностей частот локальных максимумов.
Учёт указанных закономерностей в соответствии с предложенными новыми действиями, условиями и порядком их выполнения, позволяет решить поставленную техническую задачу: обеспечить определение частоты модуляции по спектру радиосигналов.
Указанные преимущества и особенности настоящего изобретения поясняются вариантом его осуществления со ссылками на прилагаемые фигуры.
На фиг.1 показана структурная схема аппаратуры радиоконтроля, реализующая предложенный способ;
на фиг.2 – спектр мощности сигнала с синусоидальной модуляцией;
на фиг.3 – зависимость относительной средней квадратичной погрешности измерений от истинного значения частоты модуляции.
Аппаратура радиоконтроля фиг.1 содержит последовательно соединённые антенну 1, радиоприёмное устройство 2, анализатор спектра 3, пороговый элемент 4, устройство определения максимума 5, через первый вход блок весовой обработки 6, устройство вычитания 7 и блок медианной оценки 8. Выход анализатора спектра 3 подключен ко второму входу блока весовой обработки 6. Выходом аппаратуры радиоконтроля является выход функционального преобразователя 7.
Проиллюстрируем предлагаемое решение на примере выполнения способа с представлением анализируемого сигнала с синусоидальной частотной модуляцией в виде дискретных отсчётов
Figure 00000001
, (1)
где
Figure 00000002
– номер временного отсчёта при общем количестве
Figure 00000003
,
Figure 00000004
– амплитуда,
Figure 00000005
– несущая (промежуточная) частота,
Figure 00000006
– девиация частоты,
Figure 00000007
– частота модуляции,
Figure 00000008
– фаза модулирующей функции,
Figure 00000009
– шум приёма.
Временные параметры выражены в единицах периода дискретизации по времени
Figure 00000010
, частотные в бинах, единицах периода дискретизации по частоте
Figure 00000011
.
Преобразование (1) может выполняться с помощью цифрового радиоприёмного устройства [4. Побережский Е.С. Цифровые радиоприёмные устройства. – М.: Радио и связь, 1987, с. 62-73], при этом частота дискретизации
Figure 00000012
устанавливается в соответствии с теоремой Котельникова.
Возможен анализ сигнала радиопередатчика как при приёме его сигнала на антенну 1 (измерение по эфиру), так и при непосредственном подключении выхода радиопередатчика к входу радиоприемного устройства 2 (контактное измерение).
Полоса частот пропускания радиоприёмного устройства 2 и спектрального анализа анализатора спектра 3 должна включать полосу частот, занимаемых сигналом.
Дискретизированный сигнал (1) преобразуют в спектр мощности в анализаторе спектра 3, при этом его умножают на весовое окно, осуществляют дискретное преобразование Фурье и определяют квадраты модулей полученных коэффициентов Фурье [5. Херрис Ф. Использование окон при гармоническом анализе методом дискретного преобразования Фурье. ТИИР, 1978, т. 66 № 1, с. 60-96]
Figure 00000013
, (2)
где
Figure 00000014
– весовое окно Блэкмана-Хэрриса,
Figure 00000015
– номер спектральной составляющей при общем количестве
Figure 00000016
.
Весовая обработка выполняется для снижения погрешностей, обусловленных эффектом Гиббса, взаимного просачивания спектральных составляющих непериодических процессов при ограничении времени приёма сигнала. С учётом этого необходимо соблюдать также следующее условие: на интервале наблюдения должно укладываться не менее трёх периодов модуляции, то есть
Figure 00000017
.
Полученный спектр мощности (2) фильтруют, исключая спектральные составляющие по уровню менее заданного порога, с помощью порогового элемента 4:
Figure 00000018
(3)
Порог определяется прямо пропорционально дисперсии шума
Figure 00000019
и обратно пропорционально числу отсчётов, рекомендуемое его значение
Figure 00000020
. (4)
На фиг.2 показан спектр мощности сигнала с частотой модуляцией.
Установлены следующие параметры: отношение амплитуды сигнала к среднему квадратическому значению шума 10, частота модуляции и девиация частоты 0,6 кГц, несущая (промежуточная) частота сигнала 50 кГц, частота дискретизации 200 кГц, соответственно период дискретизации по времени 5 мкс, число временных отсчётов 1024, соответственно время наблюдения 5,1 мс, а период дискретизации по частоте 195 Гц и на интервале наблюдения укладывается 3,1 периода модуляции.
Пунктиром отмечен порог фильтрации в соответствии с формулой (4). Значения спектра и порога нормированы на максимальное значение спектра.
В соответствии с фиг.2 при пороговой фильтрации устраняются шумовые спектральные элементы ниже порогового уровня.
В устройстве 5 определения максимума определяют частоту каждого локального максимума отфильтрованного спектра. На фиг.2 таких максимумов пять, их значения отмечены жирными точками. Частоты максимумов определяют в порядке возрастания:
Figure 00000021
, где x – знак транспонирования. В общем случае фиксации
Figure 00000022
максимумов определяют совокупность всех их частот:
Figure 00000023
, где
Figure 00000024
– номер максимума.
Полученные значения частот уточняют в блоке 6 весовой обработки путём усреднения частот каждого локального максимума и соседних с ним дискретных составляющих с весом равным значениям спектра мощности этих составляющих
Figure 00000025
. (5)
Частоты максимумов и весовые коэффициенты
Figure 00000026
поступают соответственно с устройства 5 и анализатора спектра 3.
В устройстве вычитания 7 по соседним уточнённым частотам локальных максимумов определяют попарные разности
Figure 00000027
, где
Figure 00000028
, а в блоке медианной оценки 8 среднее медианное значение этих разностей
Figure 00000029
с переводом измеренной частоты модуляции, выраженной в бинах, в систему единиц СИ умножением на величину бина:
Figure 00000030
.
Точность измерения частоты модуляции предлагаемым способом иллюстрируется фиг.3, где приведена зависимость относительной средней квадратичной погрешности измерения частоты модуляции
Figure 00000031
от её истинного значения. Сплошная линия для объёма выборки
Figure 00000032
, пунктирная при увеличении её в два раза. Установлена девиация частоты 10 кГц, фаза модулирующей функции равновероятна в пределах области её определения
Figure 00000033
, число статистических испытаний в каждой точке, выделенной на фиг.3 кружком,
Figure 00000034
.
В соответствии с фиг.3 измерение частоты модуляции обеспечивается в достаточно широких пределах от 0,3 (0,6) кГц до 40 кГц с относительной погрешностью менее 1%. Ограничение верхнего измеряемого значения частоты модуляции обусловлено расширением спектра сигнала за пределы полосы частот спектрального анализа. Ограничение снизу вытекает из условия трёхкратного превышения интервала наблюдения периода модуляции. Пределы измерений могут быть расширены увеличением частоты дискретизации и времени наблюдения.
Таким образом, предложенное техническое решение обеспечивает определение частоты модуляции по спектру радиосигналов.

Claims (1)

  1. Способ измерения частоты модуляции, включающий приём частотно-модулированного радиосигнала и преобразование его в спектр мощности, отличающийся тем, что спектр мощности фильтруют, исключая дискретные составляющие по уровню менее заданного порога, определяют, по мере увеличения, частоту каждого локального максимума отфильтрованного спектра и уточняют её путём усреднения частот локального максимума и соседних с ним дискретных составляющих с весом, равным значениям спектра мощности этих составляющих, затем по соседним уточнённым частотам определяют попарные разности и их среднее медианное значение.
RU2020139965A 2020-12-04 2020-12-04 Способ измерения частоты модуляции RU2760744C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2020139965A RU2760744C1 (ru) 2020-12-04 2020-12-04 Способ измерения частоты модуляции

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2020139965A RU2760744C1 (ru) 2020-12-04 2020-12-04 Способ измерения частоты модуляции

Publications (1)

Publication Number Publication Date
RU2760744C1 true RU2760744C1 (ru) 2021-11-30

Family

ID=79174023

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2020139965A RU2760744C1 (ru) 2020-12-04 2020-12-04 Способ измерения частоты модуляции

Country Status (1)

Country Link
RU (1) RU2760744C1 (ru)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8416975D0 (en) * 1983-08-26 1984-08-08 Gen Electric Determining frequency of time varying signals
EP0229443A1 (en) * 1985-09-16 1987-07-22 Westinghouse Brake And Signal Holdings Limited Method and apparatus for identifying the components of a signal
RU2144733C1 (ru) * 1993-06-07 2000-01-20 Алькатель Н.В. Пакет канала сигнализации для системы связи с опорным сигналом, модулированным по закону, зависящему от времени
RU2668342C2 (ru) * 2017-03-10 2018-09-28 Акционерное общество "ИРКОС" Способ измерения частотного сдвига между радиосигналами

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8416975D0 (en) * 1983-08-26 1984-08-08 Gen Electric Determining frequency of time varying signals
EP0229443A1 (en) * 1985-09-16 1987-07-22 Westinghouse Brake And Signal Holdings Limited Method and apparatus for identifying the components of a signal
RU2144733C1 (ru) * 1993-06-07 2000-01-20 Алькатель Н.В. Пакет канала сигнализации для системы связи с опорным сигналом, модулированным по закону, зависящему от времени
RU2668342C2 (ru) * 2017-03-10 2018-09-28 Акционерное общество "ИРКОС" Способ измерения частотного сдвига между радиосигналами

Similar Documents

Publication Publication Date Title
US7400387B2 (en) Method for monitoring an optical transmission line, corresponding measuring device and optical transmitter
KR100978652B1 (ko) 스펙트럼 검출 알고리듬 및 방법
US6101218A (en) Leakage power measurement apparatus and method
JPH08248070A (ja) 周波数スペクトル分析装置
CN101334434A (zh) 利用小波变换提取电磁泄漏信号的电磁环境测试系统
KR102169078B1 (ko) 광대역 측정 신호들을 측정하기 위한 방법 및 측정 장치
Bao et al. A histogram-based segmentation method for wideband spectrum sensing in cognitive radios
CN101551419A (zh) 一种相位噪声检测方法及系统
CN102809686A (zh) 扫频源瞬时频率检测方法与系统
US8175558B2 (en) Radio communication apparatus and method
CN104251675B (zh) 一种多目标微变形实时遥测方法与系统
RU2760744C1 (ru) Способ измерения частоты модуляции
CN109001729A (zh) 太赫兹成像中调频连续波线性度实时校准方法及其系统
KR102122758B1 (ko) 생체신호 측정용 레이더에서 레이더 수집신호의 랜덤노이즈 제거방법 및 그 장치
CN102694593A (zh) 一种光无源器件的谱特性的测试方法
CN106353776A (zh) 一种非理想信道下测量零值无偏的频域抗干扰系统及其方法
RU2758342C1 (ru) Спектральный способ измерения девиации частоты
CN114217301B (zh) 一种高精度边坡监测雷达目标检测及测距方法
CN109521269A (zh) 一种幅度调制信号数字化测频方法
RU2394371C1 (ru) Устройство для определения оптимальных рабочих частот ионосферного радиоканала
CN107783084A (zh) 恒虚警检测及数据处理平台
CN113341219B (zh) 恒频交流供电系统频率调制幅度测量方法及装置
CN101404546A (zh) 基于信号矩特征的频谱占用度测量方法
CN107783089B (zh) 单元平均选小门限检测平台
CN110375939A (zh) 一种航天器结构振动试验健康状态评估的改进方法