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

Метода термодинамика(практика)

.pdf
Скачиваний:
256
Добавлен:
23.02.2015
Размер:
1.14 Mб
Скачать

20

 

 

 

 

 

N 1

 

(N 1)!

 

 

 

 

 

 

N 1

 

 

 

 

 

 

 

 

 

= Np

 

 

 

pν (1 p)( N 1) −ν = NpWN 1(ν) = Np .

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ν=0 [(N 1) − ν]! ν!

 

 

 

 

 

 

ν=0

 

 

 

 

 

1. Для случая N >> n и p <<1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (n) = N(N 1) ... (N n +1)

pn (1p)N n

n! =

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ( pN)n (1 1 N ) ... (1(n 1) N) (1p)N

n! ( pN )n (1p)N n!.

Поскольку

 

Np =

 

 

 

 

постоянная величина,

в

том числе, при

условии

 

n

 

lim Np =

 

, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N → ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (n) = ( pN)n (1 pN N )N

n!.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует распределение Пуассона:

 

 

 

 

 

 

 

 

 

 

 

 

lim WN (n) =W (n) = ( pN )n epN n! = (n )n e

 

 

n!.

 

 

 

 

 

 

n

 

 

 

 

 

 

N → ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Условие n >>1 и

N >>1

означает, что

 

и n можно считать большими, а

n

разность n =

 

n <<

 

 

.

Прологарифмируем

распределение

Пуассона,

n

n

используя приближенную формулу Стирлинга ln n! nln n n :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lnW (n) = nln n n [nln n n]= −(n − ∆n )ln (1− ∆n) n − ∆n

 

 

 

 

{

 

 

 

 

 

 

 

 

)

(

 

 

 

 

 

 

 

 

)

 

}

 

 

n

(

(1− ∆n)

n

−∆n n

(n)2 2n2 ...

+ ∆n n

 

 

 

= −(n)2

2n .

 

 

 

 

 

Рассмотренный случай с учетом условия нормировки и приводит к распределению Гаусса:

 

 

 

 

 

 

 

(n

 

)2

 

W (n) =

1

 

 

exp

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2π(n)

 

 

 

 

2(n)

 

 

 

 

 

 

 

 

 

 

 

 

П р и м е р 3.

Распределение случайной величины x

задано

функцией

f (x) = Aexp((x 1)2 ). Найти: среднее значение

 

и

вероятность,

что

x

x xo =1.1.

 

 

 

 

 

 

Р е ш е н и е :

После определения коэффициента

нормировки

A =1

π

находим среднее значение:

21

x =(1 π) [(x 1) +1]e( x1)2 dx =1.

2.После замены переменной y = x 1 выражение для искомой вероятности примет вид

 

1

y0

y2

 

1

 

0

 

y2

y0

y2

 

w(x x0 ) w( y y0 ) =

 

e

 

dy =

 

 

e

 

dy + e

 

dy .

 

 

 

 

 

 

π − ∞

 

 

π

− ∞

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первое слагаемое в скобках сводится к интегралу Пуассона и равно 0.5 π . Второе – интеграл ошибок. Но поскольку y0 = 0.1 <1, то можно подынтегральную функцию разложить в ряд, а затем почленно проинтегрировать

y0

1

 

y0

ey2 dy =

 

y2ndy =

n!

0

n=0

0

n=0

y2n+1

0 .

(2n +1) n!

Вычисления дают w(x x0 ) 0.5 + 0.1 π.

П р и м е р 4. Плотность вероятности случайной величины задана функцией Гаусса (2πσ2 )1 exp (x x0 )2 2σ2 , где x0 – ее среднее значение, а

σ– среднее квадратичное отклонение. Каково должно быть значениеσ, чтобы вероятность попадания ее в интервал xx x+ ( x = x0 + D , D > ∆ > 0 ) достигала максимума?

Р е ш е н и е . Искомая вероятность при замене переменной t = (x x0 )σ описывается формулой:

 

 

) = 1

x +

 

 

2σ2 )dx = 1

( D+∆) σ

