RU2517316C1 - Adaptive digital predictor - Google Patents

Adaptive digital predictor Download PDF

Info

Publication number
RU2517316C1
RU2517316C1 RU2012150888/08A RU2012150888A RU2517316C1 RU 2517316 C1 RU2517316 C1 RU 2517316C1 RU 2012150888/08 A RU2012150888/08 A RU 2012150888/08A RU 2012150888 A RU2012150888 A RU 2012150888A RU 2517316 C1 RU2517316 C1 RU 2517316C1
Authority
RU
Russia
Prior art keywords
subunit
adder
input
output
block
Prior art date
Application number
RU2012150888/08A
Other languages
Russian (ru)
Inventor
Камиль Хабибович Гильфанов
Роберт Накибович Хизбуллин
Андрей Геннадиевич Магданов
Original Assignee
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (ФГБОУ ВПО "КГЭУ")
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (ФГБОУ ВПО "КГЭУ") filed Critical Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (ФГБОУ ВПО "КГЭУ")
Priority to RU2012150888/08A priority Critical patent/RU2517316C1/en
Application granted granted Critical
Publication of RU2517316C1 publication Critical patent/RU2517316C1/en

Links

Images

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

FIELD: information technology.
SUBSTANCE: in the adaptive digital predictor, both inputs of terms of the first adder of the quadratic prediction subunit are respectively connected to data outputs of the first and second second-derivative calculation subunits, and the output of the first adder is connected to the first input of the second adder, the second input of which is connected to the output of a linear prediction subunit.
EFFECT: reduced hardware costs when predicting an input process.
5 dwg

Description

Изобретение относится к автоматике и вычислительной технике и может быть использовано для прогнозирования стационарных и нестационарных случайных процессов, повышения качества и точности управления в цифровых системах контроля и наведения различных (в т.ч. баллистических) объектов.The invention relates to automation and computer technology and can be used to predict stationary and non-stationary random processes, improve the quality and accuracy of control in digital control systems and guidance of various (including ballistic) objects.

Известно устройство для адаптивной экстраполяции (прогноза) по авт. св. СССР №1572281, МПК G06F 15/353, 1988), содержащий блок сглаживания, блок экстраполяции (прогноза), содержащий три вычитателя, регистр задания времени прогноза, выходной сумматор расчета квадратичного прогноза и блок оценки первой производной. Устройство функционально ограничено и имеет сложную схему коммутации узлов и блоков.A device for adaptive extrapolation (forecast) according to ed. St. USSR No. 1572281, IPC G06F 15/353, 1988), comprising a smoothing unit, an extrapolation (forecast) unit, comprising three subtractors, a prediction time setting register, an output adder for calculating a quadratic forecast, and a first derivative estimation unit. The device is functionally limited and has a complex circuit switching nodes and blocks.

Наиболее близким по технической сущности к заявляемому устройству является, выбранное в качестве прототипа, адаптивное цифровое прогнозирующее устройство (по патенту РФ №2451328, МПК G06F 15/00, 20.05.2012, бюл. №14), содержащий блок сглаживания и блок прогноза, в состав которого входят: три вычитателя, узел управления динамикой прогноза, два субблока квадратичного и линейного прогнозов, два субблока расчета второй призводной в (n-1)-й и (n-2)-й расчетных точках предыстории входного процесса и блок адаптации. Это устройство имеет относительно большой объем оборудования.The closest in technical essence to the claimed device is the adaptive digital predictive device (selected by the prototype of the Russian Federation No. 2451328, IPC G06F 15/00, 05/20/2012, bull. No. 14), selected as a prototype, containing a smoothing unit and a forecast unit, in the structure of which includes: three subtractors, a forecast dynamics control unit, two subblocks of quadratic and linear forecasts, two subblocks of the second lead calculation at the (n-1) th and (n-2) th calculation points of the input process history and an adaptation block. This device has a relatively large amount of equipment.

Техническая задача для предлагаемого устройства заключается в упрощении устройства, т.е. снижении аппаратурных затрат при его реализации.The technical problem for the proposed device is to simplify the device, i.e. reduction of hardware costs during its implementation.

