Измерение случайных процессов

Измерение случайных процессов

Измерение случайных процессов
Реферат на тему .

Содержание

Общие сведения об измерениях. . . . . . . . . . . . . . . . . . . стр 3.
Измерения математического ожидания и дисперсии случайного процесса. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . стр 9.
Измерение функций распределения вероятности. . . . стр 11.
Измерения корреляционной функции. . . . . . . . . . . . . . стр 13.
Анализ спектра мощности. . . . . . . . . . . . . . . . . . . . . . . стр 14.
Приложения . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . стр 16.
Список литературы. . . . . . . . . . . . . . . . . . . . . . . . . . . . . стр 17.

ИЗМЕРЕНИЯ ВЕРОЯТНОСТНЫХ ХАРАКТЕРИСТИК СЛУЧАЙНЫХ ПРОЦЕССОВ
1. ОБЩИЕ СВЕДЕНИЯ
Измерения вероятностных характеристик случайных процес­сов (статистические измерения) составляют один из наиболее быстро развивающихся разделов измерительной техники. В на­стоящее время область распространения статистических методов исследования и обработки сигналов измерительной информации практически безгранична. Связь, навигация, управление, диагно­стика (техническая, медицинская), исследование среды и многие другие области немыслимы без знания и использования свойств сигналов и помех, описываемых их вероятностными характери­стиками.
Потребность в изучении свойств случайных процессов приве­ла к развитию соответствующих методов и средств (преимуще­ственно электрических). Появление анализаторов функций рас­пределения вероятностей, коррелометров, измерителей математи­ческого ожидания, дисперсиометров и других видов измерителей вероятностных характеристик открыло новые возможности в об­ласти создания современной информационной и управляющей техники.
Рассмотрим необходимые исходные определения и общие сведения о статистических измерениях.
В теории статистических измерений используют следующие понятия и их аналоги, заимствованные из теории случайных функций (аналоги из математической статистики) реализация случайного процесса (выборочная функция), мгновенное значе­ние (выборочное значение), совокупность мгновенных значений (выборка), вероятностная характеристика (предел выборочного среднего).
Введем следующие обозначения Х (t) — случайный процесс;
i-порядковый номер реализации случайного процесса Х (t);
xi(tj) —мгновенное значение процесса Х (t), соответствующее значению (i-й реализации в j-й момент времени. Случайным назы­вают процесс Х (t), мгновенные значения которого xi (tj) суть случайные величины.
На рис.1 представлена в качестве примера совокупность реализации случайного процесса, воспроизводящих зависимости некоторого параметра Х от времени t.
В теории случайных процессов их полное описание произво­дится с помощью систем вероятностных характеристик многомерных функций распределения вероятности, моментных функ­ций, характеристических функций и т. п. В теории статистиче­ских измерений исследуемый случайный процесс представляется своими реализациями, причем полное представление осуществля­ется с помощью так называемого ансамбля, т. е. бесконечной совокупностью реализаций. Ансамбль — математическая аб­стракция, модель рассматриваемого процесса, но конкретные реализации, используемые в измерительном эксперименте, пред­ставляют собой физические объекты или явления и входят в ан­самбль как его неотъемлемая часть.
Если случайный процесс представлен ансамблем реализации xi (t), i=1, 2, …, со, то вероятностная характеристика в может быть определена усреднением по совокупности, т.е.
N
[X (t)]=lim 1/N  g[xi(t)], (1)
N  i =1
где g [Xi (t)]— некоторое преобразование, лежащее в основе оп­ределения вероятностной характеристики . Так, например, при определении дисперсии g [Xi (t)]= xi (t). При этом полагаем, что процесс характеризуется нулевым математическим ожиданием.
Вместо усреднения по совокупности может быть использовано усреднение по времени с использованием k-й реализации xk (t) и тогда
T
 [X(t)]= lim 1/T  g[xi(t)]dt. (2)
T
Например, при определении математического ожидания

T
M [X (t)]= lim 1/T  xk (t) dt. (3)
T  0
В общем случае результаты усреднения по совокупности (1) и по времени (2) неодинаковы. Предел выборочного среднего по совокупности (1) представляет собой вероятност­ную характеристику, выражающую зависимость вероятностных свойств процесса от текущего времени. Предел выборочного среднего по времени (2) представляет собой вероятностную характеристику, выражающую зависимость вероятностных свойств процесса от номера реализации.
Наличие и отсутствие зависимости вероятностных характери­стик от времени или от номера реализации определяет такие фундаментальные свойства процесса, как стационарность и эрго­дичность. Стационарным, называется процесс, вероятностные ха­рактеристики которого не зависят от времени; соответственно эргодическим называется процесс, вероятностные характеристи­ки которого не зависят от номера реализации.
Следовательно, стационарный неэргодический случайный процесс — это такой процесс, у которого эквивалентны времен­ные сечения (вероятностные характеристики не зависят от теку­щего времени), но не эквивалентны реализации (вероятностные характеристики зависят от номера реализации). Нестационар­ный эргодический процесс — это процесс, у которого эквивалент­ны реализации (вероятностные характеристики не зависят от номера реализации), но не эквивалентны временные сечения (вероятностные характеристики зависят от текущего времени). Классифицируя случайные процессы на основе этих призна­ков (стационарность и эргодичность), получаем следующие четы­ре класса процессов стационарные эргодические, стационарные неэргодические, нестационарные эргодические, нестационарные неэргодические.
Учет и использование описанных свойств случайных процес­сов играет большую роль при планировании эксперимента по определению их вероятностных характеристик.
Поскольку измерение представляет собой процедуру нахож­дения величины опытным путем с помощью специальных техни­ческих средств, реализующих алгоритм, включающий в себя операцию сравнения с известной величиной, в статических изме­рениях должна применяться мера, воспроизводящая известную величину.
Типовые алгоритмы измерений вероятностных характеристик случайных процессов, различающиеся способом применения ме­ры в процессе измерений, представляются в следующем виде

* [X (t)]= KSdg [X (t)]; (4)

* [X (t)]= Sd Kg [X (t)]; (5)

* [X (t)]= Sd gK [X (t)]; (6)

где Sd—оператор усреднения; К—оператор сравнения;
* [X (t)]—результат измерения характеристики  [X (t)].
Данные алгоритмы различаются порядком выполнения опе­раций. Операция сравнения с образцовой мерой (К) может быть заключительной [см. (4)], выполняться после реализации оператора g, но до усреднения [см. (5)] и, наконец, быть началь­ной [см. (6)]. Соответствующие обобщенные структурные схе­мы средств измерений значений вероятностных характеристик представлены на рис. 2.
На этих рисунках для обозначения блоков, реализующих операторы, входящие в выражения (4) — (6), используют­ся те же обозначения. Так, g — устройство, выполняющее пре­образование, лежащее в основе определения вероятностной ха­рактеристики ; Sd — устройство усреднения (сумматор или ин­тегратор); К— компаратор (сравнивающее устройство), а М—мера, с помощью которой формируется известная величина (., g., x.)

Представленное на рис. 2, а средство измерений реализует следующую процедуру на вход поступает совокупность реализа­ций {xi (t)} (при использовании усреднения по времени — одна реализация xi, (t)-, на выходе узла g имеем совокупность преоб­разованных реализации {g[xi (t)]}; после усреднения получаем величину Sd {g[xi (t)]}, которая поступает на компаратор, осуще­ствляющий сравнение с известной величиной о, в результате чего получаем значение измеряемой вероятностной характеристики *[X(t)].
Отличие процедуры, реализуемой средством измерений, пред­ставленным на рис. 2, б, заключается в том, что после формиро­вания совокупности {g [xi (t)]} она поступает не на усреднитель, а на компаратор, который выполняет сравнение с известной вели­чиной go; на выходе компаратора формируется числовой массив {g* [xi (ti)]} и усреднение выполняется в числовой форме. На выхо­де усреднителя Sd имеем результат измерения * [X (t)].
Средство измерений (рис. 2, в) основано на формировании массива числовых эквивалентов мгновенных значений реализа­ции случайного процесса Х (t), после чего преобразование g и ус­реднение выполняются в числовой форме. Это устройство эквива­лентно последовательному соединению аналого-цифрового пре­образователя (АЦП) и вычислительного устройства (процессо­ра). На выходе АЦП формируется массив мгновенных значений, а процессор по определенной программе обеспечивает реализа­цию операторов g и Sd,
Погрешность результата измерения вероятностной характе­ристики случайного процесса
* [X(t)]=*[X(t)]-  [ X(t)]. (7)
Для статистических измерений характерно обязательное на­личие составляющей методической погрешности, обусловленной конечностью объема выборочных данных о мгновенных значени­ях реализации случайного процесса, ибо при проведении физиче­ского эксперимента принципиально не может быть использован бесконечный ансамбль реализации или бесконечный временной интервал. Соотношение (7) определяет результирующую по­грешность, включающую в себя как методическую, так и инстру­ментальную составляющие. В дальнейшем будут приводиться соотношения только для определения специфической для стати­стических измерений методической погрешности, обусловленной конечностью числа реализации и временного интервала.
2. ИЗМЕРЕНИЯ МАТЕМАТИЧЕСКОГО ОЖИДАНИЯ ИДИСПЕРСИИ СЛУЧАЙНОГО ПРОЦЕССА
Математическое ожидание и дисперсия случайного процес­са — основные числовые вероятностные характеристики, измере­ние которых играет большую роль в практике научных исследова­ний, управления технологическими процессами и испытаний.
При измерении математического ожидания результатом из­мерения является среднее по времени или по совокупности мгно­венных значений реализации исследуемого случайного процесса. Усреднение по времени применяется на практике существенно чаще, чем усреднение по совокупности, поскольку работать с од­ной реализацией удобнее и проще, чем с совокупностью. На рис. 3 приведена структурная схема устройства, реали­зующего алгоритм
t
M* [X (t)]= 1/T xk (t) dt.
t-T

На рисунке Д—преобразователь измеряемой величины в электрический сигнал (датчик); НП — нормирующий преобра­зователь, превращающий входной сигнал в стандартный по виду и диапазону значений; И — интегратор; УС — устройство сопря­жения, обеспечивающее согласование выхода интегратора со входами цифрового вольтметра и регистрирующего прибора;
ЦИП — цифровой прибор (например, цифровой вольтметр);
РП—регистрирующий прибор (самопишущий прибор).
Для оценки среднего квадратичeского значения погрешности, обусловленной конечностью объема выборочных данных,
можно пользоваться следующими соотношениями
1/2
 =[2D[X(t)]  k/T]
M
при усреднении по времени T и

1/2
 =[D[X(t)]/N]
M

при усреднении по совокупности N. Здесь D[X (t)]—дисперсия процесса X(t), а  k — интервал корреляции. Дисперсия случайного процесса характеризует математиче­ское ожидание квадрата отклонений мгновенных значений реали­зации случайного процесса от математического ожидания. Таким образом,
T 2
D[X(t)]= lim 1/T  [xk (t)-[X(t)]] dt
T 0
или
N 2
D[X(t)]= lim 1/N  [xi(t)-[X(t)]] dt
N i=1
Возможны различные варианты построения устройств для измерения дисперсии случайного процесса — дисперсиометров. На рис. 4 приведена структурная схема средства измерений дисперсии случайного процесса, т. е. работающего согласно вы­ражению
t t 2
D* [X(t)]=1/T  [xk (t)- 1/T1  xk (t)dt] dt
t-T t-T1
На рисунке НП — нормирующий преобразователь; И1 и И2 — интеграторы; ВУ— вычитающее устройство; КУ— квадратирующее устройство; УС — устройство сопряжения; ЦИП — цифро­вой прибор; РП — регистрирующий прибор.
Средняя квадратическая погрешность из-за конечности объема выборочных данных о мгновенных значениях Х (t) может быть определена с помощью соотношений
2 1/2
 =[2D[X (t)]  k/T]
M

, где D[X2 (t)]— дисперсия Х (t); T—время усред­нения.
При усреднении по совокупности N реализаций
2 1/2
 =[D[X (t)] /N]
D

3. ИЗМЕРЕНИЕ ФУНКЦИЙ РАСПРЕДЕЛЕНИЯ ВЕРОЯТНОСТЕЙ
Одномерная интегральная функция распределения вероятно­сти F (X) равна вероятности того, что мгновенное значение про­извольной реализации в произвольный момент времени меньше установленного уровня, т. е. Xi (ti)  X. Функция F (X) определя­ется как предел выборочного среднего

F (X)= lim Sd [ [x (t) ,X]],
d

1 при x (t)  X
Где [x(t) ,X]=
0 при x (t) > X

Поскольку интегральные F (X) и дифференциальные w (X) функции распределения вероятности связаны между собой со­отношениями
X
w (X) =(dF (X))/dX ; F (X)=  w (X) dX
-
справедливо выражение
w (X) = lim ((F(X+X)-F (X))/X)= lim ((Sd [[x(t) ,X]])/X)
X X

1 при X < x (t)  X+X
где  [x(t) ,X]=
0 при x (t)  X, x (t) > X+X

В качестве примера рассмотрим средство измерений для определения интегральной функции распределения вероятности уровня электрического сигнала. Схема средства измерений, реа­лизующего алгоритм
t
F* (X)=1/T   [xk(t) ,X]dt ,
t-T

показана на рис. 5, где ПУ — пороговое устройство, формиру­ющее сигнал X k (t}—X; ФУ—формирующее устройство; И—интегратор, на выходе которого получается сигнал F* (X) при установленных значениях Х и Т; УС — устройство сопряжения;
ЦИП — цифровой прибор; РП — регистрирующий прибор.
Средняя квадратическая погрешность из-за конечности объема выборки определяется для F {X) с помощью соотношения

2 1/2
 =[2(F — F )  k/T]
F

при усреднении по времени и с помощью соотношения
2 1/2
 =[2(F — F )/N]
F
при усреднении по совокупно­сти. Для (X) соответствующие соотношения имеют вид
2 1/2
 =[2(w — w X)  k/T]
w

2 1/2
и  =[(w — w X)/N]
w

В приведенных соотношениях F и w — истинные значения измеряемых функ­ций при данном X.

4. ИЗМЕРЕНИЯ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ
Для случайного процесса с нулевым математическим ожида­нием корреляционная функция равна
Rx (s,) = lim Sd[xi (t) xi-s (t-)],
d
где  и s — соответственно сдвиг во времени и в пространстве реализации перемножаемых мгновенных значений.
В практических задачах большую роль играют стационарные случайные процессы, т. е. процессы с постоянными вероятностны­ми характеристиками, не зависящими от текущего времени. Сре­ди случайных процессов можно выделить эргодические процессы, для которых
t
Rx () = lim 1/T x (t) x (t-)dt,
T 0

Большое значение корреляционного анализа в различных областях науки и техники привело к созданию множества измери­тельных приборов для измерений корреляционных функций — коррелометров.
Типовая структура коррелометра, в котором используется усреднение по времени, представлена на рис. 6. При этом реализуется следующий алгоритм
t
R*x () = 1/T xk (t) xk (t-)dt,
t-T

Как видно, после нормирующего преобразователя НП сигнал поступает в устройство временной задержки УЗ и на перемножа­ющее устройство ПУ, осуществляющее перемножение мгновен­ных значений, сдвинутых по времени на интервал т. Далее с по­мощью интегратора И выполняется усреднение, после которого результирующий сигнал через УС подается на цифровой прибор ЦИП или регистрирующий прибор РП.
Средние квадратические погрешности, обусловленные ко­нечностью объема выборочных данных о мгновенных значениях реализации процесса Х (t), оцениваются с помощью соотноше­ний

1/2
 ={2D[xk (t) xk (t-)]  k/T}
R

при усреднении по времени Т и
1/2
 ={D[xk (t) xk (t-)]/N}
R

при усреднении по совокупности.

5. АНАЛИЗ СПЕКТРА МОЩНОСТИ
Спектр мощности характеризует ее частотное распределение, и он может быть определен в соответствии со следующими форму­лами
2
Sx(w) = lim 1/T  xiT (w) 
T
Где
t -jwt’
XiT (w) =  xi (t’) e dt’
t-T
На рис. 7 изображена схема анализатора спектра мощно­сти случайного процесса Х (t).
С выхода нормирующего преобразователя НП i-я реализация случайного процесса xi (t) поступает на блок Ф, выполняющий преобразование Фурье, после чего узлом Кв производится возве­дение в квадрат и нормирование с учетом интервала усреднения Т. С помощью устройства сопряжения УС сформированный сиг­нал поступает на ЦИП и регистратор РП.
В настоящее время отечественной промышленностью серийно выпускаются анализаторы случайных процессов. К ним относят­ся многофункциональный статистический преобразователь Ф790, корреллометр Ф7016, комплекс измерителей характеристик случайных сигналов Х6-4/а, многофункциональные измерители ве­роятностных характеристик Ф36 и Ф37, анализаторы спектра Ф4326, Ф4327, Ф7058 и др. С помощью этих приборов и устройств можно измерять математические ожидания и дисперсии, а также значения функций распределения вероятности, корреляционных и спектральных функций с последующим восстановлением вида самих функций. Перечисленные анализаторы рассчитаны в ос­новном на унифицированный входной сигнал и позволяют изме­рить от 256 до 4096 ординат анализируемой функции. Погреш­ность измерения не превышает ±5 %.
Кроме того, для определения вероятностных характеристик случайных сигналов могут использоваться электроизмеритель­ные приборы, предназначенные для измерения среднего и дей­ствующего значений сигнала. Для определения среднего значе­ния применяют магнитоэлектрические приборы и цифровые ин­тегрирующие приборы. Для определения среднего квадратического отклонения используют приборы, показания которых определяются действующим значением сигнала (термоэлектри­ческие, электростатические и др.).
Корреляционные устройства получили применение в различ­ных областях науки и техники для измерения различных величин. В качестве примера можно указать корреляционное устройство для измерения скорости прокатки. Эти устройства измеряют кор­реляционную функцию, зависящую от т, которая, в свою очередь, зависит от скорости прокатки.

Список литературы

1.Метрология и электроизмерительные приборы. Душин М .Е.М. Энергоатомиздат,1986.

2.Метрология, стандартизация и измерения в технике связи. Под ред. Б.П. Хромого
М. Радио и связь, 1986.

3.Основы метрологии и стандартизации. Голубева В. П. М . Вектор, 1996.