Научная статья на тему 'Измерение плотности вероятности эргодического случайного процесса с линейной корреляционной функцией'

Измерение плотности вероятности эргодического случайного процесса с линейной корреляционной функцией Текст научной статьи по специальности «Математика»

CC BY
240
38
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭРГОДИЧЕСКИЙ СЛУЧАЙНЫЙ ПРОЦЕСС / ЛИНЕЙНАЯ КОРРЕЛЯЦИОННАЯ ФУНКЦИЯ / АЛГОРИТМЫ ИЗМЕРЕНИЯ / ЭКСТРАПОЛЯЦИЯ И ИНТЕРПОЛЯЦИЯ / ERGODIC RANDOM PROCESS / LINEAR CORRELATION FUNCTION / MEASUREMENT ALGORITHMS / EXTRAPOLATION AND INTERPOLATION

Аннотация научной статьи по математике, автор научной работы — Заико Александр Иванович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Measuring probability density of ergodic random having linear correlation function

The paper proposes algorithms and error estimations for measuring probability density of random process having linearly decreasing correlation function, considering reconstruction of its trajectory between samples.

Текст научной работы на тему «Измерение плотности вероятности эргодического случайного процесса с линейной корреляционной функцией»

ISSN 1992-6502 (Print)_

2017. Т. 21, № 4 (78). С. 121-128

Вестник УГАТУ

ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru

УДК 519.216

Измерение плотности вероятности эргодического случайного процесса

с линейной корреляционной функцией

А. И. Заико

zaiko@ugatu.ac.ru

ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)

Поступила в редакцию 22.09.2016

Аннотация. Приведены алгоритмы и погрешности измерения плотности вероятности случайного процесса с равномерным законом распределения плотности вероятности

Ч/ /■» Ч/ Ч/ | ч/

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

Ключевые слова: эргодический случайный процесс; линейная корреляционная функция; алгоритмы измерения; экстраполяция и интерполяция.

ВВЕДЕНИЕ

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

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

В статье рассматривается процедура цифрового измерения одномерной плотности вероятности на примере оригинального случайного процесса с равномерным законом распределения и линейно убывающей корреляционной функцией при равномерных распределениях погрешности квантования по уровню и шагах дискретизации во времени [2, 3]. Характеристики погрешностей измерений плотностей вероятностей получены с применением комплексного подхода к определению погрешностей. Он позволяет избежать некорректного суммирования погрешностей дискретизации и квантования, а также погрешностей восстановления реализации процесса между отсчетами, и учесть их взаимное влияние друг на друга [4, 5].

При цифровых измерениях реализация случайного процесса равномерно дис-

кретизируются во времени с шагом Т0 и квантуется по уровню с шириной кванта 2Ак. Получаются дискретные отсчеты ха, где г - номер отсчета, датируемого моментом времени ti, а 1 - номер кванта, соответствующий уровню квантования х1. Номера

отсчетов

принимают

значения

г = -п, ..., -1, 0,1, ..., п, где 2п +1 - количество отсчетов, а 2пт0 - длительность реализации.

Количество уровней квантования обозначим через Ь, а номера уровней квантования 1 = 1, 2, ..., Ь [6].

Оценка X]) одномерной плотности

вероятности при интерполяции реализации процесса между отсчетами имеет вид [1, 6]

1 П—1Т0

ЫХ])и = — 'Х-пк, ..,ха,.., ХпГЩ, (1)

0 г=-п ог

где м>1[х|X; Х-пк, ..,ха,.., хпг] - одномерная

условная плотность вероятности процесса в момент времени 0 < X < Т0 между соседними

отсчетами после получения всех отсчетов х-пк, ..,х,1,.., хпг реализации; к, г = 1, 2, ...,Ь -номера уровней квантования.

При экстраполяции в будущее реализация х(() восстанавливается на г интервале дискретизации по предшествующим отсчетам х-пк,...,хй и продляется еще на один

2п+1 интервал tn < t < tn + Т0. Тогда оценка (н'1[Х]) одномерной плотности вероятности при экстраполяции

МХ ]) э =

(2п + 1)Т '

£ I *1[ X X;

, ]<&. (2)

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

м^1[Х| X; х-пк, ..., хй ] при экстраполяции

приобретает вид щ[ х| X; ха ]. Обозначим

