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

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

.pdf
Скачиваний:
9
Добавлен:
24.10.2023
Размер:
9.53 Mб
Скачать

Для определения

(.to и

jib

используя

формулу

(80),

запишем

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 = , H 0 j V « e d e

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

- Сх = - И 1 j' е~ав

dQ - цо j &е-ав

dB = - - ^ - - - ^ Но-

 

 

 

 

о

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но = а;

Ці = c i a

1

 

 

 

 

 

 

Определив

Но и Hi, подставляем

импульсную

переходную

функцию в ог­

раничивающие условия

(79)

и

находим

еще два

уравнения

 

 

 

 

А Т2

+ £ 0

+

D 0 ;

Сга-1

 

А Т"

 

 

АлТ3

 

 

а = А0Т+-^—

 

 

= - ^ —

+ — і — + Ц,Г /

Решая систему четырех уравнении, получим

 

 

 

 

 

 

 

 

 

 

4 а 2 ( а 2 7 2

+

ЗаТ + 3)

 

 

 

х а — 1) 6а 2

 

 

 

( а 2 Г 2

+

6аГ +

12) (аГ +

2)

 

( а 2 Г 2 + 6аГ +

12)

 

 

Лі

 

 

 

ба»

 

 

 

 

(Сісс — 1) 12g2

 

 

 

 

~

а 2

Г 2

+ 6 а Г + 1 2

+

Т (а-Т2 -\- 6аТ +

12)

'

 

 

 

 

 

 

 

2 (2д2Г2

+

9схГ +

12) а

 

( d « — 1) 6 (аГ +

2)

_

:

( а 2 Г 2

+

6аТ +

12) (аГ +

2)

~

Т {а-Т-

+

баГ + 12)

'

 

 

 

 

2 Г (аТ +

3)

 

 

, ( C T K —

1) 6 («Г +

12)

 

~~ (а 2

Г 2

+

6аГ +

12) (аТ +

2)

~ г

Г (а 2 Г 2

+

6аГ +

12)

Подставляя

А0,

Ah

Е0

и D0

в функцию

k*(t),

находим

оптимальную им­

пульсную переходную функцию эквивалентного корректирующего устройства.

Теперь предположим, что g{()=go;

Я т ( т ) = 0 ;

Rn(r) =№б(т); Н(р)=Нв(р)

=

= 1, а импульсная переходная

функция

 

 

 

 

 

ч ч

А + Р

 

Р

ч

 

принадлежит

неминимально-фазовой

заданной

части,

так как соответствую-

 

 

 

1 — fis

 

 

 

щая ей передаточная функция W0(s)

=

.

 

 

 

В данном

случае

 

1 +

as

 

 

 

 

 

 

 

 

 

ІУ2 (1+В2 со2 )

и поэтому

k* (т) = Ао + 73ie*"T + £ , e ^ T , 0 < т < Т,

где

Подставляя k*(x) в интегральное уравнение (87), получим два соот­ ношения для определения А В , ВІ и В2:

А0а + —

В1е%'г

 

 

В.,ех*т

 

= 0;

 

 

 

 

 

 

я.

 

 

 

д-

= 0,

 

+ — z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— Xi

 

—• — Хо

 

 

 

 

 

 

а

 

 

+

К

 

а

 

 

а

 

 

 

 

 

 

 

 

 

 

а

 

 

 

Используя выражение (80), запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

J L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

Mo J

а +

р

 

а

(5 •6(8)

£ІЄ =

 

ц 0

Г а + Р

 

 

= И.о

 

и поэтому с учетом формулы (79) третье уравнение для

определения А 0 ,

S f

и В 2 принимает форму

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 =

А0Т

-

 

(е^т

-

1) +

PS2 (е^7 "

- 1).

 

 

 

 

Решая совместно систему трех уравнений, находим

 

 

 

 

 

 

 

 

 

 

 

Р [(Р — а ) 2 е ? - ' г

— ( а +

 

Р ) 2 е ^ г

]

 

 

 

 

 

 

U P - - а ) 2

е*"7"

-- ( а

+ р ) 2 е ^ г ] ( р Г -Ь

4 р 2 ) _ 4 р 2 (ps -- а 2 )

 

 

 

 

 

 

 

2

-

а 2 ) [(Р +

а) ех>т

-

р

 

+

а]

 

 

 

 

 

 

К Р -

- а ) 2

Й ? " Г

