книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем
.pdf190 |
Т Е П Л О В Ы Е Ф Л У К Т У А Ц И О Н Н Ы Е |
К О Л Е БА Н И Я |
[ГЛ . V |
определяются уравнениями (2.31), (2.34): |
|
||
|
А — А а — а?А — A SA |
= 0, |
(5.9) |
|
М + аМ + Мат — 5 |
= 0. |
(5.10) |
Введем следующие обозначения для квадратных матриц
моментов (ковариационных матриц) размерности п |
х п: |
||
= |М1М*1|. |
п |
|
|
|
muiiqic] j = |
|
|
м ря= м\р= |м [ 2 |
|
||
|
/=1 |
п |
(5.11) |
|
|
||
|
|
- 1 2 ПЧ|М lqiqK) |
|
п |
п |
1=1 |
|
mk,5sj | = mM qqm, |
|
||
р р — м [21 яц,д, 2 |
|
||
1=1 |
s= l |
|
|
где |
|
|
|
Матрицу M представим в виде блочной матрицы:
|
|
М - |
М |
Р Р |
М |
Р Я |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
М |
я |
р |
М |
я я |
|
|
|
Уравнение (5.10) можно записать в форме |
|
|
||||||||
йрр м рд + |
г т |
М |
Р Р |
М |
Р Я |
|
+ |
|
|
|
|
с |
|
|
|
|
|
|
|
|
|
Й яя |
— тГ1 0 |
|
|
|
|
|
|
|
|
|
|
М |
Я Р |
М |
Я Я |
|
|
|
|
||
|
|
+ |
М |
Р Р |
|
М |
Р Я |
т ~ 1 г т - |
т ~ х |
0 |
|
|
|
с |
0 |
0 0 |
|||||
|
|
|
М |
Я Р |
|
М |
Я Я |
|||
|
|
|
|
|
|
|
||||
или в раскрытой до матриц-блоков форме: |
|
|
||||||||
М рр + гт~хМрр + |
сМ9р + МррПГ1^ |
+ МрдС = |
S v, |
|
||||||
M vq + гтГШрд + |
cMqq — М ррШ 1 = |
0, |
i' |
(Ь.12) |
||||||
# |
|
|
|
|
|
|
|
|
||
М д р — П Г Ш р р + М д р Г Г Г Ь Т + М д д С = 0 ,
Mqq — тГШрд — МдрПГ1 = 0-
Для случая системы, находящейся в тепловом равно весии, когда справедлива формула (5.4), эти уравне-
§ 5.1] |
Ф Л У К Т У А Ц И О Н Н Ы Е Т Е П Л О В Ы Е К О Л Е Б А Н И Я |
191 |
ния в установившемся режиме (МРР = M pq = M q„ = 0)
имеют исключительно простое решение [5.5], [5.6]:
Мд, = кТ с-1, Мрр = АТтге, МР(, = |
= 0. (5.13) |
При этом
= кТт ~\ |
(5.14) |
Средние кинетическая и потенциальная энергии флукт} а ционных колебаний такой системы равны
4 - 2 |
лЧкМ[д4?к] = -L nkT , |
n |
(5.15) |
4 " s |
с« м [ м * ] = 4 _ п А г » |
что находится в полном согласии с положениями кинети ческой теории. Согласно (5.13), (5.14) матрицы коэффи циентов распределения равновесной системы равны
Aqq — |
■М~* — |
^ я я |
1 |
пг. |
|
ш яя ~ |
кТ С’ А яя |
кТ |
|
||
Поэтому плотность равновесного распределения может
быть представлена |
в виде |
|
|
Р = d exp [4" |
+ |
4" $*Am 9) = |
|
|
= |
d exp |
(gTcq + 9Tmg)] . (5.16) |
Величина |
|
|
|
4 (ЯТсЯ + Ятт )
есть энергия системы (флуктуационных колебаний), и распределение (5.16) есть распределение Гиббса [1.16].
В ряде случаев в динамической системе нет потенциаль ных сил. В этих случаях с = 0 и уравнения (5.12) прини мают форму
1Мрр -|- ттп М рртп V =
lifpq + rm~xM vq — МррПГ1 = 0,
(5.17)
М д р + М д р Т П ~ V — Г Г Г Ш р р = 0 ,
lifgg — Щ-Шрд — МдрПГ1 = 0.
192 Т Е П Л О В Ы Е Ф Л У К Т У А Н И О Н Н Ы Е К О Л Е БА Н И Я [ГЛ. V
Здесь нет установившегося |
состояния |
в смысле M qq = |
|||
= const, |
однако есть |
установившееся |
состояние, при |
||
К О Т О р О М |
М рр |
C O n b t , |
М p q, |
М qp === C O n S t , М qq === C O H S t, |
|
т. е. вторые моменты обобщенных координат линейно из меняются во времени. Соответствующее решение для си стемы, находящейся в термодинамическом равновесии (справедлива формула (5.4)), имеет вид
М рр = |
кТгп, |
Мрд — кТтпг'1, |
M qр = |
кТ (гт)-1 пг, M qq = кТ [г-1 + (г1)"1]. |
|
В этом можно убедиться прямой подстановкой выражений
(5.4), (5.18) в (5.17).
Для неравновесной системы, температуры различных составных частей которой различны, а также для пере ходных режимов формулы для коэффициентов распределе ния сложнее указанных. Однако и для этих случаев мож но записать общие выражения для А — —М~х и получить
некоторые общие выводы.
Матрицу фундаментальной системы решений, удовлет воряющую уравнению
w + aw = 0, w (tr, t') = 1
(1 — единичная матрица), представим в блочной форме:
w(t, t') =
«V |
(‘. *') w gg (*'• |
<) |
где |
|
|
U > p p ( t ' , t ' ) = 1, W pg { t ' , t ' ) = 0, |
||
W q p (<',<') =0, |
W q q ( t ' , t ’) = |
1. |
В последних выражениях через 1 обозначены единичные матрицы размерности п X п.
Уравнение w + |
aw — 0 |
в |
раскрытой до |
матриц-бло |
||||
ков форме имеет |
вид |
|
|
|
|
|
|
|
Wpp + |
rmr'wpp + |
cwgp = |
0, |
wpq + r in xWpg + |
cwqq = |
0, |
||
|
W g p — m~xWpp = |
0 , |
|
wqq — |
m~xwvq = |
0 . |
||
Отсюда |
вытекает |
|
|
|
|
|
|
|
m ib gp + rWqp + CWgp = |
0 , |
w qv ( f , t') = |
0 , |
|
|
|||
w qv ( t \ t') = m~\ |
Wpp = mWgp. |
8 5.1] Ф Л У К Т У А Ц И О Н Н Ы Е Т Е П Л О В Ы Е К О Л Е Б А Н И Я 193
Решение уравнения (5.10) запишем в виде (2.35)i
I
М (t) = w (t, 0) М (0) ш' (/, 0) + 1 w (t, О Sw (t, 1') dt'. (5.20)
О
Необходимо различать понятия термодинамического рав новесия и равновесного или установившегося состояния в смысле уровней шумов. Если элементы системы имеют различные температуры, то термодинамическое состоя ние неравновесное. Для поддержания такого состояния (поддержания различных температур) необходимы, вооб
ще говоря, |
подвод |
и |
отвод |
тепла. |
|
|
|
Допустим, что за счет подвода и отвода тепла поддер |
|||||||
живается |
неизменное |
неравновесное термодинамическое |
|||||
состояние, |
т. е. различные |
значения |
Тih = |
const. Па |
|||
раметры |
системы |
считаем |
постоянными: |
m = const, |
|||
с = const, |
г |
— const. |
Такую систему |
будем называть |
|||
неравновесной (в термодинамическом смысле) стационар ной. Весовые функции стационарной системы суть функ
ции разности аргументов:
w (t, t') = w (t — t').
Если невозмущенное движение системы асимптотиче ски устойчиво, т. е. ш (£) -> 0 при t -*■ оо, то как в равно
весной, так и неравновесной (в термодинамическом смы сле) стационарных системах с течением времени устано вятся определенные уровни флуктуационных тепловых колебаний:
|
( |
М (оо) = lim М (<) = |
lim § w(t — t') Swт (t — f) dt' = |
t— OQ |
t—+OQg |
|
OO |
|
= Jw (t) Swt(t)dx, |
|
0 |
При блочном представлении матриц имеем
"ррЛ , ря |
W V V (т ) ШРЯ(т> |
Л1яг Л1яя |
*9Я (Т) |
S a |
0 |
w l v |
№ |
w qT p (X ) |
|
|
|||
0 |
0 |
w p q |
( V |
U g q ( Т ) |
|
|
т/27 А. А. Красовский
1 9 4 ТЕПЛОВЫЕ ФЛУКТУАЦИОННЫЕ КОЛЕБАНИЯ 1ГЛ. V
Отсюда |
следует: |
оо |
|
|
|
Мрр = |
J |
Щ р (t) SvwlP(т) dx, |
(5.21) |
|
|
0 |
|
|
|
|
00 |
|
|
|
М qq — |
J* м ? д р ( т ) S^W qp ( т ) dx. |
(5.22) |
|
|
|
0 |
|
|
Первая |
из этих формул |
с учетом выражений |
|
|
|
М рр = п г М - m , ц>рР = т й ? д р |
|
||
дает |
|
00 |
|
|
|
|
|
|
|
|
М qq = |
J ^< 7P ( ^ ) Sqi0qp(x) dx. |
(5.23) |
|
|
|
о |
|
|
Обозначая элементы матрицы wqp (т) просто через wth (т),
на основании (5.22), (5.23), (5.5) записываем:
п |
оо |
М [<мЛ = к 2 |
(т•iii r + Тн-v^v) | Щу (Т) юлi (т) dx, (5.24) |
V, Ц=1 |
О |
П00
М [м ,] = А |
2 ( ^ Г У(Л+ |
r ^ v ) |
(r)wJlx(x)dx, |
(5.25) |
|
|
V,H-=1 |
|
о |
|
|
где wij (2) |
удовлетворяет системе |
уравнений |
|
||
|
п |
п |
|
п |
|
|
2 |
+ 2 |
+ 2 Civ“’vj =0 |
(5.26) |
|
|
v = l |
V = J |
V = 1 |
|
|
при начальных |
условиях |
|
|
|
|
|
|
*.i (0) = |
, Wy} (0) = 0, |
(5.27) |
|
| m | — определитель матрицы m, тп)ч — алгебраическое дополнение элемента mj4 этого определителя.
Широко известны [5.8] формулы для вычисления ин тегральных квадратичных оценок вида
00 |
оо |
|
J u>iv (Т) Wp. (т) dx, |
J Щу (т) щX (т) dx. |
|
и |
|
о |
Если преобразования Лапласа функций wty (t), Wj (t), яв
ляющиеся дробно-рациональными функциями параметра
t 5.1] |
Ф Л У К Т У А Н И О Н Н Ы Е Т Е П Л О В Ы Е К О Л Е Б А Н И Я |
195 |
||||||
преобразования s, представлены в виде |
|
|
||||||
|
|
|
bf + |
bft + .- . + bn*71-1 |
||||
W h (s) = ^ wu it) exp (— st) dt = |
oo + |
ais + |
. .. + |
» |
||||
|
|
|
|
&ns |
||||
|
Wu (t) Wfr (t) dt = |
n |
2 |
, |
AeB?w, |
(5.28) |
||
|
|
|
_ |
|
|
|
||
|
00 |
— 02 |
04 • |
• |
• |
|
|
|
|
0 |
01 |
-*• Оз • |
• |
• |
|
|
|
|
0 |
*— Оо |
02 • |
• |
• |
|
|
|
|
|
• |
• . |
* . |
|
|
|
|
|
|
* |
• • |
* ап-\ |
|
|
||
Дq — алгебраическое |
дополнение |
элемента |
n-й |
строки, |
||||
<7~го столбца определителя Д; |
|
|
|
|
|
|||
|
Blw = b?b¥, В ^ } = K b f - |
Ь?ЪУ - bfb |
|
|||||
|
В?* = ь № - |
Ъ?Ъ? - |
|
+ |
b5%w' + |
Ь ^ , |
||
в Т 1= №
Обычно основной интерес представляют дисперсии обоб щенных координат и скоростей
М [Qi] — к |
(T ’V[o.rvjji -f- Т |л^Гр.7) j" u>;v (t)iDiy.(t)dt, |
|
|
v.H=1 |
о |
(5.29) |
|
n |
oo |
||
|
|||
M [0i] = к 2 |
(Tv|j,r4\i -f- Tpvr\iv) J iou (t) Wip (<) dt. |
|
|
v,H-=l |
|
|
Формулы (5.29), даже без использования выражений для интегральных квадратичных оценок (5.28), позволяют сделать следующий общий вывод. Допустим, что темпера туры различных элементов системы находятся в границах
^raln |
Tib <С ^maxi |
i, к = 1 , 2 , . . . , Т1‘, |
тогда уровень тепловых шумов в любой точке данной не равновесной системы больше уровня шумов равновесной системы с температурой Гт щ и меньше уровня шумов равновесной системы с температурой Гщах.
196 |
Т Е П Л О В Ы Е Ф Л У К Т У А Ц И О Н Н Ы Е КОЛЕБАН ИЯ |
[ГЛ. V |
Для доказательства обратим внимание на то, что в пас сивных системах интегралы
оо |
оо |
|
j u>iv (t) wiV. (t) dt, Jiolv(t) wiXx(t) dt |
(5.30) |
|
о |
о |
|
всегда положительны. |
Действительно, всегда |
можно за |
дать такую воображаемую систему, в которой все вели
чины |
TVvjx 4- Tav^v, |
кроме |
одной, |
равны нулю |
или |
за счет равенства нулю активных |
сопротивлений |
rvlj., |
|||
или за |
счет равенства |
нулю |
абсолютных температур |
|
|
(третье начало термодинамики в подобном воображаемом эксперименте можно во внимание не принимать).
Для подобной |
системы согласно |
(5.29) |
||
|
к(71v^rvji |
во |
|
|
М [<??] = |
T^Av^V'') J |
(t) Wty. (0 |
||
|
|
|
О |
|
|
|
|
оо |
|
М [g?] = |
/с |
|
+ T ^ r^ ) j |
wu (t) wiv.(() dt. |
|
|
|
0 |
|
Величины M [gf], |
М [gf], |
|
по своей природе |
|
и условиям положительны, поэтому положительны и ин тегралы (5.30).
Таким образом, дисперсии (5.29) в любой точке нерав новесной системы представляют собой в отношении тем ператур элементов 7\ц линейные формы с положитель ными коэффициентами. Такая форма принимает наимень шее значение на нижней границе аргументов (темпера тур) и наибольшее значение на верхней границе аргумен тов. Положение доказано.
Для равновесных систем справедливы весьма простые формулы (5.13), (5.14) вторых моментов флуктуационных колебаний. Эти формулы в сочетании с доказанным поло жением позволяют для неравновесной системы занисать следующие неравенства:
кТпн,, < М [gf,] < кТшах
(5.31)
max
6 5.1] |
Ф Л У К Т У А Ц И О Н Н Ы Е Т Е П Л О В Ы Е К О Л Е БА Н И Я |
197 |
где |
| т | — определитель матрицы пг, ?nli — алгебраи |
|
ческое дополнение диагонального элемента тц этого
определителя, | с | — определитель матрицы с, с*1— алге браическое дополнение диагонального элемента Сц этого
определителя.
Перейдем к рассмотрению конкретных примеров. Тепловые шумы в много
обмоточном |
трансформаторе. |
Rr rn |
т г22 |
Несколько |
обмоток с коэф |
|
|
фициентами |
индуктивности |
|
|
L ti (г = 1 , 2 , . . . , п) и взаи-
моиндуктивности Мтис (U
=1, 2, . . ., п) подсоединены
кактивным сопротивлениям
B t (г = 1 , 2 , . . ., п) с темпе-
ратурамиГ|= const (рис. 5.1),
находящимися в пределах
T’mln < T i |
T m a x • |
Ферромагнитный сердечник отсутствует, так что источ никами шумов служат толь ко активные сопротивления.
I |
-______ I |
1 J |
G |
^ П - / ~ ГП - 1 . п - I |
Я 'п ~ гпп |
Рис. 5.1.
Уравнения цепей обмоток
+ Rii + |
Мгц - ^ п |
+ |
Мтц, |
+ |
• • • + |
М тт -£■ = ф ;, |
|||
|
|
|
i — 1, 2, . . ., л, |
|
|
||||
можно записать |
в виде (5.2): |
|
|
|
|||||
где |
|
|
тд + |
|
rq = <р, |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
L Tl |
|
л / Т12 . |
. . |
|
|
h |
||
т = |
м тг1 |
L T‘Z2 |
* |
' |
" |
М Тчп |
. Я = |
h |
|
• |
|||||||||
|
М Тп\ |
М Тп2 • • |
• |
L Тп |
|
in U |
|||
|
Я I |
0 |
. . . |
0 |
|
|
|
|
|
г = |
О /?а |
. . . |
б |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
о |
о |
. . . д_ |
|
|
|
|
||
198 |
Т Е П Л О В Ы Е Ф Л У К Т У А Ц И О Н Н Ы Е КОЛЕБА Н И Я |
[ГЛ. V |
В отношении токов (матрица q) данная система устойчи
ва, и в ней устанавливается равновесное (в смысле шумов) состояние. Согласно (5.29) для этого состояния
(5.32)
а матрица w (t) — || юц (t) || удовлетворяет уравнению
miv + rw = 0, w (0) = тг1.
Дисперсии токов согласно (5.31) находятся в пределах
кТш1п^ < М [ijj < kTmaх^ , |
(5.33) |
|
L T l |
М |
т п |
|
|
||
т п | = |
|
L |
T 2 2 |
. |
. |
. |
^ T ln |
• |
• |
• |
М Г2П |
|
|
М Т п 1 М т п 2 • • • L T n |
|
|
|||
— алгебраические |
дополнения |
диагональных членов |
|||||
этого |
определителя. |
|
|
|
|
|
|
Для системы с одинаковыми температурами сопротив |
|||||||
лений |
|
|
|
„н |
|
|
|
|
|
|
|
|
|
||
|
|
М [ф = кТ I m | |
|
|
|||
В частности, при п = |
3 |
выражения в развернутой форме |
|||||
имеют |
вид |
|
|
|
|
|
|
|
М [ф |
кТ___________ 1 |
~ ^гз_________ |
|
|||
|
|
1 + |
2(112(113(123 - |
( ф - й-13 - |
(1^3 |
’ |
|
|
|
|
|||||
|
М [ф |
кТ___________ 1 |
—И^хз________ |
|
|||
|
L T% |
1 + |
2(112(113(123 - |
( Ф - ( ф - |
( ф |
’ |
|
|
|
||||||
|
М [ф |
кТ___________ 1 |
— t*fa________ |
|
|||
|
L T3 |
1 + |
2(112(113(123 — |
( l * 2 — (1*3 — |
(Аая |
|
|
где |
|
|
|||||
м Т12 |
|
|
М Т13 |
|
|
|
|
|
|
|
|
|
М Т23 |
||
Pl2 — V |
^13 — У L TlL T3 |
У LT%LTa |
|||||
|
|
|
|
|
Н-23 ; |
|
|
— коэффициенты связи.
§ 5.1] Ф Л У К Т У А Н И О Н Н Ы Е Т Е П Л О В Ы Е К О Л Е Б А Н И Я 199
При наличии лишь двух |
обмоток |
|
|
|
|
Р13 — О, P '2 3 — О , |
|
|
|
кТ |
i |
|
кТ |
i |
М[«!] = иТ1 |
1-иЬ |
МГ£1] = |
L T2 |
i - p i . ’ |
При стремлении коэффициента связи в последних форму лах к единице, что соответствует исчезновению потоков рассеяния, дисперсии шумовых токов неограниченно воз растают. Так и должно быть, ибо при отсутствии потоков рассеяния эквивалентная схема двух связанных обмоток вырождается в цепь с чисто активными сопротивлениями, дисперсии шумовых токов в которой по формуле Найкви ста (без квантового множителя) бесконечно велики.
Рассмотрим теплообмен за счет тепловых шумов. Допустим, что обмотки и сопротивления имеют абсолют ную тепловую изоляцию, так что обмен тепловой энергией может происходить только за счет шумовых токов. Оче видно, что формулами (5.29), а значит, и (5.32) можно пользоваться не только при постоянных температурах сопротивлений, но и при медленно (в сравнении с элек трическими переходными процессами) изменяющихся тем пературах.
Энергия, выделяемая в единицу времени в сопротив лении R t за счет шумовых токов, создаваемых всеми дру
гими сопротивлениями, согласно (5.32) равна
2/сД4 [ S |
З Д J wfj (0 dt - |
TiR, j ivl (0 dt] . |
3 = 1 |
о |
о |
Энергия в единицу времени, расходуемая сопротивлением на создание шумовых токов в остальных элементах схе мы, равна
П ОО |
00 |
2kRiTl [ 2 R\ j w)i (t) dt — R i j w\ (t) d t\ .
3=1 |
о |
0 |
Разность этих энергий равна произведению теплоем кости ci i-ro сопротивления (точнее, г-й теплоизолирован ной цепи) на производную температуры T ti
ПСО
CiTi = 2kRt [ 2 (Tj - |
ТО Rj j К |
(t) dt] . |
(5.34) |
i=i |
о |
J |
|
