- •3 В.К. Игнатьев. Статистическая радиофизика
- •Литература:
- •Основы теории вероятностей
- •1.1.Предмет статистической радиофизики
- •1.2.Случайные события. Вероятность
- •1.3.Случайные величины. Распределение вероятностей
- •1.4. Закон больших чисел. Аксиома измерений
- •1.5.Совместные распределения. Условные функции распределения
- •1.6. Характеристическая функция. Семиинварианты
- •1.7.Центральная предельная теорема
- •Случайный импульсный процесс
- •2.1.Функции случайной величины
- •2.2.Пуассоновский импульсный процесс
- •Случайные функции
- •3.1.Понятие случайной функции. Задание случайной функции
- •3.2.Сходимость, дифференцируемость, непрерывность и стационарность случайных функций
- •3.3.Моменты случайных функций
- •3.4.Эргодические случайные процессы
- •Корреляционная теория случайных процессов
- •4.1.Функция автокорреляции
- •4.2.Спектральная плотность интенсивности случайных процессов
- •4.3. Случайные последовательности
- •Воздействие случайного процесса на линейную систему
- •5.1. Спектральные характеристики процесса на выходе линейной системы
- •5.2.Распределение вероятностей на выходе линейной системы
- •5.3.Узкополосный гауссов процесс
- •5.4.Спектр колебаний с флуктуирующей частотой
- •5.5. Спектральное оценивание
- •Нелинейные преобразования случайных процессов
- •6.1.Нелинейное безинерциальное преобразование
- •6.2.Корреляционные функции на выходе умножителя частоты
- •Марковские процессы
- •7.1. Процесс без последействия
- •7.2. Уравнение Смолуховского
- •7.3.Марковский процесс с дискретными состояниями
- •7.4.Двумерные случайные блуждания
- •7.5.Уравнение Фоккера – Планка – Колмогорова
- •Стохастические дифференциальные уравнения
- •8.1. Постановка задачи
- •8.2. Случайные функции с независимыми приращениями
- •8.3. Усреднение точного решения стохастического дифференциального уравнения
- •8.4. Уравнение для средних
- •8.5.Уравнение Лиувилля
- •Случайные поля
- •9.1.Функция автокорреляции и спектр случайного поля
- •9.2.Электромагнитная волна в статистически неоднородной среде
- •9.3.Метод медленно меняющихся амплитуд
- •9.4.Плоская случайная волна в диспергирующей среде
- •Флуктуации в электрических цепях
- •10.1. Тепловой шум в линейных диссипативных системах
- •10.2.Дробовой шум
- •10.3.Фликкер-шум
- •10.4. Шумы электронно-дырочного перехода
- •10.5. Шум биполярного транзистора
- •10.6.Шумы полевых транзисторов
- •10.7.Шумы усилителей
- •Флуктуации в лазерных системах
- •11.1.Корреляционная функция одномодового лазерного излучения
- •11.2. Корреляционная функция многомодового лазера
- •11.3.Флуктуации в одномодовом лазере
- •Содержание
9.4.Плоская случайная волна в диспергирующей среде
Пусть источник поля находится в сечении z= 0 и создаёт волну с комплексной амплитудой:
. (9.0)
Случайный процесс A0(t) будем считать стационарным с корреляционной функцией вида:
B0() = A0(t)A0*(t + ) = I0exp(–2/k2).
Решение линейного дифференциального уравнения (9.14)с граничными условиями (9.15)можно записать в виде свёртки:
, (9.0)
где – функция Грина.
Умножим теперь уравнение (9.14),записанное дляA(z,t) на A*(z,t1) и сложим с уравнением(9.14), записанным дляA*(z,t1) и умноженным на A(z,t), в результате получим:
.
При усреднении уравнения вследствие стационарности процесса временные производные пропадут и получим B(z,)/z= 0, то есть корреляционная функция, а следовательно, и спектр стационарного процесса неизменны:
B(z, ) = B(0, ) = B0(), G(z, ) = G(0, ) = G0().
Рассмотрим теперь модель источника сигнала в виде импульса, заполненного шумом: A0(t) =F(t)(t), гдеF(t)–огибающая,(t) – шумовое заполнение (несущая), стационарный случайный процесс, причём
(t) = 0, , R(0) = 1, max[F(t)] = F0.
Эта задача имеет большой практический интерес для передачи информации по волоконно-оптическому каналу связи, когда нужно выяснить как меняются форма и длительность импульса при его распространении по волокну. Заметим, что форма и длительность импульса на входе волокна в сечении z= 0 зависит только от огибающейF(t),I0(t) =A0(t)A0*(t)= |F(t)|2. Время же корреляции импульса определяется временем корреляции случайного заполнения:
.
Из формулы (9.16)видно, что амплитуда импульсаA(z,t) зависит от
= t – z/u
,
где обозначено .Область интегрирования здесь порядкаи0,поэтому если в дальней зоне, то можно положить, и интеграл примет вид преобразования Фурье:
. (9.0)
Тогда
. (9.0)
Toесть в среде с дисперсией первого порядка импульс в дальней зоне превращается в спектрон, интенсивность которого повторяет спектрG0() входного импульса в линейном масштабе, не меняя его формы. Так как спектр импульса при распространении не меняется, форма спектрона остается устойчивой, по мере ростаzимпульс становится более длинным и низким. При этом длительность импульсаи(z) в дальней зоне определяется шириной полосыимпульса на входе. Но1/к0+ 1/и0= (1 +)/и0, где = и0/к0 –коэффициент некогерентности. Следовательно, некогерентность сильно влияет на расходимость импульса. В частности, при сильной некогерентности>> 1,автокорреляцию в дальней зоне можно записать с учётом соотношения(9.17)в виде:
Таким образом, импульс интенсивности в дальней зоне повторяет спектр шумового заполнения . При этом можно показать, что коэффициент некогерентности в дальней зоне имеет ту же величину, что и на входе, хотя форма импульса и его корреляционная функция совсем другие.
Флуктуации в электрических цепях
10.1. Тепловой шум в линейных диссипативных системах
Даже в отсутствие внешнего воздействия на зажимах линейного двухполюсника присутствует флуктуирующее (случайное) напряжение, а по замкнутой цепи течёт случайный ток, вызванный броуновским движением электронов. Поэтому, в силу теорем Нортона и Тевенинна, реальный двухполюсник следует рассматривать как идеальный (нешумящий) двухполюсник, последовательно с которым включен генератор шумового напряжения enили параллельно –генератор шумового токаin .При этом, естественно,
, (10.0)
где z(t) – импульсная характеристика двухполюсника. Как показано в п. 5.1,спектральные интенсивности шумового токаGi(), который можно рассматривать как входное воздействие, и шумового напряженияGe() связаны соотношением вида (5.3):
Ge() = |Z()|2Gi(), (10.0)
где – полное сопротивление двухполюсника. Найдём среднюю мощность, рассеиваемую двухполюсником, воспользовавшись соотношениями (4.4) и (10.1):
(10.0)
Рассмотрим теперь два двухполюсника Z1=R1+iX1иZ2=R2+iX2, соединенные через третий чисто реактивный двухполюсникZ3=iX3и находящиеся в термостате при температуре Т.При термодинамическом равновесии должен выполняться принцип детального равновесия, то есть мощностьР21,создаваемая генераторомe1во втором двухполюсникеZ2, должна равняться мощности P12,создаваемой генераторомe2в первом двухполюсникеZ1.Воспользовавшись соотношениями (10.2)и (10.3), получим:
Поскольку равенство Р21 = P12 должно выполняться при произвольной функцииZ3(),тоGe1()R2() =Ge2()R1(), или
Ge1()/R1() = Ge2()/R2() = u(, T). (10.0)
Таким образом, отношение спектральной плотности шумового напряжения на зажимах двухполюсника к активному сопротивлению двухполюсника при термодинамическом равновесии является универсальной функцией частоты и температуры и не зависит от вида двухполюсника. Заметим, что из соотношения (10.4) следует, что при Rк() = 0 иGeк() = 0,то есть в чисто реактивных (недиссипативных) двухполюсниках, нет тепловых шумов, поэтому в реактивный двухполюсникZ3() не включен генератор шумового напряжения.
Для того чтобы найти вид функции u(,T),рассмотрим высокодобротный последовательныйRLC-контур, который описывается уравнением
,
где q –заряд на конденсаторе. Будем считать, чтоkT>>ħи выполняется принцип равнораспределения (эргодическая гипотеза) –на каждую степень свободы системы приходится энергияkT/2.Это значит, что
. (10.0)
Но в силу соотношений (4.5)и (10.2),получаем:
.
Здесь мы учли, что
контур высокодобротный, то есть
,Z()
=R+iL–i/(C).Аналогично легко получить
.
Подставляя эти выражения в соотношение (10.5),получаем с учётом того, что резонансная частота контура0может быть выбрана произвольной:
,
или
Ge() = kTR()/ = kTRe[Z()]/ – (10.0)
теорема Найквиста. Сравнивая её с соотношением (10.4),видим, что
u(, T) = kT/.
Для физических частот f=/(2) > 0 можно записать теорему Найквиста (10.6) с учетом формулы (4.6) в виде:
Ge(f) = 4kTR(f) = 4kTRe[Z(f)].
Воспользовавшись соотношением (10.2),можно записать теорему Найквиста для спектральной плотности интенсивности шумового тока:
Gi() = kTRe[Y()]/ (10.0)
где Y() = 1/Z()–полная проводимость двухполюсника. Если рассмотреть сопротивлениеR, шунтированное емкостьюC,что всегда имеет место на практике, то
.
Тогда
.