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

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

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

20

ВВЕДЕНИЕ

щающих воздействий), или,

как мы условимся говорить, чем выше

динамическая точность, тем, вообще говоря, лучше система автомати­

ческого

управления.

 

 

 

 

Однако, как выбрать воздействия, применительно

к

которым

вести расчет, если они, по существу, заранее полностью

непредска­

зуемы?

Здесь

возможны два

различных подхода к решению задачи.

Один из них

заключается в том, что

из всей возможной

совокуп­

ности

воздействий, которые

могут

быть приложены

к

системе

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

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

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

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

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

В

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

с

точки

зрения математической

статистики обработка как первой,

так и второй последовательности

данных

не имеет принципиального

различия.

Существенным здесь является

лишь быстрота их измене­

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

ческих машинах и,

наконец, подвергнуты анализу с целью получе­

ния окончательных выводов, причем между

получением данных и

выводами допустим

значительный промежуток

времени, в то же время

ВВЕДЕНИЕ

21

радиолокационная установка, определяющая координаты какого-либо

движущегося

объекта,

должна

производить

все указанные действия

несоизмеримо

быстрее

и

часто

без всякого

вмешательства

чело­

века.

 

 

 

 

 

 

 

 

 

Метод, который следует

применять при

расчете

системы

авто­

матического

управления

на

динамическую

точность,

существенно

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

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

А'1ожно также предположить, что относительно воздействий известно лишь то, что они принадлежат к определенному классу функций, например к классу однозначных, непрерывных функций, огра­

ниченных по

модулю.

 

 

 

 

Очевидно,

что возможны

также и другие

предположения и что

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

охватывать

все

эти

возможные случаи.

 

 

 

 

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

хотя развитие методов

анализа и

син­

теза

систем

автоматического

управления при

наличии воздействий

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

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

переходного процесса, величина перерегулирования и т. д., теряет смысл?

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

Например, ясно, что применение статистических методов лило

целесообразно, если воздействия представляют

собой импульсы или

внезапные изменения входа, хотя и

имеющие

случайную величину

и возникающие в случайные моменты

времени,

но разделенные про­

межутками, достаточными для того, чтобы система успевала до

появления каждого следующего импульса приходить в установившееся состояние.

22

ВВЕДЕНИЙ

Нет

необходимости применять статистические методы и в тех

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

из

участков несколькими типовыми функциями, так чтобы система

на

протяжении каждого из участков успевала достигнуть

установив­

шегося состояния.

 

 

Наоборот, если воздействия имеют вид непрерывно,

хаотическим

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

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

продукты

и т. д.),

и проведем аналогию между

ним и

процессом

в системе

связи (см.

рис. 1). Здесь информацией,

которая

нас инте­

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

готовой

продукции, по результатам измерения которых

можно судить

о том,

удовлетворяются или нет заданные технические

требования ‘).

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

процесса кодирования является здесь процесс превращения

исход­

ных технических требований в изменение действительных

физиче­

ских характеристик сырья, из которого изготовляется

готовый

продукт.

 

Это осуществляется, например, при помощи изменения скоростей

подачи, уставок, режима тепловой или химической обработки

и т. д.,

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

собственным шумам в виде мертвых

ходов, неточности уставок,

дрейфа,

износа частей, колебаний мощности источников питания

и т. д.

Каждой независимой настройке

соответствует некоторая вели-

J) R o t h s te in J.,

Information theory and quality control, IRE Convem

tion Record, r. 4, 1956,

стр. 3—11.

 

 

ВВЕДЕНИЕ

 

23

чина

или координата x t.

Если имеется п настроек,

то

определяемое

ими

состояние системы

изображается точкой X

{хи

х 2......... х п)

