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

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

 

Область техники

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

Уровень техники

Известен способ рангового прогнозирования на основе пропускной способности в MIMO-схеме, реализованный в патенте РФ №2378767, заявл. 22.12.2005, опубл. 10.01.2010, бюллетень №1, где осуществляется прогнозирование ранга канала для следующего временного интервала передачи, т.е осуществляется прогноз соотношения сигнал/помеха на одном из многих входах MIMO-схемы (схема со многими входами и выходами). Недостатком данного способа является отсутствие учета замираний и нестационарности помеховых процессов в канале связи при прогнозировании.

Известен также способ для точного прогнозирования отношения сигнала к помехе и коэффициента шума для улучшения функционирования системы связи, реализованный в патенте РФ №2255420, заявл. 12.09.2000, опубл. 27.06.2005, бюллетень №18, заключающийся в том, что осуществляется прогнозирование соотношения сигнала к помехе при различных скоростях передачи при условиях быстрого и медленного замираний на входе системы связи. Недостатком данного способа является тот факт, что прогнозирование соотношения сигнала к помехе осуществляется при постоянных коэффициентах фильтрации, не учитывающих нестационарность ионосферных процессов в декаметровых каналах связи.

Наиболее близким по технической сущности к заявляемому изобретению (прототипом) является способ прогнозирования уровней помех в коротковолновом диапазоне на отдельных частотах, реализованный в патенте РФ №2423786, заявл. 02.06.2010, опубл. 10.07.2011, бюллетень №19, заключающийся в том, что на каждой из анализируемых частот измеряют уровни помех, вычисляют их усредненные значения, по результатам которых прогнозируют уровни помех на малые интервалы времени, после чего оценивают величину ошибки прогноза, причем прогнозирование осуществляется при различных коэффициентах адаптации и по результатам предварительного прогноза выбираются коэффициенты адаптации для каждой частоты, обеспечивающие минимальные среднеквадратичные ошибки прогноза.

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

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

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

Контрольный период делится на R контрольных временных отрезков (фиг. 1) в каждом из которых K≥2 раз измеряют уровни помех, вычисляют их усредненное значение и среднеквадратическое отклонение (СКО), прогнозируют ожидаемый уровень помех и сумму уровня помех и СКО на анализируемых частотах и оценивают величину ошибки прогноза этих показателей.

Для этого r×K раз измеряют уровни помех в контрольных временных интервалах (фиг. 2), где - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT_r и r=1, 2, …, R и вычисляют усредненное значение уровня помех и СКО σ_r для каждой длительности контрольного временного интервала ΔT_r

Затем вычисляют предварительные значения уровней прогнозируемых помех и предварительные значения прогнозируемой суммы усредненного значения и СКО методом адаптивного экспоненциального сглаживания значения измеренного в момент времени Т_И и усредненного значения соответственно по формулам:

Причем экспоненциальное сглаживание выполняют последовательно для каждого контрольного временного интервала ΔT_r, принимая значение коэффициента экспоненциального сглаживания α_i с шагом Δα=0,1 в интервале Вычисленные предварительно спрогнозированные уровни помех при всех заданных α_i и для каждой длительности контрольного интервала ΔT_r запоминают, после чего измеряют дополнительно уровень помех X_i,доп и (X_i,доп+σ_к) в момент времени ΔT_r+Δt и вычисляют среднеквадратичные ошибки δ²_i,r(α_i) между и X_i,доп, a также δ²_i,r(α_i+σ_r) между и (X_i,доп+σ_r), по формулам:

Выделяют наименьшее значение среднеквадратичной ошибки из всех δ²_i,r(α_i) и δ²_i,r(α_i+σ_r), соответствующее ему значение и начальный момент времени контрольного интервала после чего в р≥2 временных интервалах, следующих за контрольным, где р=1, 2, …, S в зависимости от выбранного параметра прогнозирования X_i,доп или (X_i,доп+σ_r), коэффициенте сглаживания и начальном моменте времени контрольного временного интервала для каждой i-ой частоты рассчитывают прогнозируемый уровень помех на момент времени (р+1)Т_И+Δt методом адаптивного экспоненциального сглаживания дополнительно измеренного уровня помех X_i,доп или (X_i,доп+σ_r) соответственно на момент времени рТ_И+Δt и спрогнозированного на момент времени рТ_И+Δt по формулам:

