классические методы корреляционно

advertisement
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ
ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАМЫШИНСКИЙ ТЕХНОЛОГИЧЕСКИЙ ИНСТИТУТ (ФИЛИАЛ)
ВОЛГОГРАДСКОГО ГОСУДАРСТВЕННОГО ТЕХНИЧЕСКОГО УНИВЕРСИТЕТА
КАФЕДРА «АВТОМАТИЗИРОВАННЫЕ СИСТЕМЫ
ОБРАБОТКИ ИНФОРМАЦИИ И УПРАВЛЕНИЯ»
КЛАССИЧЕСКИЕ МЕТОДЫ
КОРРЕЛЯЦИОННО-СПЕКТРАЛЬНОГО
АНАЛИЗА ВРЕМЕННЫХ РЯДОВ
Методические указания к выполнению
лабораторной работы
РПК «Политехник»
Волгоград
2004
УДК 517.518
К 47
КЛАССИЧЕСКИЕ МЕТОДЫ КОРРЕЛЯЦИОННО-СПЕКТРАЛЬНОГО АНАЛИЗА ВРЕМЕННЫХ РЯДОВ. Методические указания к выполнению лабораторной работы / Сост. Р. Н. Каримов; Волгоград. гос.
техн. ун-т − Волгоград, 2004. 17 с.
Изучаются основные классические методы непараметрического оценивания корреляционной функции и спектральной плотности стационарных временных рядов.
Предназначены для студентов специальности «Автоматизированные
системы обработки информации и управления».
Илл. 3. Библиогр.: 6 назв.
Рецензент В. А. Марченко
Печатается по решению редакционно-издательского совета
Волгоградского государственного технического университета
© Волгоградский
государственный
технический
университет, 2004
ВВЕДЕНИЕ
Во многих отраслях науки и техники перед исследователем возникает задача, как на основе экспериментальных данных случайных процессов (СП) с непрерывным или дискретным временем, полученных на конечном интервале времени, сформировать максимально достоверное
представление о динамической модели этого процесса. В случае дискретного времени наблюдения СП представляют собой случайную последовательность или временной ряд (СВР).
Важнейшими динамическими характеристиками СВР являются автокорреляционная последовательность (АКП) и спектральная плотность
мощности (СПМ), которые широко используются при решении многих
задач АСУТП. Например, оценки АКП нужны в задачах контроля, при
планировании эксперимента, идентификации объектов управления и действующих в них возмущений.
В работе рассматриваются классические (непараметрические) оценки АКП и СПМ стационарных СВР с нормальным законом распределения вероятностей. Выполнение работы позволит получить навыки обработки экспериментальных данных в рамках корреляционной теории стационарных СП.
Цель работы: изучение непараметрических методов получения
оценок АКП и СПМ стационарного СВР.
ЗАДАНИЕ
1.
Получить у преподавателя параметры формирования авторегрессионной модели СВР.
2.
Составить программу для получения реализации заданной модели ряда авторегрессии.
3.
Вычислить теоретическую АКП и СПМ заданной модели.
4.
Выбрать необходимый объем выборки наблюдений, удовлетворяющий требуемую точность оценки.
5.
Составить программу нахождения точечных и интервальных
классических оценок АКП и СПМ.
6.
Построить таблицу значений точечных и интервальных оценок
АКП и СПМ и по этим данным построить графики точечных и интервальных оценок.
ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ
1. ОЦЕНКА АВТОКОРРЕЛЯЦИОННОЙ ФУНКЦИИ
Пусть
{x[nt ]}
1.1. Точечная оценка АКП
задана реализация стационарного эргодического СВР
с равномерным шагом дискретности t  ti  ti 1 [5]. Если
имеется N отсчётов ряда, то оценка АКП такого ряда определяется по
формуле:
Rˆ xx [kt ] 
1
N k
( x[nt  k]  mˆ x )( x[nt ]  mˆ x ). (1)
( N  k )  n 1
Здесь
1 N k
mˆ x   n 1 x[nt ]
N
(2)
оценка среднего СВР,  – шаг дискретности корреляционного сдвига,
определённым образом связанный с t . В дальнейшем будем полагать,
что t    1.1
Оценка Rˆ xx [ k ] является случайной функцией. Наиболее полной характеристикой оценки является многомерная функция распределения.
Однако, ввиду сложности последней, приходится ограничиваться вычислениями лишь простейших числовых характеристик распределения оценки, таких, как математическое ожидание M {Rˆ xx [ k ]} , дисперсия
D{Rˆ xx [ k ]} , которые в свою очередь, зависят от ненаблюдаемой истинной
АКП
R x [k ] и величин N, t .
При этом нужно различать два случая оценивания: а) среднее
временного ряда
x[k ]
известно, б) среднее
mx
временного ряда
mx
x[k ]
неизвестно.
а. Cреднее значение
mx
ряда
x[n] известно. Тогда Rˆ x [k ] явля-
ется несмещённой оценкой истинной АКП [1]
M {Rˆ xx [k ]} 
1
N k
M {~
x [n  k ]~
x [n]}  Rˆ xx [k ] ,

