Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Солодовников В.В. Статистическая динамика линейных систем автоматического управления

.pdf
Скачиваний:
11
Добавлен:
30.10.2023
Размер:
19.04 Mб
Скачать

2 Ю

ПРОХОЖДЕНИЕ СЛУЧАЙНОГО

СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ [ГЛ. V

кроме того, мы располагаем

записью процесса

е ((), коррелирован­

ного с воздействием на входе

и

не

коррелированного с возмущением

/ (0 внутри объекта.

 

 

 

 

В

этом

случае

характеристики

объекта легко определяются из

уравнения

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rxe(-) = j

Rse(x — 8) k (8) rf8

(5.155)

 

 

 

 

6

 

 

 

 

Входящие

в это

уравнение

корреляционные

функции Rxe (т) и

R iC(^) вычисляются

обычным

способом по имеющимся реализациям

х (t),

е (t) и e (t).

 

 

 

 

 

Для решения уравнения (5.155) можно применить частотный метод,

описанный

в §

7.

 

 

 

 

 

В силу физического условия независимости внутренних помех от

внешних, воздействий в качестве е можно

взять

величину ср, посту­

пающую на вход системы.

 

 

 

 

Тогда уравнение (5.155) перепишется в

виде

 

 

 

СО

 

 

 

Я*?(т) =

J ДеДт — 8)A(8)rf8.

 

(5.156)

 

о

 

 

 

В случае, если непосредственно не удается

получить

реализа­

цию <р(0, можно использовать равенство

 

 

 

<Р(0 = e (/) -|- z (0

 

 

(5.157)

и определять k(t) из уравнения

 

 

 

СО

 

 

 

 

Я* СО+ /?«W = /

IЛ .(х -» ) + /г „ (т -8 )] A(9)d0.

(5.158)

О

9]

НАЛИЧИЕ В СИСТЕМЕ ОБРАТНЫХ СВЯЗЕЙ И ВНУТРЕННИХ ПОМЕХ

241

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

спомехами внутри объекта.

Врассматриваемом случае в качестве величины е, не коррели­

рованной с

помехой / 3 (см.

рис. 5.31), на основании предположения

о независимости внутренних

шумов от внешних воздействий можно

взять одно

из воздействий:

cpL(^), <р2(0- ?з (0- Пусть e(f) = cp1(/);

е (/), х (() — записи процессов на входе и выходе исследуемого звена. Тогда динамические характеристики звена 3 определяются из

уравнения, совершенно идентичного

уравнению (5.156)

 

£ .« ,,('0 = fСО jR«Ti ( t

— 8) ft (9) db.

(5.159)

о

 

 

J6 Зак. 1083. В. В. Солодовников

Г Л А В А VI

СТАТИСТИЧЕСКИЕ МЕТОДЫ СИНТЕЗА ОПТИМАЛЬНЫХ СИСТЕМ

1. Введение

Системы автоматического управления могут рассматриваться как системы, служащие для преобразования совокупности входных воз­ действий {/?1г(0} в совокупость сигналов на выходе (хД/)} в соот­ ветствии с некоторым оператором Н :

{хДО }==«{,»,(/)}.

 

(6.1)

Н называют идеальным преобразующим оператором.

 

Однако обычно вследствие наличия возмущающих

воздействий и

ограничений, накладываемых условиями физической

реализуемости,

преобразование воздействий

в соответствии

с (6.1) оказывается не­

возможным, и задача синтеза

состоит в том,

чтобы определить такой

физически осуществимый оптимальный оператор Ф, который обеспе­ чивает наименьшую возможную при рассматриваемых условиях дина­ мическую ошибку преобразования сигналов по сравнению с идеальным законом преобразования (6.1).

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

1)

Что следует понимать под заданным законом или идеальным

оператором преобразования управляющего воздействия?

2)

Какой

информацией о воздействиях необходимо

располагать?

3)

Каким

образом характеризовать динамическую

точность вос­

произведения системой заданного закона преобразования управляющих воздействий?

2. Преобразующий оператор

Если ограничиться для простоты линейными преобразующими операторами H(s), то можно сказать, что система должна возможно более точно воспроизводить на своем выходе не само управляющее воздействие т. (t), а некоторый сигнал /г (t), связанный с сигналом т (t) заданным функциональным преобразованием

h{t) = H (jи) т

(6.2)

или, точнее,

L [ h ( t ) ] = H ( s ) L [ n i ( t ) ] .

(6.3)

21

ПРЕОБРАЗУЮЩИЙ ОПЕРАТОР

2 4 3

Особенно часто на практике возникает задача фильтрации, заклю­

чающаяся в

возможно более точном воспроизведении управляющего

воздействия

или полезного

сигнала т (t) и в возможно более полном

подавлении

возмущающего

воздействия или помехи «(/).

Иногда эта

задача называется также задачей сглаживания, так как во многих случаях помеха n{t) является гораздо более быстро изменяющейся функцией, чем полезный сигнал, и подавление помех n(t), по суще­ ству, сводится к сглаживанию кривой, характеризующей входной

сигнал m(t)-\-n(t).

 

 

 

 

 

 

 

 

Другой практической задачей, обычно

связанной с

предыдущей,

является

задача предсказания будущих значений входного

сигнала

/«. (£+ т)

при т > 0

на

основании

предыдущей истории его

измене­

ния m(t).

Эту

задачу обычно

приходится

решать совместно

с пре­

дыдущей, т. е. определять будущее значение входного сигнала

 

при т > 0, зная

историю

изменения полезного входного сигнала m (t)

с наложенными

на него

помехами п((), т. е.

 

 

 

 

 

 

<р(/) =

7н(0 +

«(0-

 

 

(6.4)

В случае задачи фильтрации или сглаживания идеальный преобра­

зующий оператор Н (s) сводится к постоянной величине h0,

так» как

при этом

требуется,

чтобы величина

на выходе воспроизводила вСе

изменения

управляющего

воздействия

или

полезного

сигнала т (t)

на входе.

 

 

 

 

 

 

 

 

 

 

Итак,

в этом случае

выражение (6.1) сводится к виду

 

 

 

 

 

h (t) =

h0m (/).

 

 

(6.5)

В случае задачи упреждения или предсказания, решаемой незави­ симо или совместно с задачей фильтрации, в идеальных условиях мы должны иметь:

h (t) = hQm (t -j- £0),

(6.6)

и, следовательно, идеальный преобразующий оператор H{s) имеет вид

H(s) = h0et°s.

(6.7)

Возможны также случаи, когда преобразующий оператор имеет вид

Н (s) = his,

|

=

(6.8)

Я(5) = /г0е - ^ , )

и тогда требуемая операция над входным сигналом сводится соответ­ ственно к дифференцированию, интегрированию или смещению во времени на промежуток t0. Преобразующий оператор H(s) может иметь и более сложный вид, например

/ / (*) = - Т 1 + А0 + М-_

(6-9)

16*

244 СТАТИСТИЧЕСКИЕ МЕТОДЫ СИНТЕЗА ОПТИМАЛЬНЫХ СИСТЕМ [ГЛ. V!

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

3. Задачи фильтрации и упреждения

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

Вслучае задачи фильтрации ошибка

е(() = in (t) х (/)

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

Обычно оказывается, что стремление уменьшить первую соста­ вляющую приводит к увеличению второй составляющей, и наоборот.

Таким образом, задача синтеза состоит

в том,

чтобы

обеспечить

оптимальное решение вопроса, при

котором значение суммы обеих

составляющих является

минимально

возможным.

Поясним

сказанное

следующим образом.

 

 

 

 

 

 

 

 

 

Предположим, что амплитудная частотная характеристика А (ш)

рассматриваемой

системы и спектральные

плотности полезного си­

 

 

 

 

гнала Sm(u>) и помехи S n(со)

 

 

 

 

имеют

вид,

изображенный

 

 

 

 

на рис.

6.1.

 

 

 

 

 

 

 

 

Так

как

у помех более

 

 

 

 

высокочастотный спектр, чем

 

 

 

 

у полезного сигнала, то ясно,

 

 

 

 

что если

мы

сузим

полосу

 

 

 

 

пропускания

нашей

системы

 

 

 

 

так, как показано на рис. 6.1

 

 

 

 

(кривая i4j((u)), то система

 

 

 

 

почти не

будет реагировать

 

 

 

 

на помехи. Однако при этом

 

 

 

 

очень

сильно

возрастает

 

 

 

 

ошибка воспроизведения по­

 

 

 

 

лезного сигнала. Если же

увеличить полосу

пропускания

(рис.

6.1,

кривая

И2(ш)),

то весьма

существенно возрастает

влияние

помех. Очевидно,

здесь

необходим

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

Известны три способа подхода к решению задачи синтеза следя­ щих систем при наличии помех.

31 ЗАДАЧИ ФИЛЬТРАЦИИ И УПРЕЖДЕНИЯ 245

Первый и наиболее простой из них применим в том случае, когда полезный сигнал имеет гораздо более низкочастотный спектр, чем

помехи,

как

это, например,

изображено

на рис.

6.2.

 

В этом случае полоса пропускания системы должна быть выбрана

достаточно

широкой

для обеспечения требуемой

точности

воспроиз­

ведения

полезного

сигнала

и в

то

же

время достаточно

узкой для

того, чтобы

она не

реагировала

на

помехи.

 

 

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

(рис. 6.3).

Согласно этому способу форма амплитудной частотной характе­ ристики KW (у'ш) разомкнутой системы должна выбираться при низких

fitiu'.Siw)

частотах, где | KW (jw) | 1 и где сконцентрирована основная энер­ гия полезного сигнала, возможно более близкой к форме амплитудной частотной характеристики полезного сигнала, а затем быстро убывать,

но, по возможности следуя за убывающей амплитудной

характери­

стикой сигнала.

 

Таким образом, удается достигнуть, с одной стороны,

равномер­

ного ослабления влияния всего основного спектра частот

полезного

сигнала на ошибку и, с другой стороны, уменьшения, насколько удастся, влияния помех благодаря быстрому убыванию модуля пере­ даточной функции KW (у'ш).

Следует, однако,

заметить, что быстрота убывания KW ничем не

ограничивается

лишь

при |/C U /|^> 1 . При тех же значениях ш, при

которых | KW |

1,

слишком быстрое убывание этой величины может

привести к уменьшению запаса устойчивости или даже к нарушению устойчивости системы.

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

246

СТАТИСТИЧЕСКИЕ МЕТОДЫ СИНТЕЗА ОПТИМАЛЬНЫХ СИСТЕМ

[ГЛ. VI

Рассмотрим теперь

несколько частных случаев задач фильтрации

н упреждения.

что на вход системы подается полезный сиг­

1.

Предположим,

нал m{t) при отсутствии помех, т. е.

л(/‘) =

0. Тогда задача идеаль­

ного

воспроизведения

этого сигнала

имеет

тривиальное

решение,

заключающееся в том, что система должна иметь передаточную функ­ цию Ф(/<о), не зависящую от частоты на всем интервале частот 0-+-шо,

соответствующем сигналу т (/):

 

 

 

Ф (jiо) =

const,

0 < to < ш0.

 

 

2. Предположим,

что полезный

сигнал представляет собой ана­

литическую функцию

g (t).

Тогда,

по крайней мере

теоретически,

задача упреждения или предсказания значения

может

быть

решена сколь угодно

точно,

так как аналитическая

функция

одно­

значно определяется, если известно ее поведение на любом интервале. 3. Предположим, что т (/) — белый шум и Тогда задача упреждения не имеет решения, так как в этом случае значение m {t-|-т) совершенно не зависит от предыдущей истории изменения сигнала т (t). 4. Пусть m(t) и п (/) будут любыми аналитическими функциями.

Тогда, очевидно, невозможно отфильтровать т (I) от

Очевидно, что никакое выделение т(() из т (/)+■ n{t) невозможно, если только не известны какие-либо отличительные особенности m(t) и n(t). Точно так же никакое предсказание невозможно, если только m(t) не обладает каким-либо известным свойством, связы­ вающим его прошлое и будущее.

4. Исходная информация о воздействиях

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

1. Известны наиболее полные в статистическом смысле сведения о входном сигнале, т. е. известные многомерные интегральные функ­

ции

распределения

вероятности F (ср)

и F (у) входного

 

сигнала

у (t) =

у (/) -|- п (() и полезного сигнала

y(t).

 

 

 

Так, например, может быть известно, что как входной сигнал ср(/),

так и полезный сигнал у(() имеют нормальные распределения,

причем

их математические

ожидания, а также

собственные и взаимные кор­

реляционные моменты являются заданными функциями времени.

2.

Мы располагаем

наименьшими сведениями

о входном

сигнале,

т. е.

известна

лишь

принадлежность

полезного

сигнала

и

помехи

к определенному классу функций.

 

 

 

 

Так,

например,

относительно полезного сигнала может

быть из­

вестно

лишь то,

что

он

принадлежит к классу сигналов, ограничен-

51

ОБЩАЯ

ФОРМУЛИРОВКА ЗАДАЧИ

СИНТЕЗА

2 4 7

пых

по модулю, т. е.

что

 

 

 

 

 

 

I у (О I ^ Д

при

всех

t

 

или что он принадлежит к классу

полиномов:

 

 

 

у{*) = Ъ с кР,

t> о,

 

/! = 0

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

3.Имеются менее полные данные, чем в первом случае, и боле полные, чем во втором.

Так, например, могут быть известны первые два момента функций распределения входного сигнала и помехи, но неизвестны сами функ­ ции распределения.

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

щуюся информацию о входном сигнале,

вне зависимости от

степени

ее полноты. К сожалению, в настоящее

время достаточно

хорошо

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

б. Динамическая точность и функция риска. Общая формулировка задачи синтеза

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

Не ставя себе здесь задачей дать сколько-нибудь полное изложе­ ние всех известных в настоящее время методов решения сформули­ рованной выше проблемы и ограничиваясь в дальнейшем в основном лишь теми результатами, которые можно довести до практических расчетов, мы все же предварительно изложим некоторые идеи, ко­ торые при дальнейшем своем развитии, по-видимому, будут иметь большое значение для статистической динамики систем автоматиче­ ского управления :).

*) С о л о д о в н и к о в В. В., Б а т к о в А. М., Б р е д и с А. А\ М а т-

в е е в П. С.,

Методы

математической статистики и теория автоматического

управления,

Труды совещания «Автоматическое управление и вычислитель­

ная техника»

(1957 г.),

Машгиз, 1958,

2 4 8

СТАТИСТИЧЕСКИЕ МЕТОДЫ СИНТЕЗА ОПТИМАЛЬНЫХ СИСТЕМ [ГЛ. VI

Наиболее широко распространенным в настоящее время показа­ телем, применяемым для характеристики динамической точности, является среднее значение квадрата ошибки

СО

 

М [ в ° - ] = f e2w ( e , t ) d e ,

(6.10)

— СО

которое при стационарных случайных воздействиях можно также

представить

в виде

 

г

 

 

 

 

 

 

 

 

i* =

lim

f B2{t)dt,

(6.11)

 

 

 

с о

_>/

 

где_

 

 

 

 

 

 

 

e(0 = A (0 — *(0-

(6.12)

Однако

в

предыдущей

главе

нами уже был

затронут вопрос

о недостатках и ограниченной области применения

этого критерия

динамической

точности.

 

 

 

Поэтому,

имея сейчас в виду дать общую формулировку проблемы

синтеза систем автоматического управления, обеспечивающих требуе­ мую динамическую точность преобразования сигнала, рассмотрим более общий критерий динамической точности, хотя в дальнейшем будем широко пользоваться критерием среднеквадратической ошибки вследствие его относительной простоты.

Этот более общий способ оценки динамической точности основан

на введении в рассмотрение так называемой функции риска.

Функция

риска применяется в теории решающих

функций, которая,

как это

мы еще увидим

в дальнейшем, имеет для

нашей проблемы исключи­

тельно большой

интерес.

 

 

Назовем функцией веса

 

 

 

W = W(y . х) = ИГ[Ф(я), у)

(6.13)

некоторую неотрицательную функцию от полезного сигнала y(t) па входе и величины х (/) — Ф(ср) на выходе системы. Функцией риска г называется математическое ожидание функции веса:

 

 

СО

о о

 

 

г[Ф,

у ] = f

f ЧГ{Ф(9), y)dF(<?ly)dF(y).

(6.14)

 

 

— ОО

— о о

 

где F (_у) — интегральная

функция распределения вероятности полез­

ного

сигнала; F (<fly) — интегральная функция распределения вероят­

ности

входного

сигнала.

 

 

В частности, при решении задачи фильтрации функция веса может

иметь

вид

 

 

 

W { Ф (<р), y } = W [ y - < b (ср)} = W ( г ) .

( 6 . 1 5 )

61

ЗНАЧЕНИЕ ТЕОРИИ ИГР И РЕШАЮЩИХ ФУНКЦИЙ

2 4 9

Частным

случаем функции веса вида (6.15) является функция веса

в виде полинома,

содержащего

четные степени ошибки:

 

 

 

tt/(e) = Tl£2 +

T2e4_)_ . . . + Т„6*«,

(6.16)

и, наконец,

еще

более частным

случаем — функция веса

 

 

 

W (е) = f ,е2,*

(6.17)

применимая в обычном критерии среднеквадратической ошибки. После сделанных выше предварительных замечаний общая проблема

синтеза динамических систем, обеспечивающих наивысшую, в приня­ том смысле, динамическую точность преобразования управляющего воздействия (полезного‘сигнала) при наличии возмущающих воздей­ ствий (помех), может быть сформулирована следующим образом.

1. Известно, что динамическая система будет находиться под влиянием управляющего y{t) и возмущающего n(t) воздействий, ко­ торые могут быть наложены друг на друга и поступать на вход си­

стемы в виде

входного

сигнала cp(f),

 

 

 

<?= <?&. «)•

2.

Задана некоторая первичная информация о статистических свой­

ствах

сигнала

y{t) и

помехи п (t) (или входного сигнала ©(/)).

которая в наилучшем случае заключается в знании всех функций

распределения вероятности,

а в

наихудшем случае — в

знании лишь

того, что сигнал у((), помеха п (/)

и входной сигнал ф (7)

принадлежат

соответственно к известным классам функций.

 

3. Требуется, пользуясь

этими данными, найти динамическую си­

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

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

выполнения

этого этапа

решения

и Представляет собой

основной

предмет дальнейшего изложения.

Второй этап решения

состоит

в реализации

найденных

оптимальных характеристик в виде действую­

щей и удовлетворяющей

заданным

техническим условиям системы.

6. Значение теории игр и решающих функций для теории автоматического управления1)

Для выполнения первого из указанных выше этапов решения

задачи

синтеза можно использовать ряд

методов.

 

В частности, большой интерес

имеет теория игр, развитая в

ра­

ботах фон Неймана и Моргенштерна2).

В этих работах показано,

что

4)

См. ссылку на стр. 247.

t er n

О., Theory of games and econo

2)

Von Ne u ma n J., Mo r g e n s

mic behavior, Princeton, 1947.

Соседние файлы в папке книги из ГПНТБ