Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 348

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

582 ГЛ . 14. О П Р Е Д Е Л Е Н И Е О П ТИ М А Л ЬН Ы Х Л И Н Е Й Н Ы Х СИСТЕМ

По заданию требуется спроектировать систему с памятью длительности 7. Точное выполнение этого условия невозможно, так как на основании общих теорем теории дифференциальных уравнений система с дробно-рациональной передаточной функцией всегда имеет бесконечную память. Однако практи­ чески можно считать память системы равной времени переходного процесса. Для определения времени переходного процесса находим по формулам (4.4.50) и (6.5.4) весовую и переходную функции спроектированной системы. Мы пре­

доставляем

читателю

проделать необходимые

выкладки самостоятельно

и приведем

только

графики весовой функции

(пунктирная кривая на

рис. 14.5.4) и переходной функции (рис. 14.5.11). Из рис. 14.5.4 и 14.5.11

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

переходная функция h (|) не отклоняется от ее установившегося значения h (оо) больше чем на 5%) около 1,55Г, что также можно считать удовлетвори­

тельным.

Для проверки точности спроектированной системы находим по фор­ муле (7.4.4) установившуюся дисперсию ее выходной переменной. Имея в виду что спектральная плотность стационарного белого шума интенсивности к

равна к/2я, получим

 

 

 

 

 

 

к

