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

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

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

456

ТОЧНОСТЬ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ [ГЛ. 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]

ПРИБЛИЖЕННОЕ ОПРЕДЕЛЕНИЕ

ПЕРЕДАТОЧНОЙ

ФУНКЦИИ

 

457

исключается при вычислении x ( t )

из выражения

(10.14).

Если

же,

кроме того, Ф(/ш, t) изменяется по форме, т. е. если

в системе

происходит изменение характера переменности параметров

(это

осо­

бенно

важно для систем, параметры которых меняются в зависимости

от входного сигнала); то Ф(/ш, t) является общим решением (10.17). Произвольные постоянные при этом определяются из условий:

^ Ф С А » . О

fift

 

/е =

0, 1,. . ., /г -1 , (10.18)

= ~ Ф ( У ш , t)

 

t=t,+

dtv

t=t:

 

 

где

tQ— момент

изменения

формы

Ф(усо, t).

Эти зюловия

вытекают

из

того практически всегда

выполненного условия, что

параметры

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

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

положении

о малой скорости изменения коэффициентов

уравне­

ния (10.17).

При этом они в большинстве практических

случаев

дают достаточные приближения к точным решениям.

 

Изложим вначале приближенный метод определения частного

решения (10.17).

В качестве первого приближения к точной

переда­

точной функции

примем выражение

 

 

Фо О . О

м о ,

t)

(10.19)

 

D

( М

t)

 

 

 

которое может быть

получено из

 

(10.17) в предположении, что

производные от Ф (y'u>,

t) в рассматриваемом промежутке изменения t

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

чески

отличается

от нуля. Кроме того, для больших значений

частоты входного

сигнала

точная передаточная функция асимптоти­

чески

стремится к

(10.19).

Это объясняется тем, что по отношению

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

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

рекуррентную

формулу, получаемую очевидным образом, из (10.17):

Ф iC/co. t) :

1

dn$ i - 1 ( 7 o > ,t)

 

4Ф/-1о. о

D О , t) L'«»(0

dt n

- b ai

cLt

 

 

 

 

( 10. 20)

l) Отметим, что (10.19) совпадает с выражением передаточной функции, получаемым по методу «замораживания» коэффициентов.

