[go: up one dir, main page]

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

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

Info

Publication number
RU2423786C1
RU2423786C1 RU2010122576/09A RU2010122576A RU2423786C1 RU 2423786 C1 RU2423786 C1 RU 2423786C1 RU 2010122576/09 A RU2010122576/09 A RU 2010122576/09A RU 2010122576 A RU2010122576 A RU 2010122576A RU 2423786 C1 RU2423786 C1 RU 2423786C1
Authority
RU
Russia
Prior art keywords
time
interference
level
predicted
time interval
Prior art date
Application number
RU2010122576/09A
Other languages
English (en)
Inventor
Владимир Феликсович Комарович (RU)
Владимир Феликсович Комарович
Роман Юрьевич Ряскин (RU)
Роман Юрьевич Ряскин
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 RU2010122576/09A priority Critical patent/RU2423786C1/ru
Application granted granted Critical
Publication of RU2423786C1 publication Critical patent/RU2423786C1/ru

Links

Images

Landscapes

  • Noise Elimination (AREA)

Abstract

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

Description

Изобретение относится к технике радиосвязи и может быть использовано для краткосрочного прогнозирования уровня помех в декаметровом диапазоне на анализируемых частотах на пунктах ионосферно-волновой и частотно-диспетчерской службы радиоцентров в условиях изменчивости помеховой обстановки.
Известен способ рангового прогнозирования на основе пропускной способности в MIMO-схеме, реализованный в патенте РФ №2378767, заявл. 22.12.2005, опубл. 10.01.2010, бюллетень №1, где осуществляется прогнозирование ранга канала для следующего временного интервала передачи, т.е осуществляется прогноз соотношения сигнал/помеха на одном из многих входах MIMO-схемы (схема со многими входами и выходами). Недостатком данного способа является отсутствие учета замираний и нестационарности помеховых процессов в канале связи при прогнозировании.
Известен также способ для точного прогнозирования отношения сигнала к помехе и коэффициента шума для улучшения функционирования системы связи, реализованный в патенте РФ №2255420, заявл. 12.09.2000, опубл. 27.06.2005, бюллетень №18, заключающийся в том, что осуществляется прогнозирование соотношения сигнала к помехе при различных скоростях передачи при условиях быстрого и медленного замираний на входе системы связи. Недостатком данного способа является тот факт, что прогнозирование соотношения сигнала к помехе осуществляется при постоянных коэффициентах фильтрации, не учитывающих нестационарность ионосферных процессов в декаметровых каналах связи.
Наиболее близким по технической сущности к заявляемому изобретению (прототипом) является способ прогнозирования уровней помех в коротковолновом диапазоне на отдельных частотах, описанный в книге Комарович В.Ф., Сосунов В.Н. Случайные радиопомехи и надежность KB связи. - М.: Связь, 1977, с.34-37. Известный способ заключается в том, что на каждой из анализируемых частот измеряют К≥2 раз уровни помех, вычисляют их усредненные значения, по результатам которых прогнозируют уровни помех на малые интервалы времени, после чего оценивают величину ошибки прогноза.
Недостатком известного способа является то, что он требует знания параметров и вида корреляционной функции случайных последовательностей уровней помех Un(t), которые определяются из статистических исследований. Однако вид корреляционной функции меняется во времени для процессов Un(t), что в конечном итоге приводит к уменьшению достоверности прогноза уровня помехи на анализируемой частоте.
Целью настоящего изобретения является разработка способа краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне, обеспечивающего повышение достоверности прогнозирования уровня помехи на анализируемой частоте на некоторое время вперед.
Поставленная цель достигается тем, что в известном способе прогнозирования уровней помех в коротковолновом диапазоне на отдельных частотах уровни помех измеряют K раз в контрольном временном интервале ΔT=TИ-T0, где T0 и TИ - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT с шагом Δt=ΔT/K, а усредненное значение уровня помех
Figure 00000001
вычисляют в пределах контрольного временного интервала, по результатам всех измерений
Figure 00000002
, где k=1, 2,…, К.
Затем методом адаптивного экспоненциального сглаживания
Figure 00000002
, измеренного в момент времени ТИ, и усредненного значения
Figure 00000001
вычисляют предварительные значения уровней прогнозируемых помех
Figure 00000003
на момент времени ΔT+Δt по формуле
Figure 00000004
причем экспоненциальное сглаживание выполняют, последовательно принимая значение коэффициента экспоненциального сглаживания αi с шагом Δα в интервале αi=0÷1, вычисленные предварительно спрогнозированные уровни помех
Figure 00000005
при всех заданных αi запоминают, после чего вычисляют среднеквадратичные ошибки
Figure 00000006
между
Figure 00000007
и дополнительно измеренным значением уровня помехи
Figure 00000008
в момент времени ΔT+Δt по формуле
Figure 00000009
Выделяют наименьшее значение среднеквадратичной ошибки
Figure 00000006
и соответствующее ему значение
Figure 00000010
, после чего в p≥2 последующих за контрольным временных интервалах, где р=1, 2…Р, при выбранном параметре сглаживания
Figure 00000011
для каждой i-ой частоты рассчитывают прогнозируемый уровень помех
Figure 00000012
на момент времени (p+1)TИ+Δt методом экспоненциального сглаживания дополнительно измеренного уровня помех
Figure 00000008
на момент времени pTИ+Δt и спрогнозированного
Figure 00000013
на момент времени pTИ+Δt
Figure 00000014
причем для временного интервала p=1 принимаем
Figure 00000015
, а после наступления временного интервала p=P вычисление прогнозируемого уровня для следующего интервала и измерение уровня
Figure 00000008
прекращают и выполняют измерения и вычисления, аналогичные контрольному временному интервалу ΔT и временным интервалам p=1…Р.
Благодаря новой совокупности существенных признаков в способе реализована возможность краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне для повышения вероятности правильного краткосрочного прогнозирования развития помеховой обстановки на анализируемых частотах.
Заявленный способ поясняется чертежами, на которых показано:
фиг.1 - временной ряд уровней помех на анализируемой частоте;
фиг.2 - результаты измерений уровней помех при различны коэффициентах адаптации на одной частоте.
В заявленном способе краткосрочного прогнозирования уровней помех в декаметровом диапазоне на каждой из N≥2 анализируемых частот производят прогноз уровней помех на заданное время вперед, для чего в контрольном временном интервале ΔT=ТИ-T0, где T0 и TИ - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT с шагом Δt=ΔT/K, как показано на фиг.1, K≥2 раз измеряют уровень помех
Figure 00000016
, где i=1, 2…, N. Затем вычисляют усредненное значение уровня помехи на каждой анализируемой частоте в пределах контрольного интервала времени
Figure 00000001
по формуле
Figure 00000017
Затем методом адаптивного экспоненциального сглаживания
Figure 00000016
, измеренного в момент времени ТИ (фиг.1), и усредненного значения
Figure 00000001
вычисляют предварительные значения уровней прогнозируемых помех
Figure 00000003
на момент времени ΔT+Δt по формуле (1), причем экспоненциальное сглаживание выполняют, последовательно принимая значение коэффициента экспоненциального сглаживания αi с шагом Δα=0,1 в интервале αi=0÷1, вычисленные предварительно спрогнозированные уровни помех
Figure 00000005
при всех заданных αi запоминают.
После предварительного прогнозирования вычисляют среднеквадратичные ошибки
Figure 00000006
между
Figure 00000005
и дополнительно измеренным значением уровня помехи
Figure 00000008
в момент времени ΔT+Δt по формуле (2).
Выделяют наименьшее значение среднеквадратичной ошибки
Figure 00000006
и соответствующее ему значение
Figure 00000010
, после чего в p≥2 последующих за контрольным временных интервалах, где р=1, 2…Р, при выбранном параметре сглаживания
Figure 00000010
для каждой i-ой частоты рассчитывают прогнозируемый уровень помех
Figure 00000012
на момент времени (p+1)TИ+Δt методом экспоненциального сглаживания дополнительно измеренного уровня помех
Figure 00000008
на момент времени pTИ+Δt и спрогнозированного
Figure 00000018
на момент времени pTИ+Δt по формуле (3), причем для временного интервала p=1 принимаем
Figure 00000015
. После наступления временного интервала p=P вычисление прогнозируемого уровня для следующего интервала и измерение уровня
Figure 00000008
прекращают и выполняют измерения и вычисления аналогичные контрольному временному интервалу ΔT и временным интервалам p=1…Р.
Таким образом, заявленный способ дает возможность осуществлять краткосрочное адаптивное прогнозирование уровней помех в декаметровом диапазоне с учетом изменяющейся помеховой обстановки.
Выявление закономерностей динамического изменения уровней помех можно производить с помощью его выравнивания и сглаживания. Экспоненциальное сглаживание - один из простейших и распространенных приемов выравнивания ряда. В его основе лежит расчет экспоненциальных средних временного ряда уровней помех (ВР) [2]. Однако в каждом конкретном случае необходимо производить адаптивный подбор коэффициента адаптивного сглаживания α, при краткосрочном прогнозирований необходимо как можно скорее отобразить изменения α и в то же время как можно лучше «чистить» ряд от случайных колебаний. С одной стороны, необходимо увеличить вес более свежих наблюдений, что может быть достигнуто повышением α, с другой стороны, для сглаживания случайных отклонений величину α нужно уменьшить. Например, результаты вычислений уровней помех при α=0,1; 0,2; 0,9 на одной частоте приведены на фиг.2. Таким образом, правильность вычисления коэффициента адаптивного сглаживания а определяет достоверность прогноза.
Процессы, протекающие в ионосфере, являются нестационарными, но на коротких интервалах времени их можно рассматривать как квазистационарные [1]. Свойства ионосферы в обычных условиях изменяются относительно медленно. Это в основном предопределяет ограниченный во времени статистически устойчивый квазистационарный характер усредненных последовательностей уровней помех. Для учета длительности квазистационарности помеховой обстановки в заявленном способе есть возможность изменения длительности контрольного временного интервала ΔT=TИ-T0, за время которого будет набираться статистика уровней помех на каждой анализируемой частоте, изменение шага Δt=ΔT/K определяет частоту выборок временного ряда, а изменение последующих за контрольным ΔT=TИ-T0 временных интервалов p, где р=1, 2…Р, - определяет длительность прогноза.
Статистическая обработка результатов реальных измерений уровней помех, проведенная на ЭВМ, учитывала реальную нестационарность процесса радиопомех. При этом была получена оценка длительности предыстории, т.е. отрезка времени ΔT, в течение которого должны вычисляться текущие средние уровни помех на каждой частоте и при котором реализуется «в среднем» максимум апостериорной вероятности. Согласно [1] оказалось, что достаточно вычислить средний уровень помехи по результатам измерений примерно за 10 мин (в зависимости от частости отсчетов, времени суток и т.д.) до момента начала прогнозирования. Таким образом, статистические исследования показали, что для повышения вероятности правильного краткосрочного прогнозирования (на 10-20 мин) необходимо непрерывно, циклично вычислять средний уровень помехи по результатам измерений на каждой частоте и производить адаптивный подбор коэффициента адаптивного сглаживания α. Отмеченное подтверждает возможность достижения сформулированного технического результата.
Литература
1. Комарович В.Ф. Случайные радиопомехи и надежность KB связи (изд. «Связь» Госкомиздата СССР, Москва, 1977, с.1-132, с.12-15).
2. Лукашин Ю.П. Адаптивные методы краткосрочного прогнозирования // Мат. Статистика для экономистов. - М.: Статистика, 1979. С.17-22, с.90-93).