W (x

x x

 

exp((x x

)2

exp(t2 2)dt .

+

2πσ2

x

0

 

2π ( D−∆) σ

 

 

 

 

Заметим,

что

lim W = 0 .

 

Чтобы найти

экстремум

W (σ) , воспользуемся

σ→0 σ→ ∞

условием W ∂σ = 0:

(D + ∆) e( D+∆)2 2σ2 = (D − ∆) e( D−∆)2 2σ2 .

22

Отсюда следует σ

(max)

=

2D ln ((D + ∆)

(D − ∆))

. Если считать интервал

 

 

 

 

 

 

узким D <1, то, представляя логарифм рядом, имеем σ ≈ D 1 (

D)2 6 .

 

 

 

 

 

 

 

П р и м е р 5.

Аппаратура состоит из

N = 2000

элементов. Вероятность

выхода из строя каждого из них w = 0.0005 . Определить вероятность выхода из строя аппаратуры.

Р е ш е н и е . Аппаратура не работает при выходе из строя хотя бы одного элемента. Искомая вероятность определяется из (2.19) при k =1:

N 1

N

WN = CNj wN j (1 w) j = CNj wN j (1w) j (1w)N =1(1w)N = 0.63 .

j=0

j=0

Здесь учтено, что распределение Бернулли удовлетворяет условию нормировки

N

 

 

 

 

 

 

 

 

 

( CNj wN j (1 w) j =1). Поскольку по

условию

задачи

 

 

N >>1, w <<1, то

j=0

 

 

 

 

 

 

 

 

 

можно воспользоваться и распределением Пуассона:

 

 

 

N 1

N

 

 

 

 

 

 

 

 

WN (Nw)N j eNw (N j) ! = (Nw)N j eNw (N j)!eNw =1eNw .

j=0

j=0

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

)

 

 

Заметим, что lim Nw =

 

и

lim (1 w)N

= lim

(1n

N)N

 

n

= eNw .

n

n

N → ∞

N → ∞

N → ∞

 

 

 

 

 

w0

w0

 

 

 

 

 

 

 

 

2.6. Задачи и упражнения для самостоятельной работы

2.1. Производится три независимых выстрела по мишени, вероятность попадания при каждом выстреле равна p = 0.4 . Определить вероятность: 1. двух и трех попаданий; 2. хотя бы одного попадания; 3. промаха при трех выстрелах.

2.2.Используя условия предыдущей задачи, найти среднее число попаданий

вмишень, их дисперсию и относительную флуктуацию.

2.3.

Плотность

вероятности

случайной

величины

x

имеет

вид

f (x) = Ae

 

x1

 

. Определить

среднее

значение,

дисперсию и

относительную

 

 

 

 

флуктуацию x .

 

 

 

 

 

 

 

 

2.4.

 

 

Система

характеризуется

распределением

 

вероятностей

f (x, y) = Axydxdy ,

когда

0 x a ,

0 y b ,

( a =1, b = 2) ,

а вне

этих

 

 

 

23

 

 

 

интервалов

f (x, y) = 0 . Найти

1)

вероятность, что система имеет

значение

величины x

в интервале x [0.5;

0.6]

при любом значении

y ; 2)

средние

значения x ,

y и их произведения.

 

 

 

 

2.5. Вероятность работы

прибора

(надежность) равна

p = 0.9 . Для

повышения надежности он дублируется такими же самыми. Сколько нужно поставить дублирующих приборов, чтобы надежность схемы была не ниже

P = 0.999 .

2.6. Плотность распределения случайных величин x и y задана законом Коши f (x, y) = A [(1 + x2 )(1 + y2 )] . Найти среднее значение x , y и вероятность того, что 1) каждая из величин x и y при произвольном значении другой будет находиться в интервале от 1 до 1; 2) две величины находятся в заданных интервалах.

2.7. Распределение случайной величины x задано распределением Гаусса.

Доказать, что x4 = 3(x2 )2 .

2.8.Найти вероятность выпадения хотя бы раз двух шестерок при 24 бросаниях пары игральных костей.

