CC BY
44
9. Чебраков Ю. В. Системный подход к прецизионным исследованиям. СПб: Изд-во Санкт-Петербург. ун-та, 1992. 167 с.
10. Шмелев В. Н. Системный подход при научном обосновании требований, включаемых в ТТЗ на разработку новых образцов бортовой радиоэлектронной аппаратуры летательных аппаратов // Мир авионики. 2000. № 3. С. 35—36.
11. Губанов В. С. Обобщенный метод наименьших квадратов. Теория и применение в астронометрии. СПб: Наука, 1997. 318 с.
Павел Павлович Парамонов Юрий Арменакович Гатчин
Борис Викторович Видин Игорь Олегович Жаринов Олег Олегович Жаринов
Сведения об авторах
д-р техн. наук, профессор; ОКБ „Электроавтоматика", Санкт-Петербург; директор; E-mail: postmaster@elavt.spb.ru
д-р техн. наук, профессор; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра проектирования компьютерных систем; E-mail: gatchin@ifmo.ru канд. техн. наук, профессор; ОКБ „Электроавтоматика", Санкт-Петербург; зам. главного конструктора; E-mail: postmaster@elavt.spb.ru канд. техн. наук, доцент; ОКБ „Электроавтоматика", Санкт-Петербург; нач. отдела; E-mail: igor_rabota@pisem.net
канд. техн. наук, доцент; Санкт-Петербургский государственный университет аэрокосмического приборостроения, кафедра вычислительных и электронных систем; E-mail: zharinov@hotbox.ru
Рекомендована кафедрой проектирования компьютерных систем СПбГУ ИТМО
Поступила в редакцию 18.01.10 г.
УДК 519.216
А. В. Денисов, А. П. Филимонов
ИСПОЛЬЗОВАНИЕ ИНФОРМАТИВНЫХ ПРИЗНАКОВ ДЛЯ ВЫЯВЛЕНИЯ МОМЕНТА НАРУШЕНИЯ СТАЦИОНАРНОСТИ СЛУЧАЙНОГО ПРОЦЕССА
Доказывается высокая информативная значимость изменения числа инверсий случайного процесса для выявления момента нарушения стационарного режима работы сложной динамической системы. Рассматриваются два критерия: параметрический, связанный с представлением дисперсий элементарных гармонических составляющих канонического разложения случайного узкополосного процесса в виде координат вектора в гильбертовом пространстве, и непараметрический, позволяющий обнаружить тренд при условии широкой априорной неопределенности процесса.
Ключевые слова: выбросы случайных процессов, оценивание сигналов, фильтрация сигналов, каноническое представление случайного процесса.
Современные технические средства позволяют существенно повысить эффективность измерений и фильтрацию случайных сигналов. Задачи обнаружения трендов и выявления моментов нарушения стационарности случайного процесса на фоне существенных помех часто возникают при измерении характеристик геофизических сигналов для определения границ месторождений полезных ископаемых (особенно с подвижных оснований), при информационно-аналитическом сопровождении геофизических работ, при обнаружении аварийных ситуаций в работе сложной динамической системы (какой является, например, плавучая буровая станция). Для решения подобных задач используется функциональная обработка
получаемой информации и логические фильтры. Подход, связанный с вычислением средне-квадратического отклонения (СКО) и спектральной плотности, недостаточно чувствителен к отдельным выбросам сигнала, которые несут полезную информацию, а распространенные интегральные фильтры порой не обладают требуемым быстродействием.
В настоящей статье предлагается использовать для математической обработки информации метод параметрической статистики, связанный с минимизацией количества членов канонического разложения, а также метод непараметрической статистики. В обоих случаях важно правильно определить вид изучаемого процесса. Использование непараметрических методов [1, 2] упрощает функциональную и структурную схему устройств обнаружения момента нарушения стационарности процесса при сохранении степени достоверности принятия решений.
Рассмотрим один из параметрических критериев, для чего установим связь между среднестатистическим числом инверсий [3, 4] и параметрами канонического представления [5—7] случайного процесса (СП), который будем считать узкополосным. Важность рассмотрения этого случая связана с квазипериодическим изменением некоторых параметров во времени при изучении различных технологических процессов. Исследуем стационарный СП а (^), который будем полагать нормальным и центрированным, а также дифференцируемым. Пусть на конечном интервале времени Дt процесс имеет N пересечений с пороговым значением а^) = 0 . В этом случае N за время М определяется как [3, 4]
^^, (1) П О
где оа и о а — соответственно СКО функции а ^) и ее производной а = —.
Л
Из формулы (1) следует, что хотя N и зависит от оа, но информативным признаком стационарности процесса является отношение указанных СКО. Процесс а(1.) также является узкополосным.
2
Сформируем для процесса а^) с нулевым средним значением и дисперсией оа его ка-
( Т Т
ноническое представление [5] в виде тригонометрического ряда на интервале I -—, —
а(г)=2
/=0
(2тс/ ^ „ . (2п/ Л/ соб I-1 1 + В/ Б1п |
(2)
Т ) ' ^ т
где Л/ и В1 — некоррелированные между собой при разных значениях / случайные величины, распределенные по нормальному закону, имеющие нулевые математические ожидания и дисперсии о2, сумма которых равна о'Оа; при выборе периода разложения Т целесообразно использовать оптимальные параметры канонического разложения [7].
В качестве математической модели процесса а(Г) зададим его корреляционную функцию К (т) в следующем виде [7]: на интервале |т|< Т/ 2 представим ее в виде произведения огибающей, характеризующейся параметром О, и гармонического заполнения с частотой в:
ч 2 Б1П От , _
К(т) "От) со*вт,
а при |т|> Т/2 положим К(т) = 0 . С ростом |т| огибающая корреляционной функции впервые обращается в нуль при т = ±тс/О. Интервал между этими двумя точками будем считать периодом Т, на котором процесс а (^ представлен в виде уравнения (2), затем продолжим это
представление на всю временную ось I, при этом функция а^) останется непрерывной. Дисперсии элементарных гармонических составляющих в этом случае определяются формулами
2 Т/2
а2 = - | К(т)ат, (3)
Т 0
2 4 Т К() (2п1 ^ а т п
а, =— КтсоБI-т Iат, Т =— .
I т 0 | т ) а
Выбор указанного периода разложения обусловлен, во-первых, быстрым уменьшением дисперсии высокочастотных составляющих процесса [7], так что в каноническом представлении достаточно учесть только несколько первых гармоник; во-вторых, задание корреляционной функции при всех значениях т с формальной точки зрения требует наличия бесконечной выборки, что практически невозможно.
В соответствии с физической сущностью многих случайных процессов в представлении
СП а($) должна отсутствовать постоянная составляющая а0 . При выборе указанного периода Т величина а2 (см. формулу (3)) отсутствует только в двух случаях: если а/Р~ 0,30 или а/Р~ 0,42. В обоих случаях можно ограничиться первой, второй и третьей гармониками, исключив высокочастотные элементарные процессы, дисперсии которых много меньше дисперсии анализируемого СП:
3
аЦ )
/=1
(2тс/ ) „ . ( 2п/ Л, СОБI 1 + В Б1П \~T~t
(4)
При а/Р- 0,42 получаем: а2 = 0,15аа, а2 = 0,45а2, а2 = 0,30а2.
Поскольку ряд (2) мажорируется сходящимся положительным числовым рядом, членами которого являются дисперсии элементарных составляющих, то его можно почленно проинтегрировать. Интегрируя соотношение (4) и полагая постоянную интегрирования равной нулю, получаем приближенное каноническое представление для а(^), согласно которому находим среднеквадратическое отклонение аа исходного СП а(/):
а — аа а2 + 1 а2 + 1 а2 аа -—, 0] +—а2 +—а3 .
а а V 1 4 2 9 3
С учетом этого соотношения формула (1) примет следующий вид:
N 2а.
П
При рассмотрении узкополосного случайного процесса изменение значения а можно использовать как критерий, характеризующий нарушение стационарного режима СП. Целесообразно, однако, в качестве информативного признака рассматривать дисперсии элементарных процессов, которые можно представить в виде координат вектора в конечномерном гильбертовом пространстве. Этот вектор полностью характеризует СП. Затем можно разделить данное пространство на множества (классы) таким образом, чтобы расстояние между ними было наибольшим, а время, необходимое для идентификации принадлежности вектора к какому-либо классу, было наименьшим.
В случае когда априорная информация о процессе отсутствует или график спектральной плотности характеризуется наличием ряда максимумов на близкорасположенных частотах, для выявления момента нарушения стационарного режима целесообразно использовать метод нелинейной фильтрации сигналов, связанный с суммированием отклонений N от нулевого значения с учетом весовых коэффициентов.
Рассмотрим этот вопрос с общих позиций. Для обнаружения тренда следует использовать непараметрическую форму задания /-критерия в виде
где £ — вектор-столбец информационных сигналов, полученных по т каналам; Ь — вектор-столбец, определяемый как стартовая пороговая точка в да-мерном пространстве; с — коэффициент, характеризующий доверительную вероятность обнаружения тренда ( он может, например, соответствовать ряду коэффициентов Стьюдента для нормального распределения СП); Е — функциональная квадратная матрица ( т х т ) непараметрического преобразования сигналов; — квадратные матрицы ( т х т ) для нелинейной обработки входного (векторного) сигнала а^); к1 — весовые коэффициенты, зависящие от стохастических характеристик сигнала а^) и измерительных каналов, с помощью коэффициентов к1 можно повысить чувствительность критерия в каждом конкретном случае, не изменяя непосредственно матрицы и время принятия решения.
На основе априорной информации строятся далее вектор-столбцы /т1п, /тах, определяющие граничные значения для каждого информационного параметра. Разность /тах- /т1п формирует новый вектор-столбец, по которому можно судить о наличии тренда и о нарушении стационарности режима работы системы.
Проверка этого критерия (в более упрощенном виде — для одного измерительного канала) была осуществлена авторами в НЦ „Вагоны" Петербургского государственного университета путей сообщения. Для обнаружения схода колесной пары грузового вагона были проведены компьютерное моделирование и натурный эксперимент. Установлено, что „интегральные критерии" (СКО и основанные на вычислении спектральной плотности) в этом случае неэффективны, а использование непараметрического критерия, подобного приведенному в формуле (5), позволяет зафиксировать время схода вагона при максимальной скорости его движения и при любой шпальной неровности не более чем за 2 с. Эта задача была решена в аналого-цифровом варианте с использованием микропроцессорного устройства, при этом предварительная фильтрация сигнала выполнялась посредством аналоговой части схемы. Был реализован алгоритм, автоматически обнаруживающий и выделяющий тренд в темпе поступления данных, т.е. в реальном масштабе времени. В простейшем случае одноканальной системы графики функций 3т1п(N) и 3тах(N) разделяют пространство (/, N) на область /т1п(N) < /(^ < 3тах(N), в которой процесс стационарен, и две области 3(N) > 3тах(N), /тт(N) > 3(N), в которых процесс становится нестационарным с возрастающим или убывающим трендом.
Таким образом, непараметрическая форма оценки реализаций СП обладает простотой использования и имеет широкую область применения, что существенно улучшает эффективность решения задачи анализа случайных процессов.
1. Дмитриев Ю. Г., Кошкин Г. М. Непараметрическое оценивание сигналов. М.: Наука, 1997. 336 с.
2. Тарасенко Ф. П. Непараметрическая статистика. Томск: Изд-во Томск. ун-та, 1976. 292 с.
3. Светлицкий В. А. Случайные колебания механических систем. М.: Машиностроение, 1976. 216 с.
4. Тихонов В. И. Выбросы случайных процессов // Успехи физ. наук. 1962. Т. 77, вып. 3. С. 449—480
5. Пугачёв В. С. Теория случайных функций и ее применение. М.: Физматгиз, 1962. 720 с.
3 = ±Ь + 2кБ +сЕ £,
(5)
V г=1
)
СПИСОК ЛИТЕРАТУРЫ
6. Райс С. Теория флуктуационных шумов // Теория передачи электрических сигналов при наличии помех: Сб. переводов / Под ред. Н. А. Железнова. М.: Изд-во иностр. лит., 1953. С. 88—238.
7. Филимонов А. П., Денисов А. В. Выбор параметров канонического разложения узкополосного случайного процесса // Изв. вузов. Приборостроение. 2007. Т. 50, № 8. С. 48—53.
Александр Владимирович Денисов
Анатолий Павлович Филимонов —
Сведения об авторах канд. физ-мат. наук, доцент; Санкт-Петербургский государственный горный институт им. Г. В. Плеханова (ТУ), кафедра высшей математики; E-mail: rectorat@spmi.ru
канд. техн. наук, доцент; Петербургский государственный университет путей сообщения, кафедра физики; E-mail: pgups200@rambler.ru
Рекомендована кафедрой высшей математики СПбГГИ (ТУ)
Поступила в редакцию 23.10.09 г.
