Динамика нестационарных систем с равномерно изменяющимися во времени коэффициентами Текст научной статьи по специальности «Математика»

ISSN Ü868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2ÜÜ2, том 12, № 2, c. 83-92

ОРИГИНАЛЬНЫЕ СТАТЬИ

УДК 629.7.05.001(02)

© С. В. Богословский, В. С. Богословский

ДИНАМИКА НЕСТАЦИОНАРНЫХ СИСТЕМ С РАВНОМЕРНО ИЗМЕНЯЮЩИМИСЯ ВО ВРЕМЕНИ КОЭФФИЦИЕНТАМИ

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

ПОСТАНОВКА ЗАДАЧИ

Наибольшее распространение линейные системы обыкновенных дифференциальных уравнений с линейно изменяющимися во времени коэффициентами получили при моделировании динамики нестационарных систем управления (НСУ) высокоскоростными динамическими объектами, называемых "нестационарными системами управления с особой точкой" [1] или "системами управления конечным положением" [2].

Математический аппарат, пригодный для анализа и синтеза НСУ, основан на использовании результатов, полученных еще в конце 18 века П. Лапласом. Он был опубликован в 1957-1958 гг. немецкими математиками Г. Людвигом и Х. Роллни-ком [1, 2] и получил дальнейшее развитие в работах профессора А.Т. Барабанова [3, 4]. Однако в этих работах отсутствуют методы получения аналитических зависимостей фазовых координат динамического объекта от параметров НСУ на всем интервале интегрирования.

Для математического описания НСУ можно воспользоваться системой линейных дифференциальных уравнений в форме [3]

бп( Р) Х1 + Qu( Р)Х 2 + ••• + 0.Ш (Р) Хп =

Х

= А( Р)^ + /1(1), т

021 (Р)Х1 + ^22 (Р)Х2 + • • ' + Q2n (Р)Хп =

= Ь2( Р) ^ + /2(0, (1)

т

Qn1 (P)X1 + Qn2 (P)X2 + • • • + Qnn (P)Xn = = Ln (p) + fn (t),

T

где Р = —----обозначение производной по време-

dt

ни; т = t - , 0 < t < (0, t0 — длительность переходного процесса; 0- (р), Ц (р) — линейные

дифференциальные операторы с постоянными коэффициентами, I = 1, 2,..., п, j = 1, 2,..., п.

В работе [3] решение системы дифференциальных уравнений (1) при нулевых начальных условиях сведено к интегральному представлению вида

Xj (t) = 2П Í {W- (P)F- (P)ept +

c-ii

1 p

+ W-l (p)F- (p)epto - — Í e^Uj (z)dz}dp,

Ф(P) І

j = 2,..., n;

c+ii

xi (t) = 2- J W (p)F- (p)ept0 x

(2)

1 p

— Í ez>( z)dz}dp,

P) J

Ф( Р)

где (р), j =1, 2, ..., п — передаточная функция

контура НСУ с разорванной нестационарной связью от входа по воздействию

1 С+1™

/а ^) =----- I вР*Еа (p)dp до выхода звена с коор-

2га 3

С—¡ж

динатой Х- ^(р) — преобразование Лапласа от функции / ^); и- (г) = ф( г) Я- (г);

R, (z) = -

; X1(t) — первая (опреде-

А( г)

ляющая) координата НСУ; Я-(г) — сумма переда-

c-ii

x

точных функций контура управления от входов г Х] (t)

по воздействию £ =------ до выхода звена с коор-

т

динатой х()

Функция р(р) может быть найдена как решение дифференциального уравнения

dр( Р) dp

= —Ф( Р) р( Р),

где

ф( р) = 2 Ц ( р)Жц( р) ;

Ж- (Р) =Агт;

А(р)

б11(Р) 012(Р) •" бш (Р)

А( ) б21(Р) б22(Р) - б2п(Р)

А( р) =

бп1(Р) бп2(Р) - бпп(Р)

А - (р) — алгебраическое дополнение элемента б- (р) определителя А(р); Ф(р) — передаточная

функция от выхода нестационарного звена 1 до

т

выхода по координате г. Передаточная функция Ф(р) называется определяющей передаточной функцией.