Поэтому в адаптивном цифровом прогнозирующем устройстве, в состав которого входят: блок сглаживания, содержащий сумматор, первый и второй реверсивные счетчики, одноканальный субблок сглаживания из последовательно соединенных сумматора и регистра, субблок задания соотношения отклонений, содержащий регистр, счетчик и элемент задержки, субблок действительных отклонений, содержащий блок инверторов, два компаратора и элемент И, субблок единичных приращений, содержащий два элемента И и инвертор, субблок управления динамической характеристикой (СУДХ), содержащий два формирователя импульсов, элемент ИЛИ, счетчик, три элемента И и триггер режима, информационный выход блока сглаживания, информационный, первый управляющий и тактирующий входы устройства; узел тактирования блока прогноза, содержащий элемент задержки, триггер, генератор импульсов, элемент И и регистр сдвига; блок прогноза, содержащий первый, второй и третий вычитатели, каждый из которых содержит блок регистровой памяти, мультиплексор, блок инверторов и сумматор, субблок квадратичного прогноза, содержащий два последовательно соединенных сумматора, субблок линейного прогноза из одного сумматора, выходные шины которого монтажно сдвинуты на один разряд в сторону младших разрядов, узел управления динамикой прогноза, содержащий регистр хранения адреса ординат (расчетных точек) предыстории входного процесса, вход которого является вторым управляющим входом устройства задающим время (интервал) прогноза, компаратор, инвертор, элемент И и счетчик адреса, информационный выход которого заведен на адресные шины мультиплексоров всех трех вычитателей, причем две входные шины узла (шина сброса в «0» счетчика адреса и шина его восстановления) соединены с выходами двух элементов И (второго и третьего) СУДХ блока сглаживания, первый субблок расчета второй производной ( y n 1 " )

Figure 00000001
во второй (n-1)-й расчетной точке предыстории входного процесса из блока инверторов и сумматора, второй субблок расчета второй производной ( y n 2 " )
Figure 00000002
в третьей (n-2)-й расчетной точке предыстории входного процесса из двух последовательно соединенных сумматоров; блок адаптации, содержащий элемент ИСКЛЮЧАЮЩЕЕ ИЛИ, на оба входа которого заведены шины знаковых разрядов выходных сумматоров обоих субблоков расчета вторых производных, инвертор, первый и второй элементы И, триггер и мультиплексор, оба информационных входа которого подключены к выходам субблоков квадратичного и линейного прогнозов, а выход мультиплексора является первым информационным выходом устройства, для решения поставленной задачи на оба входа слагаемых первого сумматора субблока квадратичного прогноза заведены соответственно информационные выходы первого и второго субблоков расчета вторых производных, а выход первого сумматора подключен к первому входу второго сумматора, второй вход которого соединен с выходом субблока линейного прогноза.Therefore, in an adaptive digital predictive device, which includes: a smoothing unit containing an adder, first and second reversible counters, a single-channel smoothing subunit from a series-connected adder and register, a deviation ratio setting subunit containing a register, counter and delay element, a valid deviation subunit containing a block of inverters, two comparators and an element And, a subunit of unit increments, containing two elements And and an inverter, a subunit of control of a dynamic nature Ikoyi (sudha) comprising two of the pulse, an OR gate, a counter, three AND gates and the trigger mode, a smoothing information output unit information, the first control inputs and the timing device; a timing unit of a forecast block comprising a delay element, a trigger, a pulse generator, an AND element, and a shift register; a prediction block containing the first, second and third subtracters, each of which contains a register memory block, a multiplexer, an inverter block and an adder, a quadratic prediction block containing two series-connected adders, a linear prediction subunit from one adder, the output buses of which are assembled offset by one bit towards the lower digits, the forecast dynamics control unit, containing a register for storing the ordinates of the ordinates (calculated points) of the history of the input process, the input of which is the second the input of the device that sets the forecast time (interval), a comparator, an inverter, an AND element, and an address counter, the information output of which is connected to the address buses of the multiplexers of all three subtractors, moreover, two input buses of the node (reset bus to “0” of the address counter and its recovery bus ) are connected to the outputs of two elements AND (second and third) SUDH smoothing unit, the first subunit of the calculation of the second derivative ( y n - one " )
Figure 00000001
in the second (n-1) -th calculation point of the history of the input process from the block of inverters and the adder, the second sub-block of the calculation of the second derivative ( y n - 2 " )
Figure 00000002
in the third (n-2) -th design point of the input process history from two series-connected adders; an adaptation unit containing an EXCLUSIVE OR element, on both inputs of which there are sign buses of the output adders of both subunits for calculating second derivatives, an inverter, the first and second I elements, a trigger and a multiplexer, both of which information inputs are connected to the outputs of the quadratic and linear prediction subunits, and the multiplexer output is the first information output of the device, for solving the problem on both inputs of the terms of the first adder of the quadratic prediction subunit, respectively about the information outputs of the first and second subunits for calculating second derivatives, and the output of the first adder is connected to the first input of the second adder, the second input of which is connected to the output of the linear prediction subunit.

Сущность изобретения поясняется чертежами, на которых изображены: на фиг.1 - блок-схема предлагаемого устройства; на фиг.2 - блок-схема блока сглаживания; на фиг.3 - блок-схема одноканального субблока сглаживания; на фиг.4 - блок-схема узла тактирования блока прогноза; на фиг.5 - блок-схема блока прогноза.The invention is illustrated by drawings, which depict: in Fig.1 - a block diagram of the proposed device; figure 2 is a block diagram of a smoothing unit; figure 3 is a block diagram of a single-channel anti-aliasing subunit; Fig. 4 is a block diagram of a timing unit of a forecast block; figure 5 is a block diagram of a forecast block.

Известны формулы операторов прогноза, полученные аналитически с помощью аппроксимирующих многочленов по четырем точкам (ординатам) предыстории входного случайного дискретного процесса по способу наименьших квадратов (Милн В.Э. Численный анализ. М.: «ИЛ», 1951, стр.212). Оператор прогноза по аппроксимирующему многочлену второй степени (квадратичному) по четырем точкам предыстории имеет вид:The predictor operator formulas are known that are obtained analytically using approximating polynomials at four points (ordinates) of the history of the input random discrete process using the least squares method (Milne V.E. Numerical analysis. M.: IL, 1951, p. 212). The prediction operator for an approximating polynomial of the second degree (quadratic) for four historical points has the form:

y п + 1 = 1 4 ( 9 y п 3 y п 1 5 y п 2 + 3 y п 3 ) , ( 1 ) [ К В 4 ]

Figure 00000003
y P + one = one four ( 9 y P - 3 y P - one - 5 y P - 2 + 3 y P - 3 ) , ( one ) [ TO AT four ]
Figure 00000003

или в первом приближении:or as a first approximation:

y п + 1 = 1 4 ( 8 y п 2 y n 1 4 y n 2 + 2 y n 3 ) , ( 2 ) [ К В 4 ]

Figure 00000004
y P + one = one four ( 8 y P - 2 y n - one - four y n - 2 + 2 y n - 3 ) , ( 2 ) [ TO AT four ]
Figure 00000004

Оператор прогноза по аппроксимирующему многочлену первой степени (линейному) по четырем точкам предыстории имеет вид:The prediction operator for an approximating polynomial of the first degree (linear) for four historical points has the form:

y п + 1 = 1 2 ( 2 y п + y п 1 y п 3 ) , ( 3 ) [ Л Н 4 ]

Figure 00000005
y P + one = one 2 ( 2 y P + y P - one - y P - 3 ) , ( 3 ) [ L N four ]
Figure 00000005

где yп - первая (текущая) расчетная точка (ордината);where y p - the first (current) calculated point (ordinate);

yп-1, yп-2, yп-3 - соответственно вторая, третья и четвертая расчетные точки (ординаты) трехуровневой предыстории входной сглаженной дискретной последовательности. В численном анализе - это система равноотстоящих точек с шагом h, в реальном масштабе времени h - интервал между точками (ординатами), т.е. время (глубина) прогноза (H).y p-1 , y p-2 , y p-3 - respectively, the second, third and fourth design points (ordinates) of the three-level history of the input smoothed discrete sequence. In numerical analysis, this is a system of equally spaced points with a step h; in real time, h is the interval between points (ordinates), i.e. forecast time (depth) (H).

По аналогии с методами исчисления конечных разностей для численного дифференцирования, интерполяции и экстраполяции обозначим:By analogy with methods for calculating finite differences for numerical differentiation, interpolation, and extrapolation, we denote:

∆y1=(2yп-yп-1) - как биразность первого уровня предыстории входной дискретной последовательности, т.е. разность между удвоенной текущей и предыдущей ординатой процесса;Δy 1 = (2y p -y p-1 ) - as the biodiversity of the first level of the history of the input discrete sequence, i.e. the difference between the doubled current and previous ordinates of the process;

∆y2=(2yп-1-yп-2) - биразность второго уровня предыстории;Δy 2 = (2y p-1 -y p-2 ) - biodiversity of the second level of history;

∆y3=(2yп-1-yп-3) - биразность третьего уровня предыстории.∆y 3 = (2y p-1 -y p-3 ) - biodiversity of the third level of history.

После преобразования уравнений (2) и (3) с целью упрощения и с учетом биразностей получим следующие эмпирические выражения для операторов квадратичного и линейного прогнозов:After transforming equations (2) and (3) in order to simplify and taking into account biodiversity, we obtain the following empirical expressions for the quadratic and linear prediction operators:

y n + 1 = Δ y 1 + Δ y 2 1 2 ( Δ y 3 + y n 1 ) , ( 4 ) [ К В 4 ]

Figure 00000006
y n + one = Δ y one + Δ y 2 - one 2 ( Δ y 3 + y n - one ) , ( four ) [ TO AT four ]
Figure 00000006

y n + 1 = 1 2 ( Δ y 1 + Δ y 31 ) , ( 5 ) [ Л Н 4 ]

Figure 00000007
y n + one = one 2 ( Δ y one + Δ y 31 ) , ( 5 ) [ L N four ]
Figure 00000007

Известны формулы численного дифференцирования для равноотстоящих точек, выраженные через значения функции в этих точках (Дж.Поллард. Справочник по вычислительным методам статистики. М.: «ФС», 1982, стр.62), в частности для четырех точек имеем:Known formulas for numerical differentiation for equally spaced points expressed in terms of the function values at these points (J. Pollard. Handbook of Computational Methods of Statistics. M .: FS, 1982, p. 62), in particular for four points we have:

y n 1 " = 1 h 2 ( y n 2 y n 1 + y n 2 ) , ( 6 )

Figure 00000008
y n - one " = one h 2 ( y n - 2 y n - one + y n - 2 ) , ( 6 )
Figure 00000008

y n 2 " = 1 h 2 ( y n 1 2 y n 2 + y n 3 ) , ( 7 )

Figure 00000009
y n - 2 " = one h 2 ( y n - one - 2 y n - 2 + y n - 3 ) , ( 7 )
Figure 00000009

С учетом биразностей формулы (6) и (7) перепишутся в следующем виде:Taking into account biodiversity, formulas (6) and (7) will be rewritten as follows:

y n 1 " = 1 h 2 ( y n 2 y n 1 + y n 2 ) = 1 h 2 ( y n Δ y 2 ) , ( 8 )

Figure 00000010
y n - one " = one h 2 ( y n - 2 y n - one + y n - 2 ) = one h 2 ( y n - Δ y 2 ) , ( 8 )
Figure 00000010

y n 2 " = 1 h 2 ( y n 1 2 y n 2 + y n 3 ) , ( 9 )

Figure 00000011
y n - 2 " = one h 2 ( y n - one - 2 y n - 2 + y n - 3 ) , ( 9 )
Figure 00000011

Предлагаемое устройство реализует оператор прогноза по формуле (5) и расчет вторых производных по уравнениям (8) и (9), причем основными элементами схемы являются сумматор и блок инверторов, а умножение коэффициентов на слагаемые выполняются соответствующими монтажными сдвигами шин последних при вводе в сумматор. Такие операции на блок-схеме (см. фиг.5) обозначены кружочком.The proposed device implements the prediction operator according to formula (5) and the calculation of the second derivatives according to equations (8) and (9), the main elements of the circuit being the adder and the inverter unit, and the multiplication of the coefficients by the terms is carried out by the corresponding mounting shifts of the latter buses when entering into the adder. Such operations on the block diagram (see figure 5) are indicated by a circle.

Вычислим сумму трех уравнений (5), (6) и (7), причем в двух последних используем только модули (без коэффициентов размерности 1 h 2

Figure 00000012
):We calculate the sum of the three equations (5), (6) and (7), and in the last two we use only modules (without dimension coefficients one h 2
Figure 00000012
):

y n + 1 [ Л Н 4 ] + [ y n 1 " ] + [ y n 2 " ] = Δ y 1 + Δ y 2 1 2 ( Δ y 3 + y n 1 )

Figure 00000013
или y n + one [ L N four ] + [ y n - one " ] + [ y n - 2 " ] = Δ y one + Δ y 2 - one 2 ( Δ y 3 + y n - one )
Figure 00000013
or

y n + 1 = y n + 1 [ Л Н 4 ] + [ y n 1 " ] + [ y n 2 " ] . ( 10 ) [ К В 4 ]

Figure 00000014
y n + one = y n + one [ L N four ] + [ y n - one " ] + [ y n - 2 " ] . ( 10 ) [ TO AT four ]
Figure 00000014

Полученный результат суммирования (10) идентичен уравнению (4) для квадратичного прогноза yn+i [КВ4], однако для реализации уравнения (10) потребуется на один сумматор меньше, т.е. решается задача упрощения устройства.The obtained summation result (10) is identical to equation (4) for the quadratic forecast yn + i [КВ4], however, to implement equation (10), one adder less is required, i.e. solves the problem of simplifying the device.

Устройство содержит (см. фиг.1) блок сглаживания 1, блок прогноза 2 и блок адаптации 3. Блок сглаживания 1 (см. патент РФ №2444123, МПК H03H 17/04, 27.02.2012, Бюл. №6) содержит (см. фиг.2) сумматор 4, субблок 5 действительных отклонений, содержащий блок инверторов 6, два компаратора 7.1 и 7.2 и элемент И 8, первый реверсивный счетчик 9, субблок 10 задания соотношения отклонений, содержащий регистр 11, счетчик 12 и элемент задержки 13, субблок 14 единичных приращений, содержащий инвертор 15 и два элемента И 16.1 и 16.2, второй реверсивный счетчик 17, субблок 18 управления динамической характеристикой, содержащий два формирователя импульсов 19.1 и 19.2, элемент ИЛИ 20, счетчик 21, три элемента И 22.1, 22.2, 22.3 и триггер режима 23; информационный вход 24 блока сглаживания и устройства, первый управляющий 25 и тактирующий 26 входы устройства и блока сглаживания; одноканальный субблок сглаживания 27 (см. авт. св. СССР №748417, кл. G06F 15/32, 1980), содержащий (см. фиг.3) последовательно соединенные сумматор 28 и регистр 29; информационный выход 30. Узел тактирования 31 блока прогноза содержит (см. фиг.4) элемент задержки 32, триггер 33, генератор импульсов 34, элемент И 35 и регистр сдвига 36. Блок прогноза 2 (см. фиг.5) содержит первый 37, второй 38 и третий 39 вычитатели, каждый из которых содержит блок регистровой памяти 40 из (А) последовательно соединенных регистров 41, мультиплексор 42, блок инверторов 43 (в предположении, что мультиплексор не имеет инверсных выходов) и сумматор 44; субблок 45 квадратичного прогноза, содержащий первый 46 и второй 47 сумматоры, выход последнего является информационным выходом субблока; субблок 48 линейного прогноза, содержащий сумматор 49; узел 50 управления динамикой прогноза, содержащий регистр 51 хранения адреса (А) ординат расчетных точек предыстории процесса, вход 52 которого является вторым управляющим входом устройства, задающим время прогноза Н=АТ (Т - цикл работы устройства, А - максимальный адрес регистра памяти 41 блока 40), компаратор 53, инвертор 54, элемент И 55 и счетчик адреса 56; субблок 57 расчета второй производной во второй (n-1)-й расчетной точке предыстории входного процесса, состоящий из блока инверторов 58, сумматора 59 и второго информационного выхода 60 устройства, субблок 61 расчета второй производной в третьей (n-2)-й расчетной точке предыстории входного процесса, содержащий первый сумматор 62, второй сумматор 63 и третий информационный выход 64 устройства. Блок адаптации 3 содержит (см. фиг.5) элемент ИСКЛЮЧАЮЩЕЕ ИЛИ 65, инвертор 66, два элемента И 67.1 и 67.2, триггер 68, мультиплексор 69 и первый информационный выход 70 устройства.The device contains (see Fig. 1) a smoothing block 1, a forecast block 2 and an adaptation block 3. The smoothing block 1 (see RF patent No. 2444123, IPC H03H 17/04, February 27, 2012, Bull. No. 6) contains (see 2) an adder 4, a subunit 5 of actual deviations, comprising a block of inverters 6, two comparators 7.1 and 7.2 and an element And 8, a first reversible counter 9, a subunit 10 for setting the deviation ratio, comprising a register 11, a counter 12 and a delay element 13, subunit 14 of unit increments, comprising an inverter 15 and two elements And 16.1 and 16.2, a second reversible counter 17, a dynamic control subunit 18 a characteristic comprising two pulse shapers 19.1 and 19.2, an OR element 20, a counter 21, three AND elements 22.1, 22.2, 22.3 and a mode 23 trigger; information input 24 of the smoothing unit and the device, the first control 25 and clocking 26 inputs of the device and the smoothing unit; single-channel smoothing subunit 27 (see ed. St. USSR No. 748417, class G06F 15/32, 1980), comprising (see FIG. 3) adder 28 and register 29 connected in series; information output 30. The timing unit 31 of the forecast block contains (see FIG. 4) a delay element 32, a trigger 33, a pulse generator 34, an And 35 element and a shift register 36. Prediction block 2 (see FIG. 5) contains the first 37, the second 38 and third 39 subtractors, each of which contains a block of register memory 40 from (A) series-connected registers 41, a multiplexer 42, a block of inverters 43 (assuming that the multiplexer does not have inverse outputs) and an adder 44; a quadratic prediction subunit 45 containing the first 46 and second 47 adders, the output of the latter is the information output of the subunit; a linear prediction subunit 48 comprising an adder 49; node 50 control the forecast dynamics, containing a register 51 storing the address (A) of the ordinates of the calculated points of the process history, the input 52 of which is the second control input of the device that sets the forecast time H = AT (T is the cycle of the device, A is the maximum address of the memory register 41 blocks 40), a comparator 53, an inverter 54, an AND element 55, and an address counter 56; a subunit 57 for calculating the second derivative in the second (n-1) -th calculation point of the input process history, consisting of a block of inverters 58, an adder 59 and a second information output 60 of the device, a subunit 61 for calculating the second derivative in the third (n-2) -th calculation point in the history of the input process, containing the first adder 62, the second adder 63 and the third information output 64 of the device. The adaptation unit 3 contains (see FIG. 5) an element EXCLUSIVE OR 65, an inverter 66, two elements AND 67.1 and 67.2, a trigger 68, a multiplexer 69 and a first information output 70 of the device.

Цикл работы устройства состоит из 5 тактов. Блок сглаживания 1 работает в двух режимах: стационарном и динамическом (переходном), причем все операции выполняются в нем за один (1-й) такт. На стационарном режиме блок сглаживает входной случайный дискретный процесс, детерминированная основа (медиана) которого может иметь постоянный, линейный или нелинейный (квадратичный) характер изменения во времени. Блок сглаживания 1 (см. фиг.2) реализует следующую модификацию оператора сигнатурного экспоненциального сглаживания:The cycle of the device consists of 5 cycles. The smoothing unit 1 operates in two modes: stationary and dynamic (transitional), and all operations are performed in it in one (1st) cycle. In stationary mode, the unit smooths the input random discrete process, the deterministic basis (median) of which can have a constant, linear or non-linear (quadratic) character of change in time. The smoothing unit 1 (see figure 2) implements the following modification of the operator of signature exponential smoothing:

y n = y n 1 + s i g n [ 1 K ( x n y n 1 ) ] , ( 11 )

Figure 00000015
y n = y n - one + s i g n [ one K ( x n - y n - one ) ] , ( eleven )
Figure 00000015

где xn и yn - входная и выходная дискреты;where x n and y n are the input and output discrete;

α=1/K - параметр сглаживания, К - параметр адаптации;α = 1 / K is the smoothing parameter, K is the adaptation parameter;

∆xn=(xn-yn-1) - текущие отклонения от медианы процесса.∆x n = (x n -y n-1 ) - current deviations from the median of the process.

В качестве критерия эффективности (точности) сглаживания выбрано соотношение d между нулевыми и действительными отклонениями ∆xn. Последние формируют текущие единичные приращения обоих знаков выходной дискреты в соответствии с сигнатурной функцией в (11):As a criterion for the efficiency (accuracy) of smoothing, the ratio d between zero and actual deviations Δx n is chosen. The latter form the current unit increments of both signs of the output discrete in accordance with the signature function in (11):

sign[∆xn/K]=0 для [∆xn-K]<0 (∆xn - нулевые отклонения),sign [∆x n / K] = 0 for [∆x n -K] <0 (∆x n are zero deviations),

sign[∆xn/K]=1 для [∆xn-К]>0 (∆xn - действительные отклонения).sign [∆x n / K] = 1 for [∆x n -K]> 0 (∆x n are real deviations).

На стационарном режиме (D=0 - признак режима) блок 1 сглаживает входную случайную последовательность дискрет до уровня, заданного соотношением d (реальный диапазон d=7÷190), которое заносится перед началом работы устройства со входа 25 в регистр 11 субблока 10 задания соотношения отклонений. Последний представляет собой управляемый делитель частоты, например при d=7 на выходе прямого переноса счетчика 12 появляется каждый седьмой тактовый импульс со входа 26, который через элемент задержки 13 перезаписывает инверсный код d из регистра 11 в счетчик 12 (для следующего цикла работы делителя) и вычитает «1» из первого реверсивного счетчика 9, содержащего код параметра адаптации K. Процесс адаптивного сглаживания заключается в следующем. Пусть (при определенном коде K в счетчике 9) возросла дисперсия входного дискретного процесса, т.е. возросло число действительных отклонений ∆xn (обоих знаков). После сравнения их с параметром адаптации K на выходе блоков компараторов 7.1 и 7.2 субблока 5 (играющего роль отрицательной обратной связи) устанавливаются логические «1» (режим работы компараторов: [∆xn>K]=«1», [∆xn<K]=«0»), поступающие на вход элемента И 8. Так как на стационарном режиме триггер режима 23 находится в состоянии «0» (D=0), то с его инверсного выхода на первый вход элемента И 8 субблока 5 также поступает логическая «1». Высокий уровень сигнала на всех входах элемента И 8 разрешает прохождение тактовых импульсов с входа 26 на суммирующий вход первого реверсивного счетчика 9 (код К в последнем увеличивается) и на вторые входы элементов И 16.1 и 16.2 субблока единичных приращений 14. Сигнал с выхода одного из них (в зависимости от знака отклонения) поступает на суммирующий (или вычитающий) вход второго реверсивного счетчика 17 результата сглаживания, т.е. реализуется сигнатурная функция (11). Процесс роста K приведет к снижению числа действительных отклонений и будет продолжаться до тех пор, пока не наступит динамическое равновесие, т.е. число импульсов, поступивших от субблока 10 на вычитающий вход счетчика 9 будет равно числу импульсов, поступивших на его суммирующий вход от субблока 5, а дисперсия выходной сглаженной дискретной последовательности останется неизменной (для d=7: на одно действительное отклонение должно приходиться семь нулевых).In stationary mode (D = 0 - a sign of mode), block 1 smoothes the input random sequence of discretes to the level specified by the relation d (real range d = 7 ÷ 190), which is entered before starting the device from input 25 to register 11 of subunit 10 of the ratio deviations. The latter is a controllable frequency divider, for example, at d = 7, the output of the direct transfer of counter 12 appears every seventh clock pulse from input 26, which through the delay element 13 overwrites the inverse code d from register 11 to counter 12 (for the next cycle of the divider) and subtracts “1” from the first reverse counter 9 containing the adaptation parameter code K. The adaptive smoothing process is as follows. Let (for a certain code K in counter 9), the variance of the input discrete process increase, i.e. the number of actual deviations ∆x n (of both signs) increased. After comparing them with the adaptation parameter K at the output of the comparator blocks 7.1 and 7.2 of subunit 5 (playing the role of negative feedback), logical “1” is set (the mode of operation of the comparators: [∆x n > K] = “1”, [∆x n < K] = "0") received at the input of element And 8. Since in the stationary mode the trigger of mode 23 is in the state "0" (D = 0), from its inverse output to the first input of element And 8 of subunit 5 logical "1". A high level of the signal at all inputs of the And 8 element allows the passage of clock pulses from input 26 to the summing input of the first reverse counter 9 (the K code in the last increases) and to the second inputs of the And 16.1 and 16.2 elements of the subunit of unit increments 14. The signal from the output of one of them (depending on the sign of the deviation) is fed to the summing (or subtracting) input of the second reverse counter 17 of the smoothing result, i.e. the signature function (11) is implemented. The growth process K will lead to a decrease in the number of actual deviations and will continue until a dynamic equilibrium sets in, i.e. the number of pulses received from the subunit 10 to the subtracting input of the counter 9 will be equal to the number of pulses received at its summing input from the subunit 5, and the variance of the output smoothed discrete sequence will remain unchanged (for d = 7: there must be seven zero deviations for one real deviation).

Переходный (динамический) режим может быть вызван запуском, ускорением, виражом, переходом с одного стационарного режима на другой и т.д., т.е. почти скачкообразным изменением процесса. Для сглаживания входной дискретной последовательности на переходных (динамических) режимах (D=l) используется одноканальный субблок сглаживания 27 (см. фиг.3), который реализует следующий оператор экспоненциального сглаживания: y n = 1 2 ( x n + y n 1 )

Figure 00000016
, т.е. с минимальной степенью сглаживания и, соответственно, с минимальным фазовым сдвигом (запаздыванием) выходной дискреты. Субблок 27 работает на обоих режимах, инициируется тактовыми импульсами со входа 26 в регистре 29, но используется только на переходном (динамическом) режиме. Для стационарного случайного процесса вероятность появления серии, например, из m=8 (восьми) отклонений от медианы (детерминированной основы) процесса одного знака подряд в соответствии с геометрическим законом распределения вероятностей равна:Transient (dynamic) mode can be caused by starting, accelerating, bending, switching from one stationary mode to another, etc., i.e. an almost abrupt change in the process. To smooth the input discrete sequence in transient (dynamic) modes (D = l), a single-channel smoothing subunit 27 is used (see Fig. 3), which implements the following exponential smoothing operator: y n = one 2 ( x n + y n - one )
Figure 00000016
, i.e. with a minimum degree of smoothing and, accordingly, with a minimum phase shift (delay) of the output discrete. Subunit 27 operates in both modes, is initiated by clock pulses from input 26 in register 29, but is used only in transition (dynamic) mode. For a stationary random process, the probability of a series appearing, for example, from m = 8 (eight) deviations from the median (deterministic basis) of a single sign in a row in accordance with the geometric law of probability distribution is equal to:

P ( x = m ) = ( 1 2 ) m = 1 / 256 0,004

Figure 00000017
P ( x = m ) = ( one 2 ) m = one / 256 0.004
Figure 00000017

т.е. настолько мала, что можно считать появление такой серии началом переходного режима. Субблок 18 фиксирует такую серию и работает следующим образом. Так как для стационарного режима наиболее вероятны отклонения разных знаков, то при смене знака в сумматоре 4 с «плюс» на «минус» и наоборот срабатывают формирователи импульсов 19.1 или 19.2 и через элемент ИЛИ 20 сбрасывают в «0» счетчик 21 и триггер режима 23 (D=0). На динамическом режиме (формирователи 19 не срабатывают) на счетчик 21 (например, 4-разрядный) непременно поступит восемь импульсов подряд с тактового входа 26. На выходе старшего разряда счетчика 21 установится логическая «1», высокий уровень сигнала которого обеспечит прохождение через первый элемент И 22.1 тактирующего импульса, который установит триггер режима 23 в «1» (D=1). Последний сигналом с инверсного выхода заблокирует работу субблока 5 действительных отклонений и, соответственно, субблока 14 единичных приращений, а высоким уровнем сигнала прямого выхода разрешит через второй элемент И 22.2 перезапись дискрет с одноканального субблока сглаживания 27 во второй реверсивный счетчик 17 результата сглаживания.those. so small that the appearance of such a series can be considered the beginning of a transitional regime. Subunit 18 captures such a series and works as follows. Since for the stationary mode the deviations of different signs are most likely, when the sign in the adder 4 changes from “plus” to “minus” and vice versa, the pulse shapers 19.1 or 19.2 are triggered and through the element OR 20 they reset counter 21 and trigger 23 to “0” (D = 0). In the dynamic mode (shapers 19 do not work), eight pulses will certainly come to the counter 21 (for example, 4-bit) from the clock input 26. At the output of the high-order bit of the counter 21, a logical “1” will be set, the high level of which will ensure passage through the first element And 22.1 clock pulse, which sets the trigger of mode 23 to "1" (D = 1). The last signal from the inverse output will block the operation of the subunit 5 of the actual deviations and, accordingly, the subunit of 14 unit increments, and by the high level of the direct output signal it will allow the discretization from the single-channel smoothing subunit 27 to the second reversal counter 17 of the smoothing result through the second element And 22.2.

По окончании переходного режима в сумматоре 4 неизбежно возникнут отклонения разных знаков, что приведет к срабатыванию формирователей импульсов 19 и, соответственно, к переключению триггера режима 23 в состояние «0» (стационарный режим сглаживания, D=0).At the end of the transition mode, in the adder 4, deviations of different signs will inevitably occur, which will lead to the operation of the pulse shapers 19 and, accordingly, to the switching of the trigger of mode 23 to the state “0” (stationary smoothing mode, D = 0).

Операции прогнозирования выполняются за три такта, соответственно 2-й, 3-й и 4-й. Формируются они серией из трех тактирующих импульсов от узла тактирования 31 (см. фиг.4). Тактирующий импульс с входа 26 обнуляет триггер 33 и записывает «1» в младший разряд сдвигового регистра 36. Тот же тактирующий импульс, задержанный элементом задержки 32, устанавливает в «1» триггер 33, разрешая тем самым прохождение импульсов от генератора 34 через элемент И 35 в регистр сдвига 36, на шинах младших разрядов которого («а», «б», «в», «г») и появляется вышеуказанная серия. Во 2-м такте производится запись ординаты текущей (первой) расчетной точки yп в первый регистр 41 блока 40 регистровой памяти первого вычитателя 37. Одновременно происходит перезапись (сдвиг) всех предшествующих ординат в соседние регистры 41 (т.е. формируется предыстория входного процесса). На адресный вход мультиплексора 42 поступает код адреса (А) ординаты предыстории со счетчика адреса 58, равный коду адреса, записанному со второго управляющего входа 54 в регистр 53 хранения адреса перед началом работы устройства и определяющий время (интервал) прогноза Н=AT. В соответствии с этим адресом ордината с выхода мультиплексора 42 (уже как вторая расчетная точка yп-1) через блок инверторов 43 поступает на вход второго слагаемого сумматора 44, на входе первого слагаемого которого стоит удвоенная ордината предшествующей расчетной точки yп. На выходе сумматора первого вычитателя 37 устанавливается биразность 1-го уровня предыстории входной дискретной последовательности. В 3-м и 4-м тактах производятся операции, аналогичные описанным выше, но уже для второго 38 и третьего 39 вычитателей, на выходах которых устанавливаются соответственно биразности 2-го и 3-го уровней предыстории. Все сумматоры в устройстве - комбинационные. По завершении 4-го такта на выходе субблока 45 в соответствии с эмпирической формулой (10) устанавливается код оценки квадратичного (нелинейного) прогноза для нестационарной входной дискретной последовательности, на выходе субблока 48 в соответствии с формулой (5) - код оценки линейного прогноза для стационарной или медленно меняющейся входной дискретной последовательности, на выходе 60 субблока 57 в соответствии с формулой (8) - код оценки второй производной во второй (n-1)-й расчетной точке предыстории входного процесса, а на выходе 64 субблока 61 в соответствии с формулой (9) - код оценки второй производной в третьей (n-2)-й расчетной точке предыстории входного процесса.Forecasting operations are performed in three cycles, respectively, 2nd, 3rd and 4th. They are formed by a series of three clock pulses from the clock node 31 (see figure 4). A clock pulse from input 26 resets trigger 33 and writes “1” to the least significant bit of the shift register 36. The same clock pulse, delayed by delay element 32, sets trigger 33 to “1”, thereby allowing pulses from generator 34 to pass through element I 35 in the shift register 36, on the tires of the least significant bits of which ("a", "b", "c", "g") the above series appears. In the 2nd step, the ordinate of the current (first) calculated point y p is recorded in the first register 41 of the block 40 of the register memory of the first subtractor 37. At the same time, all previous ordinates are overwritten (shifted) to neighboring registers 41 (i.e., the background of the input process is formed ) The address input of the multiplexer 42 receives the address code of the history ordinate (A) from the address counter 58, which is equal to the address code recorded from the second control input 54 in the address storage register 53 before the device starts operation and determines the forecast time (interval) Н = AT. In accordance with this address, the ordinate from the output of the multiplexer 42 (already as the second calculation point y p-1 ) through the inverter block 43 is fed to the input of the second term of the adder 44, the first term of which is doubled the ordinate of the previous calculation point y p . The output of the adder of the first subtractor 37 establishes the biodiversity of the 1st level of the history of the input discrete sequence. In the 3rd and 4th steps, operations are performed similar to those described above, but for the second 38 and third 39 subtractors, the outputs of which are set respectively biodities of the 2nd and 3rd levels of history. All adders in the device are combination. At the end of the 4th cycle, at the output of subunit 45, in accordance with the empirical formula (10), a quadratic (non-linear) prediction estimation code for a non-stationary discrete input sequence is set, and at the output of subunit 48, in accordance with formula (5), a linear prediction estimation code for a stationary or a slowly varying input discrete sequence, the output 60 of subunit 57 in accordance with formula (8) is the code for evaluating the second derivative at the second (n-1) -th design point of the input process history, and the output is 64 subunits 61 according to formula (9) - the second derivative of code evaluation in the third (n-2) -th design point history input process.

В 5-м такте сигнал с выхода («г») регистра сдвига 36 узла тактирования 31 в зависимости от сочетания знаков вторых производных устанавливает триггер 68 блока адаптации 3 в состояние «1» (квадратичный прогноз) или в «0» (линейный прогноз). Прямой выход триггера 68, как адресный вход мультиплексора 69, обеспечивает выбор и передачу на его единственный выход 70 мультиплексора и устройства соответствующего кода оценки прогноза.In the 5th step, the signal from the output (“g”) of the shift register 36 of the clock unit 31, depending on the combination of the signs of the second derivatives, sets the trigger 68 of the adaptation unit 3 to the state “1” (quadratic forecast) or to “0” (linear forecast) . The direct output of the trigger 68, as the address input of the multiplexer 69, provides selection and transmission to its only output 70 of the multiplexer and the device of the corresponding forecast estimation code.

Узел управления динамикой прогноза 50 предназначен для исключения операции прогноза на динамических (переходных) режимах (D=l) работы устройства путем сброса в «0» счетчика адреса 56 тактирующим сигналом (U0) из субблока 18 управления динамической характеристикой блока сглаживания. Нулевой адрес счетчика 56 на адресных шинах всех трех мультиплексоров 42 вычитателей обеспечит на все время переходного режима расчет и установку на выходах обоих субблоков прогноза 45 и 48 кода текущей дискреты yп минимально сглаженного входного процесса. С переходом устройства на стационарный режим работы (D=0) триггер режима 23 субблока 18 разрешит прохождение тактирующего импульса со входа 26 через открытый элемент И 55 на счетный вход (Ua) счетчика адреса 56. С ростом адреса в последнем на выходах обоих субблоков прогноза 45 и 48 устанавливаются коды (yn+1) прогнозируемого входного процесса с использованием информации из трехуровневой предыстории только для нового стационарного режима. Рост кода адреса h в счетчике 56 (h=aT, а=1, 2, 3, … А), т.е. восстановление заданного времени прогноза H будет продолжаться с каждым циклом до тех пор, пока он не станет равным заданному в регистре хранения 51 h=H. Компаратор 53 (режим работы: [H=h]→«1», [H≠h]→«0») в этом случае через ключ 54 и элемент И 55 закроет счетный вход счетчика адреса 56.The forecast dynamics control unit 50 is designed to exclude the forecast operation in dynamic (transient) modes (D = l) of the device operation by resetting the address counter 56 to a “0” by a clock signal (U 0 ) from the sub-block 18 of the dynamic characteristic control of the smoothing block. The zero address of the counter 56 on the address buses of all three multiplexers 42 of the subtractor will ensure the calculation and installation of the current discrete code y n of a minimally smoothed input process at the outputs of both prediction subunits 45 and 48. With the transition of the device to the stationary mode of operation (D = 0), the trigger of mode 23 of subunit 18 will allow the clock pulse to pass from input 26 through open element And 55 to the counting input (U a ) of address counter 56. As the address increases in the latter at the outputs of both prediction subunits 45 and 48, codes (y n + 1 ) of the predicted input process are set using information from a three-level history only for the new stationary mode. The growth of the address code h in the counter 56 (h = aT, a = 1, 2, 3, ... A), i.e. the restoration of the specified forecast time H will continue with each cycle until it becomes equal to the value set in the storage register 51 h = H. The comparator 53 (operating mode: [H = h] → “1”, [H ≠ h] → “0”) in this case, through the key 54 and the AND 55 element, it closes the counting input of the address counter 56.

Использование обоих операторов прогноза для контроля, слежения или управления параметрами в технических системах базируется на фундаментальных физических законах: законе инерции, законах сохранения энергии и движения, инерционности процессов нагрева/охлаждения и др., что позволяет (опираясь на предысторию процесса) рассчитывать на высокую степень надежности прогноза. О точности прогноза можно судить только по завершении события и если в период времени (интервала) прогноза не было форс-мажорных обстоятельств: удар, скачок, взрыв и т.п.The use of both prediction operators for monitoring, tracking, or controlling parameters in technical systems is based on fundamental physical laws: the law of inertia, the laws of conservation of energy and motion, the inertia of heating / cooling processes, etc., which allows (based on the background of the process) to rely on a high degree reliability of the forecast. The accuracy of the forecast can be judged only at the end of the event and if during the period of time (interval) of the forecast there were no force majeure circumstances: blow, jump, explosion, etc.

Claims (1)

Адаптивное цифровое прогнозирующее устройство, в состав которого входят: блок сглаживания, содержащий сумматор, первый и второй реверсивные счетчики, одноканальный субблок сглаживания из последовательно соединенных сумматора и регистра, субблок задания соотношения отклонений, содержащий регистр, счетчик и элемент задержки, субблок действительных отклонений, содержащий блок инверторов, два компаратора и элемент И, субблок единичных приращений, содержащий два элемента И и инвертор, субблок управления динамической характеристикой (СУДХ), содержащий два формирователя импульсов, элемент ИЛИ, счетчик, три элемента И и триггер режима, информационный выход блока сглаживания, информационный, первый управляющий и тактирующий входы устройства; узел тактирования блока прогноза, содержащий элемент задержки, триггер, генератор импульсов, элемент И и регистр сдвига; блок прогноза, содержащий первый, второй и третий вычитатели, каждый из которых содержит блок регистровой памяти, мультиплексор, блок инверторов и сумматор, субблок квадратичного прогноза, содержащий два последовательно соединенных сумматора, субблок линейного прогноза из одного сумматора, выходные шины которого монтажно сдвинуты на один разряд в сторону младших разрядов, узел управления динамикой прогноза, содержащий регистр хранения адреса ординат (расчетных точек) предыстории входного процесса, вход которого является вторым управляющим входом устройства задающим время (интервал) прогноза, компаратор, инвертор, элемент И и счетчик адреса, информационный выход которого заведен на адресные шины мультиплексоров всех трех вычитателей, причем две входные шины узла (шина сброса в «0» счетчика адреса и шина его восстановления) соединены с выходами двух элементов И (второго и третьего) СУДХ блока сглаживания, первый субблок расчета второй производной
Figure 00000018
во второй (n-1)-й расчетной точке предыстории входного процесса из блока инверторов и сумматора, второй субблок расчета второй производной
Figure 00000019
в третьей (n-2)-й расчетной точке предыстории входного процесса из двух последовательно соединенных сумматоров; блок адаптации, содержащий элемент ИСКЛЮЧАЮЩЕЕ ИЛИ, на оба входа которого заведены шины знаковых разрядов выходных сумматоров обоих субблоков расчета вторых производных, инвертор, первый и второй элементы И, триггер и мультиплексор, оба информационных входа которого подключены к выходам субблоков квадратичного и линейного прогнозов, а выход мультиплексора является первым информационным выходом устройства, отличающееся тем, что на оба входа слагаемых первого сумматора субблока квадратичного прогноза заведены соответственно информационные выходы первого и второго субблоков расчета вторых производных, а выход первого сумматора подключен к первому входу второго сумматора, второй вход которого соединен с выходом субблока линейного прогноза.
An adaptive digital predictive device, which includes: a smoothing unit containing an adder, first and second reverse counters, a single-channel smoothing subunit from a series-connected adder and register, a deviation ratio setting subunit containing a register, a counter and a delay element, a valid deviation subunit containing a block of inverters, two comparators and an element And, a subunit of unit increments containing two elements And and an inverter, a subunit of control of the dynamic characteristic (COURT ) Containing two pulses driver, an OR gate, a counter, three AND gates and the trigger mode, a smoothing information output unit information, the first control inputs and the timing device; a timing unit of a forecast block comprising a delay element, a trigger, a pulse generator, an AND element, and a shift register; a prediction block containing the first, second and third subtracters, each of which contains a register memory block, a multiplexer, an inverter block and an adder, a quadratic prediction block containing two series-connected adders, a linear prediction subunit from one adder, the output buses of which are assembled offset by one bit towards the lower digits, the forecast dynamics control unit, containing a register for storing the ordinates of the ordinates (calculated points) of the history of the input process, the input of which is the second the input of the device that sets the forecast time (interval), a comparator, an inverter, an AND element, and an address counter, the information output of which is connected to the address buses of the multiplexers of all three subtractors, moreover, two input buses of the node (reset bus to “0” of the address counter and its recovery bus ) are connected to the outputs of two elements AND (second and third) SUDH smoothing unit, the first subunit of the calculation of the second derivative
Figure 00000018
in the second (n-1) -th calculation point of the history of the input process from the block of inverters and the adder, the second sub-block of the calculation of the second derivative
Figure 00000019
in the third (n-2) -th design point of the input process history from two series-connected adders; an adaptation unit containing an EXCLUSIVE OR element, on both inputs of which there are sign buses of the output adders of both subunits for calculating second derivatives, an inverter, the first and second I elements, a trigger and a multiplexer, both of which information inputs are connected to the outputs of the quadratic and linear prediction subunits, and the multiplexer output is the first information output of the device, characterized in that the information on the two inputs of the terms of the first adder of the quadratic prediction subunit tional outputs of the first and second subunits calculate the second derivatives, and the output of the first adder is connected to a first input of the second adder, a second input coupled to an output of the linear prediction subunit.
RU2012150888/08A 2012-11-27 2012-11-27 Adaptive digital predictor RU2517316C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2012150888/08A RU2517316C1 (en) 2012-11-27 2012-11-27 Adaptive digital predictor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2012150888/08A RU2517316C1 (en) 2012-11-27 2012-11-27 Adaptive digital predictor

Publications (1)

Publication Number Publication Date
RU2517316C1 true RU2517316C1 (en) 2014-05-27

Family

ID=50779445

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2012150888/08A RU2517316C1 (en) 2012-11-27 2012-11-27 Adaptive digital predictor

Country Status (1)

Country Link
RU (1) RU2517316C1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2622851C1 (en) * 2016-08-01 2017-06-20 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital predictive device
RU2622852C1 (en) * 2016-08-01 2017-06-20 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital smoothing and predictive device
RU2680217C1 (en) * 2017-12-25 2019-02-18 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Digital predictor
RU2707417C1 (en) * 2019-05-14 2019-11-26 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital predictive device
RU2713872C1 (en) * 2019-05-14 2020-02-07 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Digital predictive device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1120288A1 (en) * 1982-08-31 1984-10-23 Войсковая Часть 73790 Device for predicting random processes
EP1510943A1 (en) * 2002-05-31 2005-03-02 Celestar Lexico-Sciences, Inc. Interaction predicting device
RU2280227C1 (en) * 2004-12-15 2006-07-20 Государственное унитарное предприятие "Конструкторское бюро приборостроения" Method for command telecontrol of rocket
RU2451328C1 (en) * 2011-05-31 2012-05-20 Государственное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (КГЭУ) Adaptive digital predictor

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1120288A1 (en) * 1982-08-31 1984-10-23 Войсковая Часть 73790 Device for predicting random processes
EP1510943A1 (en) * 2002-05-31 2005-03-02 Celestar Lexico-Sciences, Inc. Interaction predicting device
RU2280227C1 (en) * 2004-12-15 2006-07-20 Государственное унитарное предприятие "Конструкторское бюро приборостроения" Method for command telecontrol of rocket
RU2451328C1 (en) * 2011-05-31 2012-05-20 Государственное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (КГЭУ) Adaptive digital predictor

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2622851C1 (en) * 2016-08-01 2017-06-20 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital predictive device
RU2622852C1 (en) * 2016-08-01 2017-06-20 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital smoothing and predictive device
RU2680217C1 (en) * 2017-12-25 2019-02-18 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Digital predictor
RU2707417C1 (en) * 2019-05-14 2019-11-26 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Adaptive digital predictive device
RU2713872C1 (en) * 2019-05-14 2020-02-07 Федеральное государственное бюджетное образовательное учреждение высшего образования "Казанский государственный энергетический университет" (ФГБОУ ВО "КГЭУ") Digital predictive device

Similar Documents

Publication Publication Date Title
RU2451328C1 (en) Adaptive digital predictor
RU2517316C1 (en) Adaptive digital predictor
RU2455682C1 (en) Digital predictor
RU2450343C1 (en) Digital predicting and differentiating device
RU2446454C1 (en) Digital predictor
CN101615794A (en) Electrical Power System Dynamic method for estimating state based on no mark transformation Kalman filter
RU2517317C1 (en) Adaptive digital predicting and differentiating device
RU2517322C1 (en) Adaptive digital predicting and differentiating device
RU2475831C1 (en) Adaptive digital predicting and differentiating device
RU2515215C1 (en) Digital predicting and differentiating device
RU2477887C1 (en) Digital predictor
RU2535467C1 (en) Adaptive digital differentiating and predicting device
CN104834216A (en) Binomial-based wireless sensor network trust management method
RU2470359C1 (en) Digital predicting and differentiating device
RU2449350C1 (en) Digital predicting and differentiating device
RU2580452C1 (en) Signature smoothing digital device
CN115062342A (en) Privacy calculation method and device, electronic equipment and readable storage medium
RU2622852C1 (en) Adaptive digital smoothing and predictive device
RU2459241C1 (en) Digital predictor
RU2680217C1 (en) Digital predictor
RU2446461C2 (en) Digital predictor
RU2626338C1 (en) Adaptive digital smoothing and predictive device
Rostam et al. A hybrid Gaussian process approach to robust economic model predictive control
RU2629641C1 (en) Digital predictor
CN102156802B (en) A kind of equally distributed fluctuation data predication method

Legal Events

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

Effective date: 20151128