Лабораторная работа 6

advertisement
Лабораторная работа 6
Тема: Расчет и отображение элементарных статистик сигналов
Цель работы: разработка средств расчета и отображения элементарных статистик
цифровых сигналов
Об элементарных статистиках сигналов
Анализируемый сигнал, реальный или модельный, можно рассматривать как конечную
реализацию некоторого случайного процесса X t  или X n . Случайный процесс в самом
общем случае описывается набором многомерных функций распределения своих
случайных отсчетов, взятых в различные моменты времени: одномерными - f1 x1 ;t1  ,
двумерными f 2 x1 , x2 ; t1 , t 2 , трехмерными f 3 x1 , x2 , x3 ; t1 , t 2 , t 3  и т.д. Информация об
особенностях одномерной ф.р. f1 x1 ;t1  хоть и недостаточна для полного описания всех
свойств процесса, но тем не менее уже весьма полезна. В случае, если процесс является
стационарным, его одномерная ф.р. не зависит от времени (номера отсчета в дискретном
сигнале). Это дает возможность оценить числовые характеристики одномерной ф.р. (т.н.
статистики первого порядка или элементарные статистики) путем соответствующего
усреднения по совокупности всех отсчетов сигнала, рассматриваемых как выборка
независимых значений случайной величины X
Наиболее часто рассчитываются следующие элементарные статистки
1 N 1
1) Среднее xˆ   xn 
N n 0
2
1 N 1
2) Дисперсия ˆ 2   x n   xˆ
N n 0


3) Среднеквадратичное отклонение  
4) Коэффициент вариации r 

x
1
5) Коэффициент асимметрии ˆ  N
1
6) Коэффициент эксцесса ˆ  N
 2
 xn   xˆ 
3
ˆ 3
 xn   xˆ 
4
ˆ 4
3
7) Минимальное значение сигнала
xmin  min xn, n  0, N  1
8) Максимальное значение сигнала
xmax  max xn , n  0, N  1
9) Квантиль порядка 0.05
x0.05
10) Квантиль порядка 0.95
x0.95
Примечание: квантиль порядка p случайной величины X, обозначаемый x p , определяется
из условия: вероятность P( X  x p )  p . Имея выборку наблюдений случайной величины X
- xn, n  0, N  1 , оценку квантиля x p производят следующим образом. Составляется
т.н. вариационный ряд выборки x0 , x1  ...  xN  1 из упорядоченных по
возрастанию значений исходной выборки xn, n  0, N  1 . В качестве оценки квантиля
порядка p берется отсчет вариационного ряда с номером m=   p * N , где  оператор
взятия целой части от аргумента в скобках, в нотации языка Си xˆ p  (int)( p * N ) .
11) Медиана – квантиль x0.5 , имеющий смысл средней точки области возможных
значений случайной величины X, такой, что вероятности встретить большее и меньшее,
чем она, значения X равны между собой и равны 0.5. В качестве оценки медианы берут
среднее по порядку значение вариационного ряда выборки xˆ 0.5  xN / 2
12) Гистограмма распределения значений сигнала G(k) , k  0, K - 1 , представляющая
собой оценку одномерной плотности распределения значений сигнала f x  .
Отображается в виде столбцовой диаграммы (см. рис.1).
Рис. 1 Пример окна с отображением статистик и гистограммы распределения случайной
величины, рассчитанной по конечной выборке
Гистограмма рассчитывается следующим образом. Область возможных значений сигнала
xmin , xmax  разбивается на K интервалов постоянной ширины h  xmax  xmin  / K , значение
Gk  есть оценка вероятности попадания значений сигнала в k-й интервал,
N (k )
рассчитываемая по правилу G k  
, где N k  - число точек сигнала, значения
N
которых попали в k-й интервал xmin  k * h, xmin  k  1 * h ., k  0, K  1 . Количество
интервалов K задается пользователем. Иногда пользователю предоставляется
возможность указать ширину интервала h, тогда K рассчитывается автоматически по
правилу K  xmax  xmin  / h .
Задание к л/р № 6
1. Реализовать в программе блок расчета элементарных статистик
Например, так.
В основном меню резервируется пункт – Статистики. При его вызове выпадает список
названий всех сигнальных каналов. Кликая мышкой на нужный канал, пользователь
инициирует процедуру расчета и визуализации элементарных статистик данного канала.
Если это первый вызов процедуры «Статистики», то открывается общее окно
«Статистики сигналов» и в нем отображаются название указанного канала, рассчитанные
значения всех 11 выше указанных статистик и гистограмма - 12. При запросах статистик
для новых каналов, в поле основного информационного окна статистик резервируется
субокно, в котором отображается информация о новом канале. Расчет статистик
осуществляется для фрагментов сигналов n1 , n2  , актуальных в данный момент (см.
Работу 3). При любом изменении в окне отображения осциллограмм размеров или
положения актуального фрагмента должен происходить перерасчет элементарных
статистик для всех каналов, отображаемых окне «Статистики сигналов».
Примечание: желательно реализовать альтернативный способ вызова процедуры расчета и
отображения статистик – из меню, выпадающего при нажатии правой кнопки мыши в
поле осциллограммы соответствующего канала либо в окне навигации каналов.
Все
Download