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

книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем

.pdf
Скачиваний:
35
Добавлен:
19.10.2023
Размер:
7.02 Mб
Скачать

190

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

К О Л Е БА Н И Я

[ГЛ . 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

 

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