в выражении (2) частость появления отсчетов хп, равных уровню квантования х1 и

определяющих плотность распределения вероятности м>1[Х\X; ха ] = м>1[Х|X; х1 ], через

(р(х1)) = п1!(2п +1) , где п1 - количество

отсчетов хй, равных уровню квантования х1

и )) = 1. Тогда оценка (2)

1 0

ИХ])э = -{Р(х1)) | ^[Х^; х1 ^ (3)

Т г

При интерполяции реализации х^) процесса восстановление происходит по двум смежным отсчетам хи, х^^к в выражении (1)

одномерная условная плотность вероятности ^1[Х| X; х+пк, ..х-- хпг] = ^1[Х| ^ ха, х(( I где 0 < X < т0; 1,к = 1,2,...,Ь.

Обозначим частость появления отсчетов, следующих друг за другом и соответствующих уровням квантования х1 и хк,

определяющих плотность вероятности мДХ| X; х1, хк ] через (р(х, хк^ = %/2п, где п1к - количество событий, заключающихся в появлении отсчетов, с уровнями квантова-

Ь Ь

ния х и хк , а р(х,хк)) = 1 Тогда

1=1 к=1

выражение (1) примет вид

1 Т0

(м^ = -(Р(х1, хк)}|Х|X; х1, хк(4)

Из выражений (1)-(4) следует, что оценка одномерной плотности вероятности существенно зависит от одномерной условной плотности вероятности мх [ X [X; х- пк,.., хг1,., хпг ],

которая определяется свойствами и характеристиками конкретной модели случайного процесса.

Рассмотрим случайный процесс Заико с равномерным законом распределения плотности вероятности м^^Х] [2, 3]. Такая модель проста, требует минимума априорной информации и позволяет получить пригодные для инженерной практики результаты. Она описывается всего тремя параметрами: нижней X и верхней X границами изме-

1

нения случайного процесса и нормированной корреляционной функцией р(Т), которую в нашем случае положим равной

р(т) =

1 -|т|/т0, 0 <|т|< т0; 0, | т|> То,

где т - временной сдвиг; т0 - интервал корреляции.

Для него плотность вероятности

щ^Х ] =

1(Хв - Хн )-1 = (2АкЬГ, Хн < X < Хв; I 0, в остальных случаях.

Погрешность квантования по уровню случайна, стационарна и независима для отсчетов и от процесса. Опишем ее для / = 1, 2,...,Ь равномерной плотностью веро-

ятности отсчета х

^[Х|х ]=|(2Ак )-1, ^ -Ак < Х < ^ +Ак;

[0, в остальных случаях.

При экстраполяции в будущее реализацию х(^) процесса восстанавливаем по предыдущему отсчету х;. При равномерных распределениях случайного процесса и погрешности квантования условная плотность вероятности щ [Х|1; х1 ], где 0 < 1 < Т0, также равномерна [1, 6]

щ![ Х\I; Х ] = (5)

= |[Хв (1; х1) - Хн (1; х1 )]-1, Хн (1; х1) < Х < Хв (1; х1); [0, в остальных случаях,

где верхняя хв (1; х,) и нижняя Хн (1; х1) границы динамического диапазона изменения случайного процесса при экстраполяции и 0< 1 <т0 [1, 6]:

Х в(1; х1)=х1 + А к +(Х в- х1- Ак) ут 0;

Хн (1; х1 ) = х1 - Ак - (х/ - Ак - Хн )УТ 0 .

Тогда выражение (5) примет вид

щ [ Х| 1; х1 ] =

1 + (Ь -1)

2А. т„

- X

К (6)

1 х1 - Ак -х - Ак - Хн) Vт0 < Х <

< х1 +А к +(Х в - х1 -А к )

1/Т0 ,0 < 1 < Т0;

0, в остальных случаях.

В зависимости от соотношения шага дискретизации Т0 и интервала корреляции т возможны два варианта экстраполяции реализации случайного процесса и получения оценки [Х]) одномерной плотности

вероятности: т0 < т0 и т0 > т0.

При Т0 < т0 экстраполяция по последнему отсчету х осуществляется с шагом то, и интервал корреляции Т накладывается на соседний шаг дискретизации. Поэтому оценка Щ[Х^ (3) для / = 1, 2,...,Ь примет вид

ИХ ]) э =

№ )>.

X <

1п-

э 2АкТ0 Ь -1

1 + (Ь -1)70/т 0 1 + (Ь -1) х/ -А к - Х '

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

х/ -А к - Х н

х/ - Ак -(х/ - Ак - Хн )Т^Т 0 < Х < х/ - Ак; 1п [1 + (ь - 1)Т„ /Т 0 ], х/ - А к < Х < х/ + А к; 1 + (Ь -1)7,/Т 0

1п

1 + (Ь -' Хв -х/ -Ак

х/ + А к < Х < х/ + А к +(Х в - х/-А к )Т„/ т 0;

0, в остальных случаях.

При т > т0 экстраполяция по отсчету %й осуществляется на интервале 0 < 1 < т0, и щ1[Х\ 1; х/ ] описывается выражением (6). При т 0 < 1 < Т0 з ави сим о сть от п о сл ед н его отсчета хг7 отсутствует и щх[Х| 1; х1 ] = щ1[Х ] (6) [6]. В результате при 7 > т0 и / = 1, 2,...,Ь

Щ[Х])э =- "

'э 2А к Т0

[(Р(х/ )>

X <

1п-

Ь

_ | Т 0 - Т 0 1

Ь - 1 1 + (Ь - Т0 ь '

У 'х1 -А к - Хн

<р(х/» щь + I,

Т 0 ь

Ь

Хн <Х<х/-Ак;

х/ - Ак < Х < х, + А

/ + А к ;

Ь - 1

(р(х1 ))ы Ь + Т0 - Т 0 1

Ь - 1 1 + (Ь - 1) Х - х/ -А к Т0 V

У Х в - х/ -А к

хг +А к < Х < Хв;

0, в остальных случаях.

Т

0

X

Абсолютная погрешность Дэ[х ] оценки (^1 [X]) при экстраполяции и комплексном подходе к ее определению для Т0 < х0 и , = 1, 2,...,Ь [2, 6]

При интерполяции реализацию х(?) случайного процесса восстановим по двум соседним отсчетам х^ и х/1+1лк. Тогда условная

плотность вероятности ^1[Х| X; х1, хк ], где 0 < X < Т0, равна [5-7]:

Дэ[X] = <X])э -X] = х

т0 (р(x¡\_1 + (Ь - 1)То/т0 1

х

Т0 Ь -1

1 + (Ь -1)

XI-Д к - X Ь

х1 -Д к - X н

с, - Дк - (х; - Дк - Xн )То/т0 < X < х - Дк; 1п [1 + (Ь - 1)Т0/ т 0 ]-±

Т Ь -1 0 0 Ь

х1-Д к < X < х1 +Д к;

{р(х1 ы 1 + (Ь - 1)Т„/т

Т0 Ь -1

0/ ^0

1 + (Ь -1)

X - X; -Дк Ь'

X в - х1 -Д к

х+

Д к < X < х1 + Д к + (X в - х1 -Д к )Т„/ т 0;

0, в остальных случаях.

При Т0 > т0 и , = 1, 2,...,Ь [3, 4]

Д э [X ] = -?

2Д кТ0

(Р(х;))

Ь -1

(рЫ)

Ь -1

<Р(х,))

1п-

1 + (Ь -1)

--,Xн <X < X, -Дк

1п Ь -

Ь

1п

X, -Д к - X Ь

X, -Дк - Xн

X, -Дк <X <X, +Дк;

Ix-xr-^; - Ь - X +Дк < X < Xв;

ь -1 1+(ь -1)

^ - X, -Д 0, в остальных случаях.

Отметим наличие пропорциональной {Р{^с1)}/2ДК мультипликативной погрешности, которая имеет разное значение при Т < т0 и при Т0 > т0. Кроме того, при Т0 > т0 появилась аддитивная погрешность (Т0 - т0 )/2ДкЬТ0 , которая при Т0 < Т0 отсут-

ствует.

[ Xх; X, ^] =

^ в(Х; X, ч)- X н (X; X, xk)]

0, в остальных случаях,

Xн (X; X,, xk) < X <

, < Xв (Х; X, ^);

где разность верхней X (X; щ, ^) и нижней (X; X, %) границ динамического диапазона изменения случайного процесса при интерполяции [5-7]:

X в X, xk)-X н (X; xlX )= = 2Д

Ь -(Ь -1) p(x)+ рТ - X>' 1 + Р(Т0)

Возможны три варианта восстановления реализации: Т0 <т0, т0 < Т0 < 2т0 и Т0 > 2т0. Результат восстановления на интервале 0 < X < Т зависит также от взаимосвязи от-

счетов X;; и X

которая учитывается

условной

X; X, xк].

плотностью

вероятности

При Т0 <т0 и ,,к = 1, 2,...,Ь

Н^ X ^ xг, xk ] =

7 -Д к +(Xk - X,) X/Т < X < < X +Д К +(xk - X ) VTо, ,

2Д,

0 < X < Т

0, в остальных случаях.

ПРи Т0 <Т0 < 2т0 ,,к = 1, 2,...,Ь

XIX; X,, Xk ] = (2Дк )-1 х

1 + (Ь -1)-

-1 Xl -Дк -(X, -ДК - X, ^ < X <

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

, < X, + ДК +(Хв " X, - ДК

1 + (Ь - 1)Т>-С0

1 + (Ь - 1) ^ Тп

X - Дк - (xk - Дк - Xн )(Т0 - т0 Vт0 + + (Xk - X,) Vто < X < X, +Дк +

Ч -Ч! x т0 < < X ' дк , + (Хв - ^ - Дк )(Т0 - Т0 VТ0 + (xk - X ^/т0

1 Xk -Дк-(xk -Дк - Xн )(Т) - X)/ Т0 < X <

< Xk + Дк +(X в - Xk -Дк )(Т) - XV т,

0, в остальных случаях.

1

т

0

0

х

Ь

1

х

0

0 < X < Т - т0

-1

х

0

Т0 - т0 < X < т0

т0 < X < Т

При Т > 2Т0 /,к = 1, 2,...,Ь

щ[ Х\ 1; х/, хк ] = (2А К )-1 :

1 + (Ь -1)А

-1 х/ - Ак -(х/ - Ак - Хн )1/т0 < Х < , < х/ + Ак +(Хв - х/ - Ак )Vто,

0 < 1 < Т0; Хн < Х < Хв, Т0 < 1 <Т0 -Т0; -1 хк -АК -(хк - АК -Хн )(Т0 - 1)/Т0 < Х < , < хк +АК +(Хв - хк - АК )(Т - 1)/^, Т0-Т0 <1 <Т0; 0, в остальных случаях.

1 + (ь -1) Т0-1

Тп

Подставив эти значения, получим оценку плотности распределения вероятности

Ш и [8].

При Т0 < т0 и /,к = 1, 2,...,Ь

ЫХ]) и =

2А„

х <

(Р(х/, хк))

Х - х/ -А к , х/ -А к < Х < х/ + А к

х/ +А к < Х < хк -А к

(Р( х1 )), {Р(х1, хк))

х и

хк -Ак - Х, хк -Ак < Х < хк +Ак.

к > /;

х/-А к < Х < х1 +А к, к =/; 'х1 -А к - Х, х/ -А к < Х < х/ +А к ' 2Ак, хк +Ак < Х < х/ -Ак

Х - хк -А к , хк -А к < Х < хк +А к.

к < /;

0, в остальных случаях.

при Т0 < Т0 < 2Т0, /,к = 1, 2,...,Ь и к>1

ЫХ]} и =