( а + р ) 2 е ^ г ] ( р г н -

 

4р2) _

4 Р 2

( Р 2

— а 2 )

 

 

В.,

=

 

 

( Р 2 - - а 2

) [р +

а — ( Р - а)

е"'-г ]

фі

-

 

 

 

а )2

е'-.Г

_ ( а + р ) 2 ^ П

(РТ

+ 4р2)

_

4

2

)

 

 

 

[ ( Р -

 

р-2

_ а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• Процесс определения

k*{x)

заканчивается подстановкой До,

В Х

и В 2

в

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

 

 

 

 

 

 

2. Пусть теперь

g(t)=0.

В

этом

случае

имеем систему

с

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

налы m(t),

n(t)

и

u(t).

Для

такой системы среднее

значение

квадрата ошибки

определяется формулой

(82), если

положить

7"= со. Можно, как и выше, показать, что

интегральное уравне­

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

импульсной

переходной

функции k*(t),

обес­

печивающей минимум среднего значения квадрата ошибки

(82),

имеет вид

 

 

 

 

 

 

 

 

?

[R'm (т -

Є) +

Rn ( т -

Є)+R'u(r-Q)] k* (6) dQ =

 

 

 

 

 

 

 

 

 

 

(88)

 

 

 

 

CO

 

 

 

 

Интегральное уравнение (88) имеет решение

 

 

 

 

 

 

•к.і

 

 

 

 

 

 

j=T

/ = о

 

 

 

+ L ( p ) L * ( p ) M - 1 ( p ) M * - 1 (p) Я«(т) + j Я«(т-Є)х(Є)<Ю

3.

Наконец,

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

 

что g(t)=0

и

и ( 0 = 0 .

 

В данном случае среднее значение квадрата ошибки опре­

деляется

формулой (86) при

 

Т=оо,

а интегральное уравнение

принимает форму

 

 

 

 

 

 

 

 

 

 

со

/

 

 

f

 

 

 

 

 

со

 

 

I

[Rm (т - 6) + #„ (т -

 

6)] k* (8) d9

=

j ' Rm

(т — 9) х (0) d6.

(89)

О

 

 

 

 

 

 

 

 

 

—со

 

 

 

Его решение

можно

записать следующим

образом:

 

 

 

A*(T) = 2 B |

e

v +S £ / б < / ) <*) +

 

 

 

 

 

1=1

 

 

у= 0

 

 

 

 

 

+

L (р) Z* (р) М - 1

(р)М* - ' (р)

j

R'm

(т -

0) х (0) d9. .

 

 

 

 

 

 

 

 

 

 

— 00

 

 

 

.Полагая,

что

(т) = е ~ й

' ,

R„ (г) = С2 б (т),

Н (р) = 1 и №0 (s)

=

 

 

1 — ps

Ь(х, =

1 ± ї - в " ' - і - в

( т ,

, находим

 

 

 

1 - f

vs'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 т I

 

 

 

 

1 — B 2 a 2

 

 

 

 

(Y2 — P2 )

v~

 

Спектральная плотность

 

 

 

 

5 Т (со)

[2а + С 2

(сі2 + ы2)1

(1 + Р2 м2 )

 

 

2 -}-со2 ) (1 + Y 2 « 2 )

 

ф

 

 

 

лфтоэтому

 

 

 

 

 

 

і .

»

і =

l'C2 cz2 + 2g

 

1 .

<7 = — 1 .

л

 

g

; л, = — - р - •

,/И (р) Л1* (р) = [2а + С 2

2

- р2 )] (1 -

В2 р2 ); L (р) L * (р) = (а 2

— р2 ) (1 - у2 р2 )

Тогда

 

 

 

 

 

 

k* (т) = Вге^х

+ В2е%'-Г + 1 (р) L * (р) М-г (р) ЛІ*"1

(р) Х

СО

 

 

 

 

 

 

X J

 

е-«(т-Є) б (0) d6 = В і Є * - т + В2 е*"\

где

 

 

 

 

 

 

Для определения Bi и Вг подставляем

k*(x)

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

уравнение

489), из которого определим

 

 

 

 

 

1 — a2 v2 . a + ЛІ

a 4- л2

_i_ ii

' - ^

:

' ' 1

'—fJ v

^

2Y 2(a2 V2 —1)

 

 

 

 

5i

 

В Г

 

 

 

 

X

 

 

1 -f-a-f

 

 

 

 

 

 

 