в п-мерном пространстве {А'}.

Готовая продукция определяется значениями независимых пере­ менных, необходимых для полной характеристики их свойств. Поэтому свойства готовой продукции можно изобразить точкой У (y L, у 2, . . ., у п1) в ///.-мерном пространстве {К}.

Технические условия 5 служат источником информации двух типов. Информация первого типа получается в результате выбора в про­ странстве настроек {X} начального объема (А’0}, определяющего технологию производства и первоначальный способ кодирования сооб­

щения в виде технических условий S в сигнал.

Информация второго типа получается в результате выбора в пространстве свойств готовой продукции (Г) объема {К0}, рас­ сматриваемого как объем, соответствующий ее приемлемым свой­ ствам.

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

Рассмотрим, например, такие устройства для накопления инфор­ мации, как книги, фотопластинки и т. д., — это все продукты про­ изводства, в которых информация заложена, «закодирована» в виде тех или иных изменений их структуры.

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

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

При этом в задачи приемника входят:

1)непосредственная задача приема сигнала, состоящая в опре­ делении точек У пространства свойств готовой продукции, соответ­ ствующих данной партии продуктов производства;

2)задача фильтрации, состоящая в определении того, принадле­ жат ли эти точки У к допустимому элементу объема {К0} или всей остальной части пространства {К}— {К0};

3)задача статистических выводов или решений, состоящая в том, чтобы на основании выборочных дискретных или непрерывных изме­ рений на интервале Т решить, принять или забраковать партию, характеризуемую полученными точками У (поскольку некоторые из них могут находиться в (К0), а некоторые в (К) — {К0}), или, нако­ нец, произвести ее сортировку на ряд групп в зависимости от зна­ чений тех или иных параметров.

Процесс сортировки есть не

что иное,

как процесс, который

мы имеем в системе связи, когда

несколько

сообщений посылаются

24

 

 

 

 

ВВЕДЕНИЕ

 

 

 

по одной

и той

же

линии

связи,

но

достигают

различных

пунктов

назначения.

 

 

 

 

 

 

 

 

 

Если

между

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

системами связи и

системами

автома­

тического

управления можно

провести

рассмотренную выше анало­

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

 

Основное различие заключается, пожалуй,

в

следующем.

В тех­

нических

системах

связи

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

получатель сообщения

не входят

в предмет рассмотрения.

В частности,

ие ставится

вопрос

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

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

Можно себе представить, что эти обратные связи вводятся при помощи так называемых управляющих машин (УЛ4).

Функции управляющей машины в рассматриваемом случае заклю­ чаются в следующем. Если требуемые показатели качества ие удо­ влетворяются, то она должна:

1) вычислить новый начальный объем (А''} в пространстве настроек, определяющий способ кодирования (другими словами, технологию

производства), на основании известных

зависимостей Y от

X и на­

блюдаемых {А'),

так, чтобы (К) было заключено в {К0};

 

2) произвести

преобразование {АТ0} в

{Аг/}, означающее

измене­

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

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

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

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

 

 

 

ВВЕДЕНИЕ

 

 

25

дальнейшего

изложения,

следует, прежде всего, указать

па

резуль­

таты,

полученные А. Н.

Колмогоровым, А. Я.

Хинчиным

и Н. Ви­

нером.

 

 

 

 

 

 

А.

Я.

Хиичин *)

впервые дал математически строгую и в то

время

физически ясную

теорию незатухающих

статистически

стацио­

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

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

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

лиза

и

синтеза

динамических

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

Ссылки

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

изложения.

 

 

 

 

 

*)

Х и н чин

А. Я., Теория

корреляции стационарных случайных про­

цессов, Успехи матем. наук, вып. 5, 1938.

 

стационарны

3)

 

К о л м о г о р о в А. Н.,

Интерполяция и экстраполяция

случайных последовательностей,

Изв. АН СССР,

сер. матем.,

№ 5, 1941.

3)

W i e n e r

N.,

Extrapolation, interpolation and

smoothing of

stationary

time series, J. Wiley,

N. Y., 1949.

 

 

 

Г Л А В А I

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

1. Введение

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

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

2. Дифференциальные уравнения линейной динамической системы

Дифференциальные

уравнения линейной

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

системы

с k степенями свободы имеют вид

 

 

 

 

 

О п ( р ) * 1 +