причем для временного интервала р=1 принимают и

После наступления временного интервала p=S вычисление прогнозируемого уровня и суммы уровня и СКО для следующего интервала и измерение уровня X_i,доп прекращают и выполняют измерения и вычисления аналогичные контрольному временному интервалу ΔT_r и временным интервалам р=1, 2, …, S

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

Заявленный способ поясняется чертежами, на которых показано:

фиг. 1 - пример формирования контрольных временных интервалов;

фиг. 2 - временной ряд уровней помех на анализируемой частоте;

фиг. 3-значения предварительно спрогнозированных уровней помех в момент времени ΔT_r+Δt на одной частоте при некоторых различных анализируемых параметрах, коэффициентах адаптивного сглаживания и длительности контрольного интервала;

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

В заявленном способе краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне на каждой i-й частоте из множества анализируемых частот производят прогноз уровней помех на некоторое время вперед, для чего в контрольных временных интервалах где и Т_И - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT_r с шагом Δt=ΔT/K, как показано на фиг. 2, r×K раз измеряют уровни помех Затем вычисляют усредненное значение уровня помех и СКО σ_r для каждой длительности контрольного временного интервала ΔТ_г по формулам:

Затем методом адаптивного экспоненциального сглаживания и измеренных в момент времени Т_И и усредненного значения и соответственно вычисляют предварительные значения уровней прогнозируемых помех и предварительные значения прогнозируемой суммы усредненного значения и СКО на момент времени ΔT_r+Δt по формулам (1), причем экспоненциальное сглаживание выполняют последовательно для каждого контрольного временного интервала ΔT_r, принимая значение коэффициента экспоненциального сглаживания α_i с шагом Δα=0,1 в интервале Вычисленные предварительно спрогнозированные уровни помех и при всех заданных а, и для каждой длительности контрольного интервала ΔT_r запоминают.

После предварительного прогнозирования вычисляют среднеквадратические ошибки δ²_i,r(α_i) и δ²_i,r(α_i+σ_r) между и дополнительно измеренными значениями в момент времени ΔT_r+Δt уровня помехи X_i,доп и суммой (X_i,доп+σ_r) соответственно формулам (2). Выделяют наименьшее значение среднеквадратичной ошибки из совокупности полученных δ²_i,r(α_i) и δ²_i,r(α_i+σ_r), соответствующее ему значение коэффициента сглаживания и начальный момент времени контрольного интервала после чего в р≥2 временных интервалах, следующих за контрольным, где р=1, 2, …, S в зависимости от выбранного параметра прогнозирования X_i,доп ^или (X_i,доп+σ_r), коэффициенте сглаживания и начальном моменте времени контрольного временного интервала для каждой i-ой частоты рассчитывают прогнозируемый уровень помех на момент времени (р+1)Т_И+Δt методом адаптивного экспоненциального сглаживания дополнительно измеренного уровня помех X_i,доп или (X_i,доп+σ_r) соответственно на момент времени рТ_И+Δt и спрогнозированного на момент времени рТ_И+Δt по формулам (3), причем для временного интервала р=1 принимают После наступления временного интервала p=S вычисление прогнозируемого уровня и суммы уровня и СКО для следующего интервала и измерение уровня X_i,доп прекращают и выполняют измерения и вычисления аналогичные контрольному временному интервалу ΔT_r и временным интервалам р=1, 2, …, S

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

Выявления закономерностей динамического изменения уровней помех можно производить с помощью его выравнивания и сглаживания. Экспоненциальное сглаживание - один из простейших и распространенных приемов выравнивания ряда. В его основе лежит расчет экспоненциальных средних временного ряда уровней помех [2]. Однако, в каждом конкретном случае необходимо производить адаптивный подбор коэффициента адаптивного сглаживания α, при краткосрочном прогнозирований необходимо как можно скорее отобразить изменения α и в тоже время как можно лучше «чистить» ряд от случайных колебаний. С одной стороны, необходимо увеличить вес более свежих наблюдений, что может быть достигнуто повышением α, с другой стороны, для сглаживания случайных отклонений величину α нужно уменьшить. Например, результаты вычислений уровней помех, а также суммы уровня помех и СКО при α=0,3; 0,7 на одной частоте приведены на фиг. 3. Таким образом, правильность вычисления коэффициента адаптивного сглаживания а определяет достоверность прогноза.

