книги из ГПНТБ / Солодовников В.В. Расчет оптимальных систем автоматического управления при наличии помех
.pdfто импульсная |
переходная |
функция на основании |
формулы |
|
(29) |
имеет |
вид |
|||||||||||||||||||
|
|
|
|
к (х) = |
Л0 + |
Aiz + Л3 т- + |
£об (т) + |
D0 6 (х - |
Т). |
|
|
|
||||||||||||||
После |
подстановки |
fe(x) |
в |
формулу |
(40) |
получим |
три |
|
уравнения |
для |
||||||||||||||||
определения |
неизвестных Л о, Ль Л2 , Еа |
и Д>: |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
Аг, |
|
Лх |
|
2Ао |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
а |
|
а 2 |
|
a J |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л 0 |
|
|
АХТ |
|
, |
Л, |
, |
А а Т 2 |
, |
2Л2 Г |
, |
2А„ |
|
D o |
= |
0; |
|
|
||||||
|
|
а |
|
|
|
а |
|
|
а- |
|
|
а |
+ |
а 2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
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(—Иг), |
|
|
|
8Ла |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Иг'= — ' |
|
|
|
|
|
|
|
|
|
|
||||||||||
Кроме того, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
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Л2 |
|
|
|
. 4 |
|
A T |
|
Л , |
|
|
|
Л 2 Т 2 |
|
2Л2 |
|
||||||
|
|
_.о |
|
|
|
|
• £0 = 0; |
|
4- |
|
4- |
|
|
4- |
a |
|
|
4- |
4- |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
a ^ |
a ^ a 2 |
|
т |
|
|
^ a' ^ |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
2Л2 Г |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
— Г " - D o = |
0; |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
+ |
|
|
+ —— |
+ Е0 + Ц, = |
|
1; |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
Л 0 Г 2 |
|
А{Г* |
|
Л2 |
Г* |
+ |
Z V = |
|
0; |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
А0Т* , Л І Г * |
|
Л Г Т 5 |
|
|
|
8Лг |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
о |
+ ——+ |
|
- 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 |
6д СО ^ |
(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„ (со) = |
2а |
„ . |
|
|||||||
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
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 |
4ш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 |