458 ТОЧНОСТЬ ЛИНЕЙНЫХ с и с т е м с п е р е м е н н ы м и п а р а м е т р а м и [ г л . X

Величины, получаемые из (10.20), позволяют

построить ряд для

функции

вида

 

 

 

СО

 

 

 

Ф ( у Ч 0 = 2

® / О . 0.

(10.21)

 

г-о

 

 

где ФДу’со. У) определяются по формулам (10.19) и (10.20).

Этот

метод определения передаточной функции

линейной системы

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

Необходимо отметить, что в большинстве случаев получаемое этим методом решение (10.21) является степенным рядом относительно малого параметра и исследование его сходимости не представляет труда.

П р и м е р . Пусть, например, дифференциальное уравнение системы имеет вид

T(t) -§■ + * = К (О/(О-

Такую систему можно рассматривать как последовательное соедине­ ние ЯС-фильтра с переменным R и усилителя с переменным коэф­ фициентом усиления К (У).

В этом случае

 

<*о0=

1,

 

 

 

 

(О = Т (У)

 

 

и, следовательно, (10.17) имеет вид

 

 

 

т со аФ Ud t:—

+

(0 >

+ 1 ] Ф и

О = к (У).

 

Отсюда (10.19)

и (10.20)

записываются как

 

 

Фо О . 0 =

 

K(t)

 

 

 

П 0У - + 1

 

 

 

 

 

 

Ф|(У'«о. У) =

 

1

^ф/-! (/». 01.

 

 

7" (У)> +

1

dt

 

так, например,

 

 

 

 

 

\К' {t)[T (О> + i ] - r (/) / < ( У ) у о л

 

®i о . 0 =

7" (0

 

_

7 (0 у'“ 4- 1 L------------- [ Г ( 0 > +

И3-------------J И т -

д-

Отсюда второе приближение будет равно

 

 

Ф (Уч), У) =

Фо(Уш. У) ~Ь ФДУ10* О —

 

 

1

 

I ™

ТЧОКЧО

, г (0 Г (0 К (0 У М

7(0У“ + М

w

7-(Оусо+1

[Г(0У« + 1]5

Г

Сравнивая Ф0(у'ш, У) с нетрудно видеть, что первое приближение будет мало отличаться от нулевого, когда К (У) и Т (У) медленно меняются, т. е. их первые производные достаточно малы.

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.

'4]

ПРИБЛИЖЕННОЕ

ОПРЕДЕЛЕНИЕ ИМПУЛЬСНОЙ

ПЕРЕХОДНОЙ ФУНКЦИИ 463

 

Выбирая в качестве приближенных значений коэффициентов урав­

нения (10.28)

их

величины

при

t = x,

т.

е. в момент приложения

воздействия о (/ — т), перепишем

(10.28)

в

виде

 

 

 

D(p,

т)/еи(/, т) =

[D(p, т) — D (р, t))kQ(t, х),

(10.32)

где

p =

d/dt,

a

D(p, т) имеет постоянные

относительно t

коэффи­

циенты,

равные 0/(х).

 

 

 

 

 

 

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

решение в виде

ряда

 

 

 

 

 

 

 

k0(t,

 

х) = k°0(t, х) + k'0(t, x)-\ -kUt . x ) + . . . .

(10.33)

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

D(p,

т) k0(/, х) = 0,

D(p,

(10.34)

x)k Q(t,l х) [D (р, x) — D( p, t )]k о-1 (/, т).

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

времени

/, вне которого kQ(t, т)

мала по сравнению с ее максималь­

ной величиной, этот метод дает

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

решению

при г = 2, 3 ‘).

 

В заключение кратко рассмотрим экспериментальный метод опре­

деления

импульсной переходной

функции линейной системы с пере­

менными параметрами с применением электронных моделирующих

устройств.

задача

определения k(t, т) или

 

Как это было отмечено,

k0(t, х)

как функции t сводится, по

существу,

к решению однородного диф­

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

(10.26)

при

начальных условиях

(10.27)

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

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

’) Z a d e h L. A., The determination of the impulsive response of variable networks, Journal of Applied Physics, t . 21, July 1950,

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).

4] ПРИБЛИЖЕННОЕ ОПРЕДЕЛЕНИЕ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ФУНКЦИИ 465

тельно, эта область совпадает с областью определения физически осуществимого решения исходного уравнения (10.35).

Поэтому задача экспериментального определения k0(t, х) как

функции

х сводится

к решению на

модели

уравнения

(10.39), где

 

 

П

 

 

 

 

D^ p,

В)х(в) = 2 - J j r [fl‘(* — 0)*(в)Ь

(Ю.40)

 

 

/=о

 

 

 

Легко

проверить,

что в случае

линейных

систем с

постоянными

параметрами (10.40) совпадает с (10.35). Поэтому импульсная пере­ ходная функция как функция ее второй переменной может быть опре­ делена непосредственно моделированием исходного уравнения (10.35).

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

Необходимо, однако, при этом помнить, что в результате реше­

ния мы получаем не k0(t,

х), a k0(t,

t — т), которая,

очевидно, при

каждом фиксированном

t,

отличается

от

k0{t,

х) ориентацией отно­

сительно

оси

т.

 

 

 

 

 

 

П р и м е р.

Пусть линейная система с

переменными

параметрами,

имеющая

импульсную

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

k (t, х),

описывается

уравнениями

 

 

 

 

 

 

 

 

 

% ■ + ( < - Н > ) ^ - + ( « + 1 ) *1 = 3/(0,

 

 

■ ^ + ( т + » )

 

+ И * , = - / У ) .

 

 

 

 

X — XTj-f- х 2,

 

 

 

где b — const, f(t) — воздействие на

входе и

x(t) — процесс на вы­

ходе системы.

Согласно (10.28) для определения импульсной переходной функ­

ции имеем

уравнения

 

 

 

 

 

 

 

 

+

(*'+ Ь)

+

(fit +

1) k, (t,

т) =

0,

- k W —

И- ( у +

ь)

1

(М +

1) k2 (t,

х) =

0,

 

 

k ( t , x) =

ft,(#> т)+

*„(#. x)

 

 

!) См.

книгу,

цитированную на стр. 464.

 

 

 

 

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

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