Согласно [1] следует, что если строить интегральные функции распределения вероятностей уровней помех на каждой частоте за время предыстории, то для работы следует выбирать ту частоту, для которой интегральная функция распределения вероятностей уровней помех сильнее смещена влево и имеет большую крутизну (т.е. соответственно и σ_r малы). Поэтому в качестве прогнозируемого параметра может выступать такой показатель как сумма усредненного значения уровней помех и среднеквадратического отклонения на периоде предыстории.

Процессы, протекающие в ионосфере, являются нестационарными, но на коротких интервалах времени их можно рассматривать как квазистационарные [1]. Свойства ионосферы в обычных условиях изменяются относительно медленно. Это в основном предопределяет ограниченный во времени статистически устойчивый квазистационарный характер усредненных последовательностей уровней помех. Для учета длительности квазистационарности помеховой обстановки в заявленном способе есть возможность изменения длительности контрольного временного интервала за время которого будет набираться статистика уровней помех на каждой анализируемой частоте, изменение шага Δt=ΔТ/K определяет частоту выборок временного ряда, а изменение последующих за контрольным временных интервалов р, где р=1, 2, …, S - определяет длительность прогноза.

Для прогнозирования большую роль играет длительность ряда динамики изучаемого явления. С одной стороны кажется очевидным, что чем длиннее период предыстории, тем надежнее прогноз [3], с другой стороны применение длинных рядов динамики может привести к излишнему и недействительному увеличению роли прошлого и влиянию на будущее в силу реальной нестационарности процессов. Для примера в таблицу (фиг. 4) сведены значения среднеквадратических ошибок прогноза значений уровня помех и суммы уровня помех и СКО, вычисленных для одной частоты при различных длительностях контрольных интервалов и с различным коэффициентом сглаживания. Результаты вычислений показывают, что минимальную ошибку прогноза обеспечивает прогнозирование параметра на длительности предыстории равной трем временным отрезкам и с коэффициентом α=0,7. Следовательно, продолжительность учета текущей предыстории и метод ее статистической обработки и параметр прогнозирования необходимо корректировать, исходя из предварительных статистических исследований.

Таким образом, статистические исследования показали, что для повышения вероятности правильного краткосрочного прогнозирования (на 10-20 мин) необходимо по результатам измерений на каждой частоте непрерывно, циклично вычислять такие показатели, как средний уровень помехи и сумму среднего уровня и СКО, и производить адаптивный подбор коэффициента адаптивного сглаживания α для различных длительностей контрольных интервалов ΔT. Отмеченное подтверждает возможность достижения сформулированного технического результата.

Литература

1. Комарович В.Ф., Сосунов В.Н. «Случайные радиопомехи и надежность КВ связи» - М.: Связь, 1977. - 136 с: ил.

2. Лукашин Ю.П. «Адаптивные методы краткосрочного прогнозирования временных рядов» - М.: Финансы и статистика, 2003 г. - 416 с: ил.

3. Ковалева Л.Н. «Многофакторное прогнозирование на основе рядов динамики» - М.: Статистика, 1980. - 102 с.