В случае, когда различны все корни характеристического уравнения А(р) = 0 и функция Ф(р) является дробно-рациональной,

т а т

Ф(Р) = 2И Р(Р) = П (Р — ) ,

=1 Р — гк

к=1

где гк — корни характеристического уравнения А( р) = 0, (к = 1,., т); а£ — коэффициенты разложения на простейшие дроби функции Ф(р) .

Таким образом, решение задачи Коши в работе [3] сведено к вычислению интегралов (2). Основную трудность при аналитическом вычислении переходных процессов по формулам (2) представляют интегралы типа

1 С+1ж

3 (}) = — I {ер0 X 1( р) х 2т 1

С—¡ж

1 Р

— I е2Тр(2) • X2 (2) • dz} • dp,

( Р )

(3)

щую схему вычисления интегралов типа (3), пригодную для вычисления переходных процессов на всем интервале интегрирования системы уравнений (1).

Метод вычисления интегралов (3) для дробнорациональных и имеющих только простые полюсы функций Х1 (р) и Х2 (г) заключается в замене пути интегрирования вдоль мнимой оси (с = 0) на путь интегрирования по замкнутому контуру в предположении, что т > 0 (рис. 1). Этот контур включает мнимую ось и полуокружность бесконечно большого радиуса в левой полуплоскости (полуокружность на рис. 1 не показана), описанную против часовой стрелки (в положительном

а)

р ( Р )

где X 1(р), X2(г) — дробно-рациональные функции своего аргумента.

Обобщая результаты работ [1-4, 6], изложим об-

+ж

Рис. 1. Контуры интегрирования в комплексной плоскости. Пояснения в тексте

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

1 р

----- I вгтф(г) • X2(г) • йг . Интеграл по такому

Ф(р) ^

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

1 р

функции------- I егтф(г) • X2 (г) • йг . Соотношения

Ф(р) £

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

^1 Cv

сТ(" _т v=m+q

| ...<••• | •••) = у

=тТ1 Р=^

или

с—гж

v=mТq

V=1 Cv

стгж т v=mТq

Г... = _2<••.Г...)т 2 *{•••}, <4)

v=mТ1 Р=V

где 2 Яе .?{•••} — сумма вычетов в полюсах

У=Ш+1 p=zv

функции X 1(р), не совпадающих с особыми точ-

ками

1 р

функции —— [ег—<г) • X2 <г) • & ; [•••

-< Р) 1 1

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

1 Рр

Гегт—г) • X2 <г) •

-<р) __ж

Далее Г•••

dz.

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

Cv

внешнего и внутреннего интегралов, взятых вдоль луча (_тс, гу), и интеграла по окружности бесконечно малого радиуса с центром в точке

г ...=г ер'« ..../

е

■+

Г . При обходе

Cv

контура интегрирования Су эта окружность обходится по часовой стрелке. Поскольку вычисление интеграла по окружности далее сводится к вычислению вычета [7], направление обхода окружности изменим на обратное и соответственно изменим знак перед интегралом по окружности. Внешний и внутренний интегралы в формуле (3) могут быть представлены интегралами по петлям, исходящим из бесконечности в левой полуплоскости и обходящим точку гу по часовой стрелке с изменением аргумента в пределах (0,_2п), если параметр экспоненты во внешнем и во внутреннем интегралах принимает положительные значения.

В нашем случае ¿0 > 0, поэтому для внешнего интеграла можно принять [6]

,р*0 _

1

—2п<Оу •г

,рч

<5)

Для внутреннего интеграла т = t — ^ < 0 , поэтому необходимо сформировать аналитическое продолжение внутреннего интеграла в область отрицательных значений т. Решением этой задачи будет выбор контура интегрирования внутреннего интеграла в виде петли в правой полуплоскости <рис. 1, б). Такая замена может быть представлена схемой

1

е

—2пОу г

т>0

<6)

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

1 т 1 . 1 .

J<t) = — — У —2та, Гept0Х1<Р)-—¿РX ГегТр<г)Х2<г)(к

2т ^=11 _ е 2лса' * -<Р) *

V — ^т

і т —2

2т

е^т X

еРоX 1<р)^—<Р _ ¿V)"15зv <p)dp Т