о12(р)л:2 + . . .

+

а,* ( / ? ) * * =

/ + ) .

 

fl2 l ( / 0 * l+ fl22(/>)*2+

•••

+

02ft(/>)** =

/ 2 (0.

(1 п

0*1(л) *1 +

0*2(Л) * 2 +

 

+О **(р)** =

/*(0- .

 

В уравнениях (1.1)

через х ь

х 2,

. ... ** обозначены обобщенные

координаты системы,

через /j(/),

/ 2(/), . . . . /* (0 — возмущающие воз­

действия и через аа (р) — операторы вида

 

 

 

 

о<; (Р) = /пир2 +

Ьар +

сп.

 

( 1.2)

21

УРАВНЕНИЯ ЛИНЕЙНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЫ

27

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

нулями) и р = .

Уравнения (1.1) вместе с начальными условиями полностью опре­ деляют поведение динамической системы.

Всякая линейная динамическая система подчиняется принципу супер­ позиции, согласно которому эффект нескольких приложенных к системе возмущающих воздействий /Д^), Л (0 . •••• Д (0 равен сумме эффек­ тов каждого из этих воздействий в отдельности. Поэтому без потери

общности

дальнейших

выводов

можно положить в (1.1)

все

воздей­

ствия, за

исключением

одного,

например / 1(0 =

/(0>

равными нулю.

Тогда

уравнения

(1.1) примут

вид

 

 

 

 

 

 

 

a u (p )x i +

aiz(P)x 2 +

••• 4- alk (р) x k =

/

(0,

}

 

 

a2i( p ) x l -{-a22( p ) x 2-dr

... -i-a2k(p)xk =

0.

I

(1.3)

 

 

 

 

 

 

 

 

 

 

}

 

........................................................................................

 

 

 

•••4" akk (p) x h —

 

 

I

 

 

aki(p) х 1 ~\~акЛР)хг~\~

 

 

J

 

Если нас интересует изменение во времени только одной какой-

нибудь из переменных х ь х 2, . . . .

x k, например

— х, то, исключая

из уравнений (1.3)

переменные

х 2, x s..........

х к,

мы

 

получим

одно

дифференциальное уравнение л-го

порядка

(л ^

2&):

 

 

 

 

dnx^ .

dn 1х ,

.

dx .

апх —

dt‘

 

 

 

 

= Ь„ dmf(t) Ч- bm_i -j/h-1fit)

- b ^ + b0f ( t ) ,

или

dtm

dr

 

dt

D( p ) x =

M(p)f(t ),

 

где

D (p ) = anp n +

an^ p n- 14-

... 4 - alP -4- a0,

 

причем

^ ( P ) ^ b mpm-\~bm. 1pm- i +

. . . 4 - 6 lp-|-(b0,

 

 

 

 

 

 

«п (Р) «12 (Р) • ■«1А (Р)

 

£>00 =

#21 (Р) ^22(/0

• • «2* (Р)

 

 

 

 

 

 

«И (Р) «А2 (р)

•• • akk (р)

(1.4)

(1.5)

( 1.6)

(1-7)-

a оператор M(p) представляет собой минор определителя (1.7), полу­ чающийся при вычеркивании того из его столбцов, который отно­ сится к интересующей нас переменной, и той из строк, которая соответствует в (1.3) уравнению с правой частью. ■ ■

28

НЕКОТОРЫЕ ВОПРОСЫ АНАЛИЗА

[ГЛ. I

Так

например, в рассматриваемом случае:

 

 

 

а22(р)

... а2к (р)

 

 

М { р ) =

а32(р)

... аш.(р)

( 1.8)

 

 

aki (р) ■■■ 4k(P)

 

Из сказанного легко видеть, что 1114^ 11.

Решение или общий интеграл уравнения (1.4) состоит из двух частей. Первая часть решения есть общий интеграл уравнения (1.4) без правой части, т. е. решение однородного дифференциального уравнения

(апРп an-iP" ~ 1 ••• - M iP - M o )* = 0. (1.9)

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

Рассмотрим вначале решение однородного уравнения (1.9). Оно определяет свободные колебания системы, т. е. колебания, которые система, выведенная из состояния покоя некоторым воздействием, совершает после того, как это воздействие исчезло.

Будем искать решение уравнения (1.9) в виде

х = Ке>(.

(1.10)

Подставляя (1.10) в (1.9), получим:

( а пк П ~ \ г а п - ^ ' П 1_Н ••• + fli^ - h ao)^eU = 0-

(1-11)

откуда следует, что выражение (1.10) действительно является реше­ нием уравнения (1.9) при условии, что величина X является корнем характеристического уравнения

D Q-) = ап>-п-{-ап_i>" 1

а, Х-| -д0= 0.

 

(1-12)

Так как уравнение (1.12) имеет п корней Xt, X,.........Хл,

то общее

решение уравнения (1.9) (в предположении,

что все корни

различны)

может быть представлено в виде

 

 

 

*со {t) = K / lt + K 2eh ‘ + .. .

 

 

(1.13)

где значения постоянных Kv К2, . . . . Кп зависят от начальных усло­ вий. Если среди корней имеются кратные корни, то корню )ч, имею­ щему кратность V , соответствует частное решение уравнения (1.19) вида

(Kn + Kat-\- ■• •

(1.14)

Из (1.13), а также из (1.14) очевидно, что функция x CB(t), пред­ ставляющая собой решение уравнения (1.9), будет удовлетворять условию

lim x CB(t) = 0,

(1.15)

/ > СО

 

3]

 

ВЫНУЖДЕННЫЕ

КОЛЕБАНИЯ

29

т. е.

свободные

колебания системы

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

мени,

и система

будет устойчивой в

том и только

том случае, если

все

корни

характеристического уравнения (1.12)

имеют отрицатель­

ные

вещественные части.

 

 

 

3.

Вынужденные

колебания.

Частотные характеристики

 

Перейдем

теперь к

отысканию частного решения уравнения (1.4)

для

того случая, когда воздействие

f (t) представляет собой гармо­

ническую функцию времени, т. е.

может быть представлено в виде

линейной комбинации функций sin at

и cos at или

в виде

 

 

 

 

] it) /o cos(ш^Ч~То)-

(1-16)

где / 0 — амплитуда, ш — угловая частота, ср0 — фаза.

Если

предположить,

что ср0 — 0,

то

 

 

 

 

 

f ( t ) = f 0cosat.

(1.17)

Гармоническое воздействие вида (1.17) может рассматриваться как сумма двух экспоненциальных воздействий, т. е. вместо (1.17) мы

можем

записать:

 

 

 

 

 

 

 

 

 

 

 

 

/

(0 =

/о cos at = А

е*0* +

А

 

(1.18)

где j =

Y

— 1 •

рассмотрение

экспоненциальных воздействий

весьма

Заметим,

что

удобно

с математической

точки

зрения в

силу того, что как

произ­

водная, так и интеграл

от экспоненциальной

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

собой по-прежнему экспоненциальную функцию.

 

Найдем эффект, создаваемый каждым из экспоненциальных воз­

действий

в отдельности.

Положив

 

 

 

 

 

 

 

 

 

 

 

/ (0 = 4г eJ<at,

 

 

(1.19)

будем искать

частное

решение уравнения (1.4)

в виде

 

 

 

 

 

 

* ib = 4 “ Ф ( »

e/W.

 

(1.20)

Подставляя (1.19) и (1.20) в (1.4), получим:

 

 

К О ) " +

• • • + a j a

 

а0] Ф ( »

 

=

 

 

 

и, следовательно,

 

 

= [bn ( M m +

• • • +

b j a -f- й0]

(1.21)

 

 

 

 

 

 

 

 

 

 

 

 

Ф

 

_

Ьт (7ю)’" +

• • • +

b ] j a

- f b0

( 1.22)

 

 

 

 

 

 

ап (У“) + ...

+ a-[ju>-j- а0

 

 

 

 

 

 

 

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