Claims (3)

1. Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне, заключающийся в том, что на каждой из N≥2 анализируемых частот К≥2 раз измеряют уровни помех
Figure 00000019
, где i=1, 2 …, N, вычисляют их усредненное значение
Figure 00000020
, прогнозируют ожидаемый уровень помех на анализируемых частотах и оценивают величину ошибки прогноза, отличающийся тем, что уровни помех измеряют К раз в контрольном временном интервале ΔТ=ТИ0, где Т0 и ТИ - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT с шагом Δt=ΔT/K, а усредненное значение уровня помех
Figure 00000020

вычисляют в пределах контрольного временного интервала, по результатам всех измерений
Figure 00000019
, где k=1, 2 …, К, затем вычисляют предварительные значения уровней прогнозируемых помех
Figure 00000021
методом адаптивного экспоненциального сглаживания значения
Figure 00000019
, измеренного в момент времени ТИ и усредненного значения
Figure 00000020
, причем экспоненциальное сглаживание выполняют последовательно принимая значение коэффициента экспоненциального сглаживания αi с шагом Δα в интервале αi=0÷1, вычисленные предварительно спрогнозированные уровни помех
Figure 00000022
при всех заданных аi запоминают, после чего измеряют дополнительно уровень помех
Figure 00000023
в момент времени ΔT+Δt и вычисляют среднеквадратичные ошибки
Figure 00000024
между
Figure 00000022
и
Figure 00000023
, выделяют наименьшее значение среднеквадратичной ошибки
Figure 00000025
, и соответствующее ему значение
Figure 00000026
, после чего в p≥2 последующих за контрольным временных интервалах, где р=1, 2…Р при выбранном параметре сглаживания
Figure 00000026
, для каждой i-й частоты рассчитывают прогнозируемый уровень помех
Figure 00000027
на момент времени (р+1)ТИ+Δt методом адаптивного экспоненциального сглаживания дополнительно измеренного уровня помех
Figure 00000023
на момент времени рТИ+Δt и спрогнозированного
Figure 00000028
на момент времени рТИ+Δt, причем для временного интервала р=1 принимают
Figure 00000029
, а после наступления временного интервала р=Р вычисление прогнозируемого уровня для следующего интервала и измерение уровня
Figure 00000023
прекращают и выполняют измерения и вычисления аналогичные контрольному временному интервалу ΔТ и временным интервалам р=1…Р.
2. Способ по п.1, отличающийся тем, что в контрольном временном интервале ΔT для прогнозирования на i-й частоте предварительных значений уровней прогнозируемых помех
Figure 00000021
в момент времени ΔT+Δt с помощью экспоненциального сглаживания значения
Figure 00000030
, измеренного в момент времени ТИ и усредненного значения
Figure 00000020
рассчитывают значение предварительно прогнозируемого уровня помех
Figure 00000021
в момент времени ΔT+Δt, последовательно для различных значений коэффициента экспоненциального сглаживания, задаваемого в интервале 0<α<1 с шагом 0,1 по формуле
Figure 00000031
,
а среднеквадратичную ошибку
Figure 00000032
прогноза помех на i-й частоте вычисляют по формуле
Figure 00000033
3. Способ по п.1, отличающийся тем, что для интервалов р=1, 2, 3,…Р для прогнозирования на i-й частоте уровня ожидаемых помех
Figure 00000034
в момент времени (p+1)TИ+Δt с помощью экспоненциального сглаживания уровней
Figure 00000035
(для временного интервала р=1
Figure 00000036
) и
Figure 00000023
в момент времени рТИ+Δt рассчитывают значение прогнозируемого уровня помех
Figure 00000027
в момент времени (р+1)ТИ+Δt по формуле
Figure 00000037
.
RU2010122576/09A 2010-06-02 2010-06-02 Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне RU2423786C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2010122576/09A RU2423786C1 (ru) 2010-06-02 2010-06-02 Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2010122576/09A RU2423786C1 (ru) 2010-06-02 2010-06-02 Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне

Publications (1)

Publication Number Publication Date
RU2423786C1 true RU2423786C1 (ru) 2011-07-10

Family

ID=44740468

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2010122576/09A RU2423786C1 (ru) 2010-06-02 2010-06-02 Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне

Country Status (1)

Country Link
RU (1) RU2423786C1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2757999C1 (ru) * 2021-03-30 2021-10-25 Федеральное государственное казенное военное образовательное учреждение высшего образования "Военный ордена Жукова университет радиоэлектроники" Министерства обороны Российской Федерации (ФГКВОУВО "Военный ордена Жукова университет радиоэлектроники" МО РФ) Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне с переменной длительностью предыстории

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6567374B1 (en) * 1998-02-18 2003-05-20 Sony International (Europe) Gmbh Data and pilot mapping in an OFDM system
RU2255420C2 (ru) * 1999-09-13 2005-06-27 Квэлкомм Инкорпорейтед Система и способ для точного прогнозирования отношения сигнала к помехе и коэффициента шума для улучшения функционирования системы связи
RU2386214C2 (ru) * 2005-06-16 2010-04-10 Квэлкомм Инкорпорейтед Устойчивое прогнозирование ранга для системы mimo

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6567374B1 (en) * 1998-02-18 2003-05-20 Sony International (Europe) Gmbh Data and pilot mapping in an OFDM system
RU2255420C2 (ru) * 1999-09-13 2005-06-27 Квэлкомм Инкорпорейтед Система и способ для точного прогнозирования отношения сигнала к помехе и коэффициента шума для улучшения функционирования системы связи
RU2386214C2 (ru) * 2005-06-16 2010-04-10 Квэлкомм Инкорпорейтед Устойчивое прогнозирование ранга для системы mimo

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
КОМАРОВИЧ В.Ф., СОСУНОВ В.Н. Случайные помехи и надежность KB связи. - М.: Связь, 1977, с.34-37. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2757999C1 (ru) * 2021-03-30 2021-10-25 Федеральное государственное казенное военное образовательное учреждение высшего образования "Военный ордена Жукова университет радиоэлектроники" Министерства обороны Российской Федерации (ФГКВОУВО "Военный ордена Жукова университет радиоэлектроники" МО РФ) Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне с переменной длительностью предыстории

