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

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

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

3 2 0

СИНТЕЗ ПО

МИНИМУМУ СРЕДНЕ КВАДРАТИЧЕСКОЙ ОШИБКИ

[ГЛ. VII

Принимая в (7.244) для сокращения письма следующие обо­

значения:

 

 

 

 

 

напишем:

^ 9 (тг-

Rij’

(Ti)

^ O^j) kj,

 

 

П

 

 

 

 

 

 

 

(7.245)

 

Rit4i =

2 R-ijkj,

t — 1,

2......... п.

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

Замена интегрального уравнения системой алгебраических урав­ нений, т. е. замена интеграла конечной суммой, реализуется в мате­

матической машине тем, что

функции,

входящие в уравнение, пред­

ставляются

ступенчатым напряжением,

состоящим

из

достаточно

большого числа

ступенек.

Таким

образом,

система

алгебраических

уравнений (7.245) решается в машине итеративным методом.

Известны

две

основные

разновидности

итеративного

процесса:

а) простая итерация и б)

итерация

Гаусса — Зейделя.

Итерация

Гаусса — Зейделя

характеризуется

тем, что «обновление» результатов

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

функции, а не после того,

как

вычислено

данное приближение

для

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

это

имеет место при простой итерации.

Обычно как в теоретических

работах,

так и при создании

при­

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

предпочтение итерации Гаусса — Зейделя,

так как

это позволяет

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

В линейной алгебре доказано положение

о

том,

что итеративный

процесс Гаусса — Зейделя для системы (7.245)

с симметричной поло­

жительно определенной матрицей сходится к решению при произ­

вольном начальном

приближении.

 

 

Итеративный

 

процесс Зейделя для системы (7.245) строится в сле­

дующей форме:

 

 

 

 

п

 

Rh<pl

 

1

 

1

 

 

 

■2 v r 4 I — 1, 2, . . . . п.

Мш)

Ж

7=1

Ж

 

 

 

 

y=/+i

(7.246)

 

 

 

 

 

 

Здесь индекс

в

скобках

наверху означает номер

итеративного

шага. В качестве

начального

приближения используется

нуль:

 

 

 

#У = 0;

/ =

1. 2......... а.

(7.247)

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

18] ПРИБОР ДЛЯ ОПРЕДЕЛЕНИЯ ОПТИМАЛЬНОЙ ИМПУЛЬСНОЙ ПЕРЕХ. ФУНКЦИИ 321

Функциональная схема машины дана на рис. 7.17.

 

Машина состоит из двух основных частей:

 

1)

вычислителя интегральных преобразований (В И П )

и

2)

устройства памяти и «компрессии» функций (УПКФ).

ВИП вычисляет интегральные преобразования в следующей форме:

 

г

 

 

F (х) = f ( x ) J k ( x s) cp (s) ds.

(7.248)

 

о

 

Все функции в правой части (7.248) предполагаются известными функциями. Так как в интегральном уравнении (7.14) под знаком ин­ теграла имеется неизвестная функция, то вместо нее под интеграл

УПКФ

Рис. 7.17.

вводится неточное, грубое решение, которое в ходе итеративного процесса приближается к точному. Один итеративный шаг, дающий одно приближение для искомой функции, требует тридцатидвухкрат­ ного интегрирования, поскольку диапазон изменения обеих незави­ симых переменных разбивается на 32 равные части. Вычисленное значение искомой функции k (t) запоминается в блоке УПКФ и за­ тем вводится под знак интеграла при вычислении следующего при­ ближения искомой функции. Таким образом, ВИП и УПКФ обра­ зуют замкнутую систему и этим обеспечивается возможность осу­ ществления итеративного процесса. По окончании итеративного процесса продолжаются вычисления по тому же алгоритму, но только записываемая в УПКФ функция имеет тот же вид, что и считываемая.

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

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

322 СИНТЕЗ ПО МИНИМУМУ СРЕДНЕКВАДРАТИЧЕСКОЙ ОШИБКИ [гл. VII

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

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

Блок памяти основан на применении запоминающего устройства, используемого во внутренних накопителях цифровых машин, так на­ зываемого потенцпалоскопа. Потенциалоскоп представляет собой электронно-лучевую трубку, у которой люминофор заменен сигналь­ ной пластиной. Элементами памяти служат заряды элементарных конденсаторов, образуемых на поверхности этой сигнальной пла­ стины. Каждый конденсатор называется ячейкой; общее число ячеек равно 1024 (32 X 32). Информация записывается в ячейках в двоич­ ном коде (единица или нуль). Каждый столбец (32 ячейки) отводится для запоминания одной ординаты записываемой функции. В соответ­

ствии

с величиной

этой ординаты в одной какой-то ячейке записы­

вается

единица, а

в остальных — нули.

Если мысленно

представить

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

то мы уви­

дим график записанной функции.

 

 

При считывании записанной функции

осуществляется

последова­

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

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

Дадим краткую характеристику других узлов функциональной схемы рис. 7.17. Генератор -правой части ГПЧ и генератор ядра ГЯ

18] ПРИБОР ДЛЯ ОПРЕДЕЛЕНИЯ ОПТИМАЛЬНОЙ ИМПУЛЬСНОЙ ПЕРЕХ. ФУНКЦИИ 323

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

Холостой

I ход Считывшше\ К х ш Сч \ X х t Сч | X a t Сч | Л х. \

АААААг Строчная развертка потенциалоскопа

Кадровая развертка потемциалошпо

Рис. 7.18.

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

расположением во времени управляющих

импульсов,

формируемых

в блоке управляющих импульсов {БУИ).

Идея

получения необходи­

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

 

 

 

 

Здесь для большей наглядности принято

п — 4,

а

не п — 32.

Цикл работы машины, за время которого вычисляется

одна орди­

ната искомой функции, состоит из двух равных по длительности

частей: холостого хода и считывания. Период

кадровой развертки

ГЯ несколько больше, чем длительность цикла.

Считываемая

часть

кривой ядра показана на рис. 7.18 утолщенной линией.

 

Таким

образом, задача создания

генератора

функции двух

аргу­

ментов —

ядра уравнения — решена

без построения какого-либо спе­

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

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

21*

18] ПРИБОР ДЛЯ ОПРЕДЕЛЕНИЯ ОПТИМАЛЬНОЙ ИМПУЛЬСНОЙ ПЕРЕХ. ФУНКЦИИ 3 2 5

Остальные узлы вычислительной машины состоят из элементов,

применяемых

в вычислительной технике. Устройство преобразования

диагональных

членов

матрицы УПМ служит для преобразования

в нуль

диагональных

членов ядра

и состоит из одного вибратора,

схемы

совпадений и антисовпадений.

В И П состоит из. время-импульс-

ного множительного

устройства, интегрирующего, суммирующего и

накопительного блоков. Блок управляющих импульсов состоит из триггеров, одновибраторов и счетчиков. Устройство автоматической установки знака решения АУЗ служит для того, чтобы на машине можно было решать уравнения, содержащие знакопеременные функ­ ции. Дело в том, что примененные в машине функциональные гене­ раторы, устройство памяти и множительное устройство предназна» чены для работы только со знакопостоянными функциями. Поэтому к знакопеременным функциям прибавляется такая постоянная вели­ чина, которая превращает их в знакопостоянные. Величина этой кон­ станты для ядра и правой части учитывается при установке масшта­ бов, а величина константы для решения, меняющаяся в течение итеративного процесса, учитывается автоматически при помощи АУЗ. Устройство АУЗ представляет собой накопитель.

Таковы вкратце основные черты рассмотренной вычислительной машины. Прибор, общий вид которого приведен на рис. 7.19, про­ изводит вычисления с погрешностью около 10% от максимального значения. Время решения, интегрального уравнения— 10-т-15 сек, не считая установки масштабов. Для сравнения укажем, что время

решения подобных уравнений на электрической клавишной

машине

составляет

8 -^ 1 0

рабочих

дней.

 

 

Следует

отметить,

что

описываемая машина может

быть

исполь­

зована не только

для

решения уравнения (7.14), но

также

для це­

лого ряда других задач. В качестве примеров укажем на следующие возможные применения:

а) решение интегральных уравнений Фредгольма и Вольтерра первого и второго рода с ядром типа свертки по эксперименталь­ ным данным;

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

в) решение интегральных уравнений для функций, заданных ана­

литически.

В

этом

случае

требуется

дополнительное оборудование

для реализации этих функций.

Так, например, для

реализации ядра

К ( ш,

t) — cos

требуется

три дополнительно решающих усилителя

и два

умножителя;

 

 

 

 

 

 

г)

определение

импульсных

переходных функций

объектов регу­

лирования

в

процессе

их

нормальной

эксплуатации

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

с методом,

изложенным

в

гл.

V.

 

 

Г Л А В А VIII

СИНТЕЗ ОПТИМАЛЬНЫХ СИСТЕМ С КОНЕЧНОЙ «ПАМЯТЬЮ»

1.Введение

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

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

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

Таким образом, при синтезе оптимальных систем обычно необхо­

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

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

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

1]

ВВЕДЕНИЕ

327

 

 

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

Остановимся на этих требованиях несколько более подробно.

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

ской осуществимости

и представлять

собой

функцию,

отличную от

нуля

при t > 0.

Эго

означает, что система реагирует

на всю про­

шлую

историю изменения сигнала и помех

на входе,

т.

е. что

система

обладает «бесконечной памятью».

 

 

 

 

 

 

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

обычно

требуется,

чтобы

система

реагировала лишь на часть входных

данных,

ограниченных

некото­

рым конечным промежутком времени

Т.

Например, может быть не­

обходимо, чтобы

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

t = 0, в момент

времени t = Т

была

бы

уже

готова к оптимальному

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

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

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

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

гочлен от t, т. е. что

 

=

<8Л>

 

<7 =0

В зависимости от информации, которая имеется об этой составляю^ щей, будем рассматривать три случая:

1) о коэффициентах kq ничего не известно '); - ;

') Z a d e h L.

A.,

R a g a z z l n i I.

R.,

An extension of Wiener's theory of

prediction, Journal

of

Applied Physics,

t .

2,. July 195.0, стр. 645.— 655;

3 2 8

СИНТЕЗ

ОПТИМАЛЬНЫХ СИСТЕМ

С

КОНЕЧНОЙ «ПАМЯТЬЮ» [ г л . VIII

2)

значения

коэффициентов

k q заданы ');

 

 

3)

коэффициенты k q случайны и

известна

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

функ­

ция 2*)

 

л

г

 

 

 

 

 

 

 

S

2 т У

W

 

 

 

 

 

/и о /= о

 

 

 

 

 

где

 

1и =

М [ к , к , 1

 

 

 

т е.

известна матрица ||-у/у.||

(/, j =

0,

1......... г) корреляционных

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

случайных величин

k t

и £..

 

 

В первом случае мы не располагаем

никакой информацией

о ко­

эффициентах kq и при расчете

оптимальной

системы будем

требо­

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

полинома от t. Заметим, что

в этом

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

к

системе

с астатизмом

(г-|-1)-го

порядка,

реализация которой

даже при

г = 1 может

быть связана

с

существенными практическими

трудно­

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

ного класса, т. е. любого полинома вида

(8.1),

является

очень же­

стким

требованием,

которое, вообще

говоря,

может

привести

к завышенным значениям среднеквадратической ошибки.

 

Во

втором случае

предполагается, что функция g(t) представляет

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

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

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

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

Предположим, что на вход рассматриваемой системы поданы управляющее y{t) и возмущающее g{t) воздействия. Управляющее воздействие является суммой двух составляющих:

_______________

'

.У(9 = £ 4 0 + «(*).

(8 -2)

1) С о л о д о в н и к о в

В. В.,

М а т в е е в П. С., Синтез корректирующих

устройств следящих систем при

наличии помех по заданным требованиям

к динамической точности, Автоматика и телемеханика,

т. XVI, № 3, 1955.

2) С е м е н о в

В. М.,

К теории экстраполирования

случайных процессов,

Сборник научных трудовВВИА

им. Жуковского, 1954,

т. 1, стр. 67—95.

2] ПОСТАНОВКА ЗАДАЧИ 3 2 9

где g(t) — функция времени, заданная своим аналитическим выраже­ нием. Для простоты предполагается, что g(t) имеет вид полинома

(8.1) с

известными

значениями

коэффициентов k r а т(() — стацио­

нарная

случайная

функция с

заданной корреляционной функцией

Rm(x) или спектральной плотностью S m(ш).

Возмущающее воздействие (помеха) п (t) также предполагается ста­ ционарной случайной функцией времени с заданной корреляционной функцией Rn (т) или спектральной плотностью S n (ш). Для простоты будем считать, что функции m{t) и п{£) взаимно не коррелированы.

Обозначим величину на

выходе

системы через x(t), а ее импульсную

переходную функцию через k(t).

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

что

 

 

А (0 =

0,

f < 0 - ,

 

(8.3)

 

А (0 =

0.

 

 

(8.4)

•Заметим, что выражение (8.3)

представляет

собой

условие физи­

ческой осуществимости,

а выражение (8.4) — условие

качества, со­

гласно которому время переходного процесса, вызванного управляю­

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

в виде

единичной ступенчатой

функции 0,

не

должно превышать величины Т.

величина

 

Так же как и

раньше,

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

y{t)

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

h{f) =

H( p) y{ 0

(8.5)

или

ОО

 

 

 

 

 

h(t) =

j" у (t — т)у. (х) d%,

(8.6)

где

 

ОО

 

 

 

 

=

-^ 7

f Я ( У « ) «*“' * » .

(8.7)

 

 

—оо

 

В действительности же на выходе мы будем иметь

величину x(t),

а не h (I).

 

 

 

Учитывая (8.3) и (8.4) для величины x(t) на выходе, можно на­

писать:

т

* ( 0 = / [£(* — т)-}-/и (* — i ) + n(t — x ) ] k ( x ) d x .

(8.8)

о

 

!) Условие (8.4) относится к импульсной переходной функции, но оче­ видно, что если оно имеет место, то переходный процесс закончится в те­ чение времени Т и в случае воздействия в виде ступенчатой функции.

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