1

А,т1

V=1 |р_ Zv=COnst\т

-V < Р )

ттq I 1 р I

У Яе5\еРоX 1<Р)—— ГегТ-<г)X2<г^ I, .,=тт р=р^ і -<ру і

г

г

е

с— г ж

г

е

Т

Т

Т

где

Cv+ — петля в положительном направлении

(рис. 1, б), Cv- — петля в отрицательном направлении (рис.1, а); |p — zv| = const + — замкнутый контур интегрирования, охватывающий только одну особую точку z = zv в положительном направлении; ev — порядок полюса функции X2 (z); q — число полюсов функции X1 (p), не совпадающих с особыми точками функции (р(z)X2(z); G = 0 в том случае, если все мнимые части особых точек различны;

(Р) = £D3V к (т)

(Р — Zv ) к

к=0

D3v к (т) =

к!

dp

■ + т

Pv (Р)X2v (Р)

Р = ZV

Для решения поставленной задачи, обобщая результаты работ [1-4], изложим методы вычисления всех интегралов, входящих в структуру <7) интеграла <3).

ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ—КОНСТАНТ

Рассмотрим интеграл

Jlv (tо) = f єР‘0Xі (р)-¡Pdp , J P(P)

(8)

где

Интегралы по окружностям |p - zv | = const учтены в формуле (7) на случай, если X2(z) имеет кратные полюсы (pv = 2, 3, ...) и/или функция X1 ( p) имеет полюсы, совпадающие с особыми точками функции (р(z)X2 (z). Для вычисления составляющих, соответствующих полюсам (p = pk) функции X1( p) , не совпадающим с особыми точками (z = zv) функции (р(z)X2(z), используются теорема Коши о вычетах и разложение в ряды 1 pk

Тейлора функции -------- I eZTq(z)X2(z)dz .

Р(pk) t

В структуру выражения, определяющего значение интеграла (7), входят интегралы-константы и интегралы-функции времени, оставшегося до конечной текущей точки (КТТ) (— Т = to — t).

Методы решения задачи Коши, изложенные в работах [1-4], не получили такого широкого распространения, как методы операционного исчисления [7], поскольку, в частности, к настоящему времени недостаточно проработаны методы вычисления интегралов-функций времени, входящих в формулы (3), что не позволяет установить зависимость параметров переходных процессов от параметров НСУ.

В настоящей работе ставится задача разработки методов получения аналитических зависимостей параметров переходных процессов от параметров НСУ не только в КТТ, но и на всем интервале интегрирования системы уравнений (1). Это позволит оптимизировать параметры всей траектории, а не только конечного положения нестационарного динамического объекта.

ф(р) = (р _ гу) “"•фу (р); X 1(р) =

= (р _ ) _р • X1У (р); ау — нецелое число; ру

— целое число, порядок полюса функции X1(р) в точке р = гу ; X1v (р) — дробно-рациональная функция, не имеющая полюса и нуля в точке

р = гу.

Разлагая в ряд по степеням (р _ гу) функцию (р) = X1У (р)—1—, получим

Pv (Р)

^ і

J1V (to) = Х-ГА^(Zv )Jlv,(to),

и v!-

(9)

где

р a, —pv +ft

J V (t0) = f eP0(P — Zv ) dP =

Cv

1

A(V0)(ZV) = Xiv (zv ) - ( )5

Pv (Zv )

D^( zv ) =

dzV

X у (Zv ) Pv (Zv )

После замены переменной p — zv = u интеграл J\v^ (t0) преобразуется к виду

J 1vM (t0) = eZV‘0 f eUt0 U°a PV ^du .

C0—

У

Выполним замену переменной и =-----------. Если

^0

аргумент переменой и изменялся от 0 до 2п, то аргумент у будет изменяться в пределах (-п , +п), что соответствует контуру интегрирования, подобному С0+ , но при обходе его в обратном направлении. Изменяя направление обхода контура, получим

C

v

к

1

d

av—PvТ^

dp.

<10)

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

Су_

типа (8) можно провести вдоль лучей, исходящих из особых точек и параллельных действительной оси. Тогда интеграл в формуле (10) с точностью до постоянного множителя совпадет с представлением Ханкеля гамма-функции в виде интеграла по контуру [6]

Г<а) = —