Рассматривая последнее равенство как тождество, получим два следую­ щих уравнения:

о ф + 1

f

B t

_(

В,

 

 

Ях

+

Др.

_ Q

1 — ay

[ a +

A.!

a +

Я2 J

'

1 +

уХх

1 -j- ул2

Отсюда находим неизвестные В] и В2 :

 

 

 

 

 

(в +

Я 1 ) ( а + А д ) ( 1 +

т ^ | )

 

_

'(« +

М ( «

+

^) . (1+Т?ч)

1

 

( « Р + 1 )

 

'

2

 

 

(ap +

l J

^ - M

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

(ap + I) (Л, — Aj)

Соответствующая импульсной переходной функции &*(т) передаточная функция равна

( у а 2 С г

+ 2а

— аС) (1 — ар)

(ys

+

1)

 

W3 (s) = —

^

7

р

м

г т

 

' .

 

(ар +

1) (Cs -f у a2 C« + 2a)

(0s +

 

1)

 

Передаточная функция всей системы

 

 

 

 

 

 

 

( / a2C2 J. 2a

— aC) (1 — aB) (I — Bs)

Ф (S ) = Ws (s) U70 (s) = 1

1

1

v

 

p

/

v

— .

 

(ap +

1) (Cs +

/ a * C 2

+ 2a

) (ps

+ 1)

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

 

 

 

1

+ Р «

 

 

Можно

показать,

что

если Wo(s)—-

,

т. е.

принадлежит к

 

 

 

1 + ys

 

 

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

 

 

 

(S) =

( / C a » + 2a

- в С )

(l-f-та)

 

 

 

~

 

 

•и

 

 

(Cs 4- / О х *

4- 2a ) (1 4- ps)

 

 

 

 

 

 

 

 

 

 

/ C*s» +

2а — аС

 

 

 

 

Cs 4- / С2а2-+

 

 

Отсюда

видна та

существенная разница, которая имеет

место в этих

двух случаях. Можно также показать, что две последние формулы могут •быть получены и по расчетной схеме (см. рис. 3).

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

На основании схемы рис. 4 и формулы

(76)

можно

пред­

ставить

ошибку

воспроизведения системы в

следующем

виде:

 

6 (0 =

Я (р) [g (і) + т (01 -]u{t-x)b

(x) dx

-

 

T

 

b

 

 

 

со

 

— x — Q)]b (0) dQ +

— J k* (T) dx { [g(t — r — Q) + m(t — x — Q) + n(t

оb

 

-

T

 

 

CO

 

CO

 

 

 

 

 

(90)

 

f k* (T)

 

[ 6 (0) d0 j1 u (f T 0 — 0 ) b (a) da.

 

 

 

 

 

b

 

b

 

 

 

 

 

 

 

Пусть

среднее

значение

ошибки равно нулю,

тогда

 

со

 

 

 

 

 

Т

 

со

T — 6)6(8)d9.

(ftl)

j

 

 

T ) x ( t ) d r =

[ k*(x)dr

f

С учетом

формулы

(45)

выражение

(91)

можно

переписать

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j* 2a - P v

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

У

 

S а Л P -

 

9) 6 (T) d0 ,

V ~ T ) * (T)d T =•• J ** Wd x I

T -

—оо v = 0

 

 

 

 

 

0

 

 

0 v = 0

 

 

 

 

(92)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или для класса функций

(47)

имеем

 

 

 

 

 

оа

п

г

 

 

 

 

 

Т

 

со

п

 

 

і

2 a

v S

&д (т) ^ (0 *(т ) dx = J /г* (т) dx j

V av X

 

—со v = 0

ц = 1

 

 

 

 

 

0

 

6 v=0

 

 

 

 

 

 

X

У,

bl(x)Vl(t-Q)b(Q)dQ.

 

 

 

 

(93)

 

 

 

 

 

n=i

 

 

 

 

 

 

 

 

 

Если

рассматривать

уравнение

(93)

как

тождество,

то оно

определяет

(п+1)г

ограничений,

накладываемых

на

импульс­

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

 

 

 

 

 

 

 

 

 

 

 

 

 

<& =

]bl{x)k*(x)dx,

 

 

 

 

(94)

 

 

 

 

 

 

 

'

 

 

 

 

 

 

 

 

v = 0,

1,

2,

. .

. ,

п;

р, = 1,

2, .

. .,

г,

 

 

которые в данном случае зависят уже от оператора воспроизве­ дения Н(р), функции Pv (t) и импульсной переходной функции заданной части системы, а определение ограничений ведется пу-

т

{x)k*(x)dx,

тем решения системы уравнений относительно f b~^

 

получаемой из тождества (93).

 

С учетом ограничений (94) ошибку воспроизведения можно

записать в форме (81), а среднее значение квадрата

ошибки —

в виде (82).

 

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

рата

ошибки

и

одновременно

 

удовлетворяющей

 

ограничиваю­

щим условиям

(94). Для этого определяем минимум

выражения

(50). Придавая импульсной переходной

функции

k*(x)

вариа­

цию

Дт)*(т),

а

затем

выполняя

операцию

 

 

 

=0,

полу-

чим интегральное

уравнение

 

 

 

 

дД

д=о

 

 

 

 

 

 

 

 

 

 

 

 

 

] [R'm

(т — 9) +

R'n (т, -

0) + R'u ( т - 9)] k* (0) d9 = 21 І ^

^ W +

0

 

 

 

 

 

 

 

 

 

 

 

V—0 JJL=1

 

 

 

+

R'u (т) +

J Я,» (т -

 

9) х (0) d6,

0 <

т <

Т.

 

(95)

 

 

 

 

 

—со

 

 

 

 

 

 

 

 

 

 

 

 

Корреляционные функции, входящие в интегральное урав­

нение (95), определяются по формулам

(84).

 

 

 

 

 

 

Решение интегрального уравнения, если спектральная плот­

ность

5 Ф (со)

имеет

вид формулы

(24),

можно

записать сле­

дующим образом:

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

k* (т) =

п

 

г

 

 

 

 

q

 

 

 

 

 

 

 

V ' у>

Al bl (т) +

2 В, e%ix

+ v E

j

6<'> (т)

+

 

 

 

 

V = 0 | X = 1

 

 

 

l'=l

 

 

;'=0

 

 

 

 

 

 

 

+ 2 £>/ б ( Л (* - Л + L (р) L * (р) М - 1 (р) М * - 1

(р) X

 

 

X

j

/?m ( T - e ) x (8)d9 + /?L(T)

0 < т <

Т.

 

(96)

Определим

неизвестные

Лц,

5,-, £3- и Dj. Для этого

по за­

данному оператору Я ( р ) , импульсной

переходной

функции

Ь(х)

находим ограничения (94).

 

 

 

 

 

 

 

 

 

 

Пользуясь формулами (19) и (84), определяем Rm'(x),

R'n (т),

Ru'(x),Rn"{x),Ru"{*)-

 

 

 

(96), записываем импульсную

пере­

Пользуясь

выражением

ходную функцию k* (т).

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

импульсную

переходную

функцию

(96)

в

огра­

ничивающие

условия

(94)

и

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

уравнение

(95),

получим 21-\-(п+

\)г

линейных

 

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

уравнения.

Решая систему

21+ (п+\)г

уравнений, находим

 

неизвестные

Л^, Bi, Ej и Dj

и

подставляем

их

в

импульсную

 

переходную

функцию k* (т)

(96).

 

 

 

 

 

 

 

 

 

 

 

 

 

В

заключение

данного

параграфа

рассмотрим

 

особенности

определения

 

оптимальной

импульсной

переходной

функции

k*(x),

если

заданная

функция

 

времени

может

быть представ­

лена в форме (55).

 

 

 

 

 

 

 

 

 

 

 

 

 

Для данного случая на основании

расчетной

 

схемы

(рис. 4)

и формулы

(76)

 

ошибку

воспроизведения запишем

в

форме

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

т

|' T V е а? cos согт£* (т) dx

о

т_

С T V е а і х sin ©j xk* (т) dx. b

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

среднее значение квадрата ошибки — в виде

выражения (82).

Далее решается обычная вариационная задача,

позволяю­

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

k*(t),

обеспечи­

вающую минимум среднего значения квадрата ошибки (82) и

одновременно

 

удовлетворяющую

ограничивающим

 

условиям

(64). В результате

ее решения получим интегральное уравнение

 

f [Rm (т -

6) +

R'n (т -

Є) +

R'a

- Є)] k*

(Є) dQ =

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

I

 

t e

1

C ° S

 

x

+ 2 2

Tv

e ~ ' s i n ( 0 - T

+

 

-Hi

S

 

 

 

®i

 

v=0

 

 

 

 

 

 

 

; = 0

v=0

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

+ Ru(x)-\-

j ' Rm(x-Q)%(Q)dQ,

 

 

0 < т < Г .

 

(97>

 

 

 

 

 

 

CO

 

 

 

 

 

 

 

 

 

 

 

Интегральное уравнение (97) имеет решение

 

 

 

 

й*(т)=

2

2 4Tv e~a i T cosu),-T

+ 2

2 ^T v e ~ e ' T s i n a ) i T

+

 

 

t=0 v =0

 

 

 

 

 

i=0

v=0

 

 

 

 

 

+ V

Bt

c V

+ 2

 

6( / )

(T) +

2

60 >

(x-T)+L(p)

 

L * (p)

X

i = l

 

 

 

/ = 0

 

 

 

 

/ = 0

 

 

 

 

 

 

 

 

X M - i

(p)M*-'(p)

 

/?« (т —Є)и(Є)гіЄ

-i-

# в (т)

 

0 < T < 7*.

Путь

определения

неизвестных

А1Ч,

А^,

Ви

Ej

и

£>j

оста­

ется прежним.

Выше было показано, каким образом определить импульсную

переходную

функцию &*(т), если на систему действует сигнал

g(t) вида

(45) или (55), а система имеет две точки приложе-

яия воздействий. Из общего решения этих задач можно полу­ чить при определенных предположениях и некоторые частные решения. При решении задач данного параграфа предполага­ лось, что передаточные функции заданной части системы имеют только неминимально-фазовые множители в числителе. В об­ щем случае передаточные функции заданной части системы мо­ гут иметь еще и запаздывание. Ниже будет рассмотрен этот •более общий случай определения оптимальной импульсной пе- [реходной функции k* (т).

3. ПЕРЕДАТОЧНАЯ ФУНКЦИЯ СИСТЕМЫ, СОДЕРЖАЩЕЙ НЕМИНИМАЛЬНО-ФАЗОВЫЕ МНОЖИТЕЛИ В ЧИСЛИТЕЛЕ И ЗАПАЗДЫВАНИЕ

Предположим, что заданная передаточная функция {т. рис. 1) системы имеет неминимально-фазовые множители в числителе и элементы запаздывания, т. е.

WQ(s) = Wo(s)e-'>s.

Согласно расчетной схеме (рис. 4) и формулам (72) — (76) выражение для ошибки можно записать следующим образом

в (0 =

Hg

(р) g(t) + Hip)m{t)-Ju(t~i)b{T-

«О dx -

 

 

b

 

— ]k*(t)dTj[g(t

x e) + m (t T Q) + n(t

T — 0)] 6(6 —

о

b

 

 

— ТІ) dQ + ] k* (T) dx j b (6 ^ ) dQ X

о0

X J u(t — T - 0 — a)b(a — xjdo.

(98)

о

Пусть среднее значение ошибки равно нулю, тогда

Hg(p)g(t) =]k*{x)dxj

g(t-x-Q)b(Q-x1)dQ,

(99)

оb

где b(x—Ті) является импульсной переходной функцией, со­ ответствующей передаточной функции Wo(s) и

Учитывая, что заданная функция времени g(t) есть поли-

т

ном степени г, запишем g{t—т) в форме g{t—т) = ^^J(*)V X

V=0

X-±-pvg(t).

v!

При этом выражение- (99) принимает вид

 

 

 

 

 

Г

 

 

 

Т

 

 

 

 

со

 

 

 

 

Hg

(р) g (0

== V

( - 1)' ±

j" p'g (/ -

т) ft* (т) dx j'

(0 - T l ) d9. (100)

 

 

 

 

t=0

 

 

 

о

 

 

 

 

о

 

 

 

 

 

Рассматривая

выражение

(100)

как

тождество,

находим

(г+1)

ограничивающих

условий, которые

будут такими

же,

как

и условия

(79).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При выполнении ограничений (79) ошибку воспроизведения

(98) перепишем в форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

є (0 =

\

m{t — т) х (х) dx — j' u(t

— т) Ь (т — тг) dx

+

 

 

 

 

 

 

 

—со

 

 

 

 

О

 

 

 

 

 

 

 

 

 

Т

 

 

 

со

 

 

 

со

 

 

 

 

 

 

 

 

 

+

\к*{т)йт

\

b{Q— xx )dej и it — х — 9 — a) b(o — т,) da —

 

 

 

 

о

 

 

 

о

 

 

 

о

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

— 1' /г* (т) dx \ \m (t — x — Q)+n

(/—х—Є)] b (9 — xt) dQ.

(101)

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

Возводя в квадрат ошибку воспроизведения

(101)

и усред­

няя, получим

среднее

значение

квадрата

ошибки

 

 

 

 

 

со

 

 

 

со

 

 

 

 

 

Т

 

со

 

 

 

г'ск

=

\

х (х) dx

f

Я т

(х — 6) х (9) d9 — 2 J" /г* (т) dx

j' Ь (a xx ) da

X

 

 

— CO

 

 

 

— C O

 

 

 

 

 

0

 

0

 

 

 

 

X

7

flm(x

+ a - 9 ) x ( 0 ) d e + Jfe*(T)dT J

bio-xJdeX

 

 

 

 

— CO

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

X f

b ( v - x t ) d v

f [ ^ ( т +

0 - v - e )

r!-/?„^ +

c r - v - 0 ) ] X

 

 

b

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

оо

 

 

со

 

 

со

 

 

 

X

k* (В) d0 - f

j /г* (x) dx

[

b (a т,) dx j

b (v —

Tj) dv

( b {ц — xx) dr,

X

 

 

 

 

 

o

 

b

 

 

 

b

 

 

 

/

 

 

 

 

X

f b ( g - T i ) d g T\Ru(T+o

+

i\—v

 

Z-Q)k*(&)

d 0 - 2 X

 

 

оb

 

7*

со

со

со

 

X

$k*(x)dx

)" b(a — xjd a

j' b(v Xj)dv j' Ru(x + a

+

 

b

b

b

b

 

 

 

+ 9 — v)b

(9 —x^dS.

 

 

После

замены

переменных,

учитывая

обозначения

(19) и

(84), среднее значение квадрата ошибки перепишем в следую­ щей форме:

со

со

Т

со п

«& = j * (Х)d T

J #m (Т 9 ) Х

(6)d 9 2 Jfe* (Т)d T

J Я« (T

—со

—оо

0

—со

- Є + T l ) х (9) dQ + )' k* (т) dx

J [#,'„ (T -

9) +

R'n (x - 6) +

R'„ (x -

b

о

 

 

 

- 9)] ft* (9) dQ -!- &

(0) - 2 f &

(T +

T , ) ft* (T) dx.

(102)

 

b

 

 

 

Решая обычную вариационную задачу по отысканию им­ пульсной переходной функции k*{x), обеспечивающей минимум среднего значения квадрата ошибки (102) и одновременно удов­ летворяющей ограничениям (79), получим интегральное урав­ нение

] [Rm (т -

9) +

R'n (х -

0) +

R'„ (х -

9)] /г* (9) dQ

=

о

 

 

 

 

 

 

 

=

^ytf

+ R'u

(Т + Т 0

+

 

 

 

1 = 0

 

 

 

 

J'

Rm(x

9 + T1 )x(9)d9,

0 < т < 7 \

(103)

Рассуждая, как и ранее, находим решение интегрального уравнения (103):

ft* (х) =

V

Л,т' + 2 5,

е?-'х + 2

 

(т) + J] D; . б( / > (т -

Т)

 

1=0

 

1 = 1

/=о

 

/=о

 

 

 

+

L(p)L*

ф)М-і

<р)М*-1(р)

j

і? т — в +

Tj) X

 

 

 

 

 

 

—CO

 

 

 

 

 

X

х (9) d9 +

(т + T l )

 

0 < т < 7 \

 

(104)

Порядок

определения

неизвестных

Л,, В І, Ej

и

Dj остается

прежним.

 

 

 

 

 

 

 

 

 

Из формул

(102) — (104) видно, что

среднее

значение

квад­

рата ошибки и оптимальная импульсная переходная функция являются функциями запаздывания х\.

Рассмотрим некоторые частные случаи, вытекающие из об­

щего решения.

 

 

 

 

 

 

 

1. Если положить u(t)

=0,

то

 

 

 

 

г'ск

= j '

х (т) dT j

Rn

(х —

9) х (9) d9 — 2 \ ft* (х) dx

X

 

—со

—со

 

 

 

0

 

со

'

 

 

т

т

(т - 9)

+

X ]

Rm (т - 9 + T j х (9) dQ +

J А* (т) dr j

 

 

+

R'„(x-Q)]

ft*(9)d9;

 

 

 

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