n

1
N k
(3)
Будем полагать, что шаг дискретности выбран в соответствии с условиями
теоремы Котельникова.
1
где ~
x [n ]  x[n ]  m x – центрированная реализация временного ряда.
Дисперсия
 2Rˆ (k )
оценки для СВР с нормальным законом рас-
xx
пределения определится по соотношению [1]:
D{Rˆ xx [k ]}   2Rˆ (k ) 
xx
2
N k
{R [n]  Rx [n  k ]Rxx [n  k ]} (4)
N  k  n 1 xx
Для эргодических процессов
lim Rxx [k ]  0 , поэтому значения
k 
АКП вычисляют только для k<<N, т. е. индекс корреляционного сдвига
во много раз меньше числа отсчётов исходных данных. При конечном k и
увеличении N значение дисперсии стремится к нулю, следовательно,
Rˆ xx [k ] является состоятельной оценкой R xx [k ] .
Рассмотрим
 2Rˆ (k ) , при k = 0 и значениях k больших или равных
xx
интервалу корреляции
N k . Тогда дисперсии оценок можно определить
по приближенным формулам [4]:
 2Rˆ (0)  (4 4xx / N ) кв
k
при k = 0,
(5)
xx
, при
 2Rˆ (k )  (2 4xx / N ) кв
k
xx
k  Nk ,
(6)
где

2
 кв
k  2 k 1 rxx [ k ]
квадратичный интервал корреляции, а
rxx2 [k ]
(7)
– нормированная АКП. В
промежуточных точках 0  k  N k значение дисперсии можно вычислить по формуле:
2
.
 2Rˆ (k )  (2 4x / N ) кв
k {1  rxx [k ]}
(8)
xx
Из сравнения формул (5) и (6) видно, что дисперсия оценки АКП с
ростом k от нуля до k  N k убывает в два раза, а далее остается постоянной при всех
k  N k . Отметим, что  2ˆ (k ) не может служить меR
xx
рой точности оценки без указания величины самой оценки Rˆ xx [ k ] . В качестве меры точности в теории оценивания используют коэффициент
изменчивости. Коэффициент изменчивости оценки АКП находится по
формуле
FRˆ
xx

 Rˆ (k )
xx
M {Rxx [k ]}

 Rˆ (k )
xx
Rxx [k ]
.
(9)
Для эргодических процессов величина
до нуля, а
 Rˆ
xx
R xx [k ] с ростом k убывает
(k ) равна постоянной, поэтому величина FRˆ xx увели-
чивается неограниченно. Отсюда следует, что значения оценки АКП при
больших k, как правило, не несут какой-либо информации о значениях
R xx [k ] . Многие авторы рекомендуют ограничиться рассмотрением оценок Rˆ xx [ k ] только до точки k  k y , в которой F ˆ (k y , N )  1 . При
R xx
FRˆ (k y , N )  1 имеем 100 % погрешность, т. е. Rˆ xx [k y ]   Rˆ (k y ) .
xx
Отметим, что полученные оценки
 Rˆ (k )
xx
ками, так как входящие в формулы величины
xx
не являются статисти-
2x  Rxx [0] , rxx [k ] не
наблюдаются. Для получения статистик в формулах (5–9), эти величины
нужно заменить оценками S x2  Rˆ xx [0] , rˆxx [ k ] , а при определении  кв
k
в формуле (7) верхний предел суммы взять равным N k .
б. Среднее арифметическое m x неизвестно. При неизвестном
среднем оценка АКП вычисляется с помощью формул (1) и (2). В этом
случае оценка АКП будет смещенной и величину смещения и дисперсии
оценки можно вычислить используя формулы (18.43) и (18.44), приведенные в [4], п. 18.2.4. Однако, если объем выборки превышает 10 интервалов корреляции, то величиной смещения можно пренебречь, а дисперсию оценки АКП можно вычислять по формулам (5), (6) и (8) для известного среднего.
1.2. Интервальная оценка АКП
Если N  10 N k , то распределение оценки Rˆ xx [ k ] можно считать
приближённо нормальным с параметрами R xx [k ] и
вероятность  попадания
 Rˆ (k ) [1]. Тогда
xx
R xx [k ] в интервал Rˆ xx [k ]   будет равна:
P{| Rxx [k ]  Rˆ xx [k ] | }   ,
P{ Rˆ xx [k ]  Rxx [k ] /  Rˆ
или
xx
  /  Rˆ }  .
10)
xx
Величина
t  {Rˆ xx [k ]  Rxx [k ]} /  Rˆ
xx
нормально распределена с параметрами N (0, 1) . Если введем обозначение t   /  ˆ , то вместо равенства (10) получим
R xx
P{ t  t }   .
Величина
(11)
t для заданной вероятности  легко находится по табли-
цам нормального закона распределения. Действительно, из интегрального нормального закона распределения в силу его симметрии имеем
  (t )  (t )  2(t )  1 .
Отсюда
t   1 ((1  ) / 2) ,
  t Rˆ ( N ).
и
xx
Таким образом, интервальная оценка АКП равна:
I R xx  ( Rˆ xx [k ]  t S Rˆ (k )) ,
(12)
xx
где S (k ) является оценкой  2 (k ) вычисляемой по формулам (5, 6, 8).
Rˆ
Rˆ
xx
xx
2. ОЦЕНКА СПЕКТРАЛЬНОЙ ПЛОТНОСТИ МОЩНОСТИ
Известны два эквивалентных классических непараметрических метода определения СПМ. Один из них основан на вычислении квадрата
модуля прямого преобразования Фурье исходной реализации и носит
название метода периодограмм. Другой метод использует преобразование Фурье предварительно вычисленной АКП и называется коррелограммным методом.
2.1. Периодограммная оценка СПМ
Рассмотрим сначала реализацию x(t ), t [T / 2, T / 2] непрерывного СП с нулевым средним. Оценку СПМ этой реализации определим прямым преобразованием Фурье по соотношению
Sˆ xx () 
T /2
 T / 2 x(t )e  jt dt
2
.
(13)
Эта оценка СПМ является несостоятельной. Для получения состоятельной оценки используют периодограмму, которую определяют усреднением СПМ отдельных сегментов. Исходную реализацию разбивают на M
сегментов, каждый из которых длины Lr  T / M . После этого по формуле (13) находят оценки СПМ
(r )
Sˆ xx
()
для каждого r-го сегмента. Оцен-
ку СПМ наблюдаемой реализации находят как среднеарифметическое
M
(r )
Sˆ xx ()  (1 / n) r 1 Sˆ xx
() .
(14)
Математическое ожидание оценки равно


MSˆ xx ()  S xx ()  2L (1  )Rxx ()e  j d  S xx ()   см ,
L
т. е. она является смещенной и смещение
Если предположить, что
 см
тем больше, чем меньше L.
(r )
Sˆ xx
() некоррелированы, то дисперсию оцен-
ки СПМ можно найти по формуле
 2Sˆ ()  (1/ m 2 )r 1 2S () ,
m
r
xx
причем
2 () ,
 2S ()  S xx
