9. ГЛАВНЫЕ КОМПОНЕНТЫ В ЧАСТОТНОЙ ОБЛАСТИ
9.1. Введение
В предыдущей главе была рассмотрена задача аппроксимации стационарного ряда при помощи пропущенного через линейный фильтр другого стационарного ряда. В этой главе мы изучим проблему аппроксимации ряда посредством фильтрации не другого, а того же самого ряда; при этом ограничимся фильтрами, имеющими ранг меньший, чем число компонент ряда.
Уточним постановку задачи. Пусть рассматриваемый векторный ряд 
с 
компонентами имеет среднее
(9.1.1)
абсолютно суммируемую автоковариационную функцию
(9.1.2)
и матрицу спектральной плотности
(9.1.3)
Допустим, необходимо передать по каналам связи из одного пункта в другой значения величин 
, но при этом в нашем распоряжении только 
каналов. Предположим, что процесс передачи ряда 
по q каналам связи можно описать как фильтрацию, в результате которой получается 
-компонентный векторный ряд
(9.1.4)
(здесь 
есть 
-матричный фильтр), и собственно передачу. Пусть на выходе каналов связи принимается ряд, который выступает в качестве оценки для 
(9.1.5)
где 
-вектор с 
компонентами, а 
есть 
-фильтр. Мы ставим себе целью так выбрать 
и фильтры 
чтобы ряд 
был близок к 
Связь между 
линейна и инвариантна во времени, она описывается передаточной функцией
(9.1.6)
если 
обозначают соответственно передаточные функции фильтров 
Теперь ясно, что поставленная задача сводится к определению г 
-матрицы А (К), ранг которой меньше 
, так, чтобы величина разности
(9.1.7)
оказалась мала.
Можно было бы рассматривать эту задачу и как отыскание способа построить такой 
-мерный ряд 
который несет значительную часть информации об исходном ряде X (t). Отметим здесь, что Bowley (1920) однажды сделал такое замечание: «Показатели применяются для того, чтобы измерить изменения величины, не доступной прямому наблюдению, но о которой нам известно, что она оказывает определенное воздействие на другие величины, непосредственно наблюдаемые, причем хотя само это воздействие вызывает либо одновременное увеличение, либо уменьшение всех наблюдаемых величин, его эффект маскируется действием многих других причин, по-разному влияющих на отдельные наблюдаемые». Возможно, что введенный ряд 
будет выступать в роли ряда показателей, описывающих воздействие на 
со стороны некоторого скрытого от наблюдателя ряда. Вводя ряд 
условимся считать, что он наилучшим среди 
-мерных рядов образом позволяет восстановить 
с помощью линейных операций, инвариантных во времени.
Правомерна и другая точка зрения. Введем ряд 
описывающий ошибку или искажение
(9.1.8)
тогда получим
(9.1.9)
Ряд 
оказывается выраженным как профильтрованный вариант ряда 
имеющего меньшую размерность, плюсошибка. Приведем пример ситуации, в которой полезна такая модель: пусть 
представляет собой ряд импульсов, вызванных q землетрясениями, одновременно происходящими в разных местах; пусть 
представляет сигналы, принятые 
сейсмографами, и
пусть с 
описывает переходные явления в земной коре, связанные с землетрясениями. Сейсмологов интересуют свойства ряда 
см., например, Ricker (1940) и Robinson (1967b).
Существенным моментом решения этих задач является аппроксимация изучаемого ряда другим, имеющим меньшую размерность. В § 9.2 мы приводим обзор некоторых аспектов классического анализа главных компонент векторных переменных.