1

2 • і • sinпa

Г <—Р)“—1 е — Р dp .

Следовательно, подставляя

г = av _ pv + и + 1 выражение

Грг 1е Рdp

’г—1е - Р dp = е -*•<г —1)г

Г <—Р)г—1 е — Р dp

Cот

Cот

,<е -п• гг — еп • гг) • Цг)

в формулу <8), получаем

Л Иц <0) =

= — е^^ t ~іа>Т^—^—1еп<а>—Р\>Т^т1)і

а —pv +ц

0

е

е гр

dp =

0

Cот

= еzvt0 t -а’Т^—ц—1 (е -п<.^—^тцт1)г — еп^—^тцт1)г

—е

X

хГ<а -Pv+Ц + 1)

2п)Л / пи I , л\ и

и, поскольку е = 1, а е = (_1) , то в соответствии с формулой (9) получаем

—avТPv—1

V ^о) = -2іе^ 0 tо sinп<av - Pv + 1) х

Х X <г ) Г<а - Pv + 1) « <t )

Х X V <Zv ) __ / ч ^оЛ <11)

где

«V <tо)=2 <-1)А

-V < гv)

ц=0

г<av - pv + 1) • t0 •

ц zv );

О(и)(г )

Г (г ) = V \ V /

V) ) .

Полученный ряд 5^ (^) является асимптотическим по ^ рядом, пригодным для вычисления интеграла Jv (/0) при больших t0. Для практического использования асимптотический ряд 51(, (t0) с определенными числовыми коэффициентами может быть просуммирован каким-либо из известных методов суммирования асимптотических рядов [8].

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

Методы вычисления интегралов — функций времени вида

для

Л2v ^ - tо) = Г егт-<г)X2 <z)dz ,

<12)

входящих в структуру (7), известны для больших значений (-т), для которых удается построить эффективные асимптотические разложения, а также для малых значений (т* - т), где т* < 0 — малая величина [1, 2]. Однако выполненные нами расчеты переходных процессов по асимптотическим формулам, приведенным в [1, 2], оказались достаточно точными лишь на небольшом интервале интегрирования системы уравнений (1).

Поэтому для построения аналитических зависимостей, пригодных для оценки переходных процессов на большем интервале, в данной работе был использован другой подход — разложение интегралов (12) в ряды Маклорена по текущему времени t. При таком подходе коэффициентами степенного ряда оказываются интегралы, в подынтегральных выражениях которых остается множитель ехр(-р0), обеспечивающий сходимость интегралов при Яе(р) > 0. При этом условие сходимости выполняется именно на контурах, ветви которых уходят в правую полуплоскость комплексной переменной и по которым осуществляется интегрирование интегралов-функций времени.

Для получения расчетных формул преобразуем формулу (12) для интегралов — функций времени

Л2v ^ — tо) = Г е^-^^2 <РМр = Герт <Р — Zv )т"" —V <Р) • <Р — Zv )тА • X2v <р)Ф =

Cvт Cvт

= е^ Г ерт Р •—V < Р + Zv ) • X 2v < р + Zv )dp.

Со

о

В формуле <13) интегрирование проводится в комплексной плоскости по петле, обходящей начало координат против часовой стрелки и вытянутой вдоль луча, выходящего из начала координат в положительном направлении под таким углом, что путь интегрирования не проходит ни через одну из особых точек функции —V <Р + Zv )• X 2v <Р + Zv ).

Обозначим

= /'

Cо

= Гертр а в —V<р + Zv)• X2v<р + Zv)ф.

<14)

Интеграл (14) является аналитической функцией времени t в промежутке (0 < t < t() < t0). Благодаря соответствующему выбору контура интегрирования подынтегральная функция в (14) является аналитической функцией, и интегралы

^0ук (т) =

О)

Г ~,^еР(‘—1"0)Р —V <Р + ^ ) • X2v <Р + ^ )dР,

к = 1, 2, 3...

где

Л 0v <Т) = Л 0v0 <-t0) + Л 0v1 <—0) • 1: + , Л 0v2<-t0) Л.

т-----------------1 + ...,

2!

<15)

= Г е Ро рк а в —V <р + Zv ) • X2v <р + Zv )dp,

Cо

к = 0, 1, 2, 3,.

Разлагая функцию — v <р + zv) • X^ <р + zv) в ряд Маклорена по положительным степеням р, получим

л,

0ук

< — о) = Г е -р0 рк ^ •Х 2v < p)dp :

<16)

где

Х2v <Р) = 2 »'

1=о 1!

d1

[—V <Zv ) • X2v <Zv )].

dzv;

Следовательно, интеграл (15) может быть вычислен по формуле

¿о*Но) = |е-р°Р^-^ •Х2v<p)dp =

сходятся равномерно во всех точках промежутка 0 < t < t0. Следовательно, интеграл J0v (т) можно где разложить в ряд Маклорена

Cо

= (е -2па,г — 1

)«2vk ^0 ),

<17)

«2і* ^о ) =

= 2 <-1)1

1 =0

^2у1 <Zv ) Г<к + 1' -av - вv + 1)

1-

,к+1— а —А,т1

* п

Для практического использования асимптотический ряд 52Л (t0) может быть просуммирован каким-либо из известных методов суммирования асимптотических рядов [8]. С учетом (14)-(17) интеграл (13) может быть вычислен по формуле

<t - tо) = (е- 1)еZvT—v <Zv ) • X2v <Zv )

Г<-ау - в + 1)

+-а — А/+1

I г>

Bv <^ tо),

<18)

1 . <-а - в + 1) г <. ) t .

1 +-----------1----------о) • — +

1!

; G2vk <*о ) = ; «2v0 <*0 ) = — V <Zv ) • X2v <Zv );

« 2v0<t0)

Для практического использования асимптоти- составляющие первой суммы при нецелых av уже

ческий ряд Bv (^ t0) может быть просуммирован г

могут быть вычислены: — в соответствии с

каким-либо из известных методов [8]. Если ряды ■' J

51(, (t0) и Bv (t,t0) оказываются сходящимся, то, как С''_

показывает приводимый ниже пример, использо- формулой (11), I — в соответствии с формулой

вание операции суммирования может привести к получению точного аналитического решения. Возвращаясь к формуле (7), заметим, что все

<18).

Для целых а вычисление интеграла <13) не-

0

C

V т

обходимо производить обычным образом — с использованием теоремы Коши о вычетах:

Jov(т) = |ертр~а_в •фv(р + Zv)• X2v(р + гv)йр =

С0

= 2п • I • Яе5-ертр~а_в ^ (р + гv)X2v (р + гv)}.

р =0

В результате вычисления получится полином по т степени а + вv _ 1) [6].

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

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

ПРИМЕР: ПЕРЕХОДНЫЙ ПРОЦЕСС, ВЫЗВАННЫЙ НАЧАЛЬНЫМ ОТКЛОНЕНИЕМ ПЕРВОЙ КООРДИНАТЫ ОТ ЗАДАННОЙ ПРЯМОЛИНЕЙНОЙ ТРАЕКТОРИИ

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

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

*1 =П-—*1<0),

*о

<19)

где п — отклонение первой координаты от требуемой прямой линии в переходном процессе; х1(0) — начальное отклонение первой координаты от заданной прямой линии.

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

с передаточной функцией Ф(р) = —- (рис. 3).

р

В этом случае (р(р) = р а1 . Для того чтобы воспользоваться формулами (2) в части, относящейся к координате , необходимо заменить ее с помощью формулы (19) на координату с нулевым начальным значением. При этом очевидно,

1

что на выходе нестационарного звена — появится

т

постоянное воздействие _ — х1(0), а к производ-

t0

ной новой фазовой координаты п необходимо добавить постоянное воздействие — х1(0) . Обозна-

—

чим Жи<р) = ^; Ж,1<р) = —; ^<р) = -

1

= *1<0).

?2< Р) =

*1<0)

Полагая далее Жа1Еа1 =

= Ж11^11 + ^21^21 в формуле (2) для координаты хх^), с учетом (19) получим

0

t

0

t

0

Хї<і) = — Г <1 —еРора1 Гегтг“«1dzdp- — *1<0) =

2т * р *0 р * *0

00

— х1<0)1—— Г еРор— 2Гегтг «1dzdp - —х^

0 с —і ж ж 0

— Х1<0) — •^_______<1 - «О Г еР*0 р — —2ЛР х Г еРт р-«1

ГерІ0р«1 2dp х Гертр «1dp +

*0 2т 1 - е ~2п—

0 Cl— Clт

*1 <0) •— 1 -«1 _ і р<0 _2 р « < рт)^р *1 <0) —

еРор • р • 53<р,т)dp ------------------------------------, <20)

0 2ті 0

0 |р| =сош1 0

жк

где 53 <Р,т) = 2°31к <т)--7 = Апо<т)~-+ В311 <т)^~~—+ -

- —1 + к + 1 1 - —1 2 - —1

+

Рис. 3. Схема управления движением материальной точки по заданной прямолинейной траектории

Вычислим коэффициент D3lk (т) при к = 0:

"Г d ^0

D310(t) =

d

dp

■ + Т

Vv (p)X2v (p)

= 1.

p=0v

Xj(0) • Т 1 — a1 2П

= 0,

0 |p| =const

x1 (0)•Т

ept0p 2 • p • S3(p,T)dp —

t,

0

I ept0 p“1 —2dp =

Q-

= —2 • i • sin[n(a1 — 1)]Г(“ — 1) • (t0) “v+1;

| ezT z -“1dz = (e -2mX1i — 1)

r(2 — a) (1 — «)

(—т)“1—1. (21)

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

1). г( _а + 1)

! [^0 ]_“'+'

I ezT z -a1dz = (

—2n(X1i

-^(t)

= (e ~2па1 — 1) Г(2 a) • (—Т)“1—11,

(1 — a1)

где

A(t) =

1 + (—“ + 1) ± +

1! 10

, (—“v + 1) • (—“v + 2) 12 ,

+---------------------------------Г" + '

2!

(

t

1—

t0

0

В этом конкретном случае ряд В1 ^) оказался суммируемым и сходящимся к степенной функции.

Объединим полученные результаты

x1 (t) = —

x1 (0) т 1 (1 — а1)

х

Вычисление D311(t) и последующих коэффициентов нецелесообразно, т. к. соответствующие им степени p в разложении S3(p,T) не меньше единицы и соответствующие интегралы по замкнутому контуру равны нулю.

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

11

— • Iep0 p — • p^S 3(p,T)dp = -;

2m , ,J 1 — a

p =const 1

.0 2m 1 — e ~2л°1

х |ept0pa —2dpх |epTp_“Jdp =

Q+

= x1 (0)

t0

0

sin n« — 1)

(1 — a1)---------------------------—- х

Г(“1 — 1)Г(2 — a)

(1 — a1)

Используя

г(“1 — 1)Г(2—“) = -

свойство n

гамма-функции [9], получаем ко-

sinn(a1 — 1)

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

X1id(t) = X1 (0)

( — т\“

t0

0

(22)

Если а1 равно целому числу, то можно воспользоваться теоремой Коши о вычетах, согласно которой

I e z 1dz = 2ni • Re 5

J z =0

C0+

При “ = 2 получаем

1

«1

= 2П •

.“1 —1

TOO

—W I ept0 — a^p =

1 e C

—1

t0

t

х

e

C

Cот

г=0

= —т

Тогда

—2п—1І

-х

хш<* ) = — ^ .

t о 2т 1 — е

х|ері0р——2<1 -—1)dpх Гертр~«Чр = х1( 0 ) • т2

Cl*

*о2

<23)

у <0 = —1

А0_

t — 10

<24)

у

= —

dt

t — *0

Интегрируя <25), для t >*0 получаем 1п у - 1п у0 = — 1п<* - *0) + 1п C ,

<25)

<26)

C = <-*о) — 1 и у = Уо

Г-тЛ«1

<27)

*1

Корректность формул (22) и (23) можно подтвердить непосредственным решением дифференциального уравнения, соответствующего выбранному виду определяющей передаточной функции

Ф( р) = — (рис. 3):

р

при начальном условии у(0) = у0.

