CC BY
458
УДК 681.5.015 + 519.246.25 ББК 22.17 + 22.18
ЦИФРОВОЕ МОДЕЛИРОВАНИЕ СТАЦИОНАРНЫХ СЛУЧАЙНЫХ ПРОЦЕССОВ С ЗАДАННОЙ КОРРЕЛЯЦИОННОЙ ФУНКЦИЕЙ НА ОСНОВЕ НЕПРЕРЫВНЫХ ДРОБЕЙ
1 2 Карташов В. Я. , Новосельцева М. А.
(ГОУВПО «Кемеровский государственный
университет», Кемерово)
На основе непрерывных дробей разработан алгоритм моделирования стационарного случайного процесса, установлено соответствие между непрерывной и дискретной моделями случайного процесса, обоснован выбор шага дискретизации.
Ключевые слова: стационарный случайный процесс, корреляционная функция, непрерывная дробь, формирующий объект, белый шум, непрерывная передаточная функция, шаг дискретизации, дискретная передаточная функция, стохастическое разностное уравнение.
1. Введение
При исследовании и проектировании систем контроля, диагностики, управления и обучения сложных объектов возникает необходимость исследования их функционирования под влиянием различных факторов. Особую сложность данная задача
1 Владимир Яковлевич Карташов, заведующий кафедрой автоматизации исследований и технической кибернетики ГОУ ВПО «Кемеровский государственный университет», доктор технических наук, профессор (kartash@kemsu.ru).
2 Марина Александровна Новосельцева, доцент кафедры автоматизации исследований и технической кибернетики ГОУ ВПО «Кемеровский государственный университет», кандидат технических наук, доцент (aanov@pochta.ru).
приобретает при наличии воздействия на объект случайных факторов. При этом проведение экспериментальных и теоретических исследований объекта зачастую затруднительно ввиду их сложности и затратности. Поэтому особое значение и актуальность приобретают работы, связанные с построением, использованием и совершенствованием моделей. Моделирование, являясь одним из инструментов научного познания, нашло широкое распространение в различных сферах жизнедеятельности.
При решении задач имитационного моделирования часто возникает необходимость в формировании дискретных последовательностей стационарных случайных процессов с заданным видом корреляционной функции. Обычно в таких ситуациях закон распределения случайного процесса не принимается во внимание [3, 17, 20].
2. Постановка задачи
Задача цифрового моделирования стационарного случайного процесса у(() с непрерывной корреляционной функцией Ryy(t) заключается в воспроизведении на ЭВМ дискретной последовательности стационарного случайного процесса у^Д() с заданной корреляционной функцией, где k = 0, ..., N Дt - шаг дискретизации [3, 5, 17, 20]. Таким образом, при имитации стационарного случайного процесса осуществляется переход от его непрерывной модели в форме корреляционной функции к дискретной, позволяющей сформировать реализацию у(М().
Точностный анализ полученной модели случайного процесса осуществляется на основе критерия цифрового моделирования [20]. В качестве критерия в данной работе будет использоваться дисперсия расхождений между модельной корреляционной функцией Яуу ), рассчитанной по получен-
ной реализации у^ДО, и истинной корреляционной функцией Ryy(t), взятой с шагом дискретизации Дt [20]:
1 m
(1) De = - X ^) -Me )2,
mk=0
где
е^ ) = Ям (^ )- Яуу (^ ),
1 т
ме = - X e(k),
mk=0
т - число значений корреляционных функций. При этом шаг дискретизации выбирается из некоторого множества допустимых шагов дискретизации Дt е ДТ.
Данная задача цифрового моделирования обычно сводится к определению характеристик формирующего фильтра (линейного преобразователя) при задаваемых характеристиках входного и выходного сигналов (рис. 1) [3, 7, 17, 20, 28]. Если линейный динамический объект с передаточной функцией G(s) находится под воздействием непрерывного стационарного случайного сигнала х(0, имеющего корреляционную функцию Я^), то на выходе данного объекта в установившемся режиме формируется непрерывный стационарный случайный сигнал у(0 с корреляционной функцией Я^) [17]. Связь между случайными процессами х(0 и у(0 выражается через передаточную функцию G(s) формирующего фильтра [5, 17].
ООО
у©
Рис. 1. Формирующий фильтр
Известно, что спектральная плотность выходного сигнала Syy(<ю) формирующего фильтра определяется в соответствии с
выражением [1, 3, 17, 20, 28]
і і2
(2) ^(«) = |£(]ш)\2 («),
где 8хх(ю) - спектральная плотность входного сигнала;
I |2
\0(]а)\ - квадрат модуля амплитудно-фазовой частотной
характеристики системы с непрерывной передаточной функцией (НПФ) G(s).
В качестве входного сигнала обычно используют белый шум - стационарный случайный процесс с корреляционной функцией До;(0 = ^(0. Если входное воздействие объекта -белый шум, то спектральная плотность реакции в установившемся режиме имеет вид [17, 21]
1 I |2
(3) 5 («) = - в(]Ы)\2 .
Р
На основе (3) определяется НПФ в(5) формирующего фильтра [5, 17]:
У & _ А^т + Ат -1 +... + Ао
,п і Е> „п-1
(4) в^) = ,
^ ^ Х(*) В^ + Вя_1^‘ +... + В0
где т, п - целые положительные числа (т < п); Ат, ..., А0, Вп, ..., В0 - постоянные коэффициенты; Х(5), Г(^) - преобразование Лапласа входного и выходного сигналов соответственно.
Основная задача цифрового моделирования реакции у(?) формирующего фильтра заключается в переходе от НПФ в(я) к дробно-рациональному выражению дискретной передаточной (ДПФ) в(2, ДО:
т
V а^
І
(5) в( 2, А? ) =
У (г, Аї) _ атг т +... + а1г 1 + а0 _ г_0
Х(2, М) ъп^ +... + V +1 1 + ^
i=1
где У(г, ДО, Х(г, ДО - ^-преобразования числовых последовательностей значений случайных входного и выходного процессов; ат, ..., а0, Ъп, ..., Ъ1 - постоянные коэффициенты.
При переходе от в(я) к в(2, ДО особо актуальной является проблема установления взаимно однозначного соответствия между непрерывной и дискретной моделями. Данная проблема в литературе остается малоисследованной, хотя именно взаимно однозначное соответствие позволяет получать достоверные результаты цифрового моделирования.
Как видно из (5), для получения в(г, ДО необходимо задать шаг дискретизации Д?. Теоретическое обоснование выбора шага дискретизации и исследование его влияния на взаимно одно-
значное соответствие между непрерывной и дискретной моделями являются важными вопросами цифрового моделирования.
Ниже будет приведен обзор существующих методов цифрового моделирования стационарных случайных процессов, т. е. алгоритмов перехода от непрерывной модели формирующего фильтра к его дискретному представлению.
3. Обзор методов цифрового моделирования стационарных случайных процессов с заданной корреляционной функцией
Для осуществления перехода от непрерывной модели корреляционной функции Яуу(?) к дискретной реализации стационарного случайного процесса у(кД?) используется ряд алгоритмов. Численные методы дают методические ошибки при получении реализаций случайных процессов, величина которых определяется выбранным шагом интегрирования [28, 30], поэтому в данной работе они рассматриваться не будут.
Одним из способов моделирования является метод скользящего суммирования [3, 7, 20, 28], в котором представление оператора линейного преобразования формирующего фильтра осуществляется в виде скользящей суммы с весовыми коэффициентами ак
N
(6) у(п) =Х акх(п - к).
к = 0
Существует ряд способов определения ак на основе значений корреляционной функции выходного сигнала, например, с помощью решения нелинейной системы уравнений [3, 20], разложения в ряд Фурье [3], дискретизации интеграла свертки [28], метода факторизации [3] и т.д. Помимо трудностей, возникающих при нахождении коэффициентов ак, данный метод имеет методическую погрешность ввиду ограниченности количества членов ряда N.
Более быстродействующим методом моделирования является рекуррентный алгоритм авторегрессии - скользящего среднего [7, 17, 20, 28]:
(7) У(к) = £ а,Х (к - О -£ ЪгУ(к -О,
г=0 г =1
где аг, Ъг - коэффициенты модели авторегрессии со скользящим средним; га и п - целые положительные числа (т < п). Данный метод не имеет методической погрешности [3, 5], а параметры моделирующего алгоритма выражаются в явном виде через параметры корреляционной функции. Однако существующие методы дискретизации, факторизации и т.п., основанные на применении (7) [3, 5, 18], наиболее эффективны, когда корреляционная функция моделируемого процесса имеет невысокий порядок. В этих случаях алгоритмы просты и не имеют методической погрешности [3, 5], а их параметры аг, Ъг выражаются в явном виде через параметры корреляционной функции. Если вид корреляционной функции усложняется (например, у НПФ формирующего фильтра есть кратные полюса), то рекуррентные алгоритмы становятся приближенными [3, 5].
Методы моделирования с использованием канонических и неканонических представлений [28] основаны на представлении модели случайного процесса в виде детерминированной функции случайных величин. Канонические представления имеют вид:
¥
(8) уи) = £ Укфк ^),
к=1
где ук - некоррелированные величины; фк - детерминированные функции, определяемые по корреляционной теории случайных процессов с помощью известных алгоритмов. Однако при моделировании данным методом в разложениях учитывается конечное число членов, поэтому метод имеет методическую ошибку, зависящую от этого числа.
В неканонических представлениях модель процесса задается в виде нелинейной функции от конечного числа случайных величин. Неканонические представления не обладают асимптотическим свойством, т. е. оценка корреляционной функции по N измерениям не стремится к истинной корреляционной функции с ростом N [28].
Другой подход к моделированию стационарных случайных процессов основан на том, что переменные преобразования Лапласа в НПФ 0(5) (4) и ^-преобразования в ДПФ 0(2, А^ (5) связаны соотношением 2 = е[3, 17, 24, 25]. Нелинейность этого преобразования обусловила многочисленные попытки его приближения при переходе от аналоговых представлений к дискретным. Для осуществления перехода от НПФ к ДПФ используются дробно-рациональные приближения, в частности, прямые и обратные разности Эйлера, билинейное преобразование и их многочисленные модификации [10, 21, 25]. Однако данные аппроксимации преобразования 2 = е5А не обладают свойствами первоначального преобразования, которые связаны с понятиями многозначности, периодичности, однолистности и т. п. [10, 21, 30].
С другой стороны, преобразование 2 = е5А в полной мере отвечает требованиям взаимно однозначного соответствия при отображении непрерывного пространства в дискретное [21]:
1) для сохранения селективных свойств непрерывной системы необходимо, чтобы ось гу из 5-плоскости отображалась в единичную окружность на 2-плоскости;
2) для того чтобы устойчивая непрерывная система переходила в устойчивую дискретную систему, необходимо, чтобы точки из левой полуплоскости 5-плоскости (х < 0) отображались вовнутрь единичной окружности в 2-плоскости (|2| < 1).
Таким образом, использование непосредственно 2-преобразования является более совершенным и рациональным подходом к построению цифровых моделей.
Например, в [5, 17] разработан достаточно эффективный метод моделирования случайных процессов на основе 2-преобразования 2 = е6'А\ который является точным при простом виде НПФ формирующего объекта (знаменатель НПФ - полином не выше 2-й степени). Если же выражение НПФ сложное, то авторы предлагают осуществить переход к дискретной модели при помощи билинейного преобразования 2 (2 -1)
At (2 + 1)
Как уже было отмечено ранее, подстановочные методы не обладают необходимой точностью моделирования и не позволяют устанавливать взаимно однозначное соответствие между непрерывной и дискретной моделями. В [5] авторы приводят пример моделирования стационарного случайного процесса со спектральной плотностью
^уу(^ ) = Т~л с 2 .
р(4 + 5а + ю )
При использовании 2-преобразования (шаг дискретизации равен 0,1) модель случайного процесса на выходе формирующего фильтра имеет вид:
у (к At) = 0,1218 х(к -1) + 1,7235у (к -1) - 0,741у (к - 2) . Билинейное преобразование искажает структуру модели, и конечно-разностное уравнение принимает другой вид: y(kAt) = 0,003х(к) + 0,006 х(к -1) + 0,003х(к - 2) +
+ 1,7235у(к -1) - 0,741у (к - 2).
Таким образом, существующие методы цифрового моделирования стационарных случайных процессов обладают рядом недостатков: отсутствие взаимно однозначного соответствия между непрерывной и получаемой дискретной моделями, наличие погрешности моделирования, ориентация на фиксированный шаг дискретизации, выбор которого теоретически не обоснован и не обеспечивает взаимно однозначного соответствия между непрерывной и дискретной моделями стационарного случайного процесса.
Задачей данной статьи является разработка и исследование способа взаимно однозначного перехода от непрерывной к точной дискретной модели в форме ДПФ, причем предлагаемый способ должен восстанавливать как порядок, так и значения параметров дробно-рационального выражения для дискретной модели формирующего объекта. Под отношением Я «быть точной цифровой моделью объекта» понимается взаимно однозначное соответствие М\ЯМ2 (М\ - объект моделирования, М2 -его цифровая модель), обладающее свойством эквивалентности и сохраняющее динамические свойства объекта [15].
4. Использование теории непрерывных дробей при цифровом моделировании стационарных случайных процессов
Любой непрерывный стационарный процесс с корреляционной функцией Яуу(г) и спектральной плотностью 1 ¥
(10) ^(«) = - {Лу(г)в-шгёг,
р —¥
которая является дробно-рациональной функцией от круговой частоты ш, можно рассматривать как реакцию непрерывного стохастического линейного динамического объекта (далее по тексту - формирующего объекта) в установившемся режиме на входное воздействие х(г) в форме белого шума [17]. Известно, что НПФ такого объекта будет иметь вид [12, 17]
т) ^(„) = Лху (э) = Ат8т + Ат—^т 1 + ... + А0
V ' П / \ ГУ п ГУ П—1 7~ч ’
Лхх (5) Впэ + Вп—1э + ... + А)
где Яхх(э), Яху(э) - преобразование Лапласа корреляционной функции входного сигнала и взаимной корреляционной функции входного и выходного сигналов соответственно.
Если х(г) - белый шум с корреляционной функцией вида
Го, г ф о,
(12) Лхх(г) = 3(г) И ’ 0’
[¥,г = 0;
то выражение (11) упрощается:
О( ч Кху(5) Кху(э) Л ( ч
(13) О(э) =—А—т =—-— = Яу (э).
ф(г)) 1 ^
Используя основное определение (11) и введя 2-преобразование 2 = еэАг, можно оценить ДПФ формирующего объекта [11]:
¥
М) ?.я,<п^'
хуV 5 ’ п=0
(14) О( г, Аг) =
Лхх (г Аг)
хх
п=0
п
= ^ (0) + К*У (А/)2_ + ••• + К*У (пА<)2 ~« + •••
^ (0) + ^ («А/) 2 _ + ••• + ^ («А/) 2 - « + •••’ где Яху(г, Д/), Яхх(2, Д/) - 2-преобразования взаимной корреляционной функции вход-выходного сигналов и корреляционной функции входного сигнала соответственно •
Для определения 2-преобразования ^-функции введена специальная функция, называемая функцией единичного отсчета, которая определяется следующим образом [24]:
* Ц « = 0,
(15) 8[«] = 8 {пМ) = <
[0, п ^ 0^
Таким образом, 2-преобразование функции единичного отсчета равно 1 Тогда выражение (14) приобретает вид
Я„, (0) + Я (А/) 2_ + ••• + Я («А/) 2_п + •••
(16) 0( 2, А/) = ху --------^ \--------------------------
1 + 0 • 2 + ••• + 0 • 2 « + •••
Модель ДПФ формирующего объекта в форме дробнорационального выражения имеет вид:
т
_ _ "V а.2_
(17) 0(2, А/) = ат2 т+ - + а2,+ а =-й_ = Ш. ,
Ъ2-« + ••• + Ъ, 2 +1 - А, _ 0(2)
« 1 1 + ^ Ъ12 ^пК ’
i=1
где Рт(2), 0«(2) - полиномы от комплексной переменной 2-преобразования; т, п - целые положительные числа - порядки этих полиномов •
Известно также другое определение ДПФ [6, 12] в виде
Аї ^ у(пА )
(18) 0( г, Аї) =...........
У (2Аї) = п=о
г п
X (2,Аї) аї £ х(пАї уп
п=0
У (0) + У (Аї) 2 - +... + у(пАї) 2~п +... х (0) + х (пАї) 2 - +... + х (пАї) 2~п +...
где Х(г, ДО, У(г, ДО - г-преобразования числовых последовательностей значений случайного входного {х(пД0} ¥=0 и случайного выходного {у(п Д )}¥ 0 сигналов.
Благодаря операторным свойствам ^-преобразования дискретную модель формирующего объекта можно представить как в форме стохастического разностного уравнения
т п
(19) у(Ш) = £ агХ((к - О Д^) ЪгУ ((к - 1)Д) ,
г=0 г=1
относящегося к классу процессов авторегрессии со скользящим средним, так и в форме детерминированного разностного уравнения
тп
(20) Ку (к Д) = £ агКхх ((к - г)Д) -£ ЪгКху ({к - / )Д/).
г =0 г =1
Выражение (19), полученное на основе корреляционной функции Яуу((), позволяет осуществить цифровое моделирование случайного процесса у(0, являющегося откликом формирующего объекта на воздействие х(0 в виде белого шума с математическим ожиданием Мх = 0 и дисперсией Бх = а2х . Таким образом, для получения конечно-разностного уравнения (19) необходимо найти ДПФ формирующего объекта 0(г, Д^. В качестве аппарата для нахождения 0(г, ДО, т. е. структурнопараметрической идентификации (БР-идентификации) формирующего объекта [11-13], будем использовать теорию непрерывных дробей [11].
Все имеющиеся способы разложения в непрерывную дробь основаны на представлении аналитической функции в виде степенного ряда либо отношения степенных рядов [11]. Поэтому определение (14) ДПФ как отношения двух степенных рядов от переменной ^-преобразования, а также свойства непрерывных дробей приводят к возможности их применения для аппроксимации ДПФ линейного стохастического объекта.
На основе заданной корреляционной функции Яуу(() стационарного случайного процесса у(0 осуществляется переход к НПФ формирующего объекта [5, 17]. Полученная таким обра-
зом НПФ будет полностью совпадать с Дху(я) - преобразованием Лапласа взаимной корреляционной функции вход-выходного сигналов формирующего объекта. Взяв обратное преобразование Лапласа от О(я), получим весовую функцию й(?) формирующего объекта.
Наиболее приемлемым и простым способом перехода к непрерывной дроби является модифицированный алгоритм В. Висковатова [10].
Зададим шаг дискретизации А?. Будем считать, что Н(0) Ф 0. На основании (16) определяется идентифицирующая матрица:
(-1) - строка (
(21)
1 0 0 . .0
И (0) И(Аї) И(2Аї) . . И(пАї)
а1(0) а1 (Аї) а1 (2 Аї) . . а1 (пАї)
а2 (0 ) а2(Аї) а2 (2Аї) . . а2 (пАї )
ат (0) ат (А ї ) ат (2Аї) . а п .. ї
(0) - сторока 1 - сторока
т - сторока
в которой (-1)-строка представляет собой импульсную функцию [24], (0)-строка содержит значения весовой функции Н(пА?) формирующего объекта в моменты времени {пД?}^ . (-1)- и (0)-
строки являются начальными условиями при построении матрицы, а элементы ат(пА?) последовательно определяются с помощью соотношения [13]:
(22) ат (пАї) =
ат
г((п + 1)Аї) ат-1((п + 1)Аї)
ат-2(0)
а
-1(0)
где а-1 (пД?) = £(пД?); а0(пД?) = Н(пД); т = 1, 2, 3, ..., а п = 0, 1, 2, ... Тогда элементы первого столбца матрицы (21) порождают частные числители правильной С-дроби [9, 11], что и позволяет получить ДПФ формирующего объекта:
й(0) а (0) ъ- = й(0)
(23) 0( г, Аї) =
1
1
1 + -
а1(о);
і а2
1+ 2
(0)г
При аппроксимации дробно-рациональных функций конечного порядка частные числители в (23) обладают тенденцией уменьшения их величины. В модифицированном алгоритме В. Висковатова при аппроксимации дробно-рациональной функции конечного порядка в матрице (21) наблюдается появление нулевой строки, номер которой позволяет идентифицировать порядок функции.
Если в некоторой 7-й строке (7 = 0, 1, 2, ...) матрицы (21) конечное число к7 первых элементов равны 0, а последующие элементы отличны от нуля, то необходимо осуществить сдвиг влево на к7 элементов до появления в нулевом столбце ненулевого элемента и далее продолжить определение других элементов матрицы (21) по правилу (22). Для 7-й строки при восстановлении правильной С-дроби (23) элемент а7(0) умножается на г-1(7 . В этом случае от правильной С-дроби (23) переходим к обобщенной С-дроби [11] следующего вида:
Ь(0) г- к° а(0) г ~к‘
(24) О( г, Д?) =
Свернув непрерывную дробь вида (23), получим ДПФ формирующего объекта и далее перейдем к дискретной модели случайного процесса вида (19).
Для того чтобы смоделировать случайный процесс на основе (19) с дисперсией, равной Оу = Дуу(0), необходимо рассчитать дисперсию Вх = <с2х белого шума х(?). Поскольку Мх = Му = 0, то для определения Вх можно воспользоваться равенством, полученным на основе (19):
^ т п Л2
(25) М (у2 (к)) = МI ^ а,х (к - 7) -£ Ъгу(к - 7) .
V 7=0 7=1 )
Проведя несложные преобразования в (25), возможно получить формулу связи Ву и Вх через коэффициенты полученной модели.
Важным вопросом при цифровом моделировании случайных процессов на основе (19) является определение такого значения куст, начиная с которого случайный процесс у(кА?) можно считать стационарным, а режим воспроизведения -
установившемся. Поэтому необходимо отбросить искаженный участок моделируемого процесса у(кА) [3], где
к ° —, (куст — 1).
Если модель формирующего объекта устойчива, то весовая функция с ростом к стремится к нулю. Для определения времени переходного процесса в литературе используют различные методы, например, в [17] его определяют через корни характеристического уравнения цифровой модели и используя логарифмические частотные характеристики.
В данной работе временем переходного процесса куст будем считать тот интервал, по истечении которого значения весовой функции не превысят заданное малое значение А [18, 26]
\к(к)| < А .
Обычно выбирают А = 0,05.
Пример 1. При исследовании различных технических систем (например, промышленных электроустановок и следящих систем [28]) возникают случайные процессы с корреляционной функцией вида
(26) Яуу (Г) = е-.
Требуется осуществить цифровое моделирование случайного процесса у(() с корреляционной функцией (26).
На основе корреляционной функции (26) определяем спектральную плотность [17]. Далее согласно (3) находится выражение для квадрата амплитудно-фазовой частотной характеристики формирующего объекта
22 0(]о)\ = -2— . со1 +1
По корням числителя и знаменателя, а также их коэффициентам [17] определяем НПФ формирующего объекта
л/2
(27) 0(8) = ±~.
5 + 1
Используя обратное преобразование Лапласа, находим весовую функцию формирующего объекта, равную
(28) ИЦ) = д/2е ‘.
Возьмем Дґ = 1, тогда идентифицирующая матрица в модифицированном алгоритме В. Висковатова примет вид:
( 1 0 0 0 0 . >
1,4142 0,5203 0,1914 0,0704 0,0259 .
- 0,3679 - 0,1353 - 0,0498 - 0,0183 -
0 0 0
V /
Последняя строка матрицы содержит нулевые элементы, следовательно, производим останов вычислительной процедуры. Тогда ДПФ формирующего объекта аппроксимируется элементами первого столбца полученной матрицы, порождающими конечную правильную С-дробь:
1,4142
(29) 0( г) =
1 - 0,3679 г-1
ДПФ приводит к дискретной математической модели случайного процесса вида (19):
(30) у(кМ) = 1,4142 х(к М) + 0,3679у((к - 1)М).
Определим длительность переходного процесса. Значение весовой функции Н(4) = 0,0259, следовательно, куст = 4.
Рассчитаем дисперсию, с которой необходимо моделировать белый шум х(7) с учетом того, что дисперсия моделируемого процесса Ву = Руу(0) = 1. Тогда
М(у2 (к)) = 1,4142 2М(х2(к)) + 0,36792М(у2(к -1)),
Д, = 1,4142 2Д + 0,36792Д,
у ■> X ■> у ?
откуда Вх = 0,4323. График процесса у(кД7) (к = 0, ..., 900 + куст) приведен на рис. 2.
Отбросим переходный процесс и построим корреляционную функцию смоделированного процесса у(кД) по формуле [1, 22]:
1 N - к
(31) яуу, (kДt) = £ у(/Д) у № + к ) Д<),
N к 1 =1
где N = 900. Графики истинной корреляционной функции (26) и модельной Яуу (кД^) приведены на рис. 3. Дисперсия расхождений (1) Д = 0,0002, т.е. сравнение истинных и модельных значений корреляционных функций позволяет сделать вывод о высокой точности цифрового моделирования.
Рис. 3. Графики корреляционных функций к Примеру 1
5. Выбор шага дискретизации при цифровом моделировании стационарных случайных процессов
При разработке методов моделирования постоянно подчеркивается значимость величины шага дискретизации моделируемого случайного процесса, однако не приводится конкретных правил по его выбору. Авторы обычно ограничиваются реко-
мендациями о его задании на основе сравнения с интервалом корреляции случайного процесса (шаг дискретизации должен быть значительно меньше интервала корреляции) [17]. В [20] авторы дополнительно устанавливают зависимость между погрешностью восстановления корреляционной функции и величиной Д7, объемом выборки и числом модельных экспериментов. Помимо этого для определения шага дискретизации предлагается анализировать составляющие корреляционной функции на предмет выявления ее характерных особенностей [17, 20] либо строго зафиксировать шаг дискретизации для каждого конкретного вида корреляционной функции случайного процесса [20]. Последнее предложение не всегда может удовлетворять условию практической задачи или совпадать с ней, что ограничивает область физической применимости данного метода. Кроме того, большинство методов дает хорошие результаты моделирования только лишь при малых шагах дискретизации (Д7 ^ 0).
Исследуем возможности выбора величины шага дискретизации при цифровом моделировании случайных процессов, а также определим, какое влияние он оказывает на восстановление ДПФ.
Рассмотрим с единых позиций влияние величины шага дискретизации Д7 на взаимно однозначное соответствие между НПФ и ДПФ. Начало таких исследований заложено в трудах [10, 30]. Следует отметить, что при переходе от НПФ к ДПФ трудна формализация определения дополнительных нулей и полюсов, появляющихся в ДПФ объекта.
Как уже говорилось ранее, наиболее распространенным переходом от НПФ к ДПФ является преобразование г =
£-плоскости в г-плоскость. Оно отвечает требованиям взаимно однозначного соответствия при отображении непрерывного пространства в дискретное. Такое преобразование представляет собой периодическую функцию, которая отображает основную
полосу на всю г-плоскость с разрезом вдоль
отрицательной вещественной полуоси (рис. 4).
Все дополнительные полосы отображаются в аналогичную г-плоскость, и для взаимной однозначности необходимо в г-плоскости рассматривать многолистную (риманову) поверхность, для которой «склеивание» отдельных листов осуществляется вдоль линии разрезов. При отображении основной полосы 5-плоскости в г-плоскость отметим следующее: интервал
I — г —, г — I в 5-плоскости отображается в единичную окруж-^ & ЛГ 0
—
ность ББАБС; прямые у = ±г— отображаются на разрез
Л(
(-да, 0] вещественной оси г-плоскости; точки 5 с положительной вещественной частью, принадлежащие основной полосе, - во внешность единичного круга г-плоскости, а точки 5 с отрицательной вещественной частью, принадлежащие основной полосе, - во внутренность единичного круга в г-плоскости (это обуславливает устойчивость линейных динамических объектов).
Наличие разреза нарушает взаимно однозначное соответствие полуполосы 5-плоскости и круга единичного радиуса г-плоскости, так как при обратном к г = е5М логарифмическом отображении в 5-плоскость [10]
1
1
(33) 5 = — 1п г + г—аге г Л? 1 1 Л?
одной точке на этом разрезе соответствуют две точки на границе указанной полосы, так как а^ г = -arg г [16].
ляются шагом дискретизации и его изменениями, рассмотрим несколько случаев выбора шага А?.
Случай 1.
Предположим, что А? выбран таким образом, что нули и полюса объекта не попали в основную полосу
В таком случае структура дискретной модели соответствует непрерывной, однако параметры модели, определяемые через нули и полюса формирующего объекта, будут соответствовать их проекциям в основную полосу [10]. Таким образом, происходит параметрическое искажение модели.
Для получения точной цифровой модели формирующего объекта необходимо уменьшить шаг дискретизации А? до такого значения, при котором нули и полюса объекта попадут в основную полосу, причем дальнейшее уменьшение А? не будет влиять на расположение нулей и полюсов, - они будут оставаться на своих местах.
Случай 2.
Шаг дискретизации А? выбран так, что нули и полюса полученной ДПФ попали на границу основной полосы. При таком шаге дискретизации ни структура, ни параметры модели не будут определены верно [12]. Происходит структурнопараметрическая подмена модели (рис. 5).
Поскольку многие свойства преобразования г = е5А? опреде-
Рис. 5. Изменение структуры корреляционной функции при увеличении шага дискретизации
Согласно достаточному условию структурно-
параметрической идентифицируемости [10], взаимно однозначное соответствие между непрерывной и дискретной моделями объекта возможно установить только при уменьшении шага А? и только для тех А?, при которых основная полоса в 5-плоскости покрывает истинные нули и полюса объекта. Так, если кроме действительных нулей и полюсов существуют пары комплексно-сопряженных нулей и полюсов НПФ, то максимальный шаг дискретизации, при котором выполняется условие структурно-параметрической идентифицируемости, выбирается из следующего соотношения [30]:
(34) Л?• тах 1т(51м5нт, 5]п^) < п.
г
Случай 3.
Выбор шага дискретизации А? соответствует случаю попадания нулей и полюсов в основную полосу (-гп / Л? < у < гп / Л?). Структурно-параметрическое соответствие между ДПФ и НПФ установлено, и при дальнейшем уменьшении А? ширина основной полосы увеличивается, однако нули и полюса в 5-плоскости остаются на прежних местах [10]; при этом соответствующие образы нулей и полюсов в г-плоскости продолжают перемещаться внутри единичного круга, стремясь к точке (1, 0) (рис. 6).
+■ ‘У
IV
Л1 -»□ Т
Л1 -»□ I
о
71
л
1
и
+
Рис. 6. Траетории движения нулей и полюсов в 2-плоскости при уменьшении А(
Случай 4.
Шаг дискретизации А7 ^ 0. При этом отсчеты корреляционной функции отражают лишь ее начало, на котором ярко не выражены ее особенности. Как показано на рис. 6, все нули и полюса объекта в 2-плоскости при А7 ^ 0 стремятся к точке (1, 0), т.е. попадают в е-окрестность точки (1, 0) и становятся неразличимы с некоторой заданной погрешностью е (порог неразличимости).
Задав некоторую е-окрестность точки (1, 0), можно оценить шаг дискретизации А7кр, при котором все нули и полюса объекта в 2-плоскости становятся неразличимы с заданной погрешностью е:
где через 8пшх обозначен нуль или полюс НПФ, который соответствует нулю или полюсу ДПФ, максимально удаленному от точки (1, 0) 2-плоскости.
При цифровом моделировании стационарных случайных процессов необходимо также учитывать условие
где ткор - интервал корреляции процесса. При невыполнении условия ^ > ткор) измерения процесса становятся некоррелированными и получить дискретную модель формирующего объекта невозможно. Согласно [4, 22], на практике интервал корреля-
(36) Кр <& < ^кор
ции определяют как такое значение аргумента нормированной корреляционной функции
Гуу Т) = куу Т)/ куу (0), начиная с которого выполняется соотношение
(37) \гуу (т)| £ 0,05 для всех т > тюр .
Таким образом, на начальном этапе цифрового моделирования необходимо выбирать Д следующим образом.
На основе НПФ определяются все нули и полюса формирующего объекта. Свойства 2-преобразования позволяют сразу определить величину шага дискретизации, при котором дискретная модель будет соответствовать непрерывной. Если среди нулей и полюсов есть комплексно-сопряженные, то необходимо выбрать шаг дискретизации исходя из условия структурнопараметрической идентифицируемости.
Если у НПФ нет комплексно сопряженных полюсов и нулей, то шаг дискретизации можно варьировать до тех пор, пока не будет достигнута требуемая точность моделирования. С одной стороны, это означает, что прообразы дискретных полюсов (нулей) формирующего объекта совпадут с требуемой точностью с истинными значениями полюсов (нулей) НПФ. А с другой стороны, корреляционная функция смоделированного процесса должна совпадать с истинной корреляционной функцией с заданной погрешностью.
Если при выбранном Д в ДПФ появляется полюса или нули, попадающие на разрез (-да, 0) вещественной оси 2-плоскости, то в НПФ соответствующие им конечные полюса (нули) отсутствют [30].
Таким образом, метод цифрового моделирования должен при фиксированном шаге дискретизации Д восстанавливать все нули и полюса дискретной модели, включая и возникающие на разрезе [30]. Этот эффект отсутствует при моделировании подстановочными методами - прямыми и обратными разностями Эйлера и с помощью билинейного преобразования. Причем для динамического объекта существует множество значений
Д е АТ, при которых осуществляется эквивалентность дискретных моделей, сохраняющая его динамические свойства.
Пример 2. В Примере 1 был выбран шаг дискретизации Д = 1. Оценим интервал корреляции случайного процесса. Поскольку нормированная корреляционная функция полностью совпадает с (26), Яуу(3) » 0,05, то ткор = 3. Выбранный шаг удовлетворяет условию А7 < ткор. Кроме того, Д = 1 соответствует Случаю 3 выбора шага дискретизации, следовательно, он может быть использован при цифровом моделировании.
Определим значение истинного полюса НПФ (27): 5п = -1. ДПФ (29) имеет полюс 2п = 0,3679, который, согласно (33), соответствует полюсу НПФ 5п = -1, совпадающему с полюсом истинной НПФ (27). Абсолютная ошибка восстановления непрерывного полюса равна нулю.
Проведем аналогичные расчеты цифрового моделирования случайного процесса у(() с корреляционной функцией (26) с другим шагом Д( = 2. Выбранный шаг удовлетворяет условию А7 < ткор и соответствует Случаю 3.
Идентифицирующая матрица примет вид:
0 0
( 1 1,4142 - 0,1353 0
0
0,1913 - 0,0183 0
0,0259
- 0,0024
0,0035
ДПФ формирующего объекта аппроксимируется элементами первого столбца матрицы, порождающими конечную правильную С-дробь:
. ч 1,4142
(38) в(г, М) =--------------г.
У ^ 1 - 0,1353 г-
ДПФ (38) имеет полюс = 0,1353, который соответствует полюсу НПФ 8п = -1, полностью совпадающему с полюсом истинной НПФ (27).
Дискретная математическая модель случайного процесса имеет вид:
(39) у(кД) = 1,4142 х(кД) + 0,1353у((к - 1)А/) .
Значение весовой функции ^2) = 0,0259, следовательно £уст = 2 определяет время переходного процесса.
Аналогично Примеру 1 рассчитаем дисперсию, с которой необходимо моделировать белый шум x(t). Тогда Dy (1 - 0,13532)
Dx =
1,4142
= 0,4908.
Смоделируем случайный процесс у^Д), k = 0, ..., 900 + куст, отбросим переходный процесс и построим корреляционную функцию смоделированного процесса у^Д) по формуле (31). Графики истинной корреляционной функции (26) и модельной
К^у (kДt) приведены на рис. 7.
Рис. 7. Графики корреляционных функций к Примеру 2 (Лt = 2)
Дисперсия расхождений (1) De = 0,0002, т.е. сравнение истинных и модельных значений корреляционных функций позволяет сделать вывод о высокой точности цифрового моделирования.
Необходимо отметить, что в [20] для цифрового моделирования случайного процесса с корреляционной функцией (26) шаг дискретизации однозначно зафиксирован и равен Дt = 0,4. Это ограничение связано с точностью моделирования. При цифровом моделировании на основе непрерывных дробей шаг дискретизации возможно варьировать в описанных ранее пределах.
Получим цифровую модель случайного процесса у(/) с корреляционной функцией (26) и шагом Дt = 0,4. Выбранный шаг удовлетворяет условию Лt < ткор и соответствует Случаю 3. Идентифицирующая матрица примет вид:
( 1 1,4142 - 0,6703
V 0
0
0,9480 - 0,4493 0
0
0,6354 - 0,3012
0
0,4259
ДПФ формирующего объекта
, ч 1,4142
(40) в(г, Дt ) =------------
4 ' 1 - 0,6703 г-1
Полюс ДПФ = 0,6703 соответствует полюсу НПФ 5П = -1, полностью совпадающему с полюсом истинной НПФ (27). Тогда дискретная математическая модель случайного процесса имеет вид:
(41) у^Д) = 1,4142 х(к Дt) + 0,6703у(^ -1)Д/).
Значение весовой функции ^9) = 0,0386. Можно принять
^ст = 10.
Дисперсия белого шума х(/) на входе равна
Dx. =
Dv (1 - 0,67032)
= 0,2753 .
1,4142 2
После моделирования у^Д/) произведем сравнение истинной корреляционной функции (26) и модельной Куу (kДt) (таблица 1).
Таблица 1. Значения истинной корреляционной функции (26) и модельной Яу (kДt) к Примеру 2 (Лt = 0,4)
Моменты времени k 0 1 2 3 4 5
Яуу(Ш) 1 0,6703 0,4493 0,3012 0,2019 0,1353
к (М/) 1,0019 0,6756 0,4569 0,3225 0,2130 0,1517
Дисперсия расхождений (1) De = 0,0002. Сравнение истинных и модельных значений корреляционной функции позволяет сделать вывод о получении точной дискретной модели случайного процесса (с точностью до вычислительных погрешностей).
В [20] модель аналогичного случайного процесса с корреляционной функцией (26) была получена аппроксимацией на основе аналитического метода Ньютона с шагом дискретизации Л/ = 0,4, при этом погрешность расхождения De = 0,0009, объем выборки который потребовался, равен 1000 значений у^Л/).
Таким образом, приведенные расчеты дискретной модели случайного процесса с корреляционной функцией (26) с различными шагами дискретизации подтверждают эффективность предлагаемого метода моделирования.
6. Исследование работоспособности алгоритма цифрового моделирования стационарных случайных процессов на основе непрерывных дробей
Проведенные выше исследования позволили сформулировать алгоритм цифрового моделирования стационарных случайных процессов с заданной корреляционной функцией. Приведем последовательность шагов алгоритма, которые необходимо сделать, чтобы получить точную цифровую модель случайного процесса.
1. Задание корреляционной функции Куу(/) моделируемого случайного процесса.
2. Определение НПФ формирующего стохастического объекта О(^), ее нулей и полюсов.
3. Нахождение весовой функции ) формирующего стохастического объекта как обратного преобразования Лапласа от
т.
4. Оценка интервала корреляции согласно (37).
5. Задание шага дискретизации Д/, удовлетворяющего условию БР'-идентифицируемости (34) и неравенству (36).
6. Построение идентифицирующей матрицы В. Висковатова (21), получение ДПФ формирующего объекта 0(г, Д/).
7. Нахождение нулей и полюсов 0(г, Д/), а также определение соответствующих им нулей и полюсов в £-плоскости. Сравнение полученных результатов с истинными нулями и полюсами НПФ формирующего стохастического объекта О(^), полученными в п. 1 (1-й критерий определения точности моделирования).
8. Расчет дисперсии Вх входного белого шума с использованием (25).
9. Моделирование стационарного случайного процесса на основе полученного конечно-разностного уравнения вида (19).
10. Отсечение переходного процесса на основе определения времени переходного процесса куст .
11. Расчет модельной корреляционной функции согласно
(31).
12. Определение точности моделирования на основе дисперсии расхождений (1). (2-й критерий определения точности моделирования).
Далее будут проведены исследования данного алгоритма моделирования на стационарных случайных процессах с различными корреляционными функциями, часто встречающимися в практических приложениях.
Пример 3. Требуется осуществить цифровое моделирование случайного процесса с корреляционной функцией вида
[3, 20]
(42) Куу(/) = е-'1 '(1 + | /1), имеющего спектральную плотность 4
(43) Б 0) =
р(0 + 1Г
Найдем НПФ формирующего объекта
(44) в(£) = —^т,
(£ +1)2
определив предварительно квадрат амплитудно-фазовой частотной характеристики. Тогда весовая функция формирующего объекта равна (рис. 8)
(45) И(ї) = 2їе~ґ.
Далее необходимо задать шаг дискретизации. Оценим интервал корреляции случайного процесса. Нормированная корреляционная функция совпадает с (42), ^(4,8) ~ 0,048. Тогда ткор = 4,8. У НПФ нет комплексно сопряженных полюсов и нулей, поэтому согласно Случаю 3 выбора шага дискретизации зададим Дґ = 0,9.
Рис. 8. Весовая функция (45)
Идентифицирующая матрица В. Висковатова примет вид:
1 0 0 0 0 0
0,7318 0,5951 0,3629 0,1967 0,1000 0,0488
- 0,8131 - 0,4959 - 0,2688 - 0,1366 - 0,0666 -
0,2033 0,1653 0,1008 0,0546 -
- 0,2033 - 0,1653 - 0,1008 -
0 0 -
4-я строка матрицы содержит нулевые элементы. Элементы первого столбца матрицы определяют конечную правильную С-дробь
0( г, М) =
0,7318
1 + -
1 +
1 - 0,2033г-Определим ДПФ формирующего объекта
Л ч 0,7318г-
(46) G( г, М ) =-----------:-----------2 •
1 - 0,8131г - + 0,1653г- 2
ДПФ (46) имеет полюса г^2 = 0,4066, которые согласно (33),
соответствует полюсам НПФ 51п2 = -1, совпадающим с полюсами истинной НПФ (44).
ДПФ (46) приводит к дискретной математической модели случайного процесса вида (19):
(47) у(к) = 0,7318 х(к -1) + 0,8131у(к -1) - 0,1653у (к - 2) . Значение весовой функции к(6) = 0,0488, поэтому можно
принять длительность переходного процесса куст = 10.
Рассчитаем дисперсию белого шума х(7) на входе М (у 2(к)) = М (0,7318 х(к -1) + 0,8131у(к -1) - 0,1653у (к - 2))2. Откуда
Ох =
Ву (1 - 0,81312- (-0,1635)2) - 2*0,8131* (-0,1653) 0,73182
= 0,9693.
На основе (47) смоделируем процесс у(кЛ7), график которого приведен на рис. 9 (к = 0, ..., 1000 + куст).
Отбросим переходный процесс и построим корреляционную функцию процесса у(£Л0 согласно (31). Графики истинной корреляционной функции (42) и модельной ЯМ (АЛ) приведены на рис. 10.
модельная - - - - истинная
Рис. 10. Графики истинной корреляционной функции (42) и модельной Яуу (АЛ) к Примеру 3
Дисперсия расхождений (1) Ое = 0,0005, т.е. сравнение истинных и модельных значений корреляционных функций позволяет сделать вывод о высокой точности цифрового моделирования. В [20] при моделировании случайного процесса с корреляционной функцией (42) погрешность расхождения составила Ое = 0,0079.
Пример 4. Случайный процесс имеет корреляционную функцию вида [3, 20, 28]
(48) Яуу ) = 0,5 е-^ соб2^.
Например, подобный вид имеют корреляционные функции, описывающие случайные процессы типа турбулентности атмосферы.
Найдем НПФ формирующего объекта
(49) 0(8) = /+^5.
£ + 2^ + 5
НПФ имеет нуль 5н = -2,2361 и два полюса 5П2 = - 1 ± 2/. Весовая функция формирующего объекта имеет вид ^ -1 ^
(50) И(ї) = е соб2ґ +
2
е БІЙ 2і.
Интервал корреляции случайного процесса равен ткор = 2. Зададим шаг дискретизации, удовлетворяющий условию 8Р-идентифицируемости, которое для данного примера примет вид:
At • |1т 5,п
< р.
Так как 1т ^ = 2, то граница основной полуполосы отстоит от оси ординат 5-плоскости на величину
р
1т 5,п
р
■ = —»1,5708. 2
При этом шаг дискретизации должен удовлетворять неравенству At < 1,5708.
Зададим At = 0,05. Идентифицирующая матрица В. Висковатова примет вид:
1 0 0 0 0 0
1 1,0052 0,9979 0,9795 0,9511 0,9142
-1,0052 - 0,9979 - 0,9795 - 0,9511 - 0,9142 -
0,0124 0,0235 0,0332 0,0416 -
- 0,9002 -1,7040 - 2,4111 -
0 0 -
Получаем конечную правильную С-дробь: 0(г, М) = - 1
1 + -
1 +
0,0124 г -1
1 - 0,9002г -
После сворачивания дроби получим ДПФ формирующего объ-
екта:
(51) 0( г, М) =
1 - 0,8879г -1
1 - 1,8929г- + 0,9048г-2 '
ДПФ (51) имеет полюса = 0,9465 ± / 0,0950, которые, согласно (33), соответствуют полюсам НПФ 51п2 = -1 ± 2/, совпадающие с полюсами истинной НПФ (49). Кроме того, ДПФ
(51) имеет нуль = 0,8879, которому в 8-плоскости соответствует значение 5н = -2,3805. Абсолютная погрешность в определении нуля составляет 0,1444. Численные эксперименты, проведенные на тестовых объектах [9, 28], показали, что погрешность в восстановлении нуля гораздо меньше влияет на результаты дискретного моделирования, чем наличие погрешности восстановления полюсов. Покажем это, построив стохастическое конечно-разностное уравнение случайного процесса:
(52) у(к) = х(к) - 0,8879х(к -1) +1,8929у (к -1) - 0,9048у(к - 2) . Рассчитаем дисперсию белого шума х(0 на входе формирующего объекта. Запишем модель (52) в общем виде:
(53) у(к) = а0х(к) + а1 х(к -1) + Ьх у(к -1) + Ь2у(к - 2).
Возьмем квадрат математического ожидания от обеих частей равенства (53):
М (у 2(к)) = М (а0 х(к) + а1 х(к -1) + Ь1 у(к -1) + Ь2 у (к - 2))2. Тогда
Бу = М (а0 х2(к) + а]2 х 2(к -1) + 2а0а1х(к) х(к -1)) +
+ 2М (а0Ь1х(к) у (к -1) + а0Ь2 х(к) у(к - 2) +
+ а1Ь1х(к -1) у(к -1) + а1Ь2 х(к -1) у(к - 2)) +
+М (Ь2 у 2(к -1) + Ь у 2(к - 2) + 2Ьр2 у(к -1) у (к - 2)) =
= а^Бх + а12^х + 2а1Ь1М (х(к -1)у (к -1)) +
+ Ь2 Бу + Ь2 Бу + 2ЬД Яуу (1).
Найдем математическое ожидание произведения М(х(к - 1)у(к - 1)):
М(х(к -1)у(к -1)) =
= М (а0х2 (к -1) + а1 х(к -1)х(к - 2) +
+ Ь1х(к -1) у (к - 2) + Ь2 х(к -1) у(к - 3)) = а0 Их.
Подставив полученное математическое ожидание, будем иметь
Оу = а0 Вх + а12Вх + 2а1Ь1а0Ох + Ь1 Иу + Ь2Иу + 2Ь1Ь2Яуу (1) ,
откуда
(1 - Ь2 - Ь22) - 2ЬЬ2 Я (1)
(54) Вх = уК 21--------/.
а0 + а1 + 2а1Ь1а0
На основе (54) получаем, что дисперсия входного процесса Их = 0,0508.
Время переходного процесса куст = 50. На рис. 11 представлен график процесса у(кДґ), полученного на основе дискретной модели (52) (объем реализации равен 1500 значений). Графики истинной корреляционной функции (48) и модельной Яуу (кДї)
приведены на рис. 12. Дисперсия расхождений (1) равна Ве = 0,0006.
Рис. 12. Графики истинной корреляционной функции и модельной (АД) к Примеру 4 (шаг дискретизации Д = 0,05)
Посмотрим, что произойдет с моделью формирующего объекта (49) при шаге дискретизации Д( = п/2. Его задание соответствует Случаю 2 выбора шага дискретизации. Построим идентифицирующую матрицу:
(
1 1
0,2079
0
0
- 0,2079
- 0,0432
0
0
0,0432
0,0090
0
- 0,0090
Тогда ДПФ объекта аппроксимируется элементами первого столбца полученной матрицы, порождающими конечную правильную С-дробь:
(55) G(2, Ді) =------------1-г .
1 + 0,2079г4
ДПФ (55) имеет полюс = -0,2079, который соответствует полюсу НПФ 5п = -1. Произошла структурно-параметрическая подмена модели.
Проследим изменения дискретной модели, взяв шаг дискретизации, не удовлетворяющий условию £Р-
идентифицируемости. Пусть Ді = 2 соответствует Случаю 1 выбора шага дискретизации. Идентифицирующая матрица В. Висковатова примет вид:
(
1
1
0,1518
- 0,0955 0,1207 0
0
- 0,1518
- 0,0085 0,0169
- 0,0213
0
0
0,0085
- 0,0013
- 0,0012 0,0016
00 0,0013 - 0,0004
0,0004 - 0,0000
-0,0001 -
0
0,0000
После сворачивания дроби получим ДПФ формирующего объекта:
(56) 0( г, М) =
1 + 0,1769г- + 0,0183г - 2
ДПФ (56) имеет полюса 21п2 =-0,0885 ± / 0.1024, которые
соответствуют полюсам НПФ ^П2 = — 1 ± 1,1416 / . Значения
полюсов не совпадают с полюсами истинной НПФ (49), а соответствуют их проекциям в основную полосу. Кроме того, ДПФ
(56) имеет нуль = -0,0252, который в 5-плоскости отсутствует. Таким образом, происходит параметрическое искажение модели.
Приведенные расчеты дискретных моделей формирующего объекта (49) и выходного случайного процесса с корреляционной функцией (48) на основе различных шагов дискретизации подтверждают ранее сделанный вывод о том, что при наличии в НПФ объекта комплексно сопряженных нулей или полюсов шаг дискретизации нужно выбирать из условия БР-идентифицируемости.
Пример 5. Корреляционная функция случайного процесса имеет вид [17]:
(57) Яуу(V) = 0,4 в~0’9^ + 0,6 е~0’8^ 008(1,251?|).
Тогда НПФ формирующего объекта
12,96( 52 + 2,115 +1,76)
(58) в(з) =
(^2 + 1,6л + 2,2)(10л + 9)
имеет два нуля
<2 =-1,0550 ± 0,8043/
и три полюса
(59)
5п2 = — 0,8000 ± 1,24901, 53п = —0,9.
Весовая функция формирующего объекта имеет вид На) = 0,5539е~ °’9' + 0,7421е~0,8‘ соб(1,2490Г ) +
+ 0,4698 б1п(1,2490^)
Интервал корреляции случайного процесса равен ткор ~ 2,4. Зададим шаг дискретизации, удовлетворяющий условию БР-идентифицируемости:
(60) At• тах 1т(51112, 5^2) < р.
Таким образом, шаг дискретизации должен удовлетворять р
неравенству Д < ---»1,5153.
1,2490
Зададим At = 0,1. Идентифицирующая матрица
В. Висковатова примет вид:
( 1 1,2960
0 0 0 0 0 ^
1,1744 1,1013 1,0229 0,9408 0,8570
- 0,9568 - 0,9061 - 0,8498 - 0,7893 - 0,7260 0,6612 0,0097 0,0179 0,0249 0,0305 0,0348
- 0,9110 -1,6873 - 2,3316 - 2,8487 0,0058 0,0160 0,0293
0
1,2300
- 0,8936 - 2,4554
0
0
Свернув конечную правильную С-дробь, получим ДПФ формирующего объекта:
(61) 0(г, А) =
1,296 - 2,3186г- +1,0439г ~2
1 - 2,7458г- + 2,5264г“2 - 0,7788г~3
ДПФ (61) имеет полюса
которые соответствует полюсам НПФ
5п2 = — 0,7999 ± 1,2490 / и 53п = —0,8997:
совпадающие с полюсами истинной НПФ (58) с точностью до погрешностей вычислений.
ДПФ (61) также имеет нули: гн2 = 0,8945 ± 0,0726 7, г3н = 0.
В 5-плоскости = 0, попадающий на разрез, отсутствует.
Значениям комплексно-сопряженных полюсов в 5-плоскости соответствуют:
5^2 = —1,0816 ± 0,8103 7.
Примем за абсолютную ошибку определения нулей величину [30]
(62) е =
51и — 5“ + 52и — 52м
и и м м
где 51 , 52 - нули (58); 51 , 52 - нули, полученные по модели
(61). Тогда
е = 2| —1,05 5 + 0,80 4 3 7 +1,08 1 6 — 0,8 1 037| =
= 2|0,0266 — 0,0067 = 2 • 0,0007 = 0,0014.
Таким образом, ошибка в определении нулей составила 0,14%.
Стохастическое конечно-разностное уравнение случайного процесса на основе (61) будет иметь вид:
у(к) = 1,2960х(к) — 2,3186х(к — 1) + 1,0439х(к — 2) +
(63)
+ 2,7458у (к — 1) — 2,5264у(к — 2) + 0,7788у(к — 3). Рассчитаем дисперсию белого шума х(0 на входе формирующего объекта. Модель (63) в общем виде: у(к) = а0 х(к) + а1 х(к — 1) + а2 х(к — 2) +
( ) + Ъ у (к — 1) + Ъ2 у(к — 2) + Ь. у (к — 3).
Возьмем квадрат математического ожидания от обеих частей равенства (64):
М (у2 (к)) = М (а0 х(к) + а1 х(к — 1) + а2 х(к — 2) +
+ Ъу(к — 1) + Ъ2 у (к — 2) + Ъ3у (к — 3))2.
Тогда
Ву = М (а2 х 2(к) + а12 х 2(к — 1) + а^ х 2(к — 2) +
+ 2а0а1х(к) х(к — 1) + 2а0а2 х(к) х(к — 2) +
+ 2а2а1х(к — 2) х(к — 1)) +
+ 2М (а0Ъ1х(к) у(к — 1) + а0Ъ2 х(к) у(к — 2) +
+ а0Ъ3 х(к) у (к — 3) + а1Ъ1х(к — 1) у (к — 1) +
+ а1Ъ2 х(к — 1) у(к — 2) + а1Ъ3 х(к — 1) у(к — 3) +
+ а2Ъ1х(к — 2) у(к — 1) + а2Ъ2 х(к — 2) у(к — 2) +
+ а2Ъ3 х(к — 2) у(к — 3)) +
+М (Ъ2 у2 (к — 1) + Ъ22 у2 (к — 2) + Ъ32 у2 (к — 3) +
+ 2Ъ1Ъ2 у (к — 1) у (к — 2) + 2Ъ1Ъ3 у(к — 1) у(к — 3) +
+ 2Ъ2Ъ3 у(к — 2) у(к — 3)) =
= а02 Ох + а12 Ох + а2 Ох + Ъ12 Оу + Ъ2 Оу + Ъ32 Оу +
+ 2Ъ1Ъ2Яуу (1) + 2Ъ1Ъ3^уу (2) + 2Ъ2Ъ3Яуу (1) +
+ 2а1Ъ1М (х(к — 1)у (к — 1)) + 2а2Ъ1М (х(к — 2) у(к — 1)) +
+ 2а2Ъ2М (х(к — 2) у(к — 2)).
Найдем математическое ожидание произведения М(х(к - 1)у(к - 1)):
М(х(к — 1)у (к — 1))=
= М (а0 х2 (к — 1) + а1 х(к — 1) х(к — 2) + а2 х(к — 1) х(к — 3) +
+ Ъ1 х(к — 1) у(к — 2) + Ъ2 х(к — 1) у (к — 3) +
+ Ъ3 х(к — 1) у(к — 4) = а0 Вх.
Математическое ожидание произведения М(х(к - 2)у(к - 1)):
М (х(к — 2)у(к — 1)) =
= М (а1х 2(к — 2) + а0 х(к — 1) х(к — 2) + а3 х(к — 2) х(к — 3) +
+ Ъ2х(к — 2)у (к — 3) + Ъ3 х(к — 2)у (к — 4) + Ъ1х(к — 2)у(к — 2) =
= а1Ох + ЪМ (х(к — 2) у(к — 2)).
Математическое ожидание произведения М(х(к - 2)у(к - 2)):
М (х(к — 2) у(к — 2)) =
= М (а0 х 2(к — 2) + а1 х(к — 3) х(к — 2) + а2 х(к — 2) х(к — 4) +
+ Ъ2 х(к — 2) у(к — 4) + Ъ3 х(к — 2) у(к — 5) +
+ Ъ1х(к — 2) у (к — 3) = а0 Ох.
Подставив полученные математические ожидания, будем иметь
Оу = а02 + а12 + а2 + +Ъ1 Оу + Ъ2 Оу + Ъ32 Оу +
+ 2аЪа0 Вх + 2Ъ^ Ду (1) + 2ЪаЪ3 Ду (2) + 2Ъ3Ъ2 Ду (1) +
+ 2а2Ъ а0 Я + 2а2Ъ1а1 + 2а2Ъ2а0 Я.
Окончательно
(65) Ох =
Оу (1 — Ъ — Ъ22 — Ъ2) — (2Ър2 + (1) — 2^^ (2)
а0 + а1 + а2 + 2а1Ъ1а0 + 2а2Ъ1 а0 + 2а2Ъ1а1 + 2а2Ъ2а0
На основе (65) получаем, что дисперсия входного процесса Д, = 0,0963.
Время переходного процесса куст = 30. На рис. 13 представлен график процесса у(кД7), (объем реализации равен 1700 значений).
Графики истинной корреляционной функции (57) и модельной Дуу (кД^) приведены на рис. 14. Дисперсия расхождений (1) равна Ое = 0,0007.
Рис. 13. График процесса (63)
Рис. 14. Графики истинной корреляционной функции и модельной ЕУу (кД) к Примеру 5
7. Выводы
С использованием теории непрерывных дробей разработан алгоритм цифрового моделирования, позволяющий построить точную дискретную модель случайного процесса с заданной корреляционной функцией.
Приведенные теоретические и модельные исследования позволяют утверждать, что предложен достаточно эффективный и простой в реализации переход к дискретному представлению стационарного случайного процесса. В отличие от существующих методов цифрового моделирования данный метод не имеет погрешности моделирования, позволяет теоретически обосновать выбор шага дискретизации, устанавливать взаимно однозначное соответствие между заданной непрерывной и получаемой дискретной моделями случайного процесса, а также может использоваться для моделирования формирующих объектов с различными функциональными свойствами.
Литература
1. БЕНДАТ ДЖ., ПИРСОЛ А. Прикладной анализ случайных процессов. - М.: Мир, 1989. - 464 с.
2. БУСЛЕНКО Н.П. Моделирование сложных систем. - М.: Наука, 1968. - 356 с.
3. БЫКОВ В.В. Цифровое моделирование в статистической радиотехнике. - М.: «Советские радио», 1971. - 328 с.
4. ВОЛГИН В.В., КАРИМОВ Р.Н. Оценка корреляционных функций в промышленных системах управления. - М.: Энергия, 1979. - 80 с.
5. ГРИДИНА Е.Г., ЛЕБЕДЕВ АН., НЕДОСЕКИН Д.Д., ЧЕРНЯВСКИЙ Е.А. Цифровое моделирование систем стационарных случайных процессов. - Л.: Энергоатомиздат, 1991. - 144 с.
6. ДЕЙЧ А.М. Методы идентификации динамических объектов.
- М.: Наука, 1985. - 240 с.
7. ЕМЕЛЬЯНОВ В.Ю. Методы моделирования стохастических систем управления: учеб. пособие. - СПб.: изд-во Балт. гос. техн. ун-та, 2004. - 168 с.
8. ЖОВИНСКИЙ АН., ЖОВИНСКИЙ В.Н. Инженерный экспресс-анализ случайных процессов. - М.: Энергия, 1979.
- 112 с.
9. КАРТАШОВ В.Я. Анализ и исследование аппроксимационных свойств непрерывных дробей при решении задачи структурно -параметрической идентификации динамических объектов // Препринт №22. - Барнаул: Изд-во Алтайского госуниверситета, 1996. - 40 с.
10. КАРТАШОВ В.Я., ИНДЕНКО ОН., АЛЕКСАНДРОВ А.В.
Влияние периода дискретизации на структурнопараметрическое взаимосоответствие между непрерывной и дискретной по времени моделями линейного динамического объекта // Препринт №15. - Барнаул: Изд-во Алтайского госуниверситета, 1996. - 36 с.
11. КАРТАШОВ В.Я. Непрерывные дроби (определения и свойства). - Кемерово: Изд-во Кемеровского госуниверситета, 1999. - 88 с.
12. КАРТАШОВ В.Я., НОВОСЕЛЬЦЕВА М.А. Идентификация стохастических объектов: учебное пособие. - Томск: Изд-во Томского государственного педагогического университета, 2008. - 104 с.
13. КАРТАШОВ В.Я., НОВОСЕЛЬЦЕВА М.А. Способ идентификации линейного объекта // Патент РФ №2146063. - 2000. -Бюл. №6.
14. КАРТАШОВ В.Я., НОВОСЕЛЬЦЕВА М.А. Структурнопараметрическая идентификация стохастических объектов с использованием непрерывных дробей // Управление большими системами. - 2008. - Вып. 21. - С. 27-48.
15. КАРТАШОВ В.Я. Эквивалентность дискретных моделей -реальность? // Промышленные контроллеры АСУ. - 2006. -№8. - С. 40-44.
16. КОРН Г., КОРН Т. Справочник по математике для научных работников и инженеров. - М.: Наука, 1974. - 832 с.
17. ЛЕБЕДЕВ АН., НЕДОСЕКИН Д.Д., СТЕКЛОВА Г.А., ЧЕРНЯВСКИЙ Е.А. Методы цифрового моделирования и идентификации стационарных случайных процессов в информационно-измерительных системах. - Л.: Энергоатомиздат. Ленингр. отд-ние, 1988. - 64 с.
18. МАКАРОВ И.М., МЕНСКИЙ Е.М. Линейные автоматические системы. - М.: Машиностроение, 1982. - 464 с.
19. МАКС Ж. Метода и техника обработки сигналов при физиче-скихизмерениях. Том 1. - М.: Мир, 1983. - 312 с.
20. Прикладной анализ случайных процессов / Под ред. С .А. Прохорова. - СНЦ РАН, 2007. - 582 с.
21. РАБИНЕР Л., ГОУЛД Б. Теория и применение цифровой обработки сигналов. - М.: Мир, 1978. - 848 с.
22. РОМАНЕНКО А.Ф. СЕРГЕЕВ Г.А. Вопросы прикладного анализа случайных процессов. - М.: Советское радио, 1968. -247 с.
23. САМАРСКИЙ А.А., МИХАЙЛОВ А.П. Математическое моделирование: Идеи. Методы. Примеры. - М.: Физматлит, 2002. - 320 с.
24. СИБЕРТ У.М. Цепи, сигналы, системы. Том 1.- М.: Мир, 1988.
- 510 с.
25. СМИТ ДЖ. М. Математическое и цифровое моделирование для инженеров и исследователей. - М.: Машиностроение, 1980.
- 271 с.
26. СОЛОДОВНИКОВ ВВ., ПЛОТНИКОВ В.Н., ЯКОВЛЕВ
А.В. Основы теории и элементы систем автоматического регулирования: учебное пособие для вузов. - М.: Машиностроение, 1985. - 536 с.
27. ФИХТЕНГОЛЬЦ Г.М. Курс дифференциального и интегрального исчисления. Том 2. - М.: Наука, 1966. - 800 с.
28. ШАЛЫГИН А.С., ПАЛАГИН Ю.И. Прикладные модели статистического моделирования. - Л.: Машиностроение, 1986. - 320 с.
29. ШЕННОН Р. Имитационное моделирование систем -искусство и наука. - М.: Мир, 1978. - 420 с.
30. ЩЕКОЧИХИНА С.Г. Разработка метода дискретного моделирования в задачах диагностики сложных объектов горной техники: дис. канд. тех. наук. - Кемерово, 1999. - 279 с.
DIGITAL DESIGN OF STATIONARY CASUAL PROCESSES WITH THE SET CROSS-CORRELATION FUNCTION ON BASIS OF THE CONTINUED FRACTIONS
Vladimir Kartashov, head of the Department of Automation of Researches and Technical Cybernetics of the Kemerovo State University, Doctor of engineering sciences, professor (kar-tash@kemsu.ru).
Marina Novoseltseva, associate professor of Department of Automation of Researches and Technical Cybernetics of the Kemerovo State University, Candidate of engineering sciences, (aanov@pochta.ru).
Abstract: On the basis of the continued fractions the algorithm of design of stationary casual process is worked out, accordance is set between continuous and discrete the models of casual process, the choice of quantization interval is reasonable.
Keywords: stationary casual process, cross-correlation function, continued fraction, forming object, white noise, continuous transmission function, quantization interval, discrete transmission function, stochastic difference equation.
Статья представлена к публикации членом редакционной коллегии В. А. Уткиным