r
следовательно,
2 () ,
2Sˆ ()  (1/ M )S xx
(15)
xx
lim  2Sˆ ()  0 ,
n
xx
т. е. оценка СПМ состоятельная – чем больше сегментов M, тем меньше
дисперсия, но тем больше смещение  см .
Отметим, что уменьшение дисперсии оценки СПМ обратно пропорционально числу сегментов. Уменьшение дисперсии с увеличением числа
сегментов будет несколько меньше в том случае, когда периодограммы
отдельных сегментов оказываются статистически зависимыми.
Модифицированная периодограмма для реализации СВР. В этом
случае оценка СПМ находится следующим образом. Для каждого сегмента находим коэффициенты дискретного преобразования Фурье (ДПФ)
с помощью быстрого преобразования Фурье (БПФ)
L 1
X r (k )   n 1 xr [n]e  j 2nk , k  0,, L  1 , r  1,, M .
(16)
Величина I r ( f k ) , называемая периодограммой, вычисляется по формуле
I r ( f k ) | X r (k ) |2 ,
(17)
где f k  k / L  частоты ДПФ. Оценка СПМ находится по формуле
1
M
(18)
Sˆ xx ( f k ) 
I ( fk ) .

r 1 r
M
Для уменьшения дисперсии этой оценки СПМ нужно увеличить количество сегментов, не уменьшая длины каждого сегмента. Такое увеличение числа сегментов можно осуществить следующим образом. Исходную последовательность длины N разбивают на M сегментов x r [n ] дли-
ной по L отсчетов, сдвинутых относительно друг друга на D отсчетов,
например, на рис. 1 D  L / 3 .
Сегменты x r [n ] связаны с исходной последовательностью x[n] соотношением
x[n]
D
D
x1[ n]
D
D
x2 [n]
D
n
x3[n]
Рис.1. Сегментация последовательности данных
xr [n]  x[n  (r  1) D],
r  1, , M .
(19)
Максимальное число сегментов будет определяться целой частью
числа M  ( N  L) / D  1. Количество сегментов выбирается с учетом
требуемой гладкости оценки СПМ и требуемой разрешающей способности по частоте. При малом L получается больше сегментов, что обеспечивает оценку с меньшей дисперсией, но с большим смещением и с
меньшим разрешением.
Окна и их основные свойства. Дальнейшее улучшение свойств
оценки возможно применением функции окна w[n] (windowing). Дискретное преобразование Фурье W ( f )  W (e j 2f ) от окна w[n] называется частотной характеристикой окна.
Рассмотрим кратко влияние функции окна на свойства оценки СПМ.
Пусть дана функция временного окна
w[n], 0  n  L  1,
w[n]  
n  L.
0,
(20)
Соответствующая частотная характеристика окна находится так
L 1
W ( f )   n  0 w(n)e  j 2 fn ,  0,5  f  0,5 .
(21)
Математическое ожидание оценки СПМ для отдельно взятого
сегмента равно [4]
M [ Sˆ xx ( f k )]  (1 / L)
1/ 2
 H ( f )S
f  1 / 2
xx
( fk  f ) ,
(22)
где H ( f ) | W ( f ) |2 . Из формулы (22) следует, что математическое ожидание искомой оценки равно свертке истинной СПМ с квадратом модуля
Фурье-преобразования от функции окна.
Рассмотрим влияние вида окна на смещение, разрешающую способность и дисперсию оценки СПМ. Формулы наиболее популярных
w[n]
1
1
2
0,5
3
4
0
L/2
n
L
Рис.2. Окна: 1 – прямоугольное, 2 – Бартлетта, 3 – Ханна, 4 – Хемминга.
окон приведены в [4], табл.18.2, а их графики во временной и частотной областях на рис. 2 и рис. 3.
1
|W(f)|
1
L=1
4
0,5
3
2
2
–3
–2
–1
f, гц
0
1
2
3
Рис.3. Спектральные окна: 1- прямоугольное, 2 – Бартлетта, 3 – Ханна, 4 - Хемминга
Прямоугольное окно. Первый боковой лепесток частотной характеритики окна (рис. 3, кривая 1) ИМЕЕТ экстремум в точке 3π/2 и составляет 2/3π величины главного максимума. Если имеется крутой
максимум СПМ на частоте ω, то в оценке появляются ложные максимумы на расстоянии  3 / 2 L,
 5 / 2L, от ω из-за совпадения