2.9.Как велика вероятность того, что при шестикратном бросании игральной кости будет выброшена а) только одна шестерка; б) по крайней мере одна шестерка.

2.10.Из двух колод вынимается наугад по карте. Определить вероятность того, что 1) обе карты окажутся масти пик; 2) хотя бы одна из карт окажется масти пик.

2.11.На полупрямой случайно ставятся три точки. Найти вероятность, что из трех отрезков, равных расстояниям этих точек от начала полупрямой, можно составить треугольник.

2.12.Найти вероятность, что лампа, работающая в момент времени t0 , не выйдет из строя до момента t0 + t при следующих условиях: 1) вероятность

зависит только от величины промежутка времени [t0 ; t0 +t] ; 2) вероятность выхода лампы из строя за промежуток времени t пропорциональнаt ;

3)события в непересекающиеся промежутки времени независимы между собой.

2.13.Определить время выхода из строя системы, содержащей n электронных ламп, если средний срок службы каждой лампы составляет τ

24

часов. Вероятность выхода из строя лампы пропорциональна времени, а коэффициент пропорциональности – q предполагается постоянным, т.е. q(t) = const .

2.14.Определить интервал времени между сменами ламп приемника, содержащего 6 ламп с гарантийным сроком службы 1000 часов. Предполагается одинаковая вероятность выхода из строя ламп в течение гарантированного срока службы и возможность эксплуатации приемника с надежностью P = 0.90 .

2.15.Определить гарантированный срок службы ламп телевизора, работающего с надежностью P = 0.90 , если потребитель меняет все 12 ламп через t = 500 часов эксплуатации.

2.16.Выход из строя радиоаппаратуры составляет n% . Для увеличения надежности поставлены m дублирующих цепей. Какова надежность системы?

2.17.Вероятность выхода из строя радиолампы до тренировки составляет

P1(в) = 0.1, после тренировки P2 (в) = 0.05 . Во сколько раз возрастает надежность радиолампы после проведения тренировок?

2.18.Идеальный газ находится в коническом сосуде высотой H . Найти функцию распределения его молекул по высоте z в отсутствие сил тяжести.

2.19.Разрыв электрической цепи наступает вследствие выхода из строя элемента k или двух элементов k1 и k2 , которые выходят из строя независимо

друг от друга соответственно с вероятностями0.3; 0.2 и 0.2 . Определить вероятность разрыва электрической цепи.

2.20.Найти значения x и x2 и x2 при равномерном распределении величины x между a и b .

2.21.Получить выражение для вероятности, найти гармонически колеблющуюся точку в интервале d x на линии колебания между крайними

значениями A и A.

2.22.Найти функцию распределения угла между полупрямыми на плоскости, из которых одна закреплена, а у другой все ориентации равновероятны.

2.23.Найти функцию распределения угла между закрепленной полупрямой

иполупрямой, все направления которой в пространстве равновероятны.

25

2.24. Плотность распределения вероятности случайных амплитуд боковой качки корабля ( x 0 ) имеет вид (закон Релея) f (x) =(xa2 ) exp(x2 2a2 ), где

a наиболее вероятное отклонение корабля. Найти среднее отклонение корабля во время качки.

2.25. Пространство заполняет газ. Вероятность встретить молекулу газа внутри бесконечно малого объема dv равнаαdv . Для любой молекулы в любой момент времени найдется какая-то молекула ближайший сосед. Расстояние до ближайшего соседа есть случайная величина. В разные моменты времени она различна. Найти функцию распределения расстояния до ближайшего соседа.

2.26.В центре основания кругового цилиндра радиуса r находится источник излучения. Найти функцию распределения случайной величины высоты попадания фотона в стенку цилиндра.

2.27.Источник излучения находится над поверхностью прозрачного вещества. Найти функцию распределения углов преломления лучей в прозрачном веществе (отношение показателей преломления в прозрачной среде

ивоздухе равно æ).

2.28.Вероятность столкновения двух молекул разреженного газа на участке dx равна adx . Найти вероятность прохождения молекулой пути x без

столкновения, среднюю длину свободного пробега. Рассмотреть возможность экспериментального определения длины свободного пробега.

2.29.Определить вероятность предсказания трех и шести чисел из 49 в “спортлото”.

2.30.Система радиолокационных станций отслеживает N = 20 независимых целей. Вероятность потери одной цели равна w = 0.05 . Оценить вероятность потери а) S – целей ( S =1, 2 ); б) хотя бы m – целей ( m =1, 2 ).

26

3.ФАЗОВОЕ ПРОСТРАНСТВО. РАСПРЕДЕЛЕНИЯ ГИББСА

3.1Уравнения движения

Если изолированная система состоит из N частиц (материальных точек), то уравнение движения j -ой точки имеет вид

 

 

 

N

 

 

 

 

 

 

 

 

 

d Pj

dt = Pj

= Fkj ,

(3.1)

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

где Pj импульс j -ой частицы,

Fkj сила, с которой k -ая частица действует

на j -ую. При

наличии в

системе

связей не

все уравнения (3.1) будут

независимыми;

существуют

ограничения,

которые

записываются

дополнительными уравнениями

 

 

 

 

 

 

 

= 0 ,

ν =1, 2, ..., m

(3.2)

 

fν r1, ..., rN

, t

 

 

 

 

 

 

 

 

относительно координат частиц rj в системе центра масс. Последние позволяют исключить m координат из общего числа 3N и ввести независимые

переменные обобщенные

координаты q1,

q2 ,...,

q3N m , через которые и

(q1, q2 ,...,

 

t ), n =1, 2,..., N .

выражаются координаты rn :

rn

= rn

q3N m ,

Для всякой системы можно ввести функцию Лагранжа

L(q1, q2 ,..., qs ; q1, q2 ,..., qs ; t )= L(q; q; t ),

где qj обобщенные координаты; q j = dq j dt обобщенные скорости; s = 3N m число степеней свободы.

Движение системы частиц в этом случае описывается уравнениями Лагранжа второго рода

d (dLdq j ) dt = dL dqj , j =1, 2, ..., s . (3.3)

Для статистической механики принята иная формулировка уравнения

движения, основанная на введении обобщенного

импульса pj

= ∂L q j и

функции Гамильтона H ( p, q) :

 

 

 

q j = ∂H pj ;

pj = −∂H qj ,

j =1, 2, ..., s .

(3.4)

В случае консервативной системы функция Гамильтона есть ее полная

27

механическая энергия

H ( p, q) =T ( p, q) +U (q) ,

(3.5)

где T ( p, q) кинетическая, а U (q) – потенциальная энергии системы.

3.2. Фазовое пространство

Под фазовым понимается воображаемое 6N -мерное пространство, осями которого являются 3N обобщенных координат и 3N обобщенных импульсов. Состояние системы из N частиц в фазовом пространстве в произвольный момент времени задается изобразительной точкой. Изменение координат и импульса любой частицы системы приводит к перемещению изобразительной точки в фазовом пространстве. Эволюцию во времени изобразительной точки описывает фазовая траектория. В силу единственности решения уравнения движения с начальными условиями траектория нигде не имеет точек пересечения. Фазовое пространство одной молекулы называется µ-пространством, для одноатомной молекулы µ-пространство шестимерное.

Элементом объема фазового пространства для частиц с непрерывным спектром значений координат и импульсов называется произведение

3

N

 

d Γ = ∏∏d p jk d qjk = d pG d qG = d Γp d Γq .

(3.6)

k =1

j=1

 

Для свободных частиц (не взаимодействуют между собой в отсутствие внешних сил) элемент объема d Γq равен объему пространстваV , где они

могут находиться: d Γq =V . При дискретном спектре значений импульса и координат элемент объема фазового пространства системы с f степенями свободы равен

d Γ = d Γp d Γq = h f ,

(3.7)

где h постоянная Планка. Здесь речь идет, кроме квантовых, и о квазиклассических коллективах. Квазиклассическими называют системы, частицы которых подчиняются классическим уравнениям движения, но не все траектории реализуются. Если такие частицы считать свободными, то