{р( х> )>.

и 2АкТ0 Ь -1

1п

1 + (Ь - 1)(Т0 -Т0 VТ0 х, -Ак-(х/-Ак -Хн)

1 + (Ь -1)

к V / к ^ н >

х - А к - Х ' X (Т0 - Т 0 V Т 0 < Х < х/ - А к

х/ - А к - Хн

1п[1 + (Ь - 1)(Т0 - Т0 VТ0], 1 + (Ь - 1)(Т0 - Т0 VТ0

х/ -А к < Х <■

/ + Ак;

1п

1 +(Ь - 1) Х - х/ - А к

х/ +А к < Х < х/ +А к +

+ (Хв -х, -Ак)(Т0 -Т0VТ

Хв - х/ - А к

0, в остальных случаях.

При т0 <Т0 < 2т0 и / = к = 1, 2,...,Ь

щ [Х]) = Т^М X \ ^ J/и т* т

2А КТ)

-А к -(х/ -А к - Хн )

Т - 2Т Т0 2Т 0

_ X (Т - Т 0 V Т 0 < Х < Т0 + (Ь - 1)(Т0 - т 0 ), < х/ + А к + (Хв - х/ - А к ) X

х(то- Т0 VТ 0;

0, в остальн^1х случаях.

При Т0 <Т0 <2Т0, /,к = 1, 2,...,Ь и к<1

щ [Х Л X

\ Н j/И 2АТ Ь -1

1п

/ ч/ Ч, хк-Ак-(хк -Ак -Хн)> 1 + (Ь - 1)(Т0 - Т0^Т0 к х(гк 'кч/т ; у н

Xi70 - Т0)!Т0 < Х <

1 + (Ь - 1) хк - А к - Х

хк - а к - Х н

[0 Т 0

< хк-А к ;

1п [1 + (Ь - 1)(Т0 - т 0 V Т 01 хк - А к < Х < хк + А

лЛ(т ч, хк +А к < Х < хк +А к "

1 + (Ь - 1)(Т0 - Т 0 )/Т 0 + (ХХв - хк -А к ^

1П1 + (Ь - 1) Х - хк - А К ' 4 7ХВ - хк - А к

0, в остальных случаях.

(Т0- т 0 VТ 0;

При Т0 > 2т0, /,к = 1, 2,...,Ь и к>1

Х ]) и =

(Р( х/ )>

и 2АкТ0 Ь - 1

1п

Ь

1 + (Ь -1) х/ - АК - Х ' ' !х1 - Ак - Хн

Х„ <Х <х, -К;

1п Ь. 1п—

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ь

1 + (Ь - 1) Х - х/ - Ак

х -Ак < Х < х1 + ак ;

х1 +Ак < Х < Хв;

Хв - х/ - ак

0, в остальных случаях.

При Т > 2т0 и / = к = 1, 2,...,Ь

Ы Х ]) и =

1 Т0 - 2т0 1

2А„

Т Ь

Х„ < Х < Х„

0, в остальн^гх случаях.

X

X

X

X

1

к

X

X

X

хк х/

Т

0

X

X

Т

0

При Т0 > 2т0, /,к = 1, 2,...,Ь и к<1 То РХк)) ^

№ ]), = 1п-

2Ак7о Ь -1 Ь

1 + (Ь - 1) Хк - Ак - Х

Хк - А к - хн

Хн <-

< Хк -Ак;

1п Ь 1п

Ь

1 + (ь - 1) Х - Хк - Ак

Хн - Хк - Ак

Хк-Ак < х < Хк + Ак;

Хк +Ак < X < Xв;

0, в остальных случаях.

Абсолютная погрешность а [х] оценки Ы [х] плотности вероятности при интерполяции реализации между отсчетами,

Т0 < т0, /,к = 1, 2,...,Ь и к > /

А и [X ] = ЫX ]) и - X ] =

(X - ,-А к)-1,

Х1-А к < х < Х1+А к;

2Ак--, Х1 +Ак <х < Хк - АК;

2А„

Хк - -Х1

(Р( Х, Хк )>

Хк- -Х1

(Р( Х/ Хк )>

7 Ь

Хк-А к < х < Хк +А к;

0, в остальных случаях. при Т0 <Т0, / = к = 1, 2,...,ь

А и X ] = ^

2А„

"(Р( х,)) - 1/ ь, Х -А к < X < х, + А к; 0. в остальных случаях.

При 7 < Т0, /,к = 1, 2,...ь и к </

А и [X х 2А„

лк /

Х1 - Хк

(Р( Х/ Хк )>

Х1 - Хк

(Р( Х/ Хк )>

Ь

Х1 -А К < X <

< Х1 +А к;

2А к--,

кЬ

Хк +А к < х < Х1-А к;

хк -а к < х <

< Хк +А к ;

0, в остальных случаях.

При т0 <т0 < 2т0, I,к = 1, 2,...,Ь и к >I А и [X ] = — х

и [ ] 2Ак

10 (Р(Х1 ^ 1 + (Ь - 1)(Т, - т0 VТ0 1 7 Ь -1 л.1т Хг - А к - X ь'

10 Ь 1 1 + (Ь -1)

Х/ - А к - X

х1 - А к - х н

к н

Х1 -Ак -(х/ -Ак - хн)(Т0 - Т0 VТ0 < х < Х1 -Ак

1

а СР(Х/))

Т0 Ь -1

1п

1 + (Ь -1) 70 Т 0

Т

0

Х/ -А к < х < Х/ +А к

Т± (р(Х/ ))]п 1 + (Ь - 1)(Т0 - Т0 VТ0 1

Т0 Ь -1 1 + (Ь -1) х - Х/ - А к Ь '

' 7 хв - Х/ - А к

Х + Ак < х < Х/ + Ак + (хв - Х/ - Ак )(Т0 - Т0 VТ0

0, в остальных случаях. При т0 < т0 < 2т0 и / = к = 1, 2,...,ь

А и [X и[ J 2А

00

х

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

20 (Р(х/ ))(Т0 - 2Т0 ) _ 1

Т0 Т 0 +(Ь - 1)(Т0 - Т 0 ) Ь '

х, -А к -(х/ -А к - X н )(Т0 - Т 0 V Т 0 < X <

< Х, +Ак +(хв - Х, -Ак )(Т0 - Т0 )/т0;

0, в остальных случаях.

При т0 < т0 < 2т0, /,к = 1, 2,...,ь и к < /

А и [х ] = — х и[ J 2А

Т0 (Р(Хк- 1)(Т0 - Т0 VТ

Т Ь -1

1 + (Ь - 1) Хк - А к - х

1

Ь''

Хк - А - х„

к н

Хк -Ак -(хк -Ак - хн)(Т0 - Т0 vТ0 < х < Хк -Ак

(Р( Хк )>

Т Ь -1

1п

1 + (Ь -1)

1

Ь'

Хк - а к < х < •

т0 (Р(Хк^ 1 + (Ь - 1)(Т0 - Т0 VТ

Т Ь -1

0 и 0 7/ и 0

1 + (Ь -1)х - Хк- А

1

Ь'

х - Хк - а к

+ Ак < х < Хк +Ак +(хв - Хк -Ак)(Т0 - Т0 VТ0

0, в остальных случаях.

х

Ь

х

1

х

Х

к

к

Хк Х1

Т

0

0

х

к

Х

к

к

х

Х

к

Х

к

к

Ь

Х/ Хк

При T0 > 2т0, l,k = 1, 2,...L и k > l

Ди [x] = —— x

P X))

ln-

T L -1 , , ir Л xi - Дк - X L н

—, X„ < X < x -Д.

1 + (l -1)-

T^PxxA in L -1,

P x,))

T L -1

ln

Xl -Дк < X < X, +Дк

1 + (L - 1)

X - xi - Дк L

--, X, +Дк <X <XB

0, в остальных случаях.

При T > 2т0 И I = k = 1, 2,...,L

Ди [x] = 0, xн <x <xb .

При T0 > 2т0, l,k = 1, 2,...,L и k < l

x

Д и [X ] = J-2Д к

X0 (P( Xk)) T, L -1

in

L

1 + (l -1) xk - Д к - x

Xk - Д к - X н

P * »in L - -1,

L

T0 L -1 X0 (P( Xk )) T0 L -1

X н < X < Xk -Д к;

Xk-Д к < X < Xk +Д к;

in

L

1 + (L -1) X - Xk - Д к X н - Xk - Д к

Xk +Дк < X < Xв;.

0, в остальных случаях.

Отметим, что для k=l при т <

0 < т0

и

т0 < Т < 2т0 имеют место мультипликативные и отсутствуют аддитивные погрешности оценки щ [х] . При Т > 2То мультипликативная пропадает и появляется аддитивная погрешность (т0 - 2т0 )/2ак т0ь.

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

ностей оценок w [X] позволяют адекватно

учесть эти изменения при экстраполяции и интерполяции, повысить точность и достоверность измерений плотностей вероятностей стационарных эргодических случайных процессов [5, 9].

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

СПИСОК ЛИТЕРАТУРЫ

1. Заико А. И. Теория точности статистических и спектральных измерений // Вестник УГАТУ. 2000. № 2. С. 175182. [ A. I. Zaiko, «Theory of accuracy of statistical and spectral measurements», (in Russian), in Vestnik UGATU, no 2, pp. 175-182, 2000. ]

2. Свид. 72200700005. Случайный процесс Заико А. И. с равномерным законом распределения. Математическая модель // Заико А. И.; зарег. ФГУП «ВНТИЦ» 28.02.07, 10 с. [ Certif. 72200700005. «A. I. Zaiko random process with a uniform distribution law. The mathematical model»; zareg. FSUE "VNTIC" 28.02.07, (in Russian), 10 p. ]

3. Заико А. И. Случайный процесс Заико с равномерным законом распределения // Вестник УГАТУ. 2008. T. 11, № 1 (28). С. 188-193. [ A. I. Zaiko, «Random process Zaiko A. I. with a uniform distribution law», (in Russian), in Vestnik UGATU, vol. 11, no 1 (28), pp. 188-193, 2008. ]

4. Заико А. И. Комплексный подход к определению погрешностей // Датчики и Системы. 2007. № 8. С. 52-59. [ A. I. Zaiko, «Integrated approach to the determination of errors», (in Russian), in Sensors and Systems, no 8, pp. 52-59, 2007. ]

5. Zaiko A. I. Random signal with uniform distribution // Measurement Techniques, 1999 Vol. 42, Juni, pp. 11-13, 1999. [ A. I. Zaiko, «Random signal with uniform distribution», in Measurement Techniques, vol. 42, Juni, pp. 11-13, 1999. ]

6. Заико А. И. Случайные процессы. Модели и измерения: учеб. пособие, М.: Изд-во МАИ, 2006, 207 с. [ A. I. Zaiko, Random process. Models and measurements: (in Russian). M.: Izd-vo MAI, 2006. ]

7. Заико А. И. Оценивание плотности вероятности эргодического случайного процесса // Матер. II междунар. научно-практич. конф. «Современные проблемы науки и образования в техническом вузе». Уфа: УГАТУ, 2015, Т. 1, С. 146-152. [ A. I. Zaiko, "Estimation of probability density in ergotic random process", (in Russian), in Proc. II International scientific-practical conference "Modern problems of science and education in a technical college", vol. 1, pp. 146-152, 2015. ]

8. Заико А. И. Измерение плотности вероятности случайного процесса с линейной корреляционной функцией // Труды междунар. научно-технич. конф. «Перспективные информационные технологии (ПИТ 2016)», Самара: Изд-во СНЦ РАН, 2016, С. 77-87. [ A. I. Zaiko, «Measurement of the probability density of a random process with a linear correlation function», (in Russian), in Proc. International scientific-technical conference "Advanced Information Technology (AIT2016)", Samara: Izd-vo SNC RAN, pp. 77-87. 2016. ]

L

z

0

T L -1

L

x

L

T

0

X в - Xl - Дк

X

0

x

ОБ АВТОРЕ

ЗАИКО Александр Иванович, проф. каф. теоретических основ электротехники. Дипл. инженер электронной техники (УАИ, 1970). Д-р техн. наук по информац. измерит. системам (ЛЭТИ, 1990). Заслуж. изобретатель РБ и РФ. Дейст. член Международ. инж. акад. и Инж. акад. РБ. Иссл. в области метрологич. обеспечения, анализа и синтеза информац. измерит. систем спец. назначения и измерения случайных процессов.

METADATA

Title: Measuring probability density of ergodic random having

linear correlation function Author: A. I. Zaiko Affiliation:

Ufa State Aviation Technical University (UGATU), Russia. Email: zaiko@ugatu.ac.ru Language: Russian.

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 21, no. 4 (78), pp. 121-128, 2017. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: The paper proposes algorithms and error estimations for measuring probability density of random process having linearly decreasing correlation function, considering reconstruction of its trajectory between samples. Key words: Ergodic random process; linear correlation function; measurement algorithms; extrapolation and interpolation. About authors:

ZAIKO, Alexander Ivanovich. Prof., Dept. of Theoretiсal Basics of Electrical Engineering. Dipl. Electronic Engineer (UGATU, 1970). Cand. (PhD) Tech. Sci. (KPtI, 1973), Dr. (Habil.) Tech. Sci. (LETI, 1990). Honored inventor of RB and RF.

i Надоели баннеры? Вы всегда можете отключить рекламу.