rpd000003197 (1012246), страница 37
Текст из файла (страница 37)
Выполняя частотную фильтрацию данных с целью сглаживания и очистки от шума, необходимо правильно выбрать полосу пропускания ФНЧ. Если ее выбрать слишком большой, то степень сглаживания будет недостаточной, а шум будет подавлен не полностью. Если полоса пропускания будет слишком узкой, то вместе с шумом могут оказаться подавленными и изменения, несущие полезную информацию. Если в технических приложениях существуют строгие критерии для определения оптимальности характеристик фильтров, то в аналитических технологиях приходится использовать в основном экспериментальные методы.
Спектральный анализ является одной из наиболее эффективных и хорошо разработанных методов обработки данных. Частотная фильтрация – только одно из его многочисленных приложений. Кроме этого он используется в корреляционном и статистическом анализе, синтезе сигналов и функций, построении моделей и т.д.
Лекция 11.doc
Лекция 11. Методы обработки и анализа стохастических сигналов.
Стохастическими называются сигналы (процессы), развитие которых во времени или пространстве подчиняется вероятностным законам. Иными словами, в каждом новом опыте, проводимом в неизменных условиях, всякий раз процесс принимает новое значение, причем заранее неизвестно, какое. Существует большое число разнообразных классификаций случайных процессов. Представляется целесообразным выделить те их разновидности, которые встречаются в экспериментальных исследованиях или полезны как достаточно удобная математическая модель, адекватная реальным процессам.
С этих позиций целесообразно разделить все стохастические сигналы на две группы [23]:
-
случайные процессы;
-
случайные поля.
Ранее в главе 1 при обсуждении классификации экспериментальных данных были приведены определения каждого из этих типов стохастических сигналов. С целью сохранения связности изложения повторим их.
Скалярный случайный процесс X(t)- есть функция неслучайного скалярного аргумента t (времени), которая при каждом фиксированном значении последнего является случайной величиной. Векторный случайный процесс – это совокупность нескольких скалярных случайных процессов, рассматриваемых совместно X(t)=( X1(t), X2(t),…, Xn(t))Т
Случайные функции, аргументом которых являются векторы, называются случайными полями, обозначаемыми как X(l1,l2,….,lk,t). Причем, если аргументы случайного поля представлены только координатами пространства l1,l2,….,lk, то такое случайное поле называется пространственным. Если же в состав аргументов случайного поля входит также и время t, то случайное поле называется пространственно-временным. Наиболее часто на практике встречаются двумерные X(l1,l2,), и трехмерные X(l1,l2,,l3) случайные поля. К случайным полям такого типа относятся двумерные (плоские) и трехмерные изображения.
11.1. Основные характеристики случайных процессов.
Описание статистических свойств случайного процесса X(t) базируется на рассмотрении множества случайных величин X1 =X(t1), X2 =X(t2) ,…., Xn =X(tn), соответствующих значениям процесса в различные моменты времени ti,i=1,..n на интервале его изменения (рис. 5.4). Считается [23], что случайный процесс исчерпывающе описан в статистическом смысле, если известны:
p(X1,t1)- одномерная плотность вероятностей значений случайного процесса в момент времени t= t1.
p(X1,t1; X2,t2)- двумерная совместная плотность распределения вероятностей значений случайного процесса в моменты времени t= t1, t= t2.
p(X1,t1;X2,t2;…. ;Xn,tn)- многомерная совместная плотность распределения вероятностей значений случайного процесса в моменты времени t= t1, t= t2,….. t=tn при n ).
Одномерная плотность распределения вероятностей p(X1,t1)-определяет вероятность попадания значений случайного процесса X(t) в момент времени t= t1 в интервал [1,1]:
Двумерная плотность p(X1,t1; X2,t2) позволяет определить вероятность того, что значения случайного процесса X(t) в моменты времени t= t1 , t= t2 будут лежать в интервалах [1,1], [2,2]:
Соответственно, многомерная плотность p(X1,t1;X2,t2;…. ;Xn,tn) позволяет определить вероятность того, что значения случайного процесса X(t) в моменты времени t= t1, t= t2,….. t=tn будут находится внутри n-мерного параллелепипеда, со сторонами [1,1], [2,2],….,[n,n]:
В отношении рассмотренных плотностей вероятностей имеют место следующие предельные свойства
lim p(X1,t1;X2,t2)= p(X1,t1),
t2-t10
lim p(X1,t1;X2,t2)= p(X1,t1) p(X2,t2).
t2-t1
Для статистического описания случайных процессов могут также использоваться условные плотности распределения вероятностей. Условная плотность вероятностей p(Xj,tj/X*i,ti) характеризует распределение вероятностей случайного процесса в момент времени t= tj ti при условии, что в момент времени t= tj наблюдалось конкретное значение (реализация) случайного процесса равное X(ti)= X*i. Соответственно, условная плотность p(Xn,tn / X*1,t1;X*2,t2;…. ;X*n-1,tn-1) характеризует плотность распределения вероятностей значений случайного процесса в момент времени t=tn при условии, что в предшествующие моменты времени t= t1, t= t2,….. t=tn-1 наблюдались конкретные значения случайного процесса равные: X(t1)= X*1, X(t2)= X*2,….., X(tn-1)= X*n-1.
Рассмотренное понятие плотности распределения вероятностей, как наиболее исчерпывающей характеристики статистических свойств случайных процессов позволяет провести их классификацию в зависимости от тех свойств, которыми обладают их совместные безусловные и условные плотности вероятностей.
1. Абсолютно случайный процесс. Процесс X(t) называется абсолютно случайным, если его значения Xi=X(ti), Xj =X(tj) независимы при сколь угодно малом t= ti - tj. Совместная плотность p(X1,t1;X2,t2;…. ;Xn,tn) распределения значений абсолютно случайного процесса X(t) в моменты времени t= t1, t= t2,….. t=tn равна произведению одномерных плотностей:
p(X1,t1;X2,t2;…. ;Xn,tn)= p(X1,t1) p(X2,t2)…. p(Xn,tn).
2. Марковский случайный процесс. Зададим на интервале [t0,tk] возможного изменения аргумента t случайного процесса X(t) временной ряд: t0 t1 …… tk. Случайный процесс X(t) называют марковским, если в отношении его справедливо соотношение для любых tj:
p(Xj,tj / X0,t0;X1,t1;…. ;Xj-1,tj-1)= p(Xj,tj / Xj-1,tj-1).
То есть специфика марковского случайного процесса проявляется в том, что для него условная плотность случайной величины Xj= X(tj) зависит только от того, каким было значение случайной величины Xj-1= X(tj-1) в предшествующий момент и не зависит от того, какими были реализации случайного процесса в предыдущие моменты tj-2,tj-3,… t0. Плотность p(Xj,tj / Xj-1,tj-1)называют переходной плотностью вероятностей марковского случайного процесса. Для марковского случайного процесса, учитывая его свойство, справедливо:
p(X0,t0;X1,t1;…. ;Xn,tn)= p(X0,t0) p(X1,t1 / X0,t0) p(X2,t2 / X1,t1)… p(Xn,tn / Xn-1,tn-1).
Иными словами, для исчерпывающего описания марковского случайного процесса достаточно задать его начальную одномерную плотность вероятностей p(X0,t0) и переходные плотности p(Xj,tj / Xj-1,tj-1).Среди марковских случайных процессов выделяют марковские случайные процессы с независимыми приращениями. Случайный процесс X(t) называют процессом с независимыми приращениями, если для любых значений t0 t1 …… tk., выбранных на интервале [t0,tk] приращения Xj= X(tj)- X(tj-1), j=1,…,k – независимые случайные величины. Очевидно, что любой случайный процесс с независимыми приращениями также является марковским, поскольку значение случайной величины X(tj) в конце каждого интервала tj= tj- tj-1 определяется предыдущим значением X(tj-1) и приращением на этом интервале Xj , не зависящим от приращений на предшествующих интервалах.
3. Гауссовский случайный процесс. В практических приложениях достаточно широкое распространение получили гауссовские случайные процессы. Процесс X(t) называется гауссовским, если совместная плотность распределения вероятностей p(X1,t1;X2,t2;…. ;Xn,tn) при любом n и любых t1,t2,,….,tn является гауссовской.
11.2. Корреляционный анализ случайных процессов.
При практической обработке экспериментальных данных для анализа случайных процессов , как и для описания случайных величин, часто используют их моментные характеристики, получение которых не сопряжено с такими трудностями, как получение на основе экспериментальных данных оценок многомерных плотностей распределения.
Простейшей моментной характеристикой случайного процесса X(t) является математическое ожидание mX( t ), представляющее собой неслучайную функцию времени t. Математическое ожидание случайного процесса mX(t) связано с одномерной функцией плотности p(X,t) следующим соотношением:
То есть математическое ожидание характеризует лишь среднее значение случайного процесса X(t) по времени. Оно не содержит в себе информацию ни о том, каков разброс значений возможных реализаций случайного процесса, ни о степени статистической взаимосвязи (корреляции) возможных реализаций случайного процесса в различные моменты времени. Предположим, что в результате проведенных экспериментальных исследований накоплены реализации Xi(t), i=1,…,N случайного процесса X(t). Тогда выборочная оценка математического ожидания m*X(t) может быть получена следующим образом: