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

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

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

то импульсная

переходная

функция на основании

формулы

 

(29)

имеет

вид

 

 

 

 

к (х) =

Л0 +

Aiz + Л3 т- +

£об (т) +

D0 6 (х -

Т).

 

 

 

После

подстановки

fe(x)

в

формулу

(40)

получим

три

 

уравнения

для

определения

неизвестных Л о, Ль Л2 , Еа

и Д>:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аг,

 

Лх

 

2Ао

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

а 2

 

a J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л 0

 

 

АХТ

 

,

Л,

,

А а Т 2

,

2 Г

,

2А„

 

D o

=

0;

 

 

 

 

а

 

 

 

а

 

 

а-

 

 

а

+

а 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

4Л,

 

 

 

 

 

 

 

 

 

 

2Л,

 

 

 

 

 

 

50 = —а

 

+ •

а 3

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

Далее из уравнения

(41)

следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

So =

( -

1)° [Роо ( -

1)° +

Рю ( - ! ) '

+

Р 2 0

( -

О3 !

 

j'

 

6(x)dx

—[x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CO

тб (T) dx — ці

 

Si

= ( -

І) 1

ІРої ( -

 

1)° +

Рп

( -

I ) 1

+

Раї ( -

І)2 ] і J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

So

=

( -

I ) 2

[p0 2 ( -

 

1)0

+

p1 8

( -

l)i

+

p M

( _

l)2 ]

 

j

x26

(x) dx — |.i2

 

 

 

S 0

=

Poo (1 —Ho).

 

Ho =

1 ;

Si =

Pu(—

m)

 

fi! =

0;

 

 

 

 

 

 

 

 

 

S 2

=

P2 2(—Иг),

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Иг'= — '

 

 

 

 

 

 

 

 

 

 

Кроме того,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0

=

1 =

j' [A0 +

Л ^ + Л 2 т 2

+ £ 0 б (т) +

D0 6

 

(x -

T)}

dx;

 

 

 

 

jix

=

0 =

f x [Л0 + Лх х +

Л3 т2 +

£ 0 б (х) +

D06

(х -

Т)] dx;

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8А->

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из =

 

 

 

~

=

f та 0

+

Лх х +

Л2 х2 +

£ 0 б (т) +

 

D„6 (х -

Г)1 dx.

 

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

определения

неизвестных

Ла, Л ь

Л2 , £ 0

и D0 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лі

 

 

2

 

 

 

. 4

 

A T

 

Л ,

 

 

 

Л 2 Т 2

 

2

 

 

 

_.о

 

 

 

 

• £0 = 0;

 

4-

 

4-

 

 

4-

a

 

 

4-

4-

 

 

 

 

 

 

 

 

 

 

 

 

 

a ^

a ^ a 2

 

т

 

 

^ a' ^

 

 

 

 

 

 

 

 

 

 

 

 

2 Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

— Г " - D o =

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

+ ——

+ Е0 + Ц, =

 

1;

 

 

 

 

 

 

 

 

 

 

 

 

Л 0 Г 2

 

А{Г*

 

Л2

Г*

+

Z V =

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А0Т* , Л І Г *

 

Л Г Т 5

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

о

+ ——+

 

- 7 — +

0 . 7 - » + — f = 0.

 

 

 

 

 

 

 

 

 

 

 

4

 

 

5

 

 

 

 

от2

 

 

 

 

 

 

 

Решая эту систему уравнений и подставляя значення /10, А\, А2, Еа н Da в k(x), найдем оптимальную импульсную переходную функцию, соответст­ вующую сформулированным выше условиям.

6. ОПТИМАЛЬНЫЕ ХАРАКТЕРИСТИКИ В СЛУЧАЕ ЭКСПОНЕНЦИАЛЬНЫХ И ГАРМОНИЧЕСКИХ ДЕТЕРМИНИРОВАННЫХ ВОЗДЕЙСТВИЙ

В § 1—5 были получены результаты, связанные с оп­ ределением оптимальных динамических характеристик, когда не­

случайная

функция

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

степе­

ни г.

 

 

g(t)

Однако

в ряде

случаев заданную функцию времени

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

Решим задачу определения оптимальной импульсной пере­ ходной функции, отличающуюся от задачи, рассмотренной в § 1,

лишь тем, что функция

g(t)

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

t, а функцию вида

 

 

g(t)

= J]

а 4 в в ' ' c o s (0/(0-

 

1=0

 

Для расчета импульсной переходной функции по-прежнему можно пользоваться схемой рис. 3, для которой преобразования Лапласа эквивалентных управляющего P(t) и возмущающего Q(t) воздействий имеют вид

P(s)

= H (s) [G (s) + M (S)] - W0

(s) U (s);

(42)

Q(s)

= G(s) + M (s) + N(s)-H(s)

[G(s) + M(s)].

 

С учетом выражений (42) и (11) ошибку воспроизведения можно записать в следующей форме:

 

оо

со

E(t)=

J [g(t — t) + m(t T)]K(t)dt—\u(t

T)b(x)dx —

 

—оо

0

_ j r g ( / _ T

) - i - m(t

— x) + n(t x)]k(x)dx

+

o-

Г

oo

 

 

 

 

 

 

 

 

+ J' k (T) d x \

u(t — x —

a)b (a)da.

(43)

 

b

b

 

 

 

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

M[&(t)]=

J

g(t — г)х(т)<2т ]g{t — x)k(x)dx

= 0

или

—оо

 

 

О

 

 

 

Т

 

 

оо

g(t x)K(x)dx

 

 

j

= §g(t

— x)k(x)dx.

(44)

—оо

0

Предположим, что в тождестве

(44)

 

 

 

 

g(() = j]a4Pv(ty,

 

 

(45)

тогда J 2 cv Л, (/ -

т) х (0 dx s

j 2J av Pv

(f -

T) /г (т) dr.

(46)

-^co v=0

 

0 v=0

 

 

 

Далее рассмотрим

класс гармонических

и

экспоненциальны*

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

ставить Pv

(t—т)

следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

 

Л , ( * - т ) =

2

 

bl(x)Wl(t),

 

 

 

(47)

где

b^

(х)

и

"Ц^ (^) 1.1=1, 2,..., г) являются линейно независи­

мыми функциями.

 

 

 

 

 

 

 

 

 

 

со

С учетом

выражения (47)

тождество

(46) перепишем в

виде

п

 

г

 

 

 

Т

п

г

 

 

 

 

 

J

Ъ

G v 5

b

l

(т ) ^

(0 « 00 ^ ^ J

2

flv

2

СО ^

(0 k СО Л .

(48)

—cov=0

ц=1

 

 

 

0 v=0

ц=1

 

 

 

 

 

Тождество (48) дает (я-Н)/" ограничений

на

импульсную

переходную

функцию, определяемых

преобразующим

операто­

ром

Н(р):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ql^T\bl(x)k(x)dx,

 

 

 

 

 

(49)

 

 

 

 

 

 

о

/г; ( . 1 = 1 , 2 , . . . ,

г;

 

 

 

 

 

v

=

0,

1, 2, . . . .

 

 

коэффициенты

<2д

имеют постоянные значения.

 

 

 

 

 

Если тождество (46) удовлетворяется, то в данном случае

будут справедливыми формулы (17) и

(18) для сигнала

ошибки

и среднего

значения квадрата

ошибки. Теперь

задача

состоит

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

k(t),

об­

ращающую

в минимум г2ск

и одновременно удовлетворяющую

условиям

(49). Для решения

этой задачи составим

выражение

v=0 ц=1

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

]

[Rm

(т- 8) +

Rn (т- 8) + К (т - 6)1 k (8) dQ =

Я ; (т) +

о

 

 

к (8) d8 + 2 2 W

 

т < Г. (51)

+

j

Rm (х-8)

0 <

—со

V=0 ц=І

Рассмотрим решение интегрального уравнения для стацио­ нарных случайных процессов с дробно-рациональной спектраль­ ной плотностью (24). Пользуясь формулой (25), сведем инте­ гральное уравнение (51) к дифференциальному

т

«>

 

 

М (р) М* (р) j G(x — Q)k (0) dQ =

j ' Rm

(T - в) x (0) dQ + R*a (т)

+

О

—со

 

 

-г- V V

yv &v

lT\

 

V=0 |x=l

 

 

 

Неоднородное дифференциальное

уравнение порядка

2k

имеет решение

 

 

 

 

f G ( t - 0 )

fe(9)d9=2

2

4^

+ 2

5 ' e

V

+

 

 

 

О

 

 

 

 

\-=0ц=1

 

1 = 1

 

 

 

 

 

 

 

 

 

 

f

/ ? я ( т - в ) х ( в ) ( Ю +

/?;(т)

0 <

г <

Г,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(52)

где

— корни характеристического уравнения

 

 

 

 

 

 

 

 

 

 

 

М(\)М*{Х)=*

0.

 

 

 

 

 

 

Применяя

к

обеим частям выражения (52) оператор

L(p)L*(p),

 

 

на основании формулы

(26) получим

 

 

 

 

 

k (*)•=

 

S

2 «

(х) +

І

В / ' т

+ І

£ / 6 ( / )

(г)

+

 

 

 

 

 

 

v=0n=l

 

 

i = l

 

/ = о

 

 

 

 

 

 

 

 

 

+

 

Ds

(x-T)

+

L<p)L*(p)М-1

(p)

X

 

 

 

X

 

й-

1

(P)

 

JRm(x-Q)K(Q)dQ

 

 

+ Rl(T)],

0 < т < Г .

(53)

M * "

 

 

 

 

 

 

 

 

 

 

-oo

 

 

 

 

J

 

 

 

 

 

Неизвестные постоянные І4^, 5,, Ej и £)j определяются по

следующей схеме.

 

 

 

 

 

 

k(t)

 

 

Подставляя

 

импульсную

переходную

функцию

из

вы­

ражения

(53)

в

интегральное уравнение (51), находим 21 ли­

нейных

алгебраических уравнений. Еще

(п+\)г

линейных

алге­

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

Решением

системы 2 1+г(п+1)

алгебраических уравнений и

подстановкой

найденных значений

Л*, Bit Ej и Dj в формулу

(53) заканчивается процесс отыскания оптимальной импульсной

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

 

Из общего решения, уравнения (53)

здесь могут быть полу­

чены те же частные случаи, что и в § 4.

 

2 Зак. 1249

33

1. Если предположить, что g(t)= О, то получим результат, полностью совпадающий со случаем, разобранным в § 4 [см.

формулы

(30) —(33)].

 

 

 

 

 

 

 

 

 

 

 

 

2. Предполагая, что u(t)—0,

 

получим систему

с одним

вхо­

дом, эквивалентные

воздействия которой имеют вид

 

 

 

 

 

 

P(t)

= H(p)[g(t)

+

m(t)};

 

 

 

 

 

Q (0 =

g (0 + m(t)

+ п (t) -Н{р)

[g(t)

+ m (t)].

 

 

Квадрат среднего значения ошибки может быть представлен

соотношением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<&= |

 

f #M(T-e)x(0)de-2p(T)dT J

 

Rm(x-B)X

X «(0)d0 +

j k (T) dx J [tfm

(T -

0)

+

# „ (T -

0)] k (0) dQ.

 

Интегральное

уравнение

и его решение принимают вид

+

j [Rm

(т -

0) +

Я „ (т -

6)] /г (0) d0 =

J Rm

(т -

0) х (0) d0

Т

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

—со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(54)

 

 

 

 

v = 0

 

ц=1

 

 

 

 

 

 

 

 

 

 

 

w = Е Ел

 

 

 

 

2ft

 

 

 

 

 

 

 

/ г

 

^ w + Е5 'е ' + Е £ / 6

w +

 

 

 

v = 0 ц = |

 

 

Г) +

i = 1

 

 

j—0

(р)

X

 

 

 

 

/=о

D,.6(" (т -

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х М *

'(р)

J

 

 

tfm(t-0)x(0)d0

0 < т < 7 \

 

Предположим, что

 

 

 

 

 

 

 

 

 

 

 

 

 

S (t)

= cos оУ;

Н (р) =

tf0;

 

m (/) = 0;

S„ (со) =

„ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

coJ +

а-

 

В этом

случае

k=0;

1=1;

q=0;

 

М{р)М*

)=2 а;

 

 

 

 

 

 

 

 

 

L(p) L* ( Р

) = а 2 - р 2 ;

 

 

 

 

 

 

Р it

т) =

cos (u0t cos со0т +

sin co0<-sin со0т;

 

 

 

 

 

 

 

 

 

я + 1

= 1;

r =

2.

 

 

 

 

 

Ha основании выражения (53)

можно записать, что

 

 

 

 

k (т) = А\ COS СО0Т +

Ай2

sin со0т + £„6 (т) -f.D0 6 (т — Т),

0 <

т < Т.

Подставляя импульсную переходную функцию k(x)

в тождество (49),

найдем

 

J #„ (cos ш0т cos соо0 + sin WQT sin coo0) 6 (0) dQ =

—CO

 

T

 

= Jj [cos со0т cos coo0 + sin co0 Tsin co„0] [A® cos co08 +

A\ sin co00 -j-

0

 

 

 

 

 

+ £ 0 б (0)+D o 6 (Q-T)]

dQ.

 

 

 

 

 

Отсюда

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

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

А®, А\,

 

Е0

И D0:

n /

Т

 

1

\

А2

 

 

 

 

 

 

 

#0 = А\

{-у

+

- ^ - s i n 2a>0TJ +

 

(1-cos2 со0Г) +

Е0

+

D 0

cosсо0 Т;

0 =

— ' - ( l

— cos2co0 T) + А\ ( —

sin2co0T ) +

D 0 sin со0 Г.

 

4co0

 

 

V 2

0

У

 

 

 

 

Подставляя

 

значение k(t)

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

уравнение

(54),

найдем

т

 

 

 

 

 

 

 

 

 

 

 

 

 

j е _ а

( [ _ 0

)

[Л° cosшо 0 +

А\ sin соо'0 + Я0 б (0) +

D0S (0 — Т)] dQ =

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 7j cos со0т +

у"s ' n

ш о т -

 

 

 

 

 

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

 

 

 

 

 

 

 

 

АI щ -

А\ а + Е0

(ofi + <og) = 0;

 

 

 

 

Л] со0 sin CO0JT — a cos со0Г — А® а sin со0Т —

-А\ Ш0 cos ш„Г + D0 (а* + cog) = 0.

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

дящие в импульсную переходную функцию k(x).

 

 

 

3. Если дополнить

второй случай, условием

g(()=0,

то получим четвер­

тый случай § 4.

 

 

 

 

 

 

7.

ОПТИМАЛЬНАЯ

ИМПУЛЬСНАЯ

ПЕРЕХОДНАЯ

ФУНКЦИЯ

В СЛУЧАЕ ВОЗДЕЙСТВИЯ ОБЩЕГО ВИДА

 

 

Найдем

оптимальные динамические

характеристики

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

более об­

щего класса,

а именно:

 

 

 

 

 

 

g(0 = S

ад,(1)еа1'cos

щі,

 

 

или

 

 

 

 

 

 

5(0 = S1=0 a/<7i(0ee''sin<o/*.

Предположим, как и ранее, что на первый вход системы по­ ступают управляющее воздействие y(t) и возмущающее воздей-

2* 35

ствие n(t). Управляющее воздействие является суммой двух со­ ставляющих:

 

 

y{t)

= g{t) +

m{t),

 

 

 

где

g(t)

— заданная функция времени;

 

 

m(t)

и n(t)

— стационарные случайные

функции

с нулевыми

 

 

средними

значениями

и с

заданными корреля­

 

 

ционными

функциями

Rm(x)

и Rn(т)

[спектраль­

 

 

ными ПЛОТНОСТЯМИ STO(to)

и S n (©)] .

 

На второй вход системы поступает сигнал

u(t).

Функция g(t) имеет вид

 

g-(0 = V af,®,

(55)

(=0

 

где

 

' Л-(0 = ? ( ( 0 ^ с о з с о , . ^

 

Задача может быть сформулирована аналогично § 1 и отли­ чается лишь тем, что функция g(t) представлена в виде фор­ мулы (55).

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

формулы (42).

 

 

 

 

 

 

 

 

 

 

 

Если

потребовать

равенства

нулю

среднего

значения ошиб­

ки воспроизведения

(43), то приходим к выражению

 

 

 

оо

 

 

 

 

 

Т

 

 

 

 

 

 

J

g(t — т)х(т)Аг =

$g{t

x)k(x)dx,

 

 

 

— О О

 

 

 

 

0

 

 

 

 

 

которое

с учетом

формулы

(55) можно переписать в виде

 

 

2

a^P^t

 

— x)k(x)dx

=

V а(

J

Р,(t — х)%(т)dx.

(57)

 

i=0

0

 

 

 

 

і=0

—оо

 

 

 

Если

в

выражении

(57) щ — произвольные

коэффициенты,

a Pi(t) —независимые

функции, то должны быть равными

чле­

ны с одинаковыми аг-, т. е.

 

 

 

 

 

 

 

 

 

JPt

(t — x)k (т) dx =

j

Pt

{t — x)% (т) dx,

(58)

 

 

О

 

 

 

 

оа

 

 

 

 

 

 

 

 

 

 

i = 0,

1, 2

 

 

 

г.

 

 

С учетом выражения

(56) уравнение

(58) можно переписать

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[ qt

(t — т) eai ( ' -

т ) cos со,. (t — x)k (т) dx

=

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

qi(t

т)еаі

( i ~ r ) cosal(t

x)K(x)

dx.

 

 

 

=

f

 

 

 

—CO

 

 

 

 

 

 

 

 

 

 

После сокращения на

еа«

имеем

 

г

^

 

 

 

|" Ш (t т

) е а'%

c o

s и і (t — t)k (х) dx =

 

о

 

 

 

 

со

 

 

cos со,- (^ — т) я (т) dx.

(59)

= j cjі (t — х) в~аіх

—со

 

 

 

 

Выражение (59) можно

свести

к виду

 

т^

J Яі

% ) е a / T

l c o s ю ^ c o s < в г т " г sin co^ sin согх] A (x) dx =

0

 

 

 

CO

q1; (^ — T) e _ a t T (cos co^ cos co/t-f sin cousin со/х) x (x) dx.

f

 

CO

 

 

Из последнего соотношения, путем приравнивания членов при

синусах

и косинусах

с различными а,-, имеем

т_

\qi(t-—x)e

а^coscottcosa>ixk(х)dx—

о

 

оо

х)е _ а 'т со5 0)^со8(0г гх(х)йх;

== j qt(t

—со

 

т

о qt (t — х) е а ' т sin сог/ sin со,хй (х) dx —

со

=j qi(t — х) е ~ а ' т sin cousin со£ хи (х) dx

—со

или

Т

_

 

 

со

 

 

 

_

 

(60)

\qi(t х)е a «T cosco,xA(T) dx =

j

 

— х)е a ' T cos со,ти (х) dx;

ЇВ

 

 

 

—оо

 

 

 

 

 

Г

 

 

 

оо

 

 

_

 

 

f<7/(/ — х)е a'tsincoiTfe(T)dT =

f

ql{t

— x)e

a «'T sinco,xx(x)dx,

(61)

0

 

 

—-CO

 

 

 

 

 

 

 

 

1 = 0, 1,

2,

3, .

.

. ,

r.

 

 

Если

теперь

в

выражения

(60)

и

(61)

подставить

вместо

функции

qi(t)

ее значение из формулы (56), получим

 

 

 

Г

 

 

 

оо

 

 

 

 

 

 

fxv e~a 'T cosco,TA(T)dT=

j xV~a 'T cosctvtK(T)dx;

 

(62)

 

О

 

 

 

—со

 

 

 

 

 

Т

 

 

 

со

 

 

 

 

 

J'xv e atzsin<uixk(x)dx

=

J Tv e_ a iT sincoi xx(x)dT;

 

(63)

 

0

 

 

—со

 

 

 

 

 

v = 0,

1,

2, 3, . . . , і;

і=0,

1,

2, . . . , г.

 

 

Из выражений

(62) и (63), полагая, что правые части пред­

ставляют

собой

постоянные

заданные

значения,

получим

(r _l_ і) (/- + 2) ограничений на импульсную переходную функцию k(t):

 

 

 

т

 

 

 

 

 

 

 

 

j1 T V e~aix

cos со; т/г (т) dx = Qv J

 

 

 

 

 

b

 

 

.

 

(64)

 

 

 

T

 

 

 

 

 

 

 

j' xv e_ t x 'T

sin со,- т/г (т) dx = Qv .

 

 

 

 

 

о

 

 

 

 

 

Если

уравнение

(57)

удовлетворяется

тождественно,

то,

ошибку

воспроизведения

можно записать в

виде выражения

(17),

а

среднее

значение

квадрата

ошибки

в форме выраже­

ния

(18).

 

 

 

 

 

 

Далее задача состоит в том, чтобы найти импульсную пере­

ходную

функцию

k(t),

обращающую

в минимум величину

г2ск

иодновременно удовлетворяющую ограничивающим условиям

(64).

Для решения этой задачи составим функционал

• / = « 4 - 2 2

v

Y ; Q t - 2 S І vvQi-

i =0

v=0

j = 0 v = 0

Решая, как и выше, вариационную задачу, приходим к не­ обходимому и достаточному условию минимума выражения (18) в виде интегрального уравнения, которому должна удовлетво­ рять импульсная переходная функция k(t):

f [Rm 6) - f

Rn

(т - 0) + R'u (т -

0)]

k (0) dQ =

 

b

 

 

 

 

 

= Ru (T) + f Rm

(x - 0) x (9) dQ +

2

v; *v e~a'T

X

- co

 

i=0

v=0

 

 

XCOSCU( -T +

2 2J YvTv e ~ ^ T s i n c o ,T.

(65)

Решение полученного интегрального уравнения проведем для класса стационарных случайных процессов с дробно-рацио­ нальной спектральной плотностью (24). Пользуясь формулой (25), сведем интегральное уравнение (65) к неоднородному диф­ ференциальному уравнению порядка 2 k:

М (р) М* (р) f G (т - 0) k (0) dQ =

j

Rm

(т - 0) x (0) dQ + Ru (x) +

'O

 

- co

 

+ 23 2

Y ; T v e _ a ' T cosco ,T +

2

2 Y^v e~a 'T sinco,T,

1=0 v=0

 

1=0

v=0

 

0 <

T <

7\

(66)

Решение неоднородного дифференциального уравнения (66) имеет вид

( G (т — 9) к (8) dQ =

2 2

T v

e _ V

cos со,- т +

 

О

 

 

 

1=0 V=0

 

 

 

 

+ 2 2 A*

e

"l' sin <B,T +

22k Bi e^ic

+

L (p) L * (p) М-*- (p)

X

r=0 v=0

 

 

 

i = I

 

 

 

 

x

M«-]

ip)

/ ? m ( T - 6 ) x ( 8 ) d e +

^ ( T )

(67)

Применяя к обеим частям решения дифференциального урав­ нения (67) оператор L(p)L*(p), с учетом выражения (26), имеем

k (т) =

2

2

А * t V e - C 8

' T c o s

<°/т

+ 2

2

^

t V е _ Е г Т s i n ш ' т +

 

 

 

; = 0 V=0

 

 

 

 

(=0 V=0

 

 

2k

е '

+

£ <Р)

(Р)

 

(р) М * - '

(/>)

J

/ ? m ( T - 0 ) x ( 0 ) d 0

+

+ 2 5/

 

1=1

 

 

 

 

 

 

 

 

 

 

 

 

+

Ru (т)

+

2

£ /6 ( /

) ( т)

+ 2

Dj8U)(x-T),

 

0 < т < 7 \

(68)

 

 

 

 

/ = 0

 

 

/ = 0

 

 

 

 

 

Неизвестные постоянные v4'v, i4'v, Ви Ej и Dj определим в следующей последовательности. Импульсную переходную функ­ цию (68) подставляем в интегральное уравнение (65). В ре­ зультате находим 2 / линейных алгебраических уравнений для определения неизвестных Alv, Alv, Bh Ej и Dj. Еще (r+1) (r+2) уравнений получим, если импульсную переходную функцию (68) подставим в ограничивающие условия (64).

Решая систему алгебраических уравнений порядка 2 / + ( г + + 1)(/-+2), находим неизвестные, входящие в импульсную пе­

реходную

функцию

(68).

Подставляя найденные

значения А1Ч,

А ,ВЬ Ej

и Dj в формулу

(68),

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

импульс­

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

 

 

 

 

Более

подробно

интересно

рассмотреть лишь

один

случай,

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

веденных выше.

 

Если предположить, что u(t)=0,

то получим систему с од­

ним основным входом, для которой среднее значение квадрата ошибки определяется по формуле (34), а интегральное уравне­ ние имеет вид

J lRm(т-Є)

+.R„(т-9)]

k(0)d0 = j

Rm(r-Q)*(9)

dQ +

 

 

 

CO

 

+ 2 2 y'vtVe

1 c o s и < т + S 2 Y v х * е

1 s i n M < T - 0

< T < T - ( 6 9 ) '

i = 0 v = 0

i'=0 v=0

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