Можно показать, что
где
Чтобы найти СПМ выходного АР процесса необходимо в (13) сделать замену
Тогда из (13) следует выражение для параметрической оценки СПМ
Выражение (14) широко используется в параметрическом методе спектрального оценивания.
В качестве параметров, полностью характеризующих спектральную оценку случайного процесса, выступают коэффициенты АР и порядок модели.
Параметрическое спектральное оценивание обладает рядом преимуществ по сравнению с традиционными методами спектрального оценивания. К ним относятся: более высокое спектральное разрешение при использовании коротких выборок, отсутствие боковых лепестков.
С помощью модели АР можно получать спектральные оценки случайных процессов со сложной формой СПМ.
Для этого может быть придется использовать модели АР большого порядка. На основе модели АР легко синтезируются оптимальные фильтры подавления, согласованные не только по частоте и полосе спектра, но и по форме спектра случайного процесса.
Достоинством формулы (14) является возможность анализировать СПМ в аналитическом виде, что невозможно сделать при использовании традиционных методов спектрального оценивания на основе преобразования Фурье.
Например, можно найти формулы для определения частоты максимумов и минимумов СПМ.
Чтобы определить положение максимума или минимума АР оценки СПМ, нужно взять производную от (14) по
При
где
3. Характеристическое уравнение модели авторегрессии
Модель АР, описываемая уравнением (1), может быть представлена в операторной форме
где оператор АР
Действие оператора сдвига z на текущий отсчет описывается следующим образом
Из условия устойчивости формирующего АР фильтра с рациональной передаточной функцией (11), следует условие стационарности АР процесса. Для проверки стационарности случайного АР процесса используется характеристическое уравнение
Если корни характеристического уравнения (19) лежат внутри единичного круга на комплексной плоскости, то процесс АР удовлетворяет условию стационарности и его корреляционная функция стационарна. Характеристическое уравнение (19) можно представить также в виде
Тогда условие стационарности заключается в том, что корни характеристического уравнения (20)
Используя (19) или (20) оператор АР (17) можно представить в виде
Из (21) следует, что уравнение АР (1) можно записать следующим образом
Сравнивая (1) и (22) найдем связь между коэффициентами АР и корнями
где первый индекс в квадратных скобках указывает на соответствующий порядок модели.
Полученные формулы оказываются весьма полезными для определения коэффициентов АР по заданным характеристикам случайного процесса.
Отметим, что корни характеристического уравнения полностью описывают модель АР.
Свойства модели зависят параметров, через которые они выражаются. Если корень действительный, то его можно представить в виде экспоненциальной функции
где
Тогда действительные корни характеристического уравнения принимают вид
Комплексные корни характеристического уравнения описываются выражениями
где
4. Генерация коррелированного случайного процесса
В задачах статистического моделирования часто возникает необходимость генерации случайного процесса с заданной корреляционной функцией или с заданной формой и характеристиками СПМ. Для этих целей эффективно использовать генератор процесса АР, показанный на рис. 2.
Генерация случайного процесса осуществляется методом порождающего случайного процесса.
Порождающий процесс в виде белого шума, обычно с гауссовой функцией распределения, пропускается через формирующий фильтр, параметры которого определяются соответствующей моделью АР.
Рисунок 2. Генератор процесса АР
Для генерации процесса нужно выбрать необходимое количество пиков СПМ. Тогда порядок модели АР равен удвоенному числу пиков. Так, для СПМ с одним пиком на ненулевой частоте, порядок модели равен 2. Для СПМ с двумя пиками порядок модели равен 4.
Затем выбирают частоту пика и его ширину полосы. Вычисленные значения корней характеристического уравнения по формулам (24), используются для нахождения коэффициентов АР.
Для этого корни подставляются в соответствующие выражения (23). Генерация процесса осуществляется с помощью рекуррентного выражения (1) с использованием порождающего белого шума a[t].
Функция распределения a[t] может быть любой, но, как правило, используют нормальное распределение с нулевым математическим ожиданием и единичной дисперсией. Белый шум с нормальным распределением получают из белого шума с равномерным распределением.