СПОСОБ ПЕРЕДАЧИ ТЕЛЕМЕТРИЧЕСКОЙ ИНФОРМАЦИИ
METHOD FOR TRANSMITTING TELEMETRIC INFORMATION

Изобретение относится к системам передачи телеметрической информации (ТМИ), искаженной при передаче помехами естественного и искусственного происхождения. Оно ориентировано на внедрение в существующую практику передачи ТМИ новых алгоритмов структурно-алгоритмических преобразований (САП) данных телеизмерений, предназначенных для обеспечения дополнительного малоизбыточного и безызбыточного помехоустойчивого кодирования передаваемых сообщений.
При этом в качестве дополнительных условий, упрощающих реализацию предлагаемого способа в существующих бортовых телеметрических системах (БРТС), выступает требование обеспечения адаптации к структурам представления данных, которые уже используются. В этом случае модернизация существующих БРТС может быть обеспечена путём перепрограммирования программируемых логических интегральных схем (ПЛИС), микроконтроллеров и микропроцессоров, составляющих основу их построения.
Известна группа изобретений, в которых для достижения высоких показателей помехозащищенности данных измерений используют дополнительное экономное помехоустойчивое кодирование. К ней относятся: [1] «Способ передачи информации и система для его осуществления», патент RU №2586605, приоритет от 22.03.2013г. и [2] «Способ передачи информации», патент RU №2609747, приоритет от 13.08.2017г.
b
b
X
N
n
m
=
m
=
n
b
n
b
n
1
j
3
j
j
1
3
3
j
2
3
j
2
n
n
(0)
(max)
В патенте [1], основу которого составляет нетрадиционное представление данных телеизмерений для представления результатов дополнительного малоизбыточного помехоустойчивого кодирования данных их образами-остатками  и , получающимися в результате арифметической операции деления их исходных значений , представленных  = 2 - разрядным двоичным позиционным кодом, на оптимально выбранные модули сравнения  2- 1 и  2+ 1. Его использование требует, в общем случае, один дополнительный разряд (+ 1)-й разряд для разрешения неоднозначности между значениями первого («нулевого» значения остатков = <00…0>, представленного последовательностью -разрядных символов «0» двоичного кода, и «последнего», самого большого по величине значения = <100…0>, отличающегося наличием самого старшего  (+ 1)-й разряда с двоичным символом «1».
X
b
,b
j
1
j
3
j
Суть математического задания на осуществление перепрограммирования существующих БРТС заключается в следующем. Исходное представление значения сообщения , выраженное с использованием образов-остатков , определяет следующая система уравнений:
X
= m
l
+ b
j
1
1
1
j
X
= m
l
+ b
,
j
3
3
3
j
(1)
m
=
-
m
=
+
-
1
3
n
n
где21и 21 оптимальные модули сравнений;
l
, l
.
X
m
m
1
3
j
1
3
- неполные частные, получаемые от деления на делители (модули сравнений  и ).
При переходе к системе сравнений, которая является эквивалентной представлению (1):
X
≡ b
(mod m
)
j
1
j
1
X
≡ b
(mod m
)
j
3
j
3
(2)
X
.
l
l
.
j
1
3
получаем сжатую математическую форму представления значения сообщения Сжатой она является потому, что в системе сравнений (2) отсутствуют значения неполных частных и
Из представления результата дополнительного кодирования в следующих двух вариантах:
C
b
>
,
b
>
>
,
j
1
j
2
3
j
2
2
(1)
= < < <(3)
C
b
>
,
b
>
>
,
j
3
j
2
1
j
2
2
(1)*
= < < <(4)
l
l
m
m
-
1
3
1
3
также следует, что помимо значений неполных частных ив канал связи не передаются и значения модулей сравнения иони являются частью ключевых данных, которыми стороны обмениваются между собой до начала работы.
X
l
l
.
m
m
b
, b
китайская
теорема об остатках
конструктивная теорема об остатках
b
, b
X
j
1
3
1
3
1
j
3
j
1
j
3
j
j
При сравнении представлений значения с использованием системы уравнений (1) и системы сравнений (2), следует, что в результатах дополнительного кодирования присутствует не явно выраженный эффект синтаксического сжатия данных. Неполные частныеи при соответствующем выборе модулей сравнений и, могут быть восстановлены на приемной стороне на основе принятых образов-остатков с использованием алгоритмов теорем об остатках, известных как  (КиТО) [3] или  (КтТО) [3], используемая в патенте [1]). Поэтому они не передаются, хотя и представляют собой, помимо образов-остатков , часть информации о значении передаваемого сообщения (1). В этом заключаются основополагающие теоретические положения сжатия данных и повышения их помехозащищённости при нетрадиционном представлении данных образами-остатками.
X
b
,b
j
1
j
3
j
В патенте [2], используется другой алгоритм дополнительного кодирования, который является эквивалентным основному теоретическому представлению значения сообщений образами-остатками в плане повышения показателей помехоустойчивости передачи информации, но отличается простотой реализации.
n
X
а
а
j
ст
j
мл
j
Его суть в упрощённом для толкования варианте заключается в том, что в исходном кодовом 2 - разрядном двоичном слове  выделяют старшее () и младшее () полуслова:
X
а
>
,
а
>
>
,
j
ст
j
2
мл
j
2
2
= < < <(5)
C
j
(1(Э))
которые затем переставляют местами с образованием дополнительно закодированного сообщения :
C
а
>
,
а
>
>
,
j
мл
j
2
ст
j
2
2
(1(Э))
= < < <(6)
m
=
.
N
N =
а
а
а
а
.
2
мл
j
ст
j
мл
j
ст
j
min
внеш
n
При таком представлении появляется необходимость использования модуля сравнения 2 Если число () разрядов, требуемых для представления значений контролируемого параметра, является нечётным, равным, например, 9, то наибольшее предпочтение имеют следующее представление: пятью, а  - четырьмя разрядами. Возможны и другие варианты разделения исходного кодового слова W на кодовые сегменты и  Выбор наиболее подходящего вариантаопределяется необходимостью соблюдения требований, которые являются заданными для данной ТМС и соответствуют формированию необходимого значения d. В качестве требований могут быть использованы возможные для ТМС, подвергаемой модернизации, варианты формирования внешних структур представления данных (S).
m
m
=
.
=
N
X
.
2
2
min
j
n
n
N
-1
N
-1
Также алгоритмы представления (5) и дополнительного кодирования (6) используют на приемной стороне, после декодирования, заключающегося в первичной обработке полученной информации с обнаружением и исправлением ошибок передачи ([5], Способ первичной обработки информации с обнаружением и исправлением ошибок передачи, патент RU № 2658795, приоритет от 30.05.2017г.). В этом случае используется операция сравнения получаемых при обработке данных по модулям сравнения (), определяемым из ряда чисел 2:2В результате этого и минимальное кодовое расстояние также будет определяться из ряда чисел 2:d2, где  – число разрядов в кодовом слове-измерении
n
X
X
,
b
b
b
n
b
b
b
b
b
b
j
j
3
j
3
j
3
j
2
3
j
2
3
j
2
3
j
2
2
3
j
2
3
j
2
1
j
n
n
(0)
(max)
(0)
(max)
(0)
(max)
n
0
0
1
1
0
1
Первые два способа [1] и [2], обладают различными недостатками. Первый [1], выбранный в качестве наиболее подходящего при передаче БМП, требует, как это было отмечено ранее, введения дополнительного (+ 1) -го символа. Этот недостаток становится несущественным при использовании в кодовом слове-измерении дополнительного символа «Контроль чётности бит (КЧБ). В предлагаемом техническом решении его также используют по своему основному предназначению: для контроля наличия ошибок в принятом на фоне помех сообщении. Однако им уже является не исходное кодовое слово-измерение а его остаток по модулю (2+1):(mod 2+1). При этом также его используют и по новому назначению: для разрешения неоднозначности «нулевого» значения остатков = <00…0>, представленного последовательностью -разрядных символов «0» двоичного кода, и «последнего», самого большого по величине значения  = <100…0>. В первом случае КЧБ представляют символом «»:  = <00…0>, а во втором - символом «»: = <100…0> ↔ <00…0>, где он выделен жирным цветом. При этом только во втором случае подстановка символа «1» противоречит логике формирования КЧБ и, следовательно, основное предназначение символа КЧБ не выполняется. Однако при этом кодовые конструкции  = <00…0> и  = <00…0> становятся контрастно выраженными и, следовательно, хорошо различимыми. При представлении второго образа остатка (mod 2-1) таких дополнительных сложностей нет.
X
а
а
j
ст
j
мл
j
Второй способ [2], основу которого составляет разделение исходного кодового слова-измерений  на старшее () и младшее () полуслова (5) и перестановку их местами (6), проблемы, связанной с нахождением дополнительного резерва в виде дополнительного символа, нет. Он относится к разряду безызбыточных способов помехоустойчивого кодирования.
внутр
внеш
внеш
N
Использование способов обеспечивает возможность формирования внутренней структуры представления данных (S) телеизмерений, которая по разрядности () двоичных слов, используемых для передачи информации, совпадает с существующей структурой представления слов-измерений и пакетов данных. Последняя представляет собой внешнюю структуру представления данных (S), которая используется в существующих БРТС в качестве телеметрических кадров, разрядности представления данных телеизмерений и сигналов синхронизации. В самом общем случае, внешняя структура представления данных (S) может включать в себя и другие данные, которые относят к числу служебных данных. К их числу относят: адресную информацию, проверочные символы избыточных помехоустойчивых кодов, защитные символы между словами, используемые для повышении точности системы синхронизации ([4], «Современная телеметрия в теории и на практике / Учебный курс», Спб.: Наука и Техника, 2007. - 672с, стр. 465).
Сущностные характеристики предлагаемого способа заключаются в том, чтобы внедрение алгоритмов САП [1,2] не требовало бы выхода за пределы ограничений, определяемые существующими правилами представления данных в БРТС, которые используются в настоящее время. В то же время их использование обеспечивало бы достижение наибольшего технического эффекта. Это позволило бы использовать существующие приёмно-регистрирующие станции (ПРС) и наземный телеметрический комплекс, в целом, без дорогостоящих аппаратурных доработок, а все вносимые новшества и изменения реализовать за счёт программных средств и разработки обновлённого ПМО.
Способы [1,2], обладая способностью обнаружения и исправления ошибок передачи за счёт появления нового инварианта в виде «группового свойства равноостаточности», имеют различные преимущества и недостатки. Ими определяются перспективные направления их реализации при модернизации существующих БРТС программными методами.
Они в наглядной форме представлены на иллюстрациях ММП и БМП, приведенных на фиг. 1, фиг. 2, фиг. 3, фиг. 4 и фиг. 5.
X
C
b
>
,
b
>
>
.
t
C
b
>
,
b
>
>
C
≡ζ
(mod
).
j
j
1
j
2
3
j
2
2
min
ζ
k
j
1
j
2
3
j
2
2
min
j
j
min
(1)
(1)
(1)
На фиг. 1 и 2 представлены в виде графиков основополагающие принципы предлагаемого помехоустойчивого кодирования, которые являются общими для предлагаемого изобретения и способов [1,2]. На фиг. 1(А) в качестве примера выбран телеметрируемый параметр (ТМП), изменение значений  которого происходит по закону косинуса. На следующей иллюстрации (фиг. 1(Б)), приведены результаты дополнительного кодирования = < < < В результате дополнительного кодирования увеличивается минимальное кодовое расстояние d, в результате чего амплитуду преобразованного ТМП увеличивают в соответствующее число раз. Но так как шкала (Ш) представления его значений ограничена разрядной сеткой представления данных, определяемой количеством разрядов N двоичного кода, то график (фиг. 1(Б)) претерпевает разрывы, длительность которых(Δ) определяют значения «равноостаточности», представленные в графическом виде на фиг. 1(В). В соответствии с групповым свойством «равноостаточности» результаты дополнительного кодирования = < < <, принятые без ошибок при делении на минимальное кодовое расстояние d, дают одно и то же значение остатка: d
t
,
C
ζ
ζ
ζ
C
C
C
C
C
l
l =
ζk
j
j
j
s
s
j
s
j
s
s
j
min
На фиг. 2 (А,Б,В) повторены соответствующие иллюстрации, приведенные на фиг. 1, но при имитации искажений значения ТМП помехами. При этом наиболее часто встречающиеся на выделенном временном интервале которых (Δ) значения остатков ζпредставляющие собой моду построенной гистограммы распределения частот встречаемости Мо(ζ), будут отождествляться с закодированными значениями , принятыми без ошибок. Другие значения остатков :≠ будут соответствовать значениям, которые были искажены при передаче по каналам в условиях помех. Они подлежат исправлению в условиях избыточности данных телеизмерений на ближайшее значение , отличающееся от  минимум кодового расстояния: | - | =  d, где1,2,3 - наименьшее число. В этом заключается смысл предлагаемого режима «мягкого» декодирования (восстановления) результатов дополнительного кодирования ТМП при приёме.
X
b
mod
b
mod
N
n
j
j
1
j
3
n
n
На фиг. 3(А) представлен случай дополнительного кодирования ТМП, относящегося к классу медленноменяющихся (ММП), а на фиг. 3(Б) и 3(В) приведены результаты графического представления значений слов-измерений образами- остатками ( 2- 1) (фиг. 3(Б)) и (2+ 1) (фиг. 3(В)) при байтовой структуре ( = 2) представления данных телеизмерений.
X
а
а
а
а
f
X
N
n
j
ст
j
мл
j
ст
j
мл
j
o
j
j
j+1
Однако для ММП наиболее подходящим становится способ [2], когда исходные кодовые слова-измерения  разделяют на старшее () и младшее () полуслова (5) с последующей перестановкой их местами (6). В этом случае значения старшего () полуслова будут меняться существенно реже по сравнению с когерентными данными младшего () полуслова, требующего такой же частоты опроса (), определяемой в соответствии с теоремой В.А.Котельникова, что и исходное слово-измерений . Но число разрядов кода, которое при этом необходимо передать будет в 2 раза меньше. Тогда необходимо для формирования стандартных слов W и W прежней разрядности  = 2, необходимо выполнить следующие операции:
а
а
f
мл
j
мл(
j
+1)
j
j+1
o
(ММП)
1) формально объединить в единую кодовую конструкцию младшие (и) полуслова двух ММП (W и W), требующих для передачи одну и ту же частоту  опроса ();
а
а
k
f
/k
ст
j
ст(
j
+1)
j
j+1
o
(ММП)
2) выполнить ту же операцию по отношению к старшим (и) полусловам тех же двух ММП (W и W), с последующей их передачей с пониженной в раз частотой опроса ().
X
f
а
а
а
j
o
ст
j
ст
j
мл
j
j
j+1
(ММП)
При таком представлении обеспечивается возможность передачи большего объёма ММП по сравнению с существующими технологиями, которые использованы в существующих БРТС. Появление такого эффекта связано с тем, что ранее представление слов-измерений  рассматривалось в качестве наименьшего информационного элемента представления данных, поэтому назначаемая частота опроса () была избыточной по отношению к входящему в его состав старшему () полуслову. Его разделение на старшее () и младшее () полуслова будет способствовать более рациональному выбору способов формирования кодовых слов ММП (W и W).
C
b
>
,
b
>
>
b
>
,
b
>
а
а
j
1
j
2
3
j
2
2
1
j
2
3
j
2
ст
j
мл
j
Однако этот приём может не дать ожидаемого выигрыша при передаче быстроменяющихся параметров (БМП). В этом случае наибольшее предпочтение приобретает способ [1], когда новые слова () будут сформированы из образов-остатков< < <, в соответствии с алгоритмами кодирования (3) и (4). В этом случае динамика изменения значений образов-остатков < <уже не будет такой контрастной, как при передаче ММП. И тогда преимущества предложения о разных частотах опроса старшего () и младшего () полуслов, в общем случае, могут быть утрачены.
а
ст
j
Эта ситуация продемонстрирована на иллюстрации (фиг. 4). Динамика изменения значений телеизмерений, графическое представление которых приведено на фиг. 4(А), существенно выше, следовательно, старшее полуслово () при использовании способа [2] необходимо передавать чаще по сравнению с ТМП, иллюстрация которого приведена на фиг. 3. Таким образом, выигрыш в возможности увеличения объемов передаваемой ТМИ будет сокращаться по мере увеличения динамики изменения значений телеизмерений. В этом заключается проблема создания адаптивных БРТС, а в перспективе и интеллектуальных их образцов. Для повышения эффективности их использования необходима предварительная работа по классификации с точки зрения выбора наиболее подходящих алгоритмов дополнительного помехоустойчивого кодирования (3, 4) или (6).
b
mod
>
b
mod
>
.
C
b
>
,
b
>
>
N
n
n
f
X
а
а
f
а
k
f
/k
1
j
2
3
j
2
j
1
j
2
3
j
2
2
o
j
ст
j
мл
j
o
ст
j
o
(1)
(БМП)
(ММП)
(ММП)
Так, при передаче БМП наибольший технический эффект можно получить при использовании способа [1]. На фиг. 4(Б) и фиг. 4(В) представлены в графическом виде изменения значений образов-остатков <(31)и<(33), соответственноНа фиг. 4 (Г)изображено представление закодированных значений = < < <, воспринимаемых как единое новое кодовое слово с разрядностью  = 2, где  = 5. Из приведенных иллюстраций следует, что для каждого из них требуется исходная частота опроса БМП (). В этом заключается принципиальное отличие предлагаемого изобретения от способа передачи ММП при использовании разделения исходного кодового слова-измерения на старшее () и младшее () полуслова, когда исходная частота его опроса () по отношению к старшему полуслова () может быть уменьшена в  раз: ().
Поэтому и сущностные характеристики представления БМП проявятся в другом: избыточность передаваемых данных не сокращается, как это было в случае с передачей ММП. Наоборот, искусственно формируют дополнительную избыточность передаваемой информации о БМП при сохранении прежней транспортной скорости передачи ТМИ. Таким образом, речь идёт об избыточности особого типа, когда она присутствует, но не находит своего выражения в увеличении общего объёма передаваемой ТМИ.
f
m
m
o
1
3
(БМП)
n
n
Сущностные характеристики предлагаемого способа заключаются в том, что значения обоих полуслов, представленных образами-остатками, будут передаваться с исходной  частотой опроса БМП (), в результате чего передаваемая информация будет полноценно продублирована без привлечения дополнительных датчиков. При этом исключается необходимость обратного восстановления значений БМП в первоначальном виде. Принимаемые данные о значениях БМП уже будут представлены в системе остаточных классов (СОК). В результате показатели эффективности их обработки в соответствии с теорией обработки данных, представленных в СОК, будут повышены. Это положение в данном изобретении проиллюстрировано на примере использования алгоритмов адаптивной нелинейной фильтрации (АНФ) (фиг. 5). На фиг. 5 представлены для сравнения значения дисперсий зашумляющей помехи при использовании алгоритмов АНФ для следующих случаев: 1) традиционного представления данных ТМП позиционным двоичным кодом и 2) когда они были на передающей стороне представлены образами-остатками (3) и при приёме алгоритмы АНФ были применены по отношению к значениям образов-остатков (фиг. 4(Б) и фиг. 4(В)), полученных с использованием модулей сравнения  = 2 -1 и  = 2 +1, соответственно.
Во втором случае значения вычисленных дисперсий зашумляющей помехи становятся в 2 - 3 раза меньше их исходного значения при существующем способе передачи БМП, что показано на иллюстрации, приведенной на фиг. 5.
Кроме того, дополнительное экономное малоизбыточное помехоустойчивое кодирование [1] рассматривают и как структурно-алгоритмическое преобразование первого этапа (САП-1), имеющее два взаимообусловленных вида: прямое структурно-алгоритмическое преобразование (ПСАП-1), определяющее операцию кодирования информации и обратное структурно-алгоритмическое преобразование (ОСАП-1), представляющее собой операцию декодирования.
T
X
F
X
j
max
j
В отличие от подобных САП-1 математических аналогов, известных, например, как прямое и обратное преобразования Фурье (ППФ и ОПФ) [4], обратное структурно-алгоритмическое преобразование (ОСАП-1) имеет два вида: универсальное ОСАП-1 (УОСАП-1), отождествляемое с понятием «жесткого» декодирования, и частное ОСАП-1 (ЧОСАП-1), которое определяют, как «мягкое» декодирование принятой и обрабатываемой информации. При этом алгоритм «жесткого» декодирования применим всегда, независимо от свойств передаваемой информации. Однако при этом в соответствии с законами Природы «за универсальность приходится расплачиваться потерями эффективности», что проявляется в практическом отсутствии возможности обнаружения и исправления ошибок передачи данных. Алгоритм ЧОСАП-1 («мягкого» декодирования) позволяет использовать естественную избыточность передаваемых цифровых данных для обнаружения и исправления ошибок передачи информации при ее приеме и обработке. Естественная избыточность цифровых данных является следствием применения теоремы В.А.Котельникова о дискретизации, в соответствии с которой интервалы (Δ) между опросами аналогового параметра или сигнала  определяют, как обратное отношение к значению удвоенного значения спектральной составляющей наибольшей частоты (2) [4]:
ΔT
= .                                                            (7)
F
T
T
F
X
, X
X
N
max
i
j
-1
j
j
+1
Поскольку вероятность появления частотной составляющей  на достаточно малом интервале времени Δτ ≥ 3Δнезначительна по величине, то выбранное значение интервалов Δ для других спектральных составляющих  спектра параметра или сигнала оказывается малым, что проявляется в корреляционной взаимосвязи соседних значений сообщений или слов-измерений (и). В телеметрии эта составляющая внутренней избыточности составляет, в среднем, 90%. При передаче других видов информации, например, речевой, акустической, потокового видео или навигационной она также значительна. Предполагается одновременное использование двух режимов декодирования принимаемых цифровых сигналов УОСАП-1 и ЧОСАП-1. В этом случае детектор ЧОСАП-1 работает под управлением УОСАП-1. Детектор УОСАП-1 на основе принятой последовательности закодированных данных, представленных  - разрядным двоичным кодом, определяет наличие их корреляционной взаимосвязи, и, если устанавливает её наличие для трёх и более следующих подряд сообщений, то подключает для их обработки, осуществляемой с целью обнаружения и исправления ошибок передачи, детектор ЧОСАП-1. Введение цифр, обозначающих этапы (i) распределённых структурно-алгоритмических преобразований (САП-i), также становится обязательным для определения того информационного сечения (i) в существующих трактах формирования и передачи данных, в котором предусмотрено дополнительное кодирование передаваемой информации.
Основной недостаток всех известных классических методов избыточного помехоустойчивого кодирования и обеспечения защиты информации заключается в том, что специфические особенности передаваемой информации не учитываются. Поэтому и их эффективность в условиях чрезвычайно быстрого роста объемов передаваемой информации и скоростей передачи не оправдывает надежд заказчиков и разработчиков систем передачи данных (СПД).
Способ [2] отличается от других аналогов, например, [1] возможностью наиболее простой технической реализации безызбыточного помехоустойчивого кодирования, благодаря чему он оказывается наиболее предпочтительным при передаче высокоскоростной информации вследствие высоких требований, относящихся к показателю оперативности осуществления структурно-алгоритмических преобразований (САП-i), прежде всего, ПСАП-i, поскольку обратные САП-i (ОСАП-i) могут быть реализованы не в реальном, а в отложенном времени, например, при организации обработки данных, в том числе, и при реализации способов первичной обработки данных ([5], патент RU № 2658795 «Способ первичной обработки информации с обнаружением и исправлением ошибок передачи», приоритет от 30.05.2017г. и [6], патент RU № 2672392 «Способ первичной обработки информации с использованием адаптивной нелинейной фильтрации данных измерений», приоритет от 27.06.2017г.).
Основу способа [5], обеспечивающего при приёме возможность обнаружения и исправления ошибок передачи, составляет использование двух видов декодирования, которые по аналогии с классической теорией помехоустойчивого кодирования называют «жёстким» и «мягким». Но принципиальное их отличие от известных информационных технологий заключается в том, что декодирования при этом подвергаются не символы двоичного кода, а выделенные из цифрового группового сигнала, основу которого составляют восстановленные при приёме его символы «1» и «0», значения кодовых конструкций или слов-измерений, каждый из которых представлен определённым их количеством (N).
Алгоритм «жёсткого» декодирования данных и сообщений использование алгоритма конструктивной теоремы об остатках (КтТО). При модулях, отличающихся друг от друга по абсолютной величине на 2 , алгоритм восстановления с использованием КтТО состоит из трёх звеньев вычислений [3,9]:
(8)
j
X
n
n
m
n -
m
n
j
3
3
В формуле (8) для упрощения записи опущен индекс «», свидетельствующий о том, что восстановленное значение представляет собой , обозначения  - означает, что Δ не делится без остатка на число , а Δ/ и ( +Δ)/читается, как: Δ и (+Δ) делится без остатка на число .
C
C
- C
|
С
C
- C
|
C
j -
X
C
C
C
, C
s
N
X
j
M
j
(
j
+1
(
j
+
s
)М
(
j
+
s
)
(
j
+
s
)+1
j
j
j
j
+1
(
j
+
s
)
(
j
+
s
)+1
j
N
N
N
В результате операции «мягкого» декодирования определяют графические фрагменты закодированных первоначальных значений БМП, заключенные между соседними максимальными значениями абсолютных разностей Δ = | ≥ 0,8 × 2, и Δ= | ≥ 0,8 × 2), равными или большими значений 0,8 × 2, где  результат дополнительного безызбыточного помехоустойчивого кодированиятого по принятому порядку счета слова-измерения или сообщения, где, определяют начало графического фрагмента, а - его окончание и начало следующего по порядку графического фрагмента, включающего в себя  выборок закодированных значений телеметрируемых параметров, а  - число разрядов, которые используют для представления слов-измерений или сообщений , обеспечивают обнаружение и исправление ошибок передачи значений телеметрируемых параметров на основе «групповых свойств равноостаточности», которые должны быть постоянными при отсутствии ошибок передачи информации в выделенных графических фрагментах телеметрируемого параметра, преобразованного при дополнительном помехоустойчивом кодировании на передающей стороне, а в результате операции «жесткого» декодирования осуществляют восстановление первоначальных результатов телеизмерений без исправления ошибок, разрядность которых известна на приемной стороне, осуществляют сглаживание или фильтрацию восстановленных в результате «жесткого» декодирования данных и по отношению к вычисленным соседним значениям телеизмерений определяют их разности, которые при выполнении операций «мягкого» декодирования с учетом разрешенных позиций для значений принятого безызбыточного помехоустойчивого кода используют в качестве допусков при выборе наиболее подходящих, кондиционных значений восстанавливаемых данных телеизмерений, определяемых в соответствии с принятым критерием достоверности, производят повторное «жесткое» декодирование данных телеизмерений, скорректированных в результате операций «мягкого» декодирования, сглаженные или отфильтрованные данные, полученные при первой операции «жесткого» декодирования сравнивают с синхронными, совпадающими по времени опроса, значениями, полученными в результате второй операции «жесткого» декодирования, результаты сравнения используют для оценивания достигаемого технического эффекта в виде оценок повышения показателей достоверности приема ой информации, а также для сравнения полученных результатов и корректировки значений сглаживания или фильтрации данных телеизмерений, полученных при первом «жестком» декодировании, определяют их близость к другим разрешенным позициям помехоустойчивого кода, сформированного в результате структурно-алгоритмических преобразований значений телеметрируемого параметра на передающей стороне, в результате чего реализуют расширенные возможности, используемые для контроля достоверности полученных результатов телеизмерений и информационной поддержки принятия решений.
Проведенный анализ позволяет перейти к описанию формулы изобретения.
N
n
b
m
b
m
X
t
jΔТ
Т
m
-
m
X
b
m
b
m
C
b
>
,
b
>
>
C
C
b
>
,
b
*>
>
ε
ε
b
>
= < b
+ ε
>
b
>
= < b
+ ε
>
ε
< ε
ε
< ε
1
j
1
3
j
3
j
j
1
3
j
1
j
1
3
j
3
j
1
j
2
3
j
2
2
2
внеш
внеш
j
j
1
j
2
3
j
2
2
1
j
3
j
1
j
2
1
j
1
j
2
3
j
2
3
j
3
j
2
1
j
1
j
3
j
3
j
n
n
*
*
*
*
*
(пв)
(пв)
Предлагаемый способ заключается в том, что на передающей стороне осуществляют сбор сигналов от датчиков измерений, включающих в себя и датчики быстроменяющихся параметров (БМП), преобразуют их в двоичный код, обеспечивают синхронизацию сформированных слов-измерений, представленных  = 2 - разрядным двоичным кодом, и формируют из них уплотненный цифровой групповой сигнал, подлежащий передаче по каналам связи, а на приемной стороне принимают полученную последовательность переданных символов двоичного кода. Он отличается тем, что на передающей стороне кодовые конструкции, сформированные на первом этапе кодирования слов-измерений или сообщений представляют образами-остатками (mod ) и (mod ), которые получают в результате операций, эквивалентных делению их дискретных значений , выбранных для передачи в моменты времени = , где Δ - интервал дискретизации по времени, определяемый в соответствии с теоремой о дискретизации  В.А.Котельникова, на установленные модули сравнения = 2 1 и = 2+ 1, результаты дополнительного кодирования исходных значений  образами-остатками (mod ) и (mod ) объединяют в новые кодовые слова  = < < <, где обозначение <  > означает представление двоичным кодом, таким образом, чтобы их разрядность и структура телеметрического кадра остались прежними, соответствующими обозначению S, для того, чтобы условия приёма и сбора телеметрической информации (ТМИ) остались неизменными, при приёме ТМИ восстанавливают переданные телеметрические кадры с установленной внешней структурой S, из которой выделяют дополнительно закодированные новые слова-измерения , применительно к данным БМП = < < <, искажённые при передаче помехами  и : <и <, после чего осуществляют их помехоустойчивое восстановление, обеспечивающее уменьшение интенсивности помехового воздействия: и и повышение, вследствие этого показателей достоверности приёма ТМИ.
C
b
>
,
b
*>
>
X
b
>
,
b
*>
,
b
>
b
*>
j
1
j
2
3
j
2
2
j
1
j
2
3
j
2
1
j
2
3
j
2,
*
*
/
*
*
Предлагаемый способ отличается также тем, что, помимо операций «жёсткого» декодирования, основу которого составляет алгоритм конструктивной теоремы об остатках, и «мягкого» декодирования переданных значений БМП на основе группового свойства равноостаточности, реализуют операции распараллеленной первичной обработки БМП с применением алгоритмов адаптивной нелинейной фильтрации, рассматривая представление значений БМП = < < <, в виде системы остаточных классов (СОК), вследствие чего обработке подвергают не исходные восстановленные значения БМП в виде  с уменьшенными ошибками передачи, а выделенные в соответствующие потоки исходных данных образы-остатки < <которые рассматривают в качестве распараллеленных исходных данных, далее обработку осуществляют в отдельности по отношению к каждому из сформированных потоков восстановленных данных: < и<а полученные при этом результаты, также являющие гомоморфными образами традиционной обработки с уменьшенными погрешностями, восстанавливают с использованием адаптивного алгоритма конструктивной теоремы об остатках, в результате чего дополнительно повышают показатели достоверности получаемой информации о БМП.
d
C
b
>
b
>
N
n
Х
=
Х
=
d
min
j
1
j
2
3
j
2
1
10
2
2
10
2
min
(тр)
Основу появления новых сущностных характеристик составляет эффект увеличения минимального кодового расстояния при использовании ПСАП-1 (дополнительного помехоустойчивого кодирования на основе представления дополнительно закодированных слов-измерений  образами-остатками < и <). Для пояснения этого эффекта необходимо взять соседние значения данных измерений при числе разрядности двоичных слов-измерений и сообщений, равной = 2 = 8, например, <115> = <01110011>и <116> = <01110100>. При традиционном способе их кодирования двоичным кодом кодовое расстояние между ними будет равно  = |115 - 116| =1. При этом в соответствии с теорией помехоустойчивого кодирования отсутствует возможность обнаружения ошибок передачи ТМИ ([7], Зюко А.Г. Помехоустойчивость и эффективность систем связи. - М.: Связь, 1972. - 360с.).
n
m
=
m
=
,
С
=
С
=
k
n
d
d
d
C
f
jΔT
C
C
– C
|
С
C
- C
|
s
C
d
C
d
модой -
1
2
1
2
10
2
2
10
min
min
min
j
j
M
j
(
j
+1
(
j
+
s
)М
(
j
+
s
)
(
j
+
s
)+1
j
min
j
j
j
min
j
j
j
j
4
4
n
(СОК)
N
N
1010
1011
При использовании предлагаемого способа и выборе модулей сравнения при = 4: 2-1 = 15 и 2+1 = 17получим: <1101>= <173> и = <1110><190>. В результате минимальное кодовое расстояние увеличилось в = 2+1 раз и стало при  = 4 равным  = |173 - 190| =17. В результате этого обеспечивается возможность обнаружения и исправления ошибок передачи информации. Увеличение минимального кодового расстояния () в 17 раз в соответствии с теорией помехоустойчивого кодирования обеспечивает возможность исправления до 4 ошибок в передаваемых словах или сообщениях. Однако условие  = 17 будет выполняться не всегда, а только на временных интервалах графических фрагментов результатов кодирования =  (), заключённых между разрывами, определяемыми на основе выполнения следующих неравенств, которые были отмечены ранее: Δ = | ≥ 0,8 × 2 и Δ= | ≥ 0,8 × 2. В результате этого получим  дополнительно закодированных значений , которые в случае отсутствия помех дадут при делении на значение  = 17 один и тот же остаток ζ:  ≡ ζ(mod ). Появившийся при этом инвариант в виде «группового свойства равноостаточности»: ζ= Const используют для контроля достоверности и целостности передаваемой ТМИ, а также для обнаружения и исправления ошибок передачи. При этом ошибку передачи εобнаруживают, если полученное значение ζбудет отличаться от наиболее встречающегося в гистограмме распределения значения, называемого  Мо(ζ).
Следующий этап подавления искажений данных телеизмерений, подвергнутых САП и искажённых помехой, связан, как уже отмечалось, связан с первичной обработкой полученной ТМИ. При этом в качестве примера рассмотрен случай использования предлагаемого модифицированного способа динамической фильтрации измерительных данных. Ориентация на динамическую фильтрацию определена с тем, что в большинстве практических задач обработка и анализ измерительных данных выполняются в отношении систем или процессов, которые отнесены к классу динамических. Это означает, что их функционирование рассматривается во временной области ([8],Кузнецов В.И. Статистическая идентификация. Адаптивная фильтрация. - Промышленные АСУ и контроллеры, 2014, № 7. - с. 20-27).
Известно [8], что теоретическую основу методов оптимальной динамической фильтрации составляет фильтр Калмана, различные реализации которого ориентированы на решения конкретных практических задач, в том числе и в условиях априорной неопределенности параметров информационной обстановки. В их числе: неопределённость параметров распределения случайных возмущений в уравнениях моделей состояний, случайной помехи в уравнениях измерений, объёма и качества измерительных данных [8].
В общем виде уравнения фильтрации в отношении динамической системы для заданных начальных условий  и матрицы ковариаций вектора  и ковариационной матрицы случайных возмущений , определяются следующей моделью состояний:
,                         (9)
где ,  - переходные матрицы состояний и возмущений .
Она дополняется моделью случайной помехи с ковариационной матрицей  в результатах измерений, которая описывается следующей моделью измерений:
,                                  (10)
где  - матрица измерений;  - случайная зашумляющая помеха.
В результате этого могут быть получены следующие зависимости, которые могут быть использованы для решения следующих задач:
- прогнозирования вектора состояния
(11)
и его ковариационной матрицы
(12)
- расчёта матричного коэффициента усиления фильтра
(13)
- определения отклонения прогнозируемой оценки от измеренного значения
; (14)
- расчёта текущей оценки  вектора состояний
(15)
и ковариационной матрицы погрешностей оценок
. (16)
Представленные уравнения фильтрации подводят к выводу о том, что основные трудности их реализации при непосредственном использовании результатов измерений, подверженных САП, без их обратных преобразований, т.е. без приведения их к исходным представлениям, связаны с необходимостью однозначной интерпретации отклонений прогнозируемой оценки вектора состояний от измеренного значения. С позиций математической интерпретации это положение следует записать в следующем аналитическом виде:
, (17)
X
m
j
откуда естественным образом возникает необходимость введения отображения, переводящего отклонения  в их требуемые значения , т.е. , которые появляются при нетрадиционном кодировании данных их образами-остатками, полученными в результате сравнений исходных значений слов-измерений по модулю .
Такое отображение может быть представлено в следующем виде:
(18)
b
j
Оно позволяет учесть особенности циклического перевода значений измеряемого параметра при вычислении его остатка  по модулю . Введенное отображение отличается простотой реализации и может быть встроено непосредственно в уравнения фильтрации.
Другой вариант преобразований, который реализуют в предлагаемом способе, рассматривается с использованием модели измерений, общий вид которой должен быть представлен с учётом выполняемых САП:
. (19)
В этом случае определение отклонения прогнозируемой оценки от измеренного значения будет описываться соотношением
,                          (20)
что не требует применения дополнительного отображения, рассмотренного ранее.
(k/k-1
b
b
n
N
b
b
(k/k-1
1
3
1
3
n
n
n
n
Второй вариант, использующий алгоритм дополнительного кодирования ПСАП-1, позволяет реализовать более лаконичные алгоритмы фильтрации, поскольку обеспечивает возможности использования закодированных измерительных данных, в том числе и без использования алгоритмов обратных САП (ОСАП-1). Результаты представления дисперсий шумовой помехи D[ε)] для случаев использования АНФ при традиционном представлении БМП (синий цвет) и при представлении образами-остатками(mod 2 - 1), (mod 2 + 1), где = /2 (красный цвет), приведены на фиг. 5. Полученные результаты позволяют сделать вывод о том, что при втором варианте обработки, когда АНФ подвергали данных представленные образами-остатками (mod 2 - 1), (mod 2 + 1) значения дисперсий шумовой помехи D[ε )] в 2-3 раза меньше. Следовательно, и точность восстановления переданных значений БМП будет выше.
В результате этого появляется возможность совместного использование различных технологий помехоустойчивого кодирования и обработки информации в СОК.
б
Большие объемы данных и быстро увеличивающиеся скорости передачи информации сопровождаются ухудшением показателей достоверности ее приема в виде вероятностей искажения бит (Р). Из-за этого она достаточно часто оказывается не пригодной для последующего использования, например, для анализа причин появления нештатных и аварийных ситуаций или результатов летных испытаний контролируемых объектов. Но в случае нештатных и аварийных ситуаций ее ценность многократно увеличивается: принятая ТМИ является уникальной, прежде всего потому, что повторное ее получение не представляется возможным. В результате причинно-следственные связи, которые привели к нештатным ситуациям и авариям, как правило, остаются нераскрытыми. Анализ такой типичной для существующей практики телеизмерений ситуации осложнен также тем, что в условиях нештатных и аварийных ситуаций значительно ухудшается качество канала связи, в результате чего существенно увеличивается число привнесенным им ошибок. Получается, что чем выше ценность информации, тем хуже ее качество. Такова объективная реальность. Возможность разрешения этого противоречия связана с обнаружением и исправлением ошибок, которые привносит радиоканал, а также с совершенствованием существующих способов обработки информации. Предлагаемый способ позволяет объединить эти два направления совершенствования существующей практики передачи и обработки информации на основе единого научно-методического подхода, связанного с нетрадиционным представлением данных, слов-измерений или сообщений их образами-остатками.
Способ обеспечивает возможность реализации наиболее простого технического решения модернизации существующей бортовой телеметрической аппаратуры (БТА), что может быть продемонстрировано на примере БРТС «ОРБИТА-IVMO» [4]. В ней для передачи быстроменяющихся параметров (БМП) используют следующие типы каналов, представленные в таблице 1.
Таблица 1 - Типы каналов БА «ОРБИТА-IVMO», обозначения каналов и номера разрядов двоичного кода в слове телеизмерений, занимаемых каналом при передаче БМП
Характеристика входного контролируемого параметра (канала)
Частота опросаf<sub>o </sub>= 512 Гц
Тип (Т) входного слова
Число БМП в канале
Количество разрядов в канале
Номер канала в слове «Р» и номера разрядов двоичного кода в слове телеизмерений, занимаемых каналами
Аналоговый, широкополосный(быстроменяющийся)
БМП
Т21Т22Т23
123
864
(со 2-го по 9-й), 1 разряд – для дополнительного разрешения неоднозначности значений 0 и 16 при сравнении по mod 2<sup>4</sup> +1 (с 1-го по 6-й разряды) – значения первого БМП(с 7 по 12 разряды) - значения второго БМП (2,3,4,1 разряды) – значения первого БМП (8,9,11,10 разряды) - значения второго БМП (5,6,7,12 разряды) – значения третьего БМП
8
6
4
В БТА «ОРБИТА-IVMO» [4] предусмотрены следующие 3 варианта формирования передаваемых значений БМП: 1) данные БМП представлены в одном канале 8-миразрядным двоичным словом (Т21); 2) данные двух БМП представлены в одном канале двумя 6-тиразрядными кодовыми словами (Т22) и 3) данные трёх БМП представлены в одном канале тремя 4-хразрядными кодовыми словами (Т23). В перечисленных случаях погрешности квантования значений БМП по амплитуде будут равны: 1) 1:2 × 100%= 0,39%; 2) 1:2 × 100%= 1,56% и 3) 1:2 × 100%= 6,25%. Это означает, что, увеличение числа передаваемых БМП в одном телеметрическом канале в БРТС может быть обеспечено только ценой увеличения погрешности телеизмерений. При этом общая погрешность восстановления БМП на приёмной стороне не должна превышать 7%. Тогда как в третьем случае только одна её составляющая, связанная с процессом квантования аналоговых значений по амплитуде и равная 6,25%, оказывается сравнимой с предельным значением погрешности.
b
b
n
N
b
b
b
n
b
b
b
b
b
b
b
1
j
3
j
3
j
1
j
3
j
2
10
2
10
1
j
3
j
3
j
2
10
1
j
2
10
3
j
2
10
3
j
2
10
1
j
2
10
n
n
4
n
n
3
3
n
3
3
3
3
1
1
1
При использовании предлагаемого способа (mod 2 - 1), (mod 2 + 1), где = /2 при передаче значений БМП байтовыми словами имеются дополнительные резервы для разрешения неоднозначности значений 0 и 16 образа-остатка (mod 2 +1). В случае использования второго варианта передачи двух БМП под представление (mod 2 - 1) и (mod 2 + 1) каждого из них имеется по  = 6/2 =3 разряда. При этом резервов для разрешения неоднозначности двух кодовых конструкций, соответствующих значениям <000>= <0> и  <000> = <8>, таких, как было в первом варианте представления значений БМП байтовыми словами, нет. В этом случае для разрешения неоднозначности используют сравнение (mod 2 - 1) и(mod 2 + 1). Когда (mod 2 + 1) = <000>= <0>, то такое же значение принимает и (mod 2 - 1) = <000>= <0>. Но в случае, требующем идентификации значения (mod 2 + 1) = <000>= <8> при отсутствии четвертого двоичного разряда, значение (mod 2 + 1) = <000>= <8> будет идентифицировано, если (mod 2 - 1) будет равно: <001>= <1>.
b
b
1
j
3
j
3
3
При этом погрешность восстановления БМП будет уменьшена на этапе обработки в СОК, например, при использовании алгоритмов АНФ применительно к каждому из остаточных процессов в отдельности, определяемых остатками (mod 2 - 1) и (mod 2 + 1).
При этом могут быть использованы и другие преимущества нетрадиционного представления данных в системе остаточных классов (СОК) ([10], Торгашев В.А. Система остаточных классов и надежность ЦВМ. М.: Сов. Радио, 1973.- 120 с.; [11], Акушский И.Я., Юдицкий Д.И. Машинная арифметика в остаточных классах. М.: Сов. Радио, 1968.-140 с.
Источники литературы
1. Способ передачи информации и система для его осуществления, патент RU №2586605, приоритет от 22.03.2013г.
2. Способ передачи информации, патент RU №2609747, приоритет от 13.08.2017г.
3. Кукушкин С.С. Теория конечных полей и информатика: В 2т. - т.1: Методы и алгоритмы, классические и нетрадиционные, основанные на использовании конструктивной теоремы об остатках. - М.: МО РФ, 2003. - 284с.
4. Современная телеметрия в теории и на практике / Учебный курс», Спб.: Наука и Техника, 2007. - 672с, стр. 465).
5. Способ первичной обработки информации с обнаружением и исправлением ошибок передачи, патент RU № 2658795, приоритет от 30.05.2017г.).
6. Способ первичной обработки информации с использованием адаптивной нелинейной фильтрации данных измерений, патент RU № 2672392, приоритет от 27.06.2017г.).
7. Зюко А.Г. Помехоустойчивость и эффективность систем связи. - М.: Связь, 1972. - 360с.).
8. Кузнецов В.И. Статистическая идентификация. Адаптивная фильтрация. - Промышленные АСУ и контроллеры, 2014, № 7. - с. 20-27).
9. Способ определения дальности до объекта с источником излучения сигналов с разными частотами, патент RU №2607639, опубл. 27.07.2016г., бюл. №21.
10. Торгашев В.А. Система остаточных классов и надежность ЦВМ. М.: Сов. Радио, 1973. 120 с.
11. Акушский И.Я., Юдицкий Д.И. Машинная арифметика в остаточных классах. М.: Сов. Радио, 1968. 140 с.