боковых экстремумов с пиком СПМ. Это может привести к неверным
выводам. Другим недостатком этого окна является его знакопеременность (первый боковой лепесток отрицательный), что приводит к отрицательным значениям СПМ. По указанным причинам прямоугольное окно применяется редко.
Ширина главного лепестка спектрального окна определяет минимальную ширину СПМ, взвешенной временным окном последовательности. Поскольку ДПФ является периодической функцией, то
наложение боковых лепестков может увеличить смещение и маскировать присутствие компонентов малой мощности. Снижение максимумов боковых лепестков окон Бартлетта и Ханна способствует обнаружению слабых компонент СПМ исследуемого сигнала. Если СПМ
содержит много мелких деталей, то для разрешения их нужно увели-
чивать длину сегментов при выбранном окне. Под разрешением понимается способность различать два локальных максимума в СПМ.
Формальное определение разрешающей способности – насколько
близкими должны быть две близко расположенные по частоте синусоиды, чтобы их спектральные отклики стали неразличимыми [4].
Разрешение примерно равно величине обратной ширине временного
окна.
Уменьшение боковых лепестков уменьшает смещение, однако
при этом увеличивается ширина главного лепестка, что ухудшает частотное разрешение. Нужен компромисс между шириной главного
лепестка и уровнем максимума бокового лепестка. Например, прямоугольное окно имеет самую малую ширину главного лепестка, но
очень большие медленно убывающие знакопеременные экстремумы
боковых лепестков. Частотные характеристики окон Бартлетта, Ханна
всюду положительны. Окно Хемминга не всюду положительно, однако первый боковой экстремум, имеющий отрицательное значение,
составляет лишь около 0,021 от высоты главного лепестка. Именно
это свойство оправдывает применение окна Хемминга.
Пусть выбран тип окна и его ширина. Для нахождения оценки
СПМ по периодограмме нужно выполнить следующие действия.
Каждый сегмент умножается на функцию этого окна и определяются
коэффициенты БПФ по формуле
L 1
X r [k ]   n  0 xr [n]w[n]e  jkn .
(23)
Периодограмма вычисляется по формуле
1
| X r [k ] |2 ,
U
L 1
U   n  0 w 2 [n] – энергия окна.
Ir ( fk ) 
где
(24)
(25)
Оценка СПМ находится как среднее значение периодограмм отдельных сегментов
1
Sˆ xx ( f k ) 
M
1
 r 1 I r ( f k )  MU  r 1| X r [k ] |2 .
M
M
(26)
Точность оценивания СПМ измеряется дисперсией. Если исследуемая последовательность нормально распределена, а СПМ является
достаточно гладкой функцией на ширине спектрального окна, то дисперсия оценки описывается выражением [27]
[ S ( f )]
D[ Sˆxx ( f k )]  xx k
M
где
2
M 1
M 1


1

2
 (v )  ,


M
v 1


{kL10 w[k ]w[k  vD]}2 .
 (v ) 
{kL10 w2 [k ]}2
(27)
(28)
2.2. Коррелограммная оценка СПМ
СПМ представляет собой дискретно-временное преобразование
Фурье корреляционной функции:

S xx ()   k   Rxx [k ] exp{  jk},   /      /  .
(29)
Оценка СПМ по конечной последовательности значений Rˆ xx [k ] равна
L 1
Sˆ xx ()   k  ( L 1) Rˆ xx [k ] exp{ jk},   /      / .
(30)
Оценка (30) является смещённой и несостоятельной. Для нахождения состоятельной оценки СПМ по оценке АКП наиболее часто употребляется
метод корреляционных окон. Предположим, что мы ограничились L<N–1
значениями АКП, где L называется точкой усечения. При увеличении L
растёт дисперсия оценки СПМ, а при уменьшении L увеличивается смещение оценки СПМ. Это противоречие – возможность уменьшения смещения только ценой большей дисперсии – занимает центральное место в
оценивании СПМ и, пожалуй, всё, что можно сделать на практике, – это
выбрать несколько различных значений L и для каждого из них вычислить СПМ. Из полученного множества выбирается та оценка, которая
наиболее удовлетворяет физической сущности изучаемого процесса.
Ограничение L равносильно умножению оценки АКП на весовую
функцию вида:
1 при 0  k  L,
w[k ]  
0 при k  L,
(31)
называемую прямоугольным корреляционным окном. Тогда оценка СПМ
примет вид:
L 1
Sˆ xx ()   N  ( L 1) w[k ]Rˆ xx [k ] exp{ jk},   /     /  . (32)
Недостатком прямоугольного окна является то, что значения оценки
СПМ при больших  могут быть отрицательными, т. е. не удовлетворять
физической сущности СПМ. Лучших результатов можно достичь, используя окно, согласованное с видом АКП. В качестве таких окон используются окна Бартлетта, Тьюки-Хеннинга, Парзена, Даниэля и др.
Наиболее часто используется простое треугольное окно Бартлетта, имеющее вид:
1 | k | / L, | k | L,
w[k ]  
| k | L.
0,
(33)
Так как АКП − чётная функция, то формулу (32) можно переписать
в виде:
L 1
Sˆ xx ()  2{Rxx [0]  2 k 1 Rˆ xx [k]w[k] cos( k)}, 0     / . (34)
В расчётах удобно полагать   1 и вычислять нормированную
оценку СПМ S xx ()  Sˆ xx () / Rˆ xx (0). Тогда
L 1
S xx ()  2[1  2 k 1 rˆx [k ]w[k ] cos( k )], 0    ,
(35)
где rˆxx [ k ] – нормированная оценка АКП.
Для численных расчётов нужно задать шаг дискретности: , т. е.
(36)
  l , l  0, 1, , F .
Так как 0     , то удобно  выбрать кратным  :
 =  /F,
(37)
где F = (2  3) L. С учётом введённых преобразований расчётная формула имеет вид:
L 1
Sˆ xx (l )  2[1  2  rˆxx [k ]W [k ] cos(lk / F )], l  0, 1, , F .
(38)
k 1
После расчёта Sˆ xx (l ) можно перейти к истинным масштабам частот и
значений СПМ:
Sˆ xx ()  Sˆ xx (l ) Rˆ xx (0),   l / F .
(39)
Рекомендуется оценку СПМ находить для трех значений параметра окна L:
N k / 2, N k , 2 N k ,
где N k равно числу сдвигов, соответствующему интервалу корреляции
случайного временного ряда.
Дисперсия оценок Sˆ xx () . Поскольку оценки с корреляционными окнами схожи с оценкой по множеству отрезков, то и диспер-
сии должны быть одинаковыми по структуре. Можно доказать, что
[5]
DSˆ xx () 
2 ()
2
S xx
L 1
2 [k ]  S xx () I ,
w

k  ( L 1)
N
N
(40)
где
L 1
I   k  ( L 1) w2 [k ]
(41)
энергия окна. Значения I, выраженные через параметр окна L, приведены в
[4], табл.18.3. Например, для окна Бартлетта, величина I равна I  (2 / 3) L .
Поэтому
DSˆ xx () 
2 ()
S xx
2
 L. ,
N
3
2.3. Интервальная оценка СПМ
В случае дискретного белого шума в [4] показано, (формула (18.69)),
что величина 2Sˆ xx () / S xx () имеет  2 -распределение с двумя степенями
свободы ( v  2 ). В случае корреляционных окон можно показать, что величина v Sˆ xx () / S xx () распределена как  2 с v степенями свободы,
v  2N / I .
где
Для
2
(42)
интервальная оценка СПМ с доверительной вероятностью
  1   (   уровень значимости) определяется по формуле


vSˆ ()
  2 / 2 (v)  xx
 12 ( / 2) (v)   .
S
(

)
x


(43)
Отсюда получаем
 v Sˆ ()
v Sˆ () 
 2 xx
 S xx ()  xx
  ,
2


/2 
 1 (  / 2 )
т. е. интервальная оценка СПМ на уровне значимости
выражением
 v Sˆ xx () v Sˆ xx () 

.
,
2
 2


1

(

/
2
)

/
2


 определяется
(44)
Чем больше число степеней свободы v , тем надёжнее оценка в
смысле уменьшения дисперсии, но нужен разумный компромисс между
числом степеней свободы
v
и смещением
 см .
Если число степеней свободы  превышает 20, то интервальную
оценку СПМ можно вычислять по приближенной формуле, считая, что
оценка СПМ распределена нормально с параметрами ( Sˆ xx (),  2 ) ,
Sˆ xx
где величина
 2Sˆ
xx
оценивается по формуле (40). Тогда приближенная
интервальная оценка будет определяться так
lS xx  {Sˆ xx ()  t S Sˆ ()}.
(45)
xx
ТРЕБОВАНИЯ К ОТЧЕТУ
Отчет должен содержать следующие разделы:
1) название и цель работы;
2) номер варианта задания и исходные данные;
3) краткие ответы на контрольные вопросы;
4) таблицы значений теоретических, точечных и интервальных
оценок АКП, СПМ и их графики;
5) выводы по работе.
1.
2.
КОНТРОЛЬНЫЕ ВОПРОСЫ
Объясните физический смысл АКП и СПМ.
Докажите, что периодограммная оценка СПМ
N 1
Sˆ xx ( f )  (1 / N )n0 x[n] exp( j 2fn)
2
идентична коррелограммной оценке СПМ
N 1
Sˆ xx ( f )  k ( N 1) Rˆ xx [k ] exp( j 2fk ) ,
если Rˆ xx [k ] смещенная АКП при максимальном числе возможных сдвигов.
3. Почему распределение оценки АКП нормального случайного
процесса отлично от нормального?
4. Выполняются ли условия эргодичности стационарного случайного процесса по отношению к АКП и СПМ?
ЛИТЕРАТУРА
1. Волгин В. В., Каримов Р. Н. Оценка корреляционных функций
в промышленных системах управления. – М.: Энергия, 1979. – 80 с.
2. Каримов Р. Н. Теоретические основы кибернетики / (СПИ). −
Саратов, 1975. – 115 с.
3. Каримов Р. Н. Статистика нестационарных случайных процессов в АСУ / (СПИ). − Саратов, 1986. – 80 с. (СПИ).
4. Каримов Р .Н. Обработка экспериментальной информации.:
Учеб. пособие. Ч.4. Анализ случайных процессов / СГТУ. Саратов, 2001.
– 104 с.
5. Дженкинс Г., Ваттс Д. Спектральный анализ и его приложения. – М.: Мир, 1971, вып. 1, – 316 с., 1972, вып. 2. – 288 с.
6. Марпл.-мл. С. Л. Цифровой спектральный анализ и его приложения. – М.: Мир, 1990. – 584.
1
1.1.
1.2.
2.
2.1.
2.2.
2.3.
СОДЕРЖАНИЕ
Введение ………………………………………………...….
Задание ………………………………….........................….
Теоретические сведения .....................…………………….
Оценка автокорреляционной функции……..………….…..
Точечная оценка АКП ..................................…………...….
Интервальная оценка АКП …...........…………………...…
Оценка спектральной плотности мощности ........................
Периодограммная оценка СПМ…………………………….
Коррелограмма оценка СПМ ................................................
Интервальная оценка СПМ …...............................................
Требования к отчету ...........................................................
Контрольные вопросы ..........................................................
Литература .............................................................................
Стр.
3
3
4
4
4
6
7
7
12
14
15
15
16
Составитель Равиль Нургалиевич Каримов
КЛАССИЧЕСКИЕ МЕТОДЫ
КОРРЕЛЯЦИОННО-СПЕКТРАЛЬНОГО
АНАЛИЗА ВРЕМЕННЫХ РЯДОВ
Методические указания к выполнению
лабораторной работы
Компьютерный набор и редакция автора
Темплан 2004 г., поз. № 207. Подписано в печать 10. 09. 2004 г.
Формат 60  84, 1/16. Бумага потребительская. Усл. печ. л. 1,06.
Усл. авт. л. 0,88. Тираж 60 экз. Заказ
Волгоградский государственный технический университет.
400131 Волгоград, просп. им. В.И. Ленина, 28.
РПК «Политехник»
Волгоградского государственного технического университета.
400131 Волгоград, ул. Советская, 35.
Отпечатано в типографии «Новый ветер», ПБОЮЛ Выдолоб Л. Ф.
Волгоградская обл., г. Камышин, ул. Ленина, 8/1.
Download