Адаптивное оптимальное в среднеквадратическом смысле дифференцирование Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 681.518.22 ББК 32.96

АДАПТИВНОЕ ОПТИМАЛЬНОЕ В СРЕДНЕКВАДРАТИЧЕСКОМ СМЫСЛЕ ДИФФЕРЕНЦИРОВАНИЕ

Гуляев С. В.1, Шубладзе А. М.2,

Малахов В. А.,3 Ольшванг В. Р.4

(Учреждение Российской академии наук Институт проблем управления РАН, Москва) Кротов А. В.5

(ОАО «Газавтоматика»)

Предлагается метод решения задачи дифференцирования, позволяющий получать оценки производных гауссовских стационарных сигналов, близкие к оптимальным в среднеквадратическом смысле, если спектральные плотности полезного сигнала и помехи известны с точностью до уровня. Решается задача с помощью специальным образом организованных нелинейных динамических систем. Найдено близкое к оптимальному решение задачи дифференцирования, когда для любых фиксированных уровней дробно-рациональных спектральных плотностей полезного сигнала и помехи параметры эквивалентной передаточной функции нелинейного динамического дифференциатора могут быть сделаны близкими к параметрам оптимального винеровского фильтра.

1 Сергей Викторович Гуляев, кандидат технических наук (svgul@inbox.ru).

2 Александр Михайлович Шубладзе, доктор технических наук, профессор (shub@ipu.ru).

3 Валерий Александрович Малахов, кандидат технических наук, (тел. (495) 334-88-81).

4 Владимир Рафаилович Ольшванг, кандидат технических наук, (тел. (495) 334-88-81).

5 Александр Васильевич Кротов, начальник управления, alex_k@gazauto. gazprom.ru

Ключевые слова: дифференцирование, адаптация, оптимальность, гауссовский шум.

1. Введение

При решении задач управления технологическими процессами часто требуется находить производные гауссовских стационарных сигналов в случае, когда интенсивность сигнала и помехи может меняться в широких пределах. Известные методы решения такой задачи, связанные с набором статистики, не всегда могут использоваться для решения практических задач управления, так как темп изменения параметров сигналов может оказаться достаточно высоким. Поэтому представляет интерес разработка таких способов оценки производных, которые могли бы решать задачу для таких случаев.

2. Постановка задачи

Постановка задачи дифференцирования следующая. Наблюдаемый сигнал

(1) 2^) = х(0 + ф(0,

где х(0 - полезный стационарный гауссовский сигнал; ф(0 -стационарная гауссовская помеха, некоррелированная с х(0. Спектральная плотность полезного сигнала

т- р

б( I Ь]Ю21 +1)

(2) /х (б, о) = ------■----,

I а}ю21 +1 1=1

где т < п, т - р > 0, т, р, п, Ь и а^ - известные параметры, б > 0 - неизвестный параметр. Спектральная плотность помехи

Я! й,а21

(3) Ъ (Я, о) = ------,

2 +1

I=1

где й < п, е, $, й, е1 - известные параметры; Я > 0 - неизвестный параметр.

76

Как следует из работ [1, 3], в этом случае передаточная функция оптимального в среднеквадратическом смысле дифференциатора

(4) И0 (о) =------1---?е-гМЛ ? (ш)/(о)е

0,4 ' 2пу(т){ ± у* (о)

где д - порядок дифференцирования оптимального фильтра, функции ц(/'ш) и ц (¡ш) находятся из уравнения

(5) у(т)у(т) = ц(7'ш)ц(-7'ш) = /(ш) + /ф(ш),

Выражение Ж0^ш) (4) можно преобразовать к виду

(6) = 1

0,4 їжу (ію) (ію У/х (ю)

(ію)9 /х (ю) у * (ію)

^ , (оУ/х (о)

где ^у ■'хУ ' - сумма вычетов функции ----*— ---, распо-

у * (¡о) ] у * (¡о)

ложенных в верхней полуплоскости комплексной плоскости. Из

(5) следует, что

(7) ц('ш) = л/я Е(О/Я, ¡ш), щ(1ю) = л/я Е(О/Я, -¡ш),

где Е - дробно-рациональная функция введенных выше параметров.

Имея в виду, что сумма вычетов дробно-рациональной функции также является дробно-рациональной функцией, передаточную функцию (6) можно представить в виде

9

I *1 (О/Я)(о )1

(8) ^ДО/Я,¡о) = ^-^-----------------------,

I * (О /Я )(о)1 + (о )п

1 =0

т. е. параметры передаточной функции И09(О/Я, ¡ш) (7) являются непрерывной функцией отношения О/Я.

3. Синтез нелинейного дифференциатора

Используем полученную передаточную функцию оптимального дифференциатора для синтеза нелинейного дифференциатора, эквивалентная передаточная функция которого

была бы близка к передаточной функции оптимального дифференциатора при любых значениях отношения О/Я. Близость этих передаточных функций при произвольных О/Я будет гарантировать высокое качество дифференцирования нелинейного дифференциатора, так как известно [1], что среднеквадратическая ошибка дифференцирования любого фильтра с передаточной функцией Иф(/'ш) определяется выражением

(9) ^ = о02+1 К (°)- И0 (¡'®)|2/(0, Я ш)dш,

где

о02 = 1 (|(°)912 / (0, -К (0, Я,

- среднеквадратическая ошибка фильтрации оптимального винеровского дифференциатора.

Нелинейные, близкие к оптимальным дифференциаторы будем синтезировать с помощью симметричных степенных функций, логику работы которых в рассматриваемом случае удобно представить в виде

(10) и а (г ) = и а (г )sign( s(t))

где ,?(^ связано с г(() передаточной функцией

(11) (¡о) = °-------------,

I *1 (О /Я)(т)1 + (¡о)п

1=0

*1 (О/Я), 1 = 0, ..., п - 1, - параметры знаменателя передаточной функции оптимального дифференциатора, а и*а (/) имеет вид

(12) иа(Г) = кда \д(Г)|“ |*(0|1-а> 0,

9(?) - линейное преобразование 2^) (1).

Расчет дифференциатора с нелинейным Í ^ -элементом (10)

наиболее удобно проводить методом статистической линеаризации, который описан в работе [2]. В связи с этим необходимо, прежде всего, найти эквивалентный коэффициент передачи

этого элемента, обеспечивающий минимальную среднеквадратическую ошибку

(12) s2 = M{[V(t)-ua(t)J } = min, где

u\t) = k. a 5(t).

J s

Условие (12) справедливо в предположении, что q(t) и s(t) имеют нулевые математические ожидания. Минимум выражения (12) ищется из условий

(13)

да2

дк. „

J s

да2 дк с

дм Л k a s(t) - ua (t)

дкі

: 0.

После дифференцирования (13) получаем

(.4) .=М<-{м})}.

1' М {2^)}

Считая, что д(0 и s(t) - независимые гауссовские сигналы, найдем из (14) а

(15) M {ua (t )s (t )} =

4 гг, „2-а „а 2а^ 2а!

■JK

2ка>а, 00

s q e

где ss - дисперсия s(t), а sq - дисперсия q(t). Преобразуя (15), получаем

(16) M {иа (t)s(t)}=-

і ¿-.а— 2-е

kqa аq аs

2ps sS q

a +1

3 - a

где Ä(a) = J xa 1e xdx. Знаменатель k. a (14)

2* оао2_а ¥ -_8

(17) М{2(t)}= 8 1Л оds = о2,

\2ко6, 0

а сам коэффициент, как следует из (14), (16) и (17),

оа ^ ^

а +1

3 - а

(18) * а=А

8 о

Следовательно, эквивалентный коэффициент передачи *■ а

1 8

(18) I а -элемента (10) зависит от отношения дисперсий, характеризующих 8^) и д(0.

Отмеченные свойства оптимального дифференциатора и

найденный эквивалентный коэффициент передачи *■ а (18)

1 8

элемента (10) дают возможность синтезировать нелинейный дифференциатор, эквивалентная передаточная функция которого была бы близка при произвольных отношениях О/Я к передаточной функции оптимального дифференциатора, рассчитанного для случая (1)-(3). Такая близость гарантирует, прежде всего, устойчивость нелинейного дифференциатора, что является основным требованием, предъявляемым к любым устройствам. В этой связи наиболее естественно и просто удовлетворить требованию устойчивости, если знаменатель эквивалентной передаточной функции дифференциатора при произвольных отношениях О/Я будет близок к передаточной функции оптимального дифференциатора, так как оптимальный дифференциатор, рассчитанный из условия минимума среднеквадратической ошибки, всегда устойчив.

Рассмотрим линейную систему, знаменатель передаточной функции которой совпадал бы при любых отношениях О/Я со знаменателем W0,q(Q/R, ш) (8):

(19) 2(() -/(0 = 8^), /(0 = Ух(Г) + У)(0,

(20) у) = *п-г (О/Я)) + у+! (t), /■ = 1, ., п,

(21) уП-ДО = ЦО/ЯШ,

(22) y0n)(t) = ^(О /Я)5^).

Нетрудно заметить, что знаменатель передаточной функции системы (19)-(22), например, от г(^) к у0(0, совпадает со знаменателем W0,q(Q/Я, ш) (8). Блок-схема дифференциатора, соответствующая этой системе, представлена на рис. 1.

Добиться близости знаменателей эквивалентной передаточной функции нелинейного дифференциатора и W0,q(Q/Я, ш) (8) можно, заменив в системе (19)-(22) элементы *1 (О /Я),

i = 1, ..., п, нелинейными элементами (10) так, чтобы их эквивалентные коэффициенты передачи при любом значении отношения О/Я были сколь угодно близки к *1 (О/Я) .

Для проведения такой замены установим связь между эквивалентным коэффициентом передачи к■ a (18) и отношением

^ s

Q/R. Пусть одним из входных сигналов элемента (10) является сигнал s(t) (10), (11), дисперсия которого в случае (1)-(3) определяется помехой, т. е.

1

Я I Л1ю2

(23) а2 = — [ — " 2р/ А-

-Лю = ЯА.

I ^2' +1

I=1

В качестве второго сигнала, входящего в логику элемента

(10), возьмем выходной сигнал q(t) низкочастотного линейного фильтра с передаточной функцией Ж^ю), на входе которого действует наблюдаемый гауссовский сигнал г(0 (1). В силу того, что фильтр низкочастотный, гауссовский сигнал q(t) будет определяться низкочастотной составляющей сигнала 2(1). Дисперсия сигнала q(t), от которой зависит к. а , на основании (1)-

(3) имеет вид

1 ¥ I |2

(24) а^2 = —/К (г'®)1 ^ О *, ®)^® = вВ + КС,

где

т-р

I Ът21 +1

1 7 I |2 _

(25) В = — / ^ (/ю) 1-------------------------ёю,

2ж •' 1 1 —■

I а}а>21 +1

і=1

I Лю21 +1

-Лю

I СгЮ21 +1

і=1

или

а я =^1ОВ + ЯС.

Выражения (18), (23) и (24) позволяют определить зависимость эквивалентного коэффициента передачи

(27) к. а =

р А

0,5а

а +1

3 - а

Б0’5а (Є /Я + С /+)

0,5а

от интересующего нас отношения О/Я для сигналов 5(1) и q(t), которые определяются дифференциальными уравнениями 5(0 = Wzs(p)z(t),

где Ж^(/«) определяется (11), а р = — и

—

g(t) = Wzg(p)z(t), р = —,

—

где ^^(р) - передаточная функция низкочастотного фильтра, входящая в выражения (24)-(26).

Теперь для совпадения знаменателей эквивалентной передаточной функции нелинейного фильтра и передаточной функции оптимального фильтра Ж0>с(О/Я, /ю) (8) достаточно каждый из коэффициентов к1, / = 0, ..., п - 1, (20)-(22) представить с

требуемой степенью точности в виде разложения по функциям

(28) (О/Я + С/В)0 5а = (О/Я + £>/, а = 2р > 0,

которые пропорциональны соответствующим эквивалентным коэффициентам передачи -элементов.

Разложения

_ N _

(29) к} (О/Я) »Iк}1 (О/Я + ^,

I =1

где (О /Я + В)Ь]1 при различных 1 - линейно независимые функции, позволяют установить связь между неопределенным пока параметром кс< (12), входящим в к, а (27), и коэффициентами

Уа 1 5

разложения кц (29) функции к1 (О/Я). Как следует из (27), (28) и (29), при выполнении равенства

кр АЬЬ1

(30) к - 1

jl

2

коэффициенты разложения кj¡ (29) совпадают с соответствующими коэффициентами перед степенными функциями в кí a (Q/R) (27). Заменяя в системе (19)-(22) каждый из кi (Q/R)

* s •*

группой элементов (10) с параметрами аj = 2в;1 и kg в соответ-

bjl

ствии с (29), (30) добьемся того, что знаменатель эквивалентной

передаточной функции нелинейного дифференциатора будет с требуемой степенью точности совпадать со знаменателем передаточной функции оптимального дифференциатора при любом отношении О/Я в постановке (1)-(3). Следовательно, степень устойчивости нелинейного дифференциатора сколь угодно близка к степени устойчивости оптимального дифференциатора.

Качество дифференцирования, определяемое величиной среднеквадратической ошибки дифференцирования, будет тем

2

выше, чем меньше \Ж</ё (О/Я, гю)-W0q (/Я,/а)1 , где

WШ(Q/R, /ю) - эквивалентная передаточная функция нелинейного дифференциатора. Наибольшая близость передаточных функций Wэк(Q/Я, /ю) и Wйq¡(Q/Я, /ю) в смысле минимума о — (9) будет в нелинейном дифференциаторе с элементами (10), описываемом следующей системой уравнений

(31) 2(0 - У(0 = 5 ^), у (t) = у (t) + У0 (t),

(32) у'(1 ) = Пп_г^) + у+^), / = 1, ., п - 1,

(33) УП (t) = к5 (t),

(34) У0п)(1 ) = Щ(!) - к5 ^),

(35) Ы1 (t) = 1 ы2ь^ ^),

I=1

где параметры управления и2^ г (10), (11), такие как Рц и к&1 ^ ,

находятся соответственно из разложения к1 (О /Я) (29) и формулы для к^^ (30), коэффициент к является пределом

(36) к = 11те/я®0 Iк01(О/Я + ^)Ам = Ите/я®0 к0 (О/Я) .

г=1

Выходными координатами нелинейного дифференциатора являются У01 ^), 1 = 0, ..., п - 1 (34). Блок-схема нелинейного дифференциатора, соответствующая системе (31)-(35), изображена на рис 2.

Введение в схему нелинейного дифференциатора элемента к (33), (34) и (36) сказалось только на числителе

(37) Ж# (в/Я,іт) = п-1

(іту (К (в/Я) - к)

х к] (в/Я)(ш)1 + (іт)п

1=о

N.

где к] =хк т , ] = 0, п - 1.

і=1

р 1 «-(ХН- 1 1 1

р Р р р р

1 1 1 1 1

— — <— — < — «—

р р р р р

Ямс. 2. Блок-схема нелинейного дифференциатора

Как показано в [3], модуль передаточной функции оптимального дифференциатора |W0 с (О/Я,/ю)| ® |(/'ю)с | , с = 0, .,

п - 1, при О/Я ^ <» и с (О/Я,/ю)| ® 0 при О/Я ^ 0. Из перво-

го предельного соотношения на основании (8) следует, что к1 (О/Я) и к1 (О/Я)(ку+1 (О/Я))- ® ¥ при О/Я ^ да. В этом

случае, имея в виду к^ (О/Я) » к1 (О/Я), получаем

85

Wэк(Q/Я, /ю) ^ (/ю)с, с = 0, ., п - 1 при О/Я ^ да, так как к} (О/Я)(к1+1(О/Я))-1 ® 0 при О/Я ^ да. В другом предельном соотношении при О/Я ^ 0, как и в оптимальном случае, WэK(Q/Я, /ю) ^ 0, так как к0‘ (О/Я) - к ® 0 . Имеет смысл отметить, что в наиболее интересном диапазоне 0 < О/я < 1, когда оптимальный дифференциатор с передаточной функцией (8) обеспечивает достаточно качественные оценки полезного сигнала и его производных, числитель передаточной функции оптимального дифференциатора с высокой степенью точности приближается функцией (/ю)с. Но в диапазоне 0 < О/Я < 1 числитель эквивалентной передаточной функции (37) нелинейного дифференциатора (31)-(35) также с высокой степенью точности приближается той же функцией (/ю)с.

Следовательно, качество дифференцирования нелинейного дифференциатора (31)-(35) близко к качеству дифференцирования оптимального дифференциатора при любом О/Я < 1.

4. Заключение

Сравнивая подстраиваемый оптимальный для любого 0 < О/Я дифференциатор (8) с нелинейным дифференциатором

(31)-(35), можно сказать, что, незначительно проигрывая в качестве дифференцирования, нелинейный дифференциатор, во-первых, реализуется существенно проще подстраиваемого оптимального дифференциатора, где в результате набора статистики определяется отношение О/Я, которое затем по нелинейным законам преобразуется, формируя переменные коэффициенты к^ (О/Я) и к^ (О/Я) (8), и, во-вторых, время адаптации к

изменяющимся параметрам О и Я (2), (3) в нелинейном дифференциаторе (31)-(35), который не требует набора статистики, существенно меньше времени адаптации в подстраиваемом оптимальном дифференциаторе.

Синтезированные в настоящей работе дифференциаторы в силу нелинейных свойств переключающих элементов адаптируются к неизвестным параметрам полезного сигнала и помехи, делая их практически всегда без получения каких-либо допол-

86

нительных оценок случайных сигналов близкими к оптимальным в среднеквадратическом смысле. Для получения подобных свойств адаптации в классических оптимальных дифференциаторах необходимо проводить статистические оценки уровней случайных сигналов и с помощью такой статистики корректировать параметры в их передаточных функциях. Очевидно, что такой процесс адаптации и долгий и сложный, поэтому практически он не применяется.

Литература

1. ЕМЕЛЬЯНОВ С. В. и др. Теория систем с переменной структурой. - М.: Наука, 1970. - 592 с.

2. КАЗАКОВ И. Е., ДОСТУПОВ Б. Г. Статистическая динамика нелинейных автоматических систем. - М.: Физмат-гиз, 1964. - 332 с.

3. ПУГАЧЕВ В. С. Теория случайных функций и ее применение к задачам автоматического управления. - М.: Физматгиз, 1962. - 883 с.

THE ADAPTIVE OPTIMAL DIFFERENTIATION BY STANDARD DEVIATION CRITERION

Sergey Gulyaev, V.A. Trapeznikov ICS of RAS, Moscow, Cand.Sc., assistant professor (svgul@inbox.ru).

Alexander Shubladze, V.A. Trapeznikov ICS of RAS, Moscow, Doctor of Science, professor (shub@ipu.ru)

Valery Malakhov, V.A. Trapeznikov ICS of RAS, Moscow, Cand.Sc., assistant professor (Moscow, Profsoyuznaya st., 65, (495)334-88-81).

Vladimir Olshvang, V.A. Trapeznikov ICS of RAS, Moscow, Cand.Sc., assistant professor (Moscow, Profsoyuznaya st., 65, (495)334-88-81).

Alexander Krotov, Head of the department of “Gazavtomatika” (Moscow, Savvinskaya em., 25), (499)580-41-22.

Abstract: The method is suggested to solve the differentiation problem. It allows building the estimates for the derivatives of Gaussian stationary signals that are close to the optimal ones by the standard deviation criterion when the spectral density of a useful signal and a noise are known to within level. The problem is solved by the specifically designed nonlinear dynamic systems. A nearly optimal solution of the differentiation problem is found when the equivalent transfer function parameters of nonlinear dynamic differentiator can be made close to the parameters of the optimal Wiener filter for any fixed level of rational spectral density of a useful signal and a noise.

Keywords: differentiation, adaptation, optimality, Gaussian

noise.

Статья представлена к публикации членом редакционной коллегии А. С. Манделем