1. Способ краткосрочного адаптивного прогнозирования уровней помех в декаметровом диапазоне, заключающийся в том, что на каждой из N≥2 анализируемых частот K≥2 раз в контрольном временном интервале где и Т_И - соответственно начальный и конечный моменты времени контрольного временного интервала ΔT с шагом Δt=ΔT/K, измеряют уровни помех X_i, где i=1, 2, …, N, вычисляют их усредненное значение в пределах контрольного временного интервала по результатам всех измерений где k=1, 2, …, K, методом экспоненциального сглаживания последовательно вычисляют предварительные значения уровней помех для коэффициентов экспоненциального сглаживания с шагом Δα_i=0,1, измеряют дополнительно уровень помех X_i,доп в момент времени Т_И+Δt, вычисляют среднеквадратические ошибки δ²_i (α_i) между и X_i,доп, выделяют наименьшее значение среднеквадратической ошибки δ²_i (α_i) и соответствующее ему значение после чего при выбранном параметре α_i для каждой i-й частоты рассчитывают прогнозируемый уровень помех, отличающийся тем, что в ходе вычислений изменяют начальный момент времени контрольного временного интервала и в пределах каждого контрольного временного интервала ΔT_r, где r=1, 2, …, R и R≥2, вычисляют усредненное значение уровня помех по результатам всех измерений и значения суммы усредненного значения и его среднеквадратического отклонения (СКО) для соответствующей длительности контрольного временного интервала ΔT_r, затем вычисляют предварительные значения уровней прогнозируемых помех и предварительные значения прогнозируемой суммы усредненного значения и СКО методом адаптивного экспоненциального сглаживания значения и измеренного в момент времени Т_И и усредненного значения и соответственно, причем экспоненциальное сглаживание выполняют последовательно для каждого контрольного временного интервала ΔT_r, принимая значение коэффициента экспоненциального сглаживания α_i с шагом Δα=0,1 в интервале вычисленные предварительно спрогнозированные уровни помех и при всех заданных α_i и для каждой длительности контрольного интервала ΔT_r запоминают, после чего измеряют дополнительно уровень помех X_i,доп и (X_i,доп+σ_r) в момент времени ΔТ+Δt и вычисляют среднеквадратичные ошибки δ²_i,r(α_i) между и X_i,доп, а также δ²_i,r(α_i+σ_r) между и (X_i,доп+σ_r), выделяют наименьшее значение среднеквадратичной ошибки из всех δ²_i,r(α_i) и δ²_i,r(α_i+σ_r), соответствующее ему значение и начальный момент времени контрольного интервала после чего в р≥2 временных интервалах, следующих за контрольным, где р=1, 2, …, S в зависимости от выбранного параметра X_i,доп или (X_i,доп+σ_r), коэффициенте сглаживания и начальном моменте времени контрольного временного интервала для каждой i-й частоты рассчитывают прогнозируемый уровень помех или на момент времени (р+1)Т_И+Δt методом адаптивного экспоненциального сглаживания дополнительно измеренного уровня помех X_i,доп или (X_i,доп+σ_r) соответственно на момент времени рТ_И+Δt и спрогнозированного на момент времени рТ_И+Δt, причем для временного интервала р=1 принимают и а после наступления временного интервала р=S вычисление прогнозируемого уровня и суммы уровня и СКО для следующего интервала и измерение уровня X_i,доп прекращают и выполняют измерения и вычисления, аналогичные контрольному временному интервалу ΔT_r и временным интервалам р=1, 2, …, S.

2. Способ по п. 1, отличающийся тем, что среднеквадратическое отклонение (СКО) уровней помех находят по формуле:

3. Способ по п. 1, отличающийся тем, что в контрольном временном интервале ΔT_r для прогнозирования на i-й частоте суммы уровней предварительно спрогнозированных помех и СКО в момент времени Т_И+Δt с помощью экспоненциального сглаживания значения измеренного в момент времени Т_И, и усредненного значения рассчитывают значение предварительно прогнозируемого уровня помех в момент времени Т_И+Δt, последовательно для различных значений коэффициента экспоненциального сглаживания, задаваемого в интервале 0<α<1 с шагом 0.1 по формуле:

а среднеквадратичную ошибку прогноза помех на i-й частоте вычисляют по формуле:

4. Способ по п. 1, отличающийся тем, что для интервалов р=1, 2, …, S для прогнозирования на i-й частоте суммы уровня ожидаемых помех и СКО в момент времени (р+1)Т_И+Δt с помощью экспоненциального сглаживания суммы (для временного интервала р-1 и (X_i,доп+σ_r), в момент времени рТ_И+Δt рассчитывают значение прогнозируемого уровня помех в момент времени (р+1)Т_И+Δt по формуле: