Понятие о случайной функции.

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

Случайной функцией называется функция, значение которой при каждом данном значении аргумента (или нескольких аргументов)

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

Аргумент случайной функции или совокупность всех ее аргументов будем обозначать буквой или буквой 5 и писать, как обычно принято, в скобках за обозначением самой функции, например Если аргумент случайной функции представляет собой совокупность скалярных переменных, то его можно рассматривать как -мерный вектор. Таким образом, аргументами случайных функций в излагаемой дальше теории могут быть произвольные скалярные или векторные величины

Случайную функцию можно также рассматривать как бесконечную (в общем случае несчетную) совокупность случайных величин, зависящую от одного или нескольких непрерывно изменяющихся параметров Каждому данному значению параметра (или параметров) соответствует одна случайная величина Вместе все случайные величины определяют случайную функцию Такая трактовка случайной функции показывает, что случайная функция как объект математического исследования значительно сложнее обычной случайной величины, а именно равноценна бесконечному (в общем случае несчетному) множеству случайных величин.

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

В физике часто приходится рассматривать случайные функции координат точки пространства. Пространство с заданным в нем распределением значений некоторой величины называется полем данной величины. Случайная функция координат точки пространства приводит

(кликните для просмотра скана)

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

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

Двумерным законом распределения случайной функции называется совместный закон распределения ее значений при двух произвольно взятых значениях аргумента Вообще -мерным законом распределения случайной функции называется закон распределения совокупности ее значений при произвольно взятых значениях аргумента Мы будем характеризовать -мерный закон распределения случайной функции ее -мерной плотностью вероятности которая в общем случае зависит от значений аргумента как от параметров.

Зная двумерную плотность вероятности случайной функции, можно определить ее одномерную плотность вероятности по формуле (15.8). В результате получим соотношение

Вообще, зная -мерную плотность вероятности случайной функции, можно определить все ее плотности вероятности чисел измерений, меньших чем пользуясь формулой (15.17). В результате

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

Если значения случайной функции при любых различных значениях аргумента являются независимыми случайными величинами, то -мерная плотность вероятности случайной функции согласно формуле (16.9) и определению независимости случайных величин (§ 16), при любом выражается через ее одномерную плотность вероятности формулой

Эта формула показывает, что исчерпывающей характеристикой случайной функции с независимыми значениями является ее одномерный закон распределения.

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

данному в § 47, условный закон распределения значения случайной функции зависит только от значения случайной величины и не зависит от значений случайных величин Поэтому, применяя последовательно общую формулу (16.17), получим для -мерной плотности вероятности марковского случайного процесса формулу

Но условная плотность вероятности на основании формулы (16.6) равна:

Формулы (48.4) и (48.5) дают:

Формулы (48.1) и (48.6) показывают, что -мерная плотность вероятности марковского случайного процесса при любом может быть определена, если известна его двумерная плотность вероятности. Следовательно, двумерный закон распределения является исчерпывающей характеристикой марковского случайного процесса.

Вторым примером случайных функций, для которых исчерпывающей характеристикой является двумерный закон распределения, могут служить нормально распределенные случайные функции. Мы будем считать, что случайная функция распределена нормально, если совокупность ее значений при любом и при любых из области изменения аргумента образует нормально распределенный случайный вектор. В § 23 мы видели, что -мерный нормальный закон распределения полностью определяется математическими ожиданиями, дисперсиями и корреляционными моментами случайных величин. Но математические ожидания и дисперсии случайных величин вполне определяются одномерным законом распределения случайной функции а их корреляционные моменты - двумерным законом распределения случайной функции Следовательно, двумерный закон распределения нормально распределенной случайной функции вполне определяет ее -мерный закон распределения при любом таким образом, является исчерпывающей ее характеристикой.

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

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

