Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.pdf

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

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

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

Добавлен: 23.10.2024

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

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

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

214 Т Е П Л О В Ы Е Ф Л У К Т У А Ц И О Н Н Ы Е К О Л Е БА Н И Я [ГЛ. У

В развернутой до блоков форме это уравнение дает

Д •• + m~lrRil +

т

+

Ицг^тГ1+ Ritcm 1

+

+ - щ \

(г +

гТ)

=

кТтГ1 (г + гТ)

 

 

Rit Hr mT1rRu +

m-'cR" -

 

-f

 

 

 

 

 

 

2 k T

И

R 0 t

(Г “Ь

й « —

(5.62)

— йёё + ЯгкгТт~г + Л.гСт'1 +

 

 

 

 

 

 

+ T S T - R,i

+

'*)«,, =

 

О,

 

 

 

И

 

 

 

 

 

л.. - й,ч - л.;+ -д4Д,; (г + D я и = о.

И

Эти уравнения для стационарной системы имеют весьма простое частное решение

R„ = kTmc^,

- кТтт~\ Я£. = Д .е = 0. (5.63)

Действительно, подставляя (5.63) в (5.62), убеждаемся, что три нижних уравнения обращаются в тождества, а верхнее уравнение приобретает вид

кТшттГ1гтп~1 4- кТшт “1/-тт “1 +

кт2

 

 

-тр,!п- тп"1(г + гт) тпГ1 +

 

откуда следует:

 

 

 

- кТггГ1 +

гт) ш-1,

гр2

 

 

 

m ,

 

m

 

1 тп

__

(5.64)

•Lш т

2

'

^ ’

Тш = Тя ( у

 

1 + 2 - £ — l)

(5.65)

(второй корень уравнения (5.64) отрицателен и не имеет физического смысла). Рассматриваемая система вполне наблюдаема по Калману. Действительно, легко прове­ рить, что при неособых матрицах тп, с, г + гт ранг ма­

трицы

||h\ аЧтг.....(ат)п~1Л^||

(5.66)


§ 5.2]

М И К РО У П РА В Л Е Н И Е

215

 

 

равен 2п. Из теории фильтра Калмана известно, что ре­

шение (5.63) в этом случае является единственным в том отношении, что к нему при t —>- оо стремятся все другие

решения, соответствующие произвольным начальным ус­ ловиям.

Итак, в установившемся состоянии корреляционные матрицы ошибок оптимального оценивания равны

Я гг= кТ и ( ^ 1 + 2

^ - - l ) c - \

 

(5.67)

Пц = кТи ( ] / i + 2

Y - — 1) m - \

Эти выражения следует рассматривать как характери­ стики предельной возможной точности контроля пассив­ ной системы при рассматриваемом неквантовом взаимодей­ ствии с измерителями. Они допускают простую интерпре­ тацию. Рассмотрим пассивную систему, аналогичную контролируемой, но с абсолютной температурой.

ТШ= ТИ ( j / l + 2 ■ £ - 1 )*

<5'68)

Согласно (5.13), (5.14) ковариационные матрицы тепловых флуктуационных колебаний в такой системе равны

M qq = кТшс~\ M qq = кТшПГ1,

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

образно рассмотреть три случая:

а) Та <^Т,

ТШ^ У 2 Т ^ ] \

 

б)

Ти = Т,

Тш = 0,732Т;

(5.69)

в)

7,и> 7 \

J

 

Наибольший интерес

представляют случаи

б) и в).

В случае б), когда шумовая температура измерителей равна температуре Т контролируемой системы, предельная


211)

Т Е П Л О В Ы Е Ф Л У К Т У А Ц Н О Н Н Ы Е К О Л Е Б А Н И Я

[ГЛ. V

достижимая точность контроля соответствует тепловым флуктуационным колебаниям при температуре 0,732г. В случае Та ^ Т предельная достижимая точность кон­ троля практически не зависит от Та и соответствует уров­

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

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

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

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

4 }

на основе (5.54), (5.58), (5.63) уравнение оптимального фильтра представляем в виде

пЩф + г ( ф + c q ф +

( г + г т) £ф =

Z. ( 5 . 7 0 )

Таким образом, стационарный оптимальный фильтр пред­ ставляет собой некоторую модель контролируемой сис­ темы, отличающуюся от последней, по существу, лишь диссипативными членами. Эти члены слабо влияют на собственные частоты. Кроме того, при Ти Т, когда

согласно предыдущему

Т Ш

< 1.

гги

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

§ 5.21

М И К Р О У П Р А В Л К HUE

217

вышать собственные частоты разомкнутой системы. Более строгое доказательство этого положения о предельной точности управления можно получить на основе теоремы разделения или стохастической эквивалентности [5.15]. Указанный порог, или барьер, управления становится особенно заметным в задачах контроля микроскопических объектов. Есть, однако, основания предполагать, что с помощью квантового взаимодействия с контролируемым объектом можно получить существенно более высокую точность контроля в сравнении с рассмотренным «класси­ ческим» случаем. К такому выводу приводят, в частности, следующие рассуждения. Пусть мы имеем систему (5.2), но с одной степенью свободы. Тогда т, г, с, q, <р — ска­

лярные величины и формулы (5.67) дают выражения для дисперсий ошибок оценивания координаты и ее произ­ водной:

к Т

R - =

кТ„

„ г 2 , кТта

М [ё21

# е е М [б2] = —^

9 9.

L

Величина

т 2М [ё2] = М [(тгеё)2] = кТшт

есть дисперсия ошибки оценивания импульса. Таким образом,

М [е2] М [(те)2] = е2 (те)2 = к^Тт — ,

С

v ? v m = k Tm

где соо = " j / " — собственная частота контролируемой

системы. Между тем соотношение неопределенности Гей­ зенберга [5.16], рассматривающее квантовое взаимодей­ ствие с контролируемым процессом, может быть записано в виде

V & V (те)2 ж ,

где h — постоянная Планка.

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


218 Т Е П Л О В Ы Е Ф Л У К Т У А Ц И О Н Н Ы Е К О Л Е Б А Н И Я [ГЛ. V

(1.2) без квантового множителя), имеет место сильное *) неравенство

kT m > h ^ - .

(5.71)

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

Предельная точность контроля координат посредством квантового излучения. Рассмотрим случай оптического контроля координат механической или электромеханиче­ ской системы посредством света с длиной волны, достаточ­ но короткой для того, чтобы его можно было считать по­ током фотонов. Изменение компонент вектора q вызывает

изменение световых потоков, улавливаемых фотоэлектрон­ ными умножителями (ФЭУ) или эквивалентными им по чувствительности фотодатчиками. Известно [5.13], [5.14], что чувствительность ФЭУ к слабым световым потокам может быть сделана настолько высокой, что с большой ве­ роятностью регистрируется каждый квант (фотон) излу­ чения. Световые потоки будем считать слабыми, но на­ столько, что число и* фотонов, поступающих в приемник

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

каждым

датчиком,

изменяется на полную величину

Еч — h\riv

(v — частота) при изменении

соответствую­

щей координаты на 6.

Величина б, конечно,

больше длины

волны К используемого света. Таким образом, крутизна

характеристики каждого фотодатчика считается

равной

Ц

фотонов

hv . Вследствие «дробового эффекта» потока

энергия, улавливаемая фотодатчиком за единицу време­ ни, флуктуирует. Согласно известному правилу дисперсия

*) При Т — 500 град, too = 2я 1000,

2я = 1,66-Ю9.


§ 5.2]

М И К Р О У П Р А В Л Е Н И Е

219

 

 

 

 

этих

флуктуаций выражается формулой

 

 

 

АЕ* = hvE4.

(5.72)

В

диапазоне частот

сигнала, существенно

меньших

nv

=

Ejhv, флуктуации

светового потока можно счи­

тать белыми шумами со спектральной плотностью (5.72). Однако нас интересует шум, приведенный к контролируе­ мой координате. Спектральная плотность приведенного шума определяется как величина (5.72), поделенная на квадрат крутизны преобразования:

Шум, свойственный самому фотоприемнику (ФЭУ), как уже отмечалось, может быть исчезающе малым, и им пре­ небрегаем. Световые потоки фотодатчиков считаем неза­ висимыми. Таким образом, вектор наблюдения в данном случае имеет вид

z = Я + £и,

(5.74)

где | и — векторный белый шум с матрицей спектральных плотностей

5к = -?-1,

(5.75)

1 — единичная матрица размерности п X п.

Флуктуа­

ции светового давления создают некоторое обратное шумо­ вое воздействие на контролируемый объект. Спектральная плотность флуктуаций светового давления имеет следую­ щий порядок величины:

 

Д Я 2

 

(5.76)

 

с2

 

 

СВ

 

 

где ссв — скорость

света.

Уравнения контролируемого

процесса запишутся в виде

 

 

mq

+ rq +

cq = ф + срд,

(5.77)

где ср — тепловые шумы с матрицей спектральных плотцоетрй (5.4), фд — шумы, создаваемые флуктуациями