Similar Documents

Publication Publication Date Title
EP3259613B1 (en) System and method for estimating time of arrival (toa)
CN112148557B (zh) 一种性能指标实时预测方法、计算机设备及存储介质
CN114298166B (zh) 一种基于无线通信网络的频谱可用性预测方法和系统
CN106054617B (zh) 网络控制系统前向通道随机时延的在线预测方法
US10241188B2 (en) Method and apparatus for obtaining time of arrival TOA when mobile terminal is located
Matisko et al. Optimality tests and adaptive Kalman filter
CN111279209A (zh) 用于改进定位模型的指纹数据预处理方法
RU2423786C1 (ru) Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне
JP2015156642A5 (ru)
CN115390112A (zh) 电力系统北斗定位信号滤波方法、装置、设备及存储介质
CN119363211B (zh) 无人机测控信号的捕获方法、相关装置及计算机存储介质
RU2757999C1 (ru) Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне с переменной длительностью предыстории
EP1692537B1 (en) Positioning method and apparatus
CN110856249B (zh) 一种基于基站的定位方法、装置及存储介质
CN115166650B (zh) 一种雷达信号识别与参数估计方法和系统
RU2394371C1 (ru) Устройство для определения оптимальных рабочих частот ионосферного радиоканала
Kartashov et al. Selection of the model for varying speed of sound for the optimal linear filter of atmosphere radio acoustic sounding systems
CN114705181A (zh) 一种基于数据自适应的星载设备健康预测方法
CN120769361B (zh) 基于地基导航系统的多径toa计算方法、系统及电子设备
Alawieh et al. Stochastic modeling of pseudolite clock errors using enhanced AR methods
EP2647249B1 (en) Method and device for noise floor estimation
CN114548459B (zh) 票务数据调控方法和系统及计算机可读存储介质
CN116418368B (zh) 考虑类型匹配偏差的电力线载波信道噪声预测方法和系统
CN103139787A (zh) 一种路测数据的干扰抑制方法及相关装置
EP4495633A1 (en) Method for estimating jamming in a global navigation satellite system receiver

Legal Events

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

Effective date: 20120603