Легко понять, что интеграл от случайной функции с некоррелированными (в частном случае независимыми) значениями представляет собой случайную функцию с некоррелированными (соответственно независимыми) приращениями на неперекрывающихся областях изменения аргумента. В § 54 будет показано, что интеграл от случайной функции с некоррелированными значениями имеет конечную дисперсию только в том случае, если дисперсия этой случайной функции бесконечна. Вследствие этого особенно важными для приложений являются случайные функции с некоррелированными значениями и бесконечной дисперсией, называемые обычно белыми шумами. Мы будем называть белым шумом любую случайную функцию с некоррелированными значениями, имеющую бесконечную дисперсию и конечную дисперсию интеграла от нее по любэй конечной области изменения аргумента. В основе этого термина лежат физические представления, связанные с быстро изменяющимися величинами, значения которых, разделенные очень малыми промежутками времени, практически независимы. Мы увидим дальше, что при разложении таких случайных функций на элементарные гармонические колебания гармоники всех частот оказываются одинаковыми по интенсивности. Эта аналогия с белым светом и послужила причиной того, что такие случайные функции называются белыми шумами. Это название удобно распространить на все случайные функции, обладающие перечисленными свойствами, независимо от физической (или математической) природы их аргументов.

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

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

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

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

Исчерпывающей характеристикой случайной функции является ее вероятностная мера, определение которой было дано в § 14 для любых случайных объектов, в том числе и для случайных функций. Вероятностную меру случайной функции можно определить, если известны ее законы распределения всех чисел измерений. Выделим сначала из множества всех возможных реализаций случайной функции X множество всех реализаций, значения которых в точках принадлежат данным числовым множествам Согласно определению вероятностной меры значение вероятностной меры случайной функции X, соответствующее множеству ее реализаций, определится формулой

Эта формула определяет вероятностную меру случайной функции X для всех множеств рассмотренного типа при любых и при любом выборе числовых множеств Поставим теперь в соответствие каждому значению аргумента случайной функции X некоторое числовое множество и рассмотрим множество А всех реализаций случайной функции значения которых при всех принадлежат соответствующим множествам Для того чтобы определить значение вероятностной меры случайной функции X для такого множества ее реализаций, поставим в соответствие каждому целому положительному разбиение области изменения аргумента случайной функции X на ячеек таким образом, чтобы размеры всех ячеек стремились к нулю при . В каждой ячейке разбиения выберем произвольную точку так, чтобы множество точек содержало все точки соответствующие предыдущим разбиениям. Обозначим через множество реализаций случайной функции X, значения которых в точках принадлежат соответственно множествам Тогда получим последовательность множеств реализаций случайной функции X, каждое из которых включает все последующие множества. Предположим, что произведение всех множеств (т. е. множество реализаций случайной функции X, принадлежащих всем множествам совпадает с исходным множеством реализаций А, если не считать некоторых исключительных реализаций, имеющих нулевую суммарную вероятность появления, при любом выборе такого множества реализаций А. Это предположение накладывает определенные ограничения на характер возможных реализаций случайной функции . А именно, необходимо, чтобы любое множество ее реализаций можно было определить с любой степенью точности, накладывая на них ограничения в конечном числе достаточно близких друг к другу точек. Полагая в формуле (48.7)

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

который и является значением вероятностной меры случайной функции X для рассматриваемого множества ее реализаций А.

Формулы (48.7) и (48.8) определяют вероятностную меру случайной функции для всех цилиндрических множеств реализаций. Этого достаточно для того, чтобы определить ее для любых множеств реализаций .

Для случайной функции можно также определить функционал распределения, который является естественным обобщением функции распределения случайной величины. В соответствии с определением функции распределения (14.13) функционалом распределения случайной функции X называется вероятность выполнения неравенства при всех значениях аргумента

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

Вероятностная мера и функционал распределения случайной функции пока не имеют большого практического значения, вследствие того, что методы вычисления интегралов типа (18.12) для произвольно заданной вероятностной меры в настоящее время еще очень мало разработаны .

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

где интеграл распространяется на всю область изменения аргумента Характеристический функционал случайной функции зависит от функции (т. е. от значений этой функции при всех значениях аргумента

Характеристический функционал является исчерпывающей характеристикой случайной функции Действительно, задавая функцию к как линейную комбинацию импульсных -функций:

получим на основании свойств -функции:

Сравнивая это выражение с (28.1), приходим к заключению, что величина представляет собой характеристическую функцию -мерного случайного вектора с составляющими Поэтому, применяя формулу (28.14), можно определить -мерную плотность вероятности случайной функции при любом значении Таким образом, если задан характеристический функционал случайной функции то его значения при частных видах функции определяют все законы распределения случайной функции.

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

где произвольные постоянные, а произвольные функции. Интеграл в показателе в формуле (48.11), очевидно, является линейным функционалом от случайной функции Сумма в показателе формулы (48.13) также является линейным функционалом от случайной функции Линейный функционал от функции можно сокращенно обозначать опуская скобки и обозначение аргумента функции х.

Обобщая определение (48.11), можно определить характеристический функционал случайной функции формулой

где А - произвольный линейный функционал. Задавая в формуле (48.15), линейный функционал А в виде интеграла или суммы, получим формулы (48.11) и (48.13) как частные случаи формулы (48.15). Формула (48.15) определяет характеристический функционал и в том случае, когда аргумент случайной функции X является вектором, одни составляющие которого представляют собой непрерывно изменяющиеся переменные, а другие составляющие являются дискретными переменными, в то время как формула (48.11) определяет характеристический функционал только в частном случае, когда все составляющие вектора являются непрерывно изменяющимися переменными.

Если характеристический функционал случайной функции X определяется формулой

где - некоторые функции, а индексы у линейных функционалов А указывают, к функциям каких аргументов они применяются, то характеристические функции всех чисел измерений случайной функции А

будут нормальными и, следовательно, случайная функция X распределена нормально. Таким образом, формула (48.16) определяет характеристический функционал нормально распределенной случайной функции. Эта формула является очевидным обобщением формулы (28.18) для характеристической функции нормально распределенного случайного вектора.

Пример 1. Найти плотности вероятности случайной функции скалярной независимой переменной с независимыми приращениями, если при ее значение равно нулю, а ее приращение на любом интервале распределено нормально и имеет математическое ожидание, равное нулю, и дисперсию

В данном случае значение случайной функции X при любом равно сумме ее значения при (равного нулю) и ее приращения на интервале Следовательно, одномерная плотность вероятности случайной функции X определяется формулой

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

o Случайной функцией называется функция X(t), значение которой при любом значении аргумента t является случайной величиной.

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

o Конкретный вид, принимаемый случайной величиной в результате опыта, называется реализацией случайной функции.

Т.к. на практике аргумент t чаще всего является временным, то случайную функцию иначе называют случайным процессом.

На рисунке изображено несколько реализаций некоторого случайного процесса.

Если зафиксировать значение аргумента t, то случайная функция X(t) превратится в случайную величину, которую называют сечением случайной функции , соответствующим моменту времени t. Будем считать распределение сечения непрерывным. Тогда Х(t) при данном t определяется плотностью распределения p(x; t).

Очевидно, p(x; t) не является исчерпывающей характеристикой случайной функции X(t), поскольку она не выражает зависимости между сечениями X(t) в разные моменты времени t. Более полную характеристику дает функция -совместная плотность распределения системы случайных величин , где t 1 и t 2 -произвольные значения аргумента t случайной функции. Еще более полную характеристику случайной функции X(t) даст совместимая плотность распределения системы трех случайных величин и т.д.

o Говорят, что случайный процесс имеет порядок n , если он полностью определяется плотностью совместимого распределения n произвольных сечений процесса, т.е. системы n случайных величин , где X(t i)-сечение процесса, отвечающее моменту времени t i , но не определяется заданием совместного распределения меньшего, чем n, числа сечений.

o Если плотность совместного распределения произвольных двух сечений процесса вполне его определяет, то такой процесс называется марковским.

Пусть имеется случайная функция X(t). Возникает задача описания ее с помощью одной или нескольких неслучайных характеристик. В качестве первой из них естественно взять функцию -математическое ожидание случайного процесса. В качестве второй берется среднее квадратическое отклонение случайного процесса .

Эти характеристики являются некоторыми функциями от t. Первая из них-это средняя траектория для всех возможных реализаций. Вторая характеризует возможный разброс реализаций случайной функции около средней траектории. Но и этих характеристик недостаточно. Важно знать зависимость величин X(t 1) и X(t 2). Эту зависимость можно характеризовать с помощью корреляционной функции или корреляционного момента.

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

У этих случайных процессов примерно одинаковые математические ожидания и средние квадратичные отклонения. Тем не менее это различные процессы. Всякая реализация для случайной функции X 1 (t) медленно меняет свои значения с изменением t, чего нельзя сказать о случайной функции X 2 (t). У первого процесса зависимость между сечениями X 1 (t) и будет больше, чем зависимость для сечений X 2 (t) и второго процесса, т.е. убывает медленнее, чем , при увеличении Δt. Во втором случае процесс быстрее «забывает» свое прошлое.

Остановимся на свойствах корреляционной функции, которые вытекают из свойств корреляционного момента пары случайных величин.

Свойство 1. Свойство симметричности .

Свойство 2. Если к случайной функции X(t) прибавить неслучайное слагаемое , то от этого корреляционная функция не изменится, т.е. .

Действительно,

Свойство 3. , где -неслучайная функция.

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

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

В отличие от числовых характеристик случайных величин, предоставляющих собой определенные числа, характеристики случайных функций представляют собой в общем случае не числа, а функции.

Математическое ожидание случайной функции определяется следующим образом. Рассмотрим сечение случайной функции при фиксированном . В этом сечении мы имеем обычную случайную величину; определим ее математическое ожидание. Очевидно, в общем случае оно зависит от , т. е. представляет собой некоторую функцию :

. (15.3.1)

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

По смыслу математическое ожидание случайной функции есть некоторая средняя функция, около которой различным образом варьируются конкретные реализации случайной функции.

На рис. 15.3.1 тонкими линиями показаны реализации случайной функции, жирной линией - ее математическое ожидание.

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

Дисперсией случайной функции называется неслучайная функция , значение которой для каждого равно дисперсии соответствующего сечения случайной функции:

. (15.3.2)

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

Очевидно, есть неотрицательная функция. Извлекая из нее квадратный корень, получим функцию - среднее квадратическое отклонение случайной функции:

. (15.3.3)

Математическое ожидание и дисперсия представляют собой весьма важные характеристики случайной функции; однако для описания основных особенностей случайной функции этих характеристик недостаточно. Чтобы убедиться в этом, рассмотрим две случайные функции и , наглядно изображенные семействами реализаций на рис. 15.3.2 и 15.3.3.

У случайных функций и примерно одинаковые математические ожидания и дисперсии; однако характер этих случайных функций резко различен. Для случайной функции (рис. 15.3.2) характерно плавное, постепенное изменение. Если, например, в точке случайная функция приняла значение, заметно превышающее среднее, то весьма вероятно, что и в точке она также примет значение больше среднего. Для случайной функции характерна ярко выраженная зависимость между ее значениями при различных . Напротив, случайная функция (рис. 15.3.3) имеет резко колебательный характер с неправильными, беспорядочными колебаниями. Для такой случайной функции характерно быстрое затухание зависимости между ее значениями по мере увеличения расстояния по между ними.

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

Пусть имеется случайная функция (рис. 15.3.4); рассмотрим два ее сечения, относящихся к различным моментам: и , т. е. две случайные величины и . Очевидно, что при близких значениях и величины и связаны тесной зависимостью: если величина приняла какое-то значение, то и величина с большой вероятностью примет значение, близкое к нему. Очевидно также, что при увеличении интервала между сечениями , зависимость величин и вообще должна убывать.

Степень зависимости величин и может быть в значительной мере охарактеризована их корреляционным моментом; очевидно, он является функцией двух аргументов и . Эта функция и называется корреляционной функцией.

Таким образом, корреляционной функцией случайной функции называется неслучайная функция двух аргументов , которая при каждой паре значений , равна корреляционному моменту соответствующих сечений случайной функции:

, (15.3.4)

, .

Вернемся к примерам случайных функций и (рис. 15.3.2 и 15.3.3). Мы видим теперь, что при одинаковых математических ожиданиях и дисперсиях случайные функции и имеют совершенно различные корреляционные функции. Корреляционная функция случайной функции медленно убывает по мере увеличения промежутка ; напротив, корреляционная функция случайной функции быстро убывает с увеличением этого промежутка.

Выясним, во что обращается корреляционная функция , когда ее аргументы совпадают. Полагая , имеем:

, (15.3.5)

т. е. при корреляционная функция обращается в дисперсию случайной функции.

Таким образом, необходимость в дисперсии как отдельной характеристике случайной функции отпадает: в качестве основных характеристик случайной функции достаточно рассматривать ее математическое ожидание и корреляционную функцию.

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

. (15.3.6)

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

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

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

Вместо корреляционной функции можно пользоваться нормированной корреляционной функцией:

, (15.3.7)

которая представляет собой коэффициент корреляции величин , . Нормированная корреляционная функция аналогична нормированной корреляционной матрице системы случайных величин. При нормированная корреляционная функция равна единице.

Пусть над случайной функцией X(t) проведено п независимых опытов (наблюдений) и в результате получено п реализаций случайной функции (рис. 15.4.1).

Рис. 15.4.1

Требуется найти оценки для характеристик случайной функции: ее математического ожидания m x (t), дисперсии D x (t) и корреляционной функции K x (t,t).

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

и зарегистрируем значения, принятые функцией X(t) в эти моменты времени. Каждому из моментов /, t 2 , ..., t m будет соответствовать п значений случайной функции.

Значения /, I, t m обычно задаются равноотстоящими; величина интервала между соседними значениями выбирается в зависимости от вида экспериментальных кривых так, чтобы по выбранным точкам можно было восстановить основной ход кривых. Часто бывает так, что интервал между соседними значениями t задается независимо от задач обработки частотой работы регистрирующего прибора (например, темпом киноаппарата).

Зарегистрированные значения X(t) заносятся в таблицу, каждая строка которой соответствует определенной реализации, а число столбцов равно числу опорных значений аргумента (табл. 15.4.1).

Таблица 15.4.1

		X 2 (?2)		x 2 U k )		X 2 {ti)		x 2 (J m)
								%i (tm)
		X„{t 2)		X„(tk)		X„ (?,)

В таблице 15.4.1 в /-Й строке помещены значения случайной функции, наблюденной в /-й реализации (/-м опыте) при значениях аргумента, / 2 , ..., t m . Символом Xj(4) обозначено значение, соответствующее /-й реализации в момент t k .

Полученный материал представляет собой не что иное, как результаты п опытов над системой т случайных величин

и обрабатывается совершенно аналогично (см. подраздел 14.3). Прежде всего находятся оценки для математических ожиданий по формуле

затем - для дисперсий

и, наконец, для корреляционных моментов

В ряде случаев бывает удобно при вычислении оценок для дисперсий и корреляционных моментов воспользоваться связью между начальными и центральными моментами и вычислять их по формулам:

При пользовании последними вариантами формул, чтобы избежать разности близких чисел, рекомендуется заранее перенести начало отсчета по оси ординат поближе к математическому ожиданию.

После того, как эти характеристики вычислены, можно, пользуясь рядом значений m x (t {),m x (t 2), m x (t m), построить зависимость m x (t) (рис. 15.4.1). Аналогично строится зависимость О х (/). Функция двух аргументов K x (t,t") воспроизводится по ее значениям в прямоугольной сетке точек. В случае надобности все эти функции аппроксимируются какими-либо аналитическими выражениями.

15.5. Методы определения характеристик преобразованных случайных функций по характеристикам исходных случайных функций

В предыдущем подразделе мы познакомились с методом непосредственного определения характеристик случайной функции из опыта. Такой метод применяется далеко не всегда. Во-первых, постановка специальных опытов, предназначенных для исследования интересующих нас случайных функций, может оказаться весьма сложной и дорогостоящей.

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

Ясно, что при этом непосредственное исследование случайных функций, характеризующих работу системы, нецелесообразно, а в ряде случаев вообще невозможно. В таких случаях в качестве основных рабочих методов применяются не прямые, а косвенные методы исследования случайных функций. Подобными косвенными методами мы уже пользовались при исследовании случайных величин: ряд глав нашего курса -10,11,12 - был посвящен нахождению законов распределения и числовых характеристик случайных величин косвенно, по законам распределения и числовым характеристикам других случайных величин, с ними связанных. Пользуясь совершенно аналогичными методами, можно определять характеристики случайных функций косвенно, по характеристикам других случайных функций, с ними связанных. Развитие таких косвенных методов и составляет главное содержание прикладной теории случайных функций.

Задача косвенного исследования случайных функций на практике обычно возникает в следующей форме.

Рис. 15.5.1

Имеется некоторая динамическая система А; под «динамической системой» мы понимаем любой прибор, прицел, счетно-решающий механизм, систему автоматического управления и т.п. Эта система может быть механической, электрической или содержать любые другие элементы. Работу системы будем представлять себе следующим образом: на вход системы непрерывно поступают какие-то входные данные; система перерабатывает их и непрерывно выдает некоторый результат. Условимся называть поступающие на вход системы данные «воздействием», а выдаваемый результат «реакцией» системы на это воздействие. В качестве воздействий могут фигурировать изменяющиеся напряжения, угловые и линейные координаты каких-либо объектов, сигналы или команды, подаваемые на систему управления, и т.п. Равным образом и реакция системы может вырабатываться в той или иной форме: в виде напряжений, угловых перемещений и т.д. Например, для прицела воздушной стрельбы воздействием является угловая координата движущейся цели, непрерывно измеряемая в процессе слежения, реакцией - угол упреждения. Рассмотрим самый простой случай: когда на вход системы А подается только одно воздействие, представляющее собой функцию времени х(/); реакция системы на это воздействие есть другая функция времени у (/). Схема работы системы А условно изображена на рис. 15.5.1. Будем говорить, что система А осуществляет над входным воздействием некоторое преобразование, в результате которого функция x(f) преобразуется в другую функцию у (/). Запишем это преобразование символически в виде:

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

При исследовании динамической системы в первую очередь решается основная задача: по заданному воздействию x(t) определить реакцию системы y(t). Однако для полного исследования системы и оценки ее технических качеств такой элементарный подход является недостаточным. В действительности воздействие х(/) никогда не поступает на вход системы в чистом виде; оно всегда искажено некоторыми случайными ошибками (возмущениями), в результате которых на систему фактически воздействует не заданная функция x(t), а случайная функция X(t) соответственно этому система вырабатывает в качестве реакции случайную функцию Y(t), также отличающуюся от теоретической реакции у (/) (рис. 15.5.2).

Рис. 15.5.2

Естественно возникает вопрос: насколько велики будут случайные искажения реакции системы при наличии случайных возмущений на ее входе? И далее: как следует выбрать параметры системы для того, чтобы эти искажения были минимальными?

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

Из двух поставленных выше задач, естественно, более простой является первая - прямая - задача. Сформулируем ее следующим образом.

На вход динамической системы А поступает случайная функция Х(1 ); система подвергает ее известному преобразованию, в результате чего на выходе системы появляется случайная функция:

Известны характеристики случайной функции X(t): математическое ожидание и корреляционная функция. Требуется найти аналогичные характеристики случайной функции Y(t). Короче, по заданным характеристикам случайной функции на входе динамической системы найти характеристики случайной функции на выходе.

Поставленная задача может быть решена совершенно точно в одном частном, но весьма важном для практики случае: когда преобразование А принадлежит к классу так называемых линейных преобразований и соответственно система А принадлежит к классу линейных систем.

Предварительные замечания. Найдем изображение Фурье от d -функции.

Очевидно, справедливо и обратное преобразование Фурье:

А также:

1. Пусть процесс представляет собой постоянную величину x(t)=A o . Как уже было выяснено ранее, корреляционная функция такого процесса равна Найдем спектральную плотность процесса путем прямого преобразования Фурье функции R(t):

Спектр процесса состоит из единственного пика типа импульсной функции, расположенной в начале координат. Таким образом, если в процессе присутствует только одна частота w =0, то это значит, что вся мощность процесса сосредоточена на этой частоте, что и подтверждает вид функции S(w). Если случайная функция содержит постоянную составляющую, т.е. среднее значение , то S(w) будет иметь разрыв непрерывности в начале координат и будет характеризоваться наличием d -функции в точке w =0.

2. Для гармонической функции X=A o sin(w 0 t+j) корреляционная функция:

Спектральная плотность равна

График S(w) будет иметь два пика типа импульсной функции, расположенных симметрично относительно начала координат при w= +w 0 и w= -w 0 . Это говорит о том, что мощность процесса сосредоточена на двух частотах +w 0 и -w 0 .

Если случайная функция имеет гармонические составляющие, то спектральная плотность имеет разрывы непрерывности в точках w = ±w 0 и характеризуется наличием двух дельта-функций, расположенных в этих точках.

Белый шум . Под белым шумом понимают случайный процесс, имеющий одинаковые значения спектральной плотности на всех частотах от -¥ до +¥ : S(w ) = Const.

Примером такого процесса при определенных допущениях являются тепловые шумы, космическое излучение и др. Корреляционная функция такого процесса равна

Таким образом R(t) представляет собой импульсную функцию, расположенную в начале координат.

Этот процесс является чисто случайным процессом, т.к. при любом t ¹0 отсутствует корреляция между последующими и предыдущими значениями случайной функции. Процесс с такой спектральной плотностью является физически нереальным, т.к. ему соответствуют бесконечно большие дисперсия и средний квадрат случайной величины:

Такому процессу соответствует бесконечно большая мощность и источник с бесконечно большой энергией.

2. Белый шум с ограниченной полосой частот. Такой процесс характеризуется спектральной плотностью вида

S(w)=C при ½w½ <w n ,

S(w) =0 при ½w½>w n .

где (-w n , w n) полоса частот для спектральной плотности.

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

Дисперсия процесса равна

5. Типовой входной сигнал следящей системы. В качестве типового сигнала принимают сигнал, график которого показан на рис.63. Скорость вращения задающего вала следящей системы сохраняет постоянное значение в течение некоторых интервалов времени t 1 , t 2 ,...

Переход от одного значения к другому совершается мгновенно. Интервалы времени подчиняются закону распределения Пуассона. Математическое ожидание

Рис.63. Типовой сигнал

График такого вида получается в первом приближении при слежении РЛС за движущейся целью. Постоянные значения скорости соответствуют движению цели по прямой. Перемена знака или величины скорости соответствует маневру цели.

Пусть m -среднее число перемен скорости за 1 с. Тогда Т=1/m будет среднее значение интервалов времени, в течение которых угловая скорость сохраняет свое постоянное значение. Применительно к РЛС это значение будет средним временем движения цели по прямой. Для определения корреляционной функции необходимо найти среднее значение произведения

При нахождении этого значения могут быть два случая.

1. Моменты времени t и t+t относятся к одному интервалу. Тогда среднее произведения угловых скоростей будет равно среднему квадрату угловой скорости или дисперсии:

2. Моменты времениt и t+t относятся к разным интервалам. Тогда среднее произведения скоростей будет равно нулю, так как величины W(t) и W(t+t) для разных интервалов можно считать независимыми величинами:

Корреляционная функция равна:

где, Р 1 - вероятность нахождения моментов времени t и t+t в одном интервале, а Р 2 =1- Р 1 вероятность нахождения их в разных интервалах.

Оценим величину Р 1 . Вероятность появления перемены скорости на малом интервале времени Dt пропорциональна этому интервалу и равна mDt или Dt/Т. Вероятность отсутствия перемены скорости для этого же интервала будет равна 1-Dt/Т. Для интервала времени t вероятность отсутствия перемены скорости т.е. вероятность нахождения моментов времени t и t+t в одном интервале постоянной скорости будет равна произведению вероятности отсутствий перемены скорости на каждом элементарном промежутке Dt, т.к. эти события независимые. Для конечного промежутка получаем, что число промежутков равно t/Dt и

Перейдя к пределу, получим