Методы и средства оперативного анализа случайных процессов - Пивоваров Ю.Н., Тарасов В.Н., Селищев Д.Н

где Pl−1 (X A ) - многочлены Чебышева П.Л.

Первые два из этих многочленов имеют вид

P0 (X A ) =1; P2 (X A ) = X A − n 2−1,

a остальные определяются по формуле

Коэффициенты α1 ,...,αl в формуле (2.27) также находятся по методу

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

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

Таким образом, рассмотрены способы определения математической зависимости между двумя составляющими объекта измерения. Точно так же решается задача и определения математической зависимости одной составляющей объекта измерения Xk от нескольких XA, XB, .Разница заключается лишь в том, что в данном случае уравнение регрессии надо искать в виде

где α1 ,...,αel как и ранее, неопределенные коэффициенты, значения которых должны быть найдены по принципу наименьших квадратов.

2.4 Математическое описание составляющих объекта измерения

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

111

математическому описанию самих составляющих, которые являются, входными сигналами ИИС.

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

В качестве обобщенной модели какой-то k-той составляющей объекта измерения можно взять модель вида

где φ(t) и Ψ(t) - некоторые детерминированные функции времени; N(t) – случайная функция времени.

Из формулы (2.31) следует, что для математического описания составляющей объекта измерения Xk(t) нужно уметь описывать детерминированные компоненты и случайную N(t).

2.4.1Методы представления детерминированных компонент составляющих объекта измерения.

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

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

k =0

где Ak –коэффициенты разложения;

φ0 (t),...,φk (t) - ортогональные координатные функции, т.е. такие, что

Здесь p(t)-весовая функция.

В качестве координатных функций могут выступать самые разнообразные функции. Так, если функция φ(t) рассматривается на конечном интервале времени от T1 до T2, то в качестве координатных функций могут быть выбраны различные ортогональные полиномы Чебышева, Лежандра и др.

Наиболее часто в качестве координатных ортогональных функций выбираются тригонометрические функции. В этом случае детерминированная функция (p(t), рассматриваемая на конечном интервале времени от T1 до T2, может быть представлена в виде ряда Фурье:

его энергии относится так называемая мощность сигнала и его действующее значение.

Мощностью сигнала φ(t) , существующего на интервале времени от T1 до T2, называется величина

Из формул (2.35), (2.39) и (2.40) видно, что все энергетические характеристики сигнала (энергия Е, мощность Р, и действующее значение) жестко связаны между собой.

Если, сигнал φ(t) представлен в виде ряда Фурье (2.33) или (2.34) то, как следует из выражений, (2.36) и (2.39), его мощность может быть определена по формуле

k =−∞

где C*k - комплексная величина, сопряженная с Ck.

В тех случаях, когда детерминированный сигнал φ(t) является

114

S(w) называется просто спектром сигнала.

Энергия сигнала φ(t) представленного в виде выражения может быть подсчитана по формуле:

a верхняя граничная частота Eb его спектра определяется из уравнения

Уравнением (2.49) целесообразно пользоваться при известной энергии сигнала, а уравнением (2.50) - при неизвестной.

Большое значение для математического описания сигналов имеет теорема Котельникова, которая утверждает, что непрерывная функция времени φ(t) не содержащая частот выше граничной wb = 2πFb , полностью

определяется отсчетами мгновенных значений φ(k∆t) в точках, отстоящих

друг от друга на интервалы ∆t = π . wb

Эта теорема позволяет представить непрерывную функцию φ(t) в виде

Если функция ц(t) с ограниченным спектром рассматривается на

конечном интервале времени Т, то точное разложение (2.51) заменяется приближенными:

Таким образом, в данном случае функция определяется в виде конечного числа n=2FbT ее отсчетов.

115

2.4.2Методы представления случайных компонент составляющих объекта измерения

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

Случайный сигнал (процесс) N (t) в общем случае может быть охарактеризован его m-мерной плотностью вероятности системы m случайных величин (N(t1),….., N(tm)), где t1,...,tm- произвольные значения аргумента t.

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

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

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

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

Для оценки интервала частот, в котором существует стационарный случайный процесс, вводят понятие эквивалентной ширины спектра мощности, которую определяют по формуле

∞∫SN (w)dw

Где SNm (w)- максимальное значение спектральной плотности.

116

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

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

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

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

3Методы оценки характеристик составляющих объекта измерения

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

анекоторые их приближенные значения, называемые оценками.

^

Под оценкой Θ характеристики Θ составляющей X (t) понимается некоторое преобразование последней

оцениваемой характеристики. В тоже время должно обеспечиваться как можно меньшее отклонение оценки от оцениваемой характеристики.

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

117

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

отклонение оценки.

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

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

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

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

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

3.1 Методы оценки одномерных моментных характеристик

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

118