(*

 

 

0,58174 (ій>)4-32,6872(гсо)2 + 256

 

(14.5.12)

*1

j

I о,273 (і(1))3+ і ,762Г2

(/ш)2 +

7,62Г(іш)+16

[2

 

 

 

— ОО

 

 

 

 

 

 

 

Производя

замену

переменных

7ш =

р, получаем

 

 

 

к

1

 

ОО

0,581 (г'р)4— 32,68 (гр)2 + 256

^

 

f]—Dw

 

Г

(14.5.13)

Т ' 2л

J

|0,2(ip)3 +

l,762(tp)2 +

7 ,6 2 (+ ) + 16|2 ^

 

 

Пользуясь

для

вычисления этого интеграла формулой

(7.7.6), находим

 

 

 

 

 

 

п ~

4’Ш

 

 

 

 

 

 

 

Ц—

m

 

 

Формула (13.6.17) дает в данном случае для дисперсии ошибки оптимальной

системы выражение

4 к


§

14.6. М ЕТОД К А Н О Н И Ч ЕС К И Х Р А ЗЛ О Ж Е Н И Й

583

Таким образом,

спроектированная система и по точности достаточно близка

к оптимальной.

Предоставляем читателю самостоятельно проверить, насколько близка к оптимальной система, которая получится, если аппроксимировать логариф­ мическую амплитудную характеристику оптимальной разомкнутой системы двумя полупрямыми с наклонами —40 дБ/дек п —20 дБ/дек, сопрягающимися при соТ = 4.

§ 14.6. Определение оптимальной линейной системы методом канонических разложений

Уравнение (14.1.1) может быть в общем случае решено методом канонических разложений. При этом решение получается в форме бесконечного ряда, конечным отрезком которого можно пользо­ ваться для приближенного определения оптимальной линейной системы.

Представим случайную функцию X (т) которая, как мы помним является суммой нерегулярной части полезного сигнала и помехи, каким-нибудь ее каноническим разложением в интервале наблю­ дения t Т ^ т ^ t:

X (t) = 2 F > v(t) (f —Г С тСО-

(14.6.1)

Как известно, это всегда можно сделать для любой случайной

функции

с

конечной дисперсией. Для этого достаточно найти

систему пар

функций аѵ (т) и хѵ (т), удовлетворяющих

условиям

(см. [54],

§

6.2 или [53],

§

60)

 

 

 

t

 

 

 

 

 

j

а\,(т)ан(т)(1т = 0У|г,

(14.6.2)

 

 

1 -Г

 

1

 

 

 

 

 

 

 

 

хѵ(ъ) = - щ

j Кх (т, o)av(a)do,

(14.6.3)

 

 

 

 

t-т

 

где для краткости через 6Ѵ(1 обозначены числа, равные единице при одинаковых индексах и нулю при различных индексах (бѵѵ = = 1, 8Ѵц = 0 при р =т^=ѵ). Предположим, что функции аѵ (т)

ихѵ (т), удовлетворяющие условиям (14.6.2) и (14.6.3), найдены,

ибудем искать решение интегрального уравнения (14.1.1) в форме

ряда

g(t,

____

 

 

т) = 2£ѵ«ѵ(т),

(14.6.4)

 

 

V

 

где

— неопределенные

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

 

Подставляя выражение (14.6.4) в уравнение (14.1.1) и выпол­ няя почленное интегрирование ряда, получим

1

 

2 £ѵ f Кх (х, o)av(r)dx=:f(t, о)

(t — Г < о < < ). (14.6.5)

V t - т


584 ГЛ . 14. О П Р Е Д Е Л Е Н И Е О П ТИ М А Л ЬН Ы Х Л И Н Е Й Н Ы Х СИСТЕМ

Но согласно

формуле '(14.6.3)

 

 

 

t

 

t

 

 

 

j Кх (х,

о) av (%)dx = j

Kx (a,

т) av(x) dr = Dvxv (о).

t - T

 

t - T

 

 

 

Подставляя

это выражение в (14.6.5), получим

 

 

2

А £ѵ*ѵ(<*)= /(*,

°)

(« —Г < о < 0 -

(14.6.6)

 

V

 

 

 

 

Таким образом, для определения неизвестных коэффициентов £ѵ достаточно представить правую часть интегрального уравнения (14.1.1) / (t, а), рассматриваемую как функция второго аргумента,

разложением

по

сопряженным

координатным функциям х ѵ (а)

в интервале наблюдения t Т ^

а ^

t.

 

Для

определения коэффициентов £ѵ умножим равенство

(14.6.6)

на

йц (ст) и проинтегрируем

результат по о в

пределах

от t Т до t.

Тогда получим

 

 

 

 

2

А &

t

 

t

 

 

j«в(°)xv(а) da —

j av-ia)f{t, а) da.

(14.6.7)

 

V

 

t-т

t-T

 

Все члены ряда в левой части этого равенства, кроме члена, для которого индекс суммирования ѵ равен р, равны нулю вследствие формулы (14.6.2). Для члена ряда, соответствующего ѵ = [х, интеграл на основании формулы (14.6.2) равен единице. Следова­ тельно, равенство (14.6.7) дает

t

f an(a)f(t, ° )do•

(14.6.8)

t-т

 

Отсюда находим

 

t

 

£ц = ту-

j аи (а)/ (** °)da-

(14.6.9)

**

t - T

 

Таким образом, неизвестные коэффициенты ряда (14.6.4)

опреде­

ляются весьма просто.

(14.6.9) коэффициентов £ѵ

(сменив

Подставляя выражение

в нем, конечно, индекс р на ѵ) в (14.6.4), получим следующую окончательную формулу для решения интегрального уравнения

(14.1.1) :

 

_____

t

 

g(t,

=

j av (<?)/(*, О) da.

(14.6.10)

V

V t - T


§ 14.6. М ЕТОД К А Н О Н И Ч Е С К И Х Р А ЗЛ О Ж Е Н И Й

585

Таким образом, метод канонических разложений случайных функций дает возможность весьма просто найти решение основного интегрального уравнения статистической теории оптимальных сис­ тем (14.1.1) в форме бесконечного ряда (14.6.10) Этот ряд определя­ ет решение интегрального уравнения (14.1.1) во всех случаях, когда его правая часть / (t, а) может быть представлена разложе-

нпем (14.6.6) по сопряженным координатным функциям х ѵ (о) слу­ чайной функции X (<) в интервале наблюдения t Т ^ а ^ L В случае, когда функция / (t, а) не может быть представлена раз­ ложением (14.6.6), интегральное уравнение (14.1.1) не имеет реше­

ния ([53], §§ 134, 135).

В предыдущих выкладках mj>i выполнили почленное интегри­ рование ряда в формуле (14.6.5) Можно доказать, что эта операция законна и что конечным отрезком ряда (14.6.10) всегда можно пользоваться для приближенного определения весовой функции оптимальной линейной системы ([53], § 135). При этом, несмотря на то, что найденная таким образом весовая функция может доволь­ но значительно отличаться от весовой функции оптимальной систе­ мы, соответствующая средняя квадратическая ошибка, как прави­ ло, лишь незначительно превосходит минимальную среднюю квад­ ратическую ошибку оптимальной системы, если взять 10 -f- 20 членов ряда (14.6.10).

Изложенный метод решения интегрального урвнения (14.1.1)

вформе бесконечного ряда гораздо проще методов, изложенных

впредыдущих параграфах. Одпако для применения этого метода

необходимо предварительно

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

аѵ (т), х ѵ (т), определяющих

каноническое разложение суммы

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

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

средней квадратической ошибки.

последовательно

Полагая в формуле

(14.6.10)

/ («, о) =

К ѵх (<, о),

ф! (о), . . ., <pN (а),


5 8 6

ГЛ . 14. О П Р Е Д Е Л Е Н И Е О П ТИ М А Л ЬН Ы Х Л И Н Е Й Н Ы Х СИСТЕМ

найдем решения всех интегральных уравнений (13.6.1) и (13.6.2):

 

 

 

 

 

 

 

____

і

 

 

 

 

 

g«o>(t, т) =

2

^

-

j <h(o)Kyx(t, a)da,

(14.6.11)

gir) (t, t) =

__

 

t

V

V

t - T

(r =l ,

N).

 

2

 

 

^

av (а) фг (o) da

(14.6.12)

 

 

 

 

t - T

 

 

 

 

 

 

 

Полагая

для краткости

 

 

 

 

 

 

 

V

V

 

 

 

 

 

 

 

 

 

Уѵій=- щ

jt

 

av(a)Kyx(t, a)da,

 

 

(14.6.13)

 

 

v

t - T

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

в-vT (t) = -Щ- £

 

av (a) cpr (a) da

(r= 1, . .. ,

N),

(14.6.14)

 

 

v

t - T

 

 

 

 

 

 

 

 

можем представить формулы (14.6.11) и (14.6.12) в виде

 

 

 

gw (t,

х) =

2

2/v (О «V (х),

 

 

(14.6.15)

 

 

g <r) (i,

X ) =

SV

«vr (0

аѵ(т)

(г = 1, . . . ,

N).

(14.6.16)

V

Подставляя эти выражения в формулы (13.6.6) и выполняя почлен­

ное интегрирование рядов,

получим

 

 

 

 

t

 

 

 

 

 

 

*ро=

J

g(0) (t, т)фр (т)гіт =

 

 

 

 

 

t - T

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

=

\ flv(t)cpv {x)dx

(p= 1, . . . ,

N),

(14.6.17)

 

t

V

t - T

 

 

 

 

bpq=

j

g(<l)(t, Т)фр(т)с2т =

 

 

 

 

 

t - T

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

= ^ j avq(t) f av (x)(pp(x)dx

(p, q= 1, ..

N).

(14.6.18)

 

 

V

t - T

 

 

 

 

Но на основании формулы (14.6.14)

 

 

 

 

 

t

 

1

 

 

 

 

 

j

av (г) фр (х) dx=

^ av (er) ФР (a) da = ZVzVp (0 •

 

t - T

t - T