книги из ГПНТБ / Солодовников В.В. Статистическая динамика линейных систем автоматического управления
.pdf456 |
ТОЧНОСТЬ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ [ГЛ. X |
эта |
задача может решаться как моделированием уравнения исходной |
системы, так и уравнения (10.17). |
|
|
Второй метод состоит в применении аналитических приближенных |
методов решения уравнения (10.7) или (10.17). Преимущество этого метода состоит в том, что при этом получается аналитическое выра жение для Ф(у'ш, t) или k(t, т), что упрощает вычисление выраже ний (10.6) и (10.14).
Рассмотрим кратко эти методы.
3. Приближенные методы определения передаточной функции линейной системы с переменными параметрами
Переходя к изложению приближенных методов определения пере даточных функций линейных систем с переменными параметрами, необходимо указать на возможность непосредственного применения
преобразования Лапласа к |
решению дифференциального |
уравне |
ния (10.1), определяющего |
систему. Как можно показать1), |
решение |
этого уравнения с помощью сравнительно несложных вычислений может быть получено в форме ряда, члены которого находятся по методу последовательных приближений. При этом вследствие того, что этот ряд при довольно общих условиях может быть представлен в виде убывающих геометрических прогрессий, обычно возможно записать решение в конечной форме. Этот метод позволяет опреде лить передаточную функцию системы, не зависящую от времени, для заданных воздействий, т. е. для различных входных сигналов будут получаться различные передаточные функции. Причина этого заключается в том, что, пользуясь этим методом, мы как бы ищем систему с постоянными параметрами, эквивалентную в смысле реакции на данный входной сигнал исследуемой переменной системе.
Поэтому, исходя из случайного характера воздействий, здесь более целесообразно рассматривать приближенные методы решения уравнения (10.17), которое приводит к передаточной функции, опре деляемой только параметрами системы и не зависящей от входных воздействий.
Итак, рассмотрим некоторые приближенные аналитические методы решения уравнения (10.17). Прежде всего решим вопрос о начальных условиях, которым должна удовлетворять искомая передаточная функция для возможности однозначного ее определения.
Можно показать2), что если функция Ф(/ш, t) на рассматривае мом промежутке не имеет скачков, то она является частным реше
нием уравнения (10.17). |
Кажущаяся произвольность в выборе решения |
*) С о л о д о в н и к о в |
В. В., Об одном применении операторного исчи |
сления к динамическим |
системам с переменными параметрами, Изв. |
АН СССР, ОТН, № 12, 1945.
2) См. работу, цитированную на стр. 453 или более подробную работу — на стр. 454.
3] |
ПРИБЛИЖЕННОЕ ОПРЕДЕЛЕНИЕ ПЕРЕДАТОЧНОЙ |
ФУНКЦИИ |
4 5 9 |
На |
рис. 10.3 приведены1) кривые реакций |
на ступенчатое |
воз |
действие 1(0. определенные по (10.14) для приближенного и точного решений. Отклонение от точного решения не превышает при этом 4%.
Рассмотрим теперь метод приближенного общего |
решения |
урав |
|||||||
нения |
(10.17). |
|
|
|
|
|
|
|
|
Обозначим |
через |
а1 постоянные, равные |
средним |
значениям |
at (t) |
||||
на рассматриваемом |
промежутке времени или, в общем случае, зна |
||||||||
чениям |
at (t) |
в какой-либо |
фиксированный |
момент |
t, и |
перепи |
|||
шем (10.17) в виде |
|
|
|
|
|
|
|
||
d"Ф О , 0 , -* |
dn *Ф ОЧ t) |
Ч- ao® (7“ . t) — |
|
|
|||||
|
|
An-1 |
■n-i |
• • • |
|
|
|||
dt n |
1 ' ' ' |
1 |
|
|
|
|
|||
|
dtl |
= M(ju>, |
|
|
01. |
(10.22) |
|||
|
|
|
|
|
|
||||
где обозначено
-dn(L>( / ы ()
D' [Ф U ш, 01 = |
К — «л (01-----jfn— |
■ |
• • • + K — «о (01 Ф О . t). |
|||||
|
|
|
|
|
|
|
|
(10.23) |
|
Применяя, как и ранее, метод последовательных приближений, |
|||||||
будем |
искать |
решение в |
виде |
ряда |
(10.21). |
|
||
|
Первое приближение опре |
|
|
|
||||
делим как решение неоднород |
|
|
|
|||||
ного дифференциального урав |
|
|
|
|||||
нения с постоянными коэффи |
|
|
|
|||||
циентами: |
|
|
|
|
|
|
||
- |
dn<I>0 0 4 t). |
| |
|
|
|
|
|
|
“л |
dtn |
' |
|
|
|
|
|
|
• .. + |
а0Ф0(уш, t) = |
M(jw, |
t), |
|
|
|
||
так |
как |
a, = const. |
(10.24) |
|
|
|
||
|
|
|
|
|
||||
|
Последующие приближения , т. е. члены ряда (10.21), опреде- |
|||||||
ляются |
как решения (10.24) с измененными правыми частями: |
|
||||||
|
|
4"ф,(уЧ 0 + |
___ + |
а0Ф(.(уш, 0 |
= О ЧФ ,_1(7а). 01- |
(10.25) |
||
Если на рассматриваемом интервале изменения t происходит изме нение формы Ф (усо, t), т. е. характера переменности коэффициентов уравнения (10.1), то необходимо искать общие решения уравне ний (10.24), (10.25). Получающиеся при этом п произвольных постоянных определяются из условий (10.18). Для этого, очевидно, необходимо знать величину передаточной функции и ее производных до (п — 1)-го порядка в момент времени, непосредственно предше ствующий моменту изменения формы Ф(уш, t).
1) См. работу, цитированную на стр. 454.
460 ТОЧНОСТЬ ЛИНЕЙНЫХ |
СИСТЕМ |
С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ [ГЛ. X |
Как показывает ряд |
задач, |
этот метод приводит к быстро схо |
дящемуся ряду (10.21) в случае, когда параметры системы являются периодическими функциями. Однако и в большинстве других случаев обычно практически достаточно определения двух-трех первых членов
этого ряда. |
|
Рассмотрим |
систему |
первого порядка, |
описываемую |
||||||||
П р и м е р . |
|||||||||||||
дифференциальным |
уравнением |
|
|
|
|
|
|
||||||
|
-^- + |
0 |
+ |
р cos 1о00 х = |
(1 +pcosco00 /( 0 - |
|
|
||||||
Дифференциальное |
уравнение, соответствующее |
передаточной функ |
|||||||||||
ции системы Ф (усо, /), имеет вид |
|
|
|
|
|
|
|||||||
аФ <j “’ |
+ |
(У'м + |
1 -Ь р cos u>00 Ф (уш. t) = |
1 + р cos ш01. |
|||||||||
Рассматривая |
интервал времени |
Т |
|
2тс/о)0, |
выбираем |
из очевид |
|||||||
ных соображений |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Оо = Д ° + |
Ь |
|
|
|
||
Тогда (10.23) |
принимает |
вид |
|
|
|
|
|
|
|||||
- Ф^ ' - |
|
(у'ш +1)Ф (уш, 0 = 0 |
+ |
р cos (о00 — рФ(уш. Осоэш^. |
|||||||||
Согласно изложенному получаем уравнения для определения при |
|||||||||||||
ближенных значений Ф (у'со, t): |
|
|
|
|
|
|
|||||||
^Ф0(уЧ_0 _|_(1 —|—уш) Ф0 (у'о), 0 = |
1 Р cos ш0£, |
|
|
||||||||||
|
(уь>. 0 |
|
^ |
j |
у*со) Ф, (уш, |
/) = |
— рФ0(у'ш, Ocosu)0C |
||||||
а'Ф; |
|
+ |
(1 + |
» |
Ф/ О . |
О = |
— рФг- 1 0 . О cos Ш<У- |
||||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
Ф0 (уш, |
0 |
= |
С е (у<0+1) 1 |
------ - |
|
|
-pu)0 sin ш0/ ■ |
|
|||||
|
|
|
|
|
|
|
y<0+ 1 |
|
|
(1 + 7 “) |
+ “о |
||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
—|—pcos u)0£ |
|
1 + J u |
|
|
|
|
|
|
|
|
|
|
|
|
УЧ>2 + »o |
||
|
|
|
|
|
|
|
|
|
|
|
(1 + |
||
ф^уш, |
|
_n„-iX+Ja)l |
1 |
|
|
p COS 10Qt |
14-y'co |
|
|||||
t) = |
Ce |
|
|
1 + j O) |
(1 -py<o)3 -f- ш3 |
||||||||
|
|
|
|
|
|
|
—|—рш0 sin u>Q( |
|
|
и т. д. |
|||
|
|
|
|
|
|
|
|
|
|
(1 + ;ш)2 + |
Ш3 |
||
4] ПРИБЛИЖЕННОЕ ОПРЕДЕЛЕНИЕ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ФУНКЦИИ 461
Здесь С определяется из начальных условий, получаемых из (10.18)
Ф(уш, ;0~) = Ф ( > , to)
по заданной функции Ф(у-ш, to) в момент, непосредственно предше ствующий моменту tQ изменения формы Ф (у‘ш, t).
Отметим, что при р 1 Ф (уи>, /) = Ф0 -(-Ф, дает хорошее при ближение к точному решению.
4. Приближенные методы определения импульсной переходной функции
Здесь удобно дать эквивалентное, но более строгое математи чески, определение импульсной переходной функции k(t. т), соот ветствующей системе (10.1).
Как известно2), неоднородное линейное дифференциальное урав нение, содержащее в правой части дельта-функцию и ее производ ные, при нулевых начальных условиях эквивалентно однородному дифференциальному уравнению с ненулевыми начальными условиями.
Поэтому |
из уравнения (10.7) |
|
прямыми вычислениями можно пока |
|||||||
зать3), что k(t, |
т) |
является |
решением однородного |
дифференциаль |
||||||
ного |
уравнения |
|
D(p, |
t)k(t, т) = 0 |
(10.26) |
|||||
|
|
|
|
|
|
|||||
при |
начальных |
условиях: |
|
|
|
|
||||
d°k (t, т) |
= |
0, |
а = |
0, |
|
1......... п — т — 2, |
|
|||
|
dttt |
/-1 |
|
|
||||||
|
|
|
|
|
|
|
|
|||
d' -v-'kit, т) |
|
|
(*) |
|
|
|
|
|||
|
dt |
|
|
|
|
|
|
(10.27) |
||
|
|
п—а—2 |
|
|
Л —р.—Г—1 |
|||||
|
|
|
|
|
||||||
- |
|
Е |
drkdt; |
Т) t . T |
2 |
|
( - 1 ) ^ +0а :+, „ +1(т) |
|||
|
г ~ п — т — \ |
|
|
|
а» о |
|
|
|
||
|
|
|
|
|
|
{х = 0, |
1........... т . |
|
||
где |
|
|
|
|
|
|
|
|
м! |
|
|
|
|
|
|
|
а |
= |
|
||
|
|
|
|
|
|
(у— Н-)1- Н-! |
|
|||
|
|
|
|
|
|
|
|
|
|
|
*) См. работу, цитированную на стр. 453. |
parameter systems, |
|||||||||
2) |
Z a d е h |
L. |
A., Initial conditions in linear varying |
|||||||
Journal |
of |
Applied |
Physics, t . 22, № 6, 1951. |
|
||||||
3) Б а т к о в |
A. M„ Решение |
|
одного класса нестационарных задач ста |
|||||||
тистической динамики САУ с применением моделей. Канд. диссертация, Мос ковский инженерно-физический институт, 1958.
462 |
ТОЧНОСТЬ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ |
ПАРАМЕТРАМИ |
[ГЛ. X |
||||||||||||
|
В частном случае, когда |
М(р, / ) = 1 , |
импульсная переходная |
||||||||||||
функция |
k0(t, т), |
соответствующая (10.1), определяется как решение |
|||||||||||||
при |
|
|
|
|
|
D(p, |
()kQ(t, х) — 0 |
|
|
|
(10.28) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d ‘ k „ ( t , |
х) |
|
С |
0, |
1 = |
0, |
1......... 11 — 2, |
(10.29) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
d t ‘ |
|
t=Тт |
1 —т т - • 1 = п — 1. |
|
|||||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
I |
a,i (t) |
|
|
|
|
|
|
|
|
Из выражения (10.6) можно установить достаточно простую связь |
||||||||||||||
k(t, |
т) |
с соответствующей |
системе |
(10.1) |
при |
М( р, |
t ) = 1 функ |
||||||||
цией kQ{t, т). |
|
|
|
|
|
|
|
|
|
|
|
||||
|
Действительно, |
рассматривая правую |
часть |
в (10.7) как воздей |
|||||||||||
ствие, |
можно записать на |
основании (10.6): |
|
|
|
|
|||||||||
|
|
|
|
|
I |
|
|
|
|
|
|
|
|
|
|
|
k{t, |
т ) = |
j ' |
k0(t, X) М (р, Х)8(Х — т) rfX, |
где |
р = |
|
||||||||
далее, |
так |
как |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
/ |
/ ( и ) — |
(“ик |
du = |
(— !)*/<*> (т). |
|
(10.30) |
|||||
то при |
/ > |
т получаем: |
|
|
|
|
|
|
|
|
|
||||
*(*. т) = |
£0(г)£0(/, |
х) — |
[ft, (х) *0 (/. |
т ) ] + . . . |
|
|
|||||||||
|
••• |
|
|
|
дт |
|
|
X)] = ЛГ(р. |
x)k0(t, |
X). |
(10.31) |
||||
|
|
|
|
|
|
|
|||||||||
где p = |
d!dz, а ЛГ(/?, т) |
является |
оператором, |
сопряженным |
опера |
||||||||||
тору Л4 (р, т) *). |
|
|
|
|
|
|
|
|
|
|
|
||||
|
Полученные выражения |
показывают, |
что |
в общем |
случае |
задача |
|||||||||
определения импульсной переходной функции сводится к решению линейного однородного дифференциального уравнения с переменными коэффициентами при ненулевых начальных условиях. При этом оче видно, что обычно более просто искать k0(t, т), а затем определять k(t, т) дифференцированием по формуле (10.31). Эта задача экви валентна, по существу, задаче определения передаточной функции и может быть приближенно решена методами, аналогичными изло женным в предыдущем параграфе.
Остановимся кратко на одном из них.
*) О сопряженных операторах см., например, С т е п а н о в В. В., Курс обыкновенных дифференциальных уравнений, ГИТТЛ, 1953.
464 ТОЧНОСТЬ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ [ГЛ. X
на основе метода сопряженной системы]), который состоит в сле дующем.
Известно*2), что импульсная переходная функция заданной системы, т. е. функция, удовлетворяющая уравнению
D(p, f)*o(Л т) = о(/ — т), t > т, р = - ^ , |
(10.35) |
рассматриваемая как функция ее второй переменной х при фиксиро ванных значениях первой, является решением дифференциального уравнения, сопряженного данному, т. е.
D*(p, x)kQ(t, |
x) = o(t — x), р = -§^. |
(Ю.36) |
где |
п |
|
|
|
|
D(p, |
= |
(10.37) |
И |
п |
|
|
|
|
D*(P, f ) x ( 0 = |
У (— 1)' - ^ - [ ^ ( 0 ^ ( 0 1 . |
(10.38) |
мdt
Однако при попытке решить (10.36) на модели возникает прин ципиальная трудность, состоящая в том, что область определения
физически осуществимого решения сопряженного |
уравнения 3) (10.36) |
||
x ^ - t , т. |
е. область, в которой значения первого |
аргумента t меньше |
|
значений |
второго аргумента т, не совпадает с областью |
физически |
|
осуществимого решения исходного уравнения (10.35) t ^ - x . |
Поэтому, |
||
чтобы определить вид импульсной переходной функции при решении сопряженного уравнения (в области физической осуществимости реше
ния исходного |
уравнения), заменим |
переменную х в (10.36) |
на t — 0. |
||||||||||||
При этом |
получим уравнение |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
Я. (А е)*о(*. |
t — |
G) = |
8(0), |
p |
= |
- f f . |
|
(10.39)' |
|||
|
Физически |
осуществимое |
решение |
этого |
уравнения |
определено |
|||||||||
в |
области |
0 ^ -0 (дельта-функция |
прикладывается |
ко входу системы |
|||||||||||
в |
момент |
0 = |
0), т. е. в области, |
где значения |
первого |
аргумента |
|||||||||
импульсной переходной функции t больше |
второго, |
t — 0. |
Следова |
||||||||||||
|
!) |
L a n i n g |
J. Н., В a 11 i n |
R. H., |
Random |
processes |
in |
automatic |
con |
||||||
trol, N. |
Y., |
1956. |
Имеется перевод: |
Л э н и н г |
Дж„ Б е т т и н |
Р. Г., |
Слу |
||||||||
чайные процессы в задачах автоматического управления, ИЛ, Москва, 1958.
2) М i 11 е г К. S„ |
Properties of Impulsive |
Responses |
and |
Green’s funct |
|
ions, IRE Transactions, Circuit theory, |
March 1955. |
гл. I. |
Отметим,- что |
||
3) О физической |
осуществимости |
решения |
см. § 7.9 |
||
при решении (10.36) аналитическими методами может быть использован
метод |
последовательных приближений, изложенный ранее. |
При этом, оче |
|
видно, |
решение определяется во всей области изменения t |
и х , |
что исклю |
чает необходимость дополнительных преобразований уравнения |
(10.36). |
||
