CC BY
43
X УПРАВЛЕНИЕ В МЕДИЦИНЕ И БИОЛОГИИ
УДК 681.51 + 519.6
РАЗЛАДКА, ГОМЕОСТАЗИС, ИЗМЕРЕНИЕ В РАМКАХ КОМПЕНСАЦИОННОГО ПРИНЦИПА РАВНОВЕСИЯ В ДИНАМИЧЕСКИХ СИСТЕМАХ Часть 3: Варианты алгоритмов
Э. П. Тихонов,
доктор техн. наук, профессор
Санкт-Петербургский государственный электротехнический университет «ЛЭТИ»
Предлагаются и исследуются на базе общего подхода, рассмотренного в первой части статьи, различные алгоритмы обнаружения разладки. Особенность алгоритмов состоит в том, что они представляют собой модификации адаптивных алгоритмов временной дискретизации и решают задачу обнаружения изменения свойств сигналов в условиях априорной неопределенности.
Ключевые слова — разладка, адаптация, измерение, адаптивные алгоритмы, сходимость.
Введение
В некоторых достаточно распространенных случаях по условиям поставленной задачи обнаружение разладки должно осуществляться либо по отдельным вероятностным характеристикам (ВХ), либо по их совокупности. В последнем случае задача существенно усложняется, так как требуется измерять соответствующие ВХ, к которым относятся не только первые два момента, а и, например, спектральные или корреляционные функции сигналов. При этом проблема неадекватности априорной информации, необходимой для решения задачи разладки, только усугубляется. Поэтому для решения поставленной задачи интересно совместить рассмотренные в первой части статьи [1] методы, использующие принцип симметрии в условиях минимума априорной информации о ВХ сигнала, и алгоритмы вероятностной адаптивной дискретизации (ВАД), исследованные в работах [2, 3]. Алгоритмы ВАД основаны на итерационной коррекции временного интервала дискретизации по результатам сравнения на каждой итерации погрешности восстановления сигнала, которая формируется в соответствии с интерполяционной формулой Лагранжа, с ее заданной величиной. Погрешность допускается априорно задавать как в абсолютном, так и в относительном виде. Однако для того, чтобы использовать рассмотренный в предыдущих ча-
стях статьи [1, 4] метод обнаружения разладки в алгоритмах ВАД, необходимо внести в эти алгоритмы соответствующие изменения. Для этого предварительно рассмотрим аналитически заданные в обобщенном виде алгоритмы ВАД.
Описание и условие сходимости адаптивного алгоритма временной дискретизации
Пусть исходный аналоговый сигнал при преобразовании в цифровой код дискретизируется с постоянным временным шагом дискретизации Дг. Эта операция выполняется во всех информационно-измерительных системах посредством аналого-цифровых преобразователей. В дальнейшем для упрощения записи будем предполагать, что Дг = 1. Известный алгоритм ВАД для поиска оптимального интервала дискретизации при восстановлении сигнала полиномом т-й степени на
т
временном сегменте Тс = ^ ту в соответствии со
1=0
статьей [2] представляется в виде
тс(к + 1) = тс(к) -
- Д(£Ме{5[фт[г, тт У(г)]}, 5о}, (1)
где Тс(й + 1) и Тс(к) — значения искомого временного сегмента на (й + 1)-м и й-м шаге (такте) итерации, й = 0, 1, 2, Д(й) — некоторая (фокусиру-
ющая) последовательность, влияющая на изменение значения искомого временного сегмента на (к + 1)-м шаге итерации в зависимости от его значения на к-м шаге итерации; 5[фт[...], у(г)] — текущая погрешность восстановления исходного сигнала у(£) функцией восстановления; фт[...] — функция восстановления, представляющая собой интерполяционный полином т-й степени, построенный по значениям сигнала в точках отсчеТ (к)
та ту (к) = с ^ ], ] = 1, 2, 3, ... т + 1; т — индекс,
значение которого определяется порядком т интерполирующего полинома; е{...} — функция меры, оценивающая вид погрешности восстановления сигнала у(£) на временном сегменте Тс; ц{е{...}, 5о} — преобразование, характеризующее изменение отклонения значения функции меры текущей погрешности восстановления сигнала у(£) полиномом т-й степени от ее заданной величины
50, геТс.
Отметим, что все временные переменные, входящие в алгоритм (1), и вытекающие из него алгоритмы, рассмотренные ниже по тексту, также кратны шагу дискретизации Дг = 1, а обозначения г и Тс сохранены для сокращенного обозначения ввиду несущественного влияния их выражения через Дг на результаты исследования.
Сходимость алгоритма (1) к аттрактору обеспечивается тем, что погрешность восстановления сигнала возрастает от нуля с ростом интервала дискретизации до некоторого максимального значения, причем на начальном участке погрешность растет монотонно. Поэтому в пределах достаточно малой погрешности восстановления (по крайней мере, до 10-20 %) существует такое значение интервала дискретизации, при котором в статистическом смысле устанавливается равенство действительной (текущей) погрешности восстановления заданной величине. По существу в этом случае речь идет о равенстве функции регрессии, которая описывает в среднем изменение погрешности восстановления в зависимости от изменения временного сегмента 0 < Тс < Тс тах, заданной величине 80. Если же построить алгоритм таким образом, что в качестве заданной погрешности восстановления использовать также текущую погрешность восстановления, но изменяющуюся противоположно, т. е. от максимального значения до нуля, то будет устанавливаться по-прежнему некоторое, но, возможно, уже другое устойчивое значение временного сегмента Тс. При этом значении временного сегмента образуется новый аттрактор, при котором обеспечивается в статистическом смысле равенство монотонно возрастающей и монотонно убывающей погрешностей восстановления сигнала (рис. 1). Заметим,
■ Рис. 1. Условия сохранения постоянства параметра Тс0 для различных по виду автокорреляционных функций сигнала
что применять интерполирующий полином выше 2-3-й степени для решения задачи разладки нецелесообразно.
Обобщенный алгоритм разладки на основе ВАД
Для обнаружения разладки алгоритм (1) в соответствии с обобщенным алгоритмом, предложенным в первой части статьи [1], преобразуем к виду
тс(к + 1) = тс(к) - Д(£Ме1[5(фт(г - Тт, Тт - ЗД),
у(г - Тт)], е2[5(фт(*, Тс(к)), у(*))]}, (2)
где ц{...} = ЦТА; чадо..), у(г))], е2[5(.), у(г))]} -
преобразование, которое устанавливается с учетом свойств функции меры или функции качества [2]; е^...] и е2[.] — преобразования со свойствами функции меры, описывающей зависимость величины отклонения сигнала у(£) и его задержанного или запомненного значения у(г - Тт) от указанных функций восстановления, т. е. определяющих погрешность восстановления сигнала; тт — временной параметр, или время памяти, на который задерживается (запоминается) относительно текущего момента времени входной сигнал у (г).
Для дальнейшей конкретизации алгоритма (1) целесообразно представить преобразование Ч{...} в виде
Ч{...} = Ч^^фтО - Тт, Тт - ЭД), уЦ - Т„))] -
- е2[5(фт(г, Тс(к)), у(г))]},
где в целях обеспечения симметрии рекомендуется брать е^...] = е2[.] = е[...].
В силу свойств симметрии и монотонного изменения усредненной погрешности восстановле-
ния в интересующих нас пределах ее изменения алгоритм (2) сходится к временному сегменту тс0, если выполняется условие
МуМе[5(фт(* - Тm, Тт - Tc0), у(* - Тт))] -
- е[5(фт(г, Тс0), у(*))]}} = 0, (3)
где Му{...} — оператор определения математического ожидания по у(г).
В случае нарушения симметрии между ВХ сигнала в прошлом и настоящем равенство (3) не выполняется, и тем самым фиксируется разладка в темпе с ее наступлением. Нарушение симметрии выражается в том, что соответствующие преобразования, выполненные для запомненных и текущих, т. е. поступающих в реальном масштабе времени, значений сигнала, представленного моделью квазистационарного случайного эргодического процесса, в среднем не совпадают. Причем стационарность сигнала и равенство (3), которые нарушаются в момент наступления разладки, снова восстанавливаются на последующем временном секторе (фрагменте) сигнала.
Пример алгоритма разладки
Представим алгоритм (1) для восстанавливающего полинома нулевой степени и постоянного шага итерации Д(п) = аг при п = 0, 1, 2,. и ц(х) = х в виде
Тс(п + 1) = Тс(п) + а* Шп) - £[п - (Тт - Тс(п))]]2 -
- К[п - Тс(п)] - £(п) ]2}, (4)
причем Тс(0) = 0.
Тогда в соответствии с (3) при отсутствии разладки после начального переходного процесса получаем
11Ш М% {к(п -Тт ) -1[п - {Тт - Тс (п))]12 -
-[| [п -Тс (п)]- ^(п)]2 } =
= 2а2 [г (Тс0)- Г (Тт -Тс0 )] = 0, (5)
где а2 и г(...) — соответственно дисперсия и нормированная автокорреляционная функция (НАКФ) сигнала.
Из уравнения (5) для стационарного фрагмента сигнала, т. е. при отсутствии разладки, следует равенство
Тс0 = Тт/2, (6)
которое и определяет в среднем аттрактор или неподвижную точку для алгоритма (4). Обратим внимание на то, что равенство (6) не зависит от НАКФ для стационарного фрагмента сигнала. При этом параметр Тт выбирается из условия 0 < Тт < тк,
где тк — интервал корреляции НАКФ. Важным является вопрос, связанный с выявлением представительности класса сигналов, для которого алгоритм (2) будет работоспособен при заданной величине Тт. В пределах допустимого класса сигналов для интервалов корреляции НАКФ должно выполняться условие 0 < Тт < тк т1п. На рис. 1 качественно показано, при каких условиях сохраняется равенство (6) для различных по виду НАКФ из заданного класса сигналов.
Оценка переходного процесса при наступлении разладки
Предположим, что, начиная с некоторого п0, скачком изменилась НАКФ, тогда после переходного периода, как это следует из условия (3), равенство (6) снова восстанавливается, но уже для другой НАКФ. Оценку переходного процесса аналитически можно осуществить при таком значении параметра Тт, когда НАКФ с погрешностью второго порядка малости можно заменить линейными членами ряда Тейлора относительно стационарного значения (6). Тогда, после усреднения алгоритма (4) и вычитания в нем справа и слева временного сегмента Тс0, получим
йс (п +1) = йс (п) + а {2а2 (1- г (Тс0)- г'(Тс0 )йс (п))-
- 2а2 (- Г1 (Тс0 ) - Г1 (Тс0 )йс (n))}, (7)
где йс (п +1) = Тс0 -Тс (п + 1) и йс (п) = Тс0 -Тс (п);
г'(Тс0) и г{(Тс0) — первые производные НАКФ до и в момент разладки для значений временного сегмента Тс0-
Из (7)следует
1-(1-2а^ у)
йс (п) = йс (0)(1-2а*У) -
2аД
п 1
1-(1-2а, у)
1-2а,у)п -—, (8)
у
где и(0) — начальное отклонение в момент наступления разладки, и(0) = 0, так как Тс (0) = Тс0 = Тт /2;
У = а1г1 (Тс0,
1+
а2г' (Тс0 )'
а1г1 (Тс0 ).
— = 1а2 (1-г (Тс0))-а2 (1-Г1 (Тс0)
Параметры у в первом приближении и — во втором приближении при разложении НАКФ и ее производной в ряд Тейлора с учетом г^(0)= 0 можно упростить, поэтому
и
Ь"(о)|т„ ( м)
У »—L^-J (р + 1)
а? |г1//( 0)Тг2 7 — » МП Л т (р-1),
где г"(0) и гД0) — вторые производные НАКФ в нуле до и в момент разладки для значений времен-
а2г"(0)
ного сегмента Тс0; р = — коэффициент, ха-
а12г1"(0)
рактеризующий разладку по НАКФ сигнала.
Представим отношение — / у приближенным равенством
— / у = р * 0,25Тт[(р - 1) / (р + 1)].
Это равенство определяет чувствительность алгоритма к разладке, поэтому (8) преобразуется к виду
йс (п)»(1 - 2а,у)п - р. (9)
Минимальная чувствительность оценивается в виде
р * 0,25тк т1д[(Р - 1) / (р + 1)]-
Из (9) следует, что динамика изменения временного сегмента Тс(п) с ростом п относительно его стационарного значения приближенно описывается формулой
Tc (n)» Tc0 +(1 - 2a,y)n - pf
0,25Tm
2 - р-1
P + 1
-2atyn. (10)
Приближенный анализ позволяет получить достаточно полную информацию о характере переходного процесса при изменении НАКФ сигнала в момент наступления разладки. Действительно, в соответствии с анализом (9) и (10) можно сделать следующие выводы:
— при переходе от одного фрагмента сигнала к другому искомый временной сегмент Тс (п) в среднем изменяется скачком, величина которого соответствует р и зависит, таким образом, от изменения ВХ сигнала;
— величина скачка р с ростом числа итераций в течение 0 < п < Тт меняется незначительно, поэтому текущее значение временного сегмента Тс(п) в момент разладки изменяется скачком и далее по экспоненте стремится к прежнему стационарному значению независимо от вида НАКФ;
— легко установить, что воздействие стационарной аддитивной помехи не приводит к смещению оценки временного сегмента Тс (п), так как не происходит нарушение условия симметрии.
и
Очевидно, что для получения более полных результатов исследования свойств алгоритма необходимо провести анализ влияния дисперсии параметра Тс(п) относительно Тс0 на ошибку обнаружения разладки. Ввиду громоздкости подобного анализа ограничимся в рамках данной статьи имитационным моделированием рассмотренного алгоритма в модификации (4), результаты которого представлены на рис. 2, а—в. Имитационное моде-
а)
4,5060 4,5050 4,5040 4,5030 4,5020 4,5010 4,5000 4,4990 4,4980 4,4970 4,4960 4,4950 4,4940
0 10 20 30 40 50 60 70 80 90 100
б)
4.506
4.504
4.502
4.500
4.498
4.496 4,494 4,492 4,490 4,488
0 10 20 30 40 50 60 70 80 90 100
в)
2.506
2.505 2,504
2.503 2,502
2.501 2,500
2.499 2,498
2.497
0 10 20 30 40 50 60 70 80 90 100
■ Рис. 2. Результаты имитационного моделирования алгоритма (4) с квадратичной функцией меры по обнаружению разладки по частоте следования для двух гармонических сигналов (для графика 1 <»! = 0,0101; <»2 = 0,01015; для графика 2 <»! = 0,01; <»2 = 0,01005): а — at = 5; б, в — at = 10
-Tm = 5,01-
[
1
S'
J m u
4 Tm = 9,01-
X -.2
1 m ^
лирование проводилось для двух гармонических сигналов, у каждого из которых скачком меняется частота следования от шг до + А<ю при i = 1, 2. При этом алгоритм обнаружения разладки (4) для гармонического сигнала имеет следующий вид:
Tc(n + 1) = Tc(n) + ott((sin(roi(n - Tm + Tc(n))) -- sin(roin))2 - (sin(roi(n - Tc(n))) - sin(roin))2), i = 1, 2.
Графики переходных процессов алгоритмов в зависимости от изменения значений исследуемых параметров (см. рис. 2) по виду совпадают с полученной аналитической зависимостью (10). В целом результаты исследования показывают, что предложенный алгоритм обнаружения разладки обладает высокой чувствительностью к изменению параметров НАКФ сигнала. Так, изменение частоты следования сигнала менее чем на 1 % приводит к легко обнаруживаемому скачку искомого параметра относительно его стационарного значения. Также дисперсия искомого параметра существенно меньше величины изменения искомого параметра в зависимости от изменения частоты гармонического сигнала.
Модифицированные алгоритмы разладки
Модификацию алгоритма (2) для равномерной меры приближения сигнала некоторой интерполяционной функцией ф(...) представим в виде
Tc (n +1) = Tc (n)-
- at {maxiy(t Tm)-фт (t-Tm, Tm-Tc (n))|-
- max|y(t) Фт (t. Tc (n))|
(11)
и применим к его исследованию результаты, полученные в работе [3] для интерполяционных полиномов любого порядка.
Для того чтобы обойти трудности анализа сходимости в среднем алгоритма (11), применим методику исследования, предложенную в работе [3], с учетом следующих известных результатов: — максимальную в пределах временного сегмента Тс погрешность интерполяции 5тах(у(г), фт(г, Тс)) полиномом т-й степени в форме Ньютона с равномерной мерой приближения представим в виде
5шах (^/(t), фт (t, т0т ))
г+1)
И1)!
max
0 £ a £ m
m+1
где |Ym(a)| = (Xm)m + 1 |[(a) (a - 1) (a - 2) ... (a - m)] =
= (xm)m + 1 a = t/xm и Tm = Tm/(m + 1);
— математическое ожидание абсолютного значения случайной величины |||, описываемой достаточно широким классом плотностей распределения вероятностей (ПРВ), можно представить следующим образом:
М|§| = Ка^,
где К — постоянная, зависящая от вида ПРВ случайной величины |; а^ — среднеквадратическое отклонение случайной величины
Для алгоритма (11) в соответствии с условием (3) получаем для определения стационарного значения временного сегмента Тс0 равенство
\(т+1)
T -T
c0/
m + 1
[(т + 1)!]
-к
[c0/
/m + 1
г+1)
[(m + 1)!]‘
0 m
%0m
Ry (m+1)](0)
R
2 (m+1)]
(0)
(12)
которое соответствует (6).
Очевидно, что при использовании полиномов нулевой степени (нулевого порядка) в алгоритмах (4) и (11) при наступлении разладки изменение временного сегмента Тс(п) будет в основном происходить только при изменении ВХ первой производной сигнала и не зависеть от постоянной составляющей сигнала. Для интерполяционного полинома первого порядка аналогичная реакция временного сегмента Тс в соответствии с алгоритмом (11) будет осуществляться уже при изменении ВХ второй производной сигнала и не зависеть от постоянной составляющей сигнала и его скорости изменения и т. д. Это утверждение непосредственно вытекает из (11) и равенства (12). Отметим, что алгоритм (11) только демонстрирует принципиальную возможность применения интерполяционных полиномов п-го порядка. Конкретный вид соответствующих алгоритмов, которые можно было бы использовать в цифровых системах, для установленного порядка интерполяционного полинома должен уточняться. Например, алгоритм, синтезированный на основе полинома первого порядка, с использованием конечных разностей представляется следующим образом:
Тс(п + 1) = Тс(п) - о*{|0,5[у(п + Тт - Тс(п))] -- у[п + 0,5 (Тт - Тс(п))]|} - о*{|0,5[у(п) +
+ у(п + Тс(п))] - у[п + 0,5 Тс(п)]|}, п = 1, 2, ....
Для достижения инвариантности параметра аг к изменению ВХ сигнала целесообразно использовать методы, предложенные еще Ньютоном и исследованные применительно к аналогичным алгоритмам в работе [5].
У
Алгоритм (11) для обнаружения разладки можно усложнить, если контролировать в нем в любых сочетаниях линейную комбинацию погрешностей восстановления для различных интерполяционных полиномов от 0-го до т-го порядка. Подобный алгоритм, например, имеет вид
Тс(п + 1) = Тс(п)-
^ 5 шах 1=0
г-Т , Тс(п) Тс(п)
+ а^|^ 5 шах 1;=0
I + 1 I + 1
Тс (п)
I + 1
где 5тах(...) — максимальная погрешность восстановления сигнала у(г) соответствующей функцией восстановления фг(...).
Инвариантность равенства (6) в случае отсутствия разладки для аттрактора сохраняется.
В рамках рассмотренного подхода разберем пример синтеза еще одной модификации алгоритма [6], которая позволяет решить задачу обнаружения разладки по изменению функции распределения вероятностей (ФРВ) или ПРВ сигнала. Алгоритм разладки первого порядка по ФРВ имеет вид
С(п + 1) = С(п) - а^^пС/М - у(п)] -
- sign[d•(n) - у(п - Тт)]},
(13)
где С(п + 1) и Су(п) — значения искомого параметра на (п + 1)-м и п-м шаге итерации фиксированной размерности, причем для всех п = 0, 1, 2, . и / = 1, 2, ..., р выполняется условие \С(п)\ < ^0 и устанавливаются начальные значения параметра для всех / в виде /0) = /(Е0/р), / = 1, 2, ..., р; р — число ординат ФРВ, по которым контролируется момент наступления разладки; Е0) — допустимый диапазон изменения искомого параметра; у(п) и у(п - Тт) — текущее и задержанное или запомненное значение сигнала; sign[...] — знаковая функция, т. е.
81£П |Х| =
1при х > 0;
- 1при х <0.
Искомый параметр С(п + 1), как уже отмечалось, представляет собой некоторую статистическую функцию, по которой принимается решение о возникновении разладки в сигнале у(п) в момент времени п0. Напомним, что решение о возникновении разладки принимается, если для указанной функции выполняется условие С(п0 + 1) > С/п, где С/п — некоторое ее заданное /-е пороговое значение, или просто •-й порог. Величина этого порога в классическом варианте определяется на основа-
нии априорных данных о виде ФРВ решающей функции, включая диапазон изменения сигнала Е0. В алгоритме (13) временной параметр Тт целесообразно выбирать из условия Тт > тк т^.
Сходимость алгоритма (13) уже доказывалась различными методами [7, 8], поэтому остановимся только на тех моментах, которые важны для рассматриваемой в данной работе задачи. Усредним правую и левую части уравнения (13), в результате получим
dj (п +1) = dj (п) -
- 2^ {) [{ (п)] - Р/(п-Тт )№ (п)]}, (14)
где Ру(п)Щп)\ и ^у(п - Тт)[Су(п)] — ФРВ сигнала в моменты времени п и п - Тт относительно усредненных значений аргументов.
Очевидно, что если ФРВ сигнала не изменяется, то Е [/(п)] совпадает с Еу(п - Тт)С/(п)] и параметр Су(п) флуктуирует относительно начальных значений С/(0) для всех /. Как только наступает изменение ФРВ, так сейчас же происходит изменение параметра Су(п), пока измененная ФРВ будет сравниваться с предыдущими значениями ФРВ на интервале Тт.
С учетом вариации ФРВ относительно первоначального ее вида в момент времени п0 в соответствии с (14) приходим к следующей формуле, характеризующей динамику отклонения усредненного искомого параметра относительно его первоначального значения:
С/(п0 + т) = С/(п0)[1 - а/Е (С/]т,
где 5Е(С/ — вариация ФРВ для указанного аргумента.
Здесь, как и прежде, возникает задача выбора порога для проверки нулевой гипотезы против установленной альтернативы. Величину порога можно выбрать, исходя из минимально допустимой вариации ФРВ.
Заключение
Основным достоинством предложенного и исследованного алгоритма является его универсальность, так как он позволяет на основе использования принципа симметрии и алгоритмов вероятностной адаптивной дискретизации решить проблему обнаружения разладки при минимуме априорной информации по ВХ наблюдаемого случайного квазистационарного сигнала. Общность представления алгоритма позволяет его использовать для решения широкого круга практических задач, например, для оперативного контроля изменения расхода жидких или газообразных продуктов при их транспортировке в трубопроводах.
Литература
1. Тихонов Э. П. Разладка, гомеостазис, измерение в рамках компенсационного принципа равновесия в динамических системах. Ч. 1: Общий анализ // Информационно-управляющие системы. 2012. № 2. С.82-88.
2. Тихонов Э. П. Вероятностные адаптивные алгоритмы дискретного представления аналоговых сигналов. Ч. 1: Исследование свойств // Информационно-управляющие системы. 2011. № 2(51). С. 8-15.
3. Тихонов Э. П. Вероятностные адаптивные алгоритмы дискретного представления аналоговых сигналов. Ч. 2: Сравнительный анализ и численные данные // Информационно-управляющие системы. 2011. № 3(52). С. 9-14.
4. Тихонов Э. П. Разладка, гомеостазис, измерение в рамках компенсационного принципа равновесия в динамических системах. Ч. 2: Исследование конкретного алгоритма // Информационно-управляю-щие системы. 2012. № 3. С. 77-83.
5. Тихонов Э. П., Федотенкова М. Л. Адаптивные алгоритмы для измерения признаков биомедицинских показателей при автоматической медицинской диагностике по электрофизиологическим по-
казателям // Биомедицинская радиоэлектроника. 2010. № 4. С. 74-77.
6. А.с. № 235412, Кл. 42т4, 7/52; МПК G 06g. Устройство для измерения функции распределения случайных сигналов / Э. П. Тихонов. № 1203876/18-24; Заявл. 18.12.67; Опубл. 16.01.69; Бюл. № 5. — 2 с.
7. Тихонов Э. П. Метод определения оценки одномерной функции распределения // Методы представления и аппаратурный анализ случайных процессов и полей: тез. докл. III Всесоюз. симп., III секц. «Синтез оптимальных операторов оценки, алгоритмы и структурные схемы измерений», Сухуми, 15-18 декабря 1970 г. / Всесоюз. НИИ электроизмерительных приборов. Л., 1970. С. 48-53.
8. Тихонов Э. П. Анализ эффективности метода оценивания с использованием опорных случайных процессов // Перспективные методы планирования и анализа экспериментов при исследовании случайных полей и процессов: тез. докл. Всесоюз. конф., Нальчик, 1-3 ноября 1982 г. М., 1982. Ч. I. С. 29-31.
УВАЖАЕМЫЕ АВТОРЫ!
Национальная электронная библиотека (НЭБ) продолжает работу по реализации проекта SCIENCE INDEX. После того как Вы зарегистрируетесь на сайте НЭБ (http://elibrary.ru/ defaultx.asp), будет создана Ваша личная страничка, содержание которой составят не только Ваши персональные данные, но и перечень всех Ваших печатных трудов, имеющихся в базе данных НЭБ, включая диссертации, патенты и тезисы к конференциям, а также сравнительные индексы цитирования: PИНЦ российский индекс научного цитирования), h (индекс Xирша) от Web of Science и h от Scopus. После создания базового варианта Вашей персональной страницы Вы получите код доступа, который позволит Вам редактировать информацию, в том числе добавлять публикации, которых нет в базе данных НЭБ, помогая создавать максимально объективную картину Вашей научной активности и цитирования Ваших трудов.
№ Б, 2Q12
79