Уравнение (24) является дифференциальным уравнением с разделяющимися переменными, т. к. может быть приведено к виду

Рис. 4. Графики переходных процессов

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

печивалось заданное значение хш<*о) в КТТ:

1к

ускорения

1к

или

*1<0)•

—1 <—1 — 1)

—1

<—т)

—1 —2

<28)

где С — постоянная интегрирования.

Для получения аналитического продолжения решения <26) дифференциального уравнения <20) в область интересующих нас значений t формально определим постоянную интегрирования С, исходя из начального условия у<0) = у0. Тогда получим

Неравенство (28) будет выполняться при любых а1> 2 и не будет выполняться при любых а1 < 2. При а1= 2 условие (28) будет выполняться, если переходный процесс начался заблаговременно, т. е. должно выполняться ограничение

*0

2 • х1 <0)

Х1к <0

| | V У

Полученное методом разделения переменных решение (27) дифференциального уравнения (24) с точностью до обозначений соответствует решениям (22) и (23), полученным другим методом — с использованием преобразования Лапласа, как для целых, так и нецелых значений параметра а1.

Графики переходных процессов при различных значениях а1 представлены на рис. 4. Располагая аналитической зависимостью первой координаты от параметров х1 (0), а1 и t0, можно выбрать зна-

Траектория снижения (22) будет соответствовать заданной прямолинейной траектории (заданной траектории равномерного изменения дистанции до КТТ) при а1 = 1 и при соответствующем выборе параметров х1 (0) = х1 (0) зад и ^ ^зад .

ЗАКЛЮЧЕНИЕ

В работе рассмотрены методы получения ана-

и

г

0

т=0

t

0

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

СПИСОК ЛИТЕРАТУРЫ

1. Ludwig G. und Rollnik H. Theorie von Regelsystemeb mit zeitlich variabler Regelstärke // Zeitschrift für angewandte Mathematik und Mechanik, Heft 11/12, 1957. Berlin: Akademie Verlag, 1957. S. 457-470.

2. Ludwig G. und Rollnik H. Erzwungene Schwingungen und Fehler bei regelsystemen mit zeitlich variabler Regelstärke // Zeitschrift für angewandte Mathematik und Mechanik, Heft 1/2, 1958. Berlin: Akademie Verlag, 1958. S. 16-20.

3. Барабанов А.Т. Теория линейных нестационарных систем с особой точкой. Устойчивость систем // Автоматика и телемеханика. 1969. № 6. С. 31-39.

4. Барабанов А.Т. Методы исследования систем с переменными коэффициентами // Методы ис-

следования нелинейных систем автоматического управления / Под ред. Р.А. Нелепина. М.: Наука, 1975. С. 132-254.

5. Федосов Е.А., Инсаров В.В., Селивохин О.С. Системы управления конечным положением в условиях противодействия среды. М.: Наука, 1989. 272 с.

6. Уиттекер Е.Т. и Ватсон Г.Н. Курс современного анализа, часть вторая: трансцендентные функции. (Пер. с англ.) Л.-М.: ГТТИ, 1934. 468 с.

7. Математическая энциклопедия, т. 4. Ок-Сло / Гл. ред. И.М. Виноградов. М.: Советская Энциклопедия, 1984. 608 с.

8. Математическая энциклопедия, т. 5. Слу-Я / Гл. ред. И.М. Виноградов. М.: Советская Энциклопедия, 1984. 624 с.

Санкт-Петербургский Государственный университет аэрокосмического приборостроения

(С. В. Богословский)

Санкт-Петербургская Академия межотраслевых наук (В. С. Богословский)

Материал поступил в редакцию 5.06.2002.

DYNAMICS OF NON-STATIONARY SYSTEMS WITH UNIFORMLY TIME-VARYING FACTORS

S. V. Bogoslovsky, V. S. Bogoslovsky*

St. Petersburg State University of Aerospace Instrumentation *St. Petersburg Academy of Interdisciplinary Sciences

The linear differential equations with floating factors are widely enough applied to simulation of the dynamics of non-stationary automatic control systems. The methods of obtaining analytical relationships between parameters of transients and those of a control system based on application of the Laplace transform to linear systems of ordinary differential equations with linearly time-varying factors are considered.