d Γp = h f V .

(3.8)

Число состояний (плотность состояний) в заданном фазовом объеме для

28

 

 

 

 

 

квантовых и квазиклассических систем равно

 

)

 

 

( p, q) d pG d qG = ∆Γ h f =

(

h f

(∂Γ ∂ε)dε.

(3.9)

1

 

3.3. Теорема Лиувилля

При движении изобразительной точки в фазовом пространстве по фазовой траектории форма элемента фазового пространства dΓ с течением времени меняется произвольным образом, однако величина его при этом остается неизменной. Это утверждение называется теоремой Лиувилля и записывается в виде df (q, p, t)dt = 0 , где f (q, p, t) функция распределения или плотность изобразительных точек.

3.4. Распределения Гиббса (РГ)

РГ отвечают на вопрос о том, какова вероятность wi найти систему в состояниях с энергией, заключенной между εi и εi + δε i ( δεi << εi , i =1, 2, 3,...). Каждому значению энергии εi отвечает своя плотность (εi ) состояний.

Для решения поставленной задачи вводится очень большое число S копий данной системы (тождественных систем) в какой-либо момент времени t . Речь идет об одном гамильтониане H ( pi , qi ) , но для разных начальных условий. Такой ансамбль систем распределен в фазовом пространстве с плотностью вероятности f ( pi , qi , t) . Выделенному интервалу энергий соответствует фазовый объем ∆Γ в определенном месте фазового пространства с n числом систем ансамбля. Поэтому вероятность нахождения системы в заданном состоянии равна

wi = lim

(n S )= f ( pi , qi , t) ∆Γ.

(3.10)

S→∞

 

 

Различают замкнутые (изолированные) и квазинезависимые, слабо

взаимодействующие с окружающей средой системы.

 

3.4.1.

М и к р о к а н о н и ч е с к о е

р а с п р е д е л е н и е .

В

термодинамике под замкнутыми условно понимают такие системы, энергия которых за все время наблюдения остается заключенной в заданных узких пределах δεi . Можно считать, что ее энергия является постоянной

H ( pi , qi ) = Ei , а изобразительные точки такого ансамбля распределены по гиперповерхности постоянной энергии. Плотность вероятности в этом случае

29

 

находится из распределения

 

f [H ( pi , qi , t)]= δ[H ( pi , qi , t) E] (E) ,

(3.11)

где δ(x) дельта-функция Дирака, которая определяется для непрерывной функции ϕ(x) при ϕ(x → ±∞) 0 операторным соотношением

ϕ(x) δ(x a) dx = ϕ(a) .

Поскольку все состояния с данной энергией εi равновероятны, вероятность, что замкнутая система находится в одном из них, пропорциональна плотности числа состояний

w(εi ) ~ (εi ) .

(3.12)

Формулы (3.11) или (3.12) суть микроканоническое распределение Гиббса.

3.4.2. К а н о н и ч е с к о е р а с п р е д е л е н и е .

Рассматриваются

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

w(εi ) = e− εi θ(εi ) e− εi θ(εi ) ,

(3.13.а)

i

 

dw(εi ) = e− ε θdΓ e− ε θdΓ = e− ε θ (∂Γ ∂ε)dε e− ε θ (∂Γ ∂ε)dε .

(3.13.б)

Первое из них относится к системам с дискретным, а второе с непрерывным спектром значений энергий. Здесь учтено условие нормировки w(εi ) =1 и

 

 

 

i

 

 

 

 

 

e− ε θ (∂Γ ∂ε)dε =1.

Модуль

статистического

распределения

или

− ∞

 

 

 

 

с т а т и с т и ч е с к а я

т е м п е р а т у р а θ = kT

является существенно

положительной величиной.

Число состояний (εi ) быстро растет с увеличением энергии и числа частиц в системе. Произведение экспоненциально убывающей с ростом энергии ( e− εi θ ) и быстро возрастающей функций приводит к резкому максимуму в распределении Гиббса. Поэтому среднее значение функций от любой физической величины макроскопических систем ( N ~ 1026 ) практически