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

книги из ГПНТБ / Климов, В. А. Некоторые прикладные методы анализа и синтеза сложных автоматических систем с использованием ЦВМ

.pdf
Скачиваний:
7
Добавлен:
22.10.2023
Размер:
16.72 Mб
Скачать

384

^ i , m a x j L — 1 • S ,

(6.-115)

Массивы (6.II4) и (6Д15) должны соответствовать диапазонам возможных значений выбираемых параметров, и, следовательно,эти 'значения должны удовлетворять условиям

X i ^ X i , m a x » i = ^ S .

( 6 . П 6 )

При осуществлений оптимизации со случайным поиском опти­ мальных значений выбираемых параметров должны в качестве вспо­ могательных использоваться алгоритмы (программы) формирования случайных чисел [23,44] . В соответствии с (6.ИЗ) каждый раз

в расчетах должно осуществляться формирование сочетания s случайных чисел.

В применявшихся алгоритмах вместо чисел (6.II3) вырабаты­

вались числа x-L, удовлетворяющие условиям

О - х^ , i = 0 f s . (6.II7)

Затем переход от х ^ к числам x-Lосуществлялся по зависимости

x i ~ х L

m a x ~ ^ m i n ) + х m i n »

 

 

L - / f S .

(6.II8)

 

 

Легко убедиться, что при использовании (6.118) числа (6.ИЗ)

удовлетворяют (6.II6), если числа

удовлетворяют (6.II7). -

Применительно к задачам оптимизации со случайным поиском

значений выбираемых параметров,

которые и используются в рас­

сматриваемых расчетных схемах,

оказалось целесообразным исполь­

зовать объединенный алгоритм определения коэффициентов харак­

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

Целесообразность составления такого алгоритма объясняется

двумя обстоятельствами.

Первое обстоятельство заключается в повышении достоверно­

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

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

385

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

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

нии этой возможности будет увеличено при данном времени счета число рассмотренных сочетаний чисел x-L и, следовательно, увели­ чена достоверность результата. .Вели же число вариантов, кото­ рое должно исследоваться, предписано заранее, то указанный выше подход позволяет со1фатить потребное время счета.

Кроме того, применение объединенного алгоритма определения коэффициентов характеристического уравнения и оценки запасов

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

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

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

Рассматриваемый объединенный алгоритм был сформирован пу­ тем включения в алгоритм вычисления коэффициентов характеристиче'ского уравнения (см.(6.5)3 процедур проверки положительности

коэффициентов характеристического уравнения и вычисления пара­ метров /п» [см. (4.140)] , с помощью которых осуществляется оцен­ ка запасов устойчивости.

Вычисления по алгоритму осуществляются до конца (для вы­ бранного сочетания случайных чисел x-L ), если оказываются по­ ложительными все коэффициенты характеристического уравнения и все параметры m s меньше некоторого конкретного числа т и с х . Это число назначается различным образом для разных расчетных схем. Вели какой-либо- коэффициент характеристического уравне­ ния a-L не будет положительным, то сразу же расчеты по объеди­ ненному алгоритму прерываются и осуществляется переход к выбо­ ру нового сочетания случайных чисел x L . Такой же переход осуществляется, если где-либо наступит условие m-t >mut .

Строка алгоритма (6.5), соответствующая вычислению коэф­ фициента А, , применительно к объединенному алгоритму запи­

сывается следующим образом:

Р, = Р J SP Pi =cl i ’

(6 .II9 )

386

 

если

q,*: 0

,

то

Cf = Р1 - q f E

 

 

 

 

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

> (6.II9)

 

чисел

XI .

 

 

 

 

 

 

 

J

Как видно, данная строка алгоритма (6.5)

дополняется толь­

ко проверкой положительности коэффициента

А,

. Это видно,

если учитывать (6.8).

 

 

 

 

 

 

 

Остальные строки указанного алгоритма дополняются еще

сравнением соответствующих параметров

mg

с числом

т исх .

Строка с номером

J

алгоритма

(6.5) применительно к объе­

диненному алгоритму получает вид

 

 

 

 

 

а)

если

^ .< 0

,

то

m gfg = f

( Ах ) ,

 

 

 

 

 

&

 

 

х = J ~ 3 + j

или

 

(6.120)

 

иначе переход

к выбору случайного

сочетания

 

 

чисел

Х[ ;

 

 

 

 

 

 

 

 

б)

если гпэ ^ т исх,

то

C^=P-q^E ; иначе переход

 

к выбору случайного сочетания чисел x L- .

 

Вели для определения последних коэффициентов характеристи­ ческого уравнения используется второе применение процедур Д.К.Фаддеева, то процедуры (6.П9) и (6.120) используются ана­

логичным образом с заменой коэффициентов А^ на козффициен- , ты A j .

Однако при этом вместо параметров т3 (4.140) должны использоваться величины , которые соответствуют более ши­

роким областям, чем параметры т 3 . Имеем

а.,

в Щ ( ^ У о 4 + 1 0 о Щ ( Щ & f

~ - <r»<*-2L a j ' . г V

aP - i l \ а У -г / \

ay - i l \ Д /» -z /J

3,i

6 а у - з

 

j ' = 4- f A .

387

Зависимости (6.I2I) получены путем гйэдстаяовки в зависимо­ стях (4.140),.если их записать применительно к уравнению (2.63), вместо коэффициентов Aj с использованием (6.13) коэф­

фициентов уравнения (6.12) и с последующим переходом к коэф­ фициентам, соответствующим первой форме записи уравнения (6.12).

Так, при получении зависимости для т ., (6.I2I)

указанные

подстановки проводились в выражении для т ^ п_у, т.е.

использо­

валось условие

 

У = тз, n-jW

 

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

гг>2 л соответствующих ^ =4тД. Для т^г и т^3использовалось усло-

вие г2’3) и ПРИ преобразованиях в их выражениях по­ лучились коэффициенты,которые вычисляются позже при применении процедур Д.К.Фаддеева.Они были приняты равными нулю и тогда по­

лучились зависимости,записанные в (6.121).Из-за использования этих зависимостей рабочие области, соответствующие параметрам

п, и получаются более широкими.

Необходимо вместе с тем иметь в виду, что при двойном при­

менении процедур Д.К.Фаддеева проверка условия /т^«т^^вместо контроля /л .с/л приводит к тому, что для некоторых параметров

Зу и с х '

/7?j последнее условие не проверяется. Поэтому необходимо этот пробел восполнять с обязательной проверкой условия /л .< mUCJC после завершения всех процедур Д.К.Фаддеева.

Теперь рассмотрим алгоритмы собственно расчетных схем оп­ тимизации, которые состоят из рассмотренного выше объединенно­

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

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

Х[ и особенностями в назначении массивов (6.II4) и (6.II5).

Для первой схемы расчетов, соответствующей чисто случайно­

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

При первом приеме числотцсх остается неизменной для всего расчета, для всех сочетаний случайных чисел cc-L . При такой

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

388

системы меньше фиксированного значения т исх

и удовлетворяют­

ся остальные требования к свойствам системы.

 

 

 

При втором приеме назначения числа тцсх

происходит его

изменение каждый раз,

когда

 

 

 

 

 

 

 

 

 

m i * m Uc x »

i = Z т

p

 

 

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

мы.

При этом тисх присваивается значение,

равное наибольшему

из

значений

т j

,

полученных при данном сочетании чисел х; .

 

При данной схеме расчетов определяются такие сочетания

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

x-L ,

которые соответствую^ последовательному

уменьшению наибольших параметров

/п^

. Это означает следующее.

 

Пусть при некотором сочетании случайных чисел л -

должно

произойти изменение ти с х

. Тогда из всех параметров

т 3 для

данного сочетания

x t-

выделяется наибольший и параметру тис

присваивается его

значение. Так делается и дальше. Йели все

значения

,

выделенные таким образом,

т.е. все т цсх .рас­

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

то это будет убывающий ряд чисел.

 

Таким образом,

при втором приеме назначения числа т исх

выделяются сочетания случайных чисел х ; ,

соответствующие умень­

шению для систем наибольших значений

т^ ,

т.е. соответствую­

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

Как при первом, так при втором приеме назначения числа т цсх для первой схемы расчетов остаются неизменными массивы (6.II4) и (6.II5), назначаемые вначале как исходные данные.

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

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

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

исполь­

зуется массив х . ,

который мы обозначим

 

 

х * , i =

1 -г S .

(6.122)

Шссивы (6.II4)

и (6.II5)

задаются по соотношениям

 

x i, m in ~ П m in x i »

(6.123)

389

 

x i, m a x = n m a x x i

(6.124)

В (6.123)

и (6.124) коэффициенты предполагаются положительными

числами,

причем пт1П‘с 1, а п тах> 1 . Выбор коэффициентов п т -

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

 

Здесь предполагается, что уравнения систем представлены

таким образом, что выбираемые параметры

jc-l могут быть только

положительными числами.

Далее особенность формирования массивов (6.123) и (6.124) состоит в том, что при изменении числа т исх одновременно из­ меняется и массив (6.122), за который принимается последнее

сочетание случайных чисел х [ , соответствующее новому числу/??^ Таким образом, вследствие изменения вместе с тисх массива

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

ции.

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

включая и такие варианты, когда требования по одним параметрам

тЪ будут использоваться как ограничения, а по другим будет

осуществляться оптишзация. Все эти и возможные другие вариан­ ты алгоритмов будут вносить особенности в содержание процедур.

§ 7. АЛГОРИТМЫ ОПРЕДЕЛЕНИЯ ПЕРЕХОДНЫХ ПРОЦЕССОВ С

ПОСЛЕДОВАТЕЛЬНЫМ ИСКЛЮЧЕНИЕМ БЬСТРОПРОТЕКАЩИХ СОСТАВЛЯЮЩИХ

Общие положения до целесообразности использования для си­ стем высоких порядков излагаемых выше алгоритмов рассмотрены в первой главе (§ 6). К материалам этой главы необходимо еще

добавить следующее.

В проблеме использования ЦВМ для определения переходных процессов затруднения возникают не только в связи с большим потребным временем счета. Эти затруднения часто связаны так­

же с выбором шага интегрирования и обеспечением устойчивости

счета {46, 68 и дрГ].

390

Перечисленные затруднения также отмечались выше в предисло­

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

ной структуры замещающих систем уравнений (см.рис.6.8 и 6.9),

Рис.6.8

по которым осуществляется интегрирование. Система уравнений,

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

Затруднения, связанные с выбором шага интегрирования,уда­

ется преодолеть путем использования результатов метода эффек-

391

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

нения в использовании ЦВМ для определения переходных процессов, связанные с большим потребным временем счета, удается преодо­

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

Рис. 6.9

392

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

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

коэффициенты уравнения системы (1 .1'■ )- Затем конструируются замещающие структурные схемы, о'которых было сказано выше и которые показаны на рис.6.8 и 6.9. Первая схема (рис.6.8) ис­ пользуется для построения процессов при скачкообразных входных

воздействиях, а вторая (рис.6.9) - при произвольных законах из­ менения воздействий. Последняя схема, следовательно, соответ­

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

ся для второй схемы т ~ п - 1 . Постоянные времени для апериоди­ ческих составляющих вычисляются по формуле

5

=

(6.125)

a n-L,

для колебательных составляющих по зависимости

Т.

-

(6.126)

 

 

a n-L-f

т .'~

(6.127)

J

~

a n-i:

В (6.125) - (6.127) под

lj ,

как и ранее, понимается сум­

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

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

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

турной схеме. Эти уравнения записываются следующим образом для принятой в качестве примера схемы (рис.6.8):

 

 

393

 

х

i

1

(

Ьт

 

l

т;

I

а п

 

 

 

 

 

l /

i

 

>

х , = У ( х , ~ х *) >

 

 

'i

 

 

 

* г = у

( * , - • * * ) * ,

*

' г

 

• •

(6.128)

• •

 

« •

т (^л-г'^л)5

Л""7

= 7^ (^я-f ~ ^-л)>

Х>, = 7^ ( х л ~ Х_х).

Здесь

X , = X .

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

чается уравнение

( а о Р П+ а , р П~’+ - ■■ + а п. , р +

ап ) х = bm f ,

(6.129)

которое отличается от ( I .I 1)

тем,

что в цравой части (6.129)

есть только одно слагаемое.

Остальные слагаемые полинома пра­

вой части уравнения (I . I*) учитываются через начальные условия

для координат системы (6.128).

Для выполнения интегрирования по системе (6.128) необхо­ димо иметь кроме значений постоянных времени [см.(6.125) +

(6.127)3 еще данные по начальным условиям для координат этой

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

JC;(+0)

=

;

 

(6.130)

 

а п -1 Г 1

 

 

а^(+0)= uJ H - Lf z

X : ( + 0 )

=

--------- i —

(6.I3I)

 

i / л \

 

P m -ij-l

 

a n - i j - Z

 

 

a n-if 1

 

Для того чтобы убедиться в справедливости формул (6.130) +

(6.I3I), достаточно сравнить начальные условия для координаты х ( х л)[значения х ( + 0) , х(+ 0 X (+0)],определяемые в соот-

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