Пара преобразований Фурье. Спектральная плотность сигнала. Спектральная плотность мощности детерминированного сигнала

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

Обобщенная формула Рэлея. Докажем важное вспомогательное положение, касающееся спектральных свойств сигналов.

Пусть два сигнала в общем случае комплекснозначные, определены своими обратными преобразованиями Фурье:

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

Здесь внутренний интеграл представляет собой, очевидно, спектральную плотность сигнала . Поэтому

Полученное соотношение представляет собой обобщенную формулу Рэлея. Легко запоминающаяся трактовка этой формулы такова: скалярное произведение двух сигналов с точностью до коэффициента пропорционально скалярному произведению их спектральных плотностей.

Обобщение понятия спектральной плотности.

Будем считать, что сигнал представляет собой абсолютно интегрируемую функцию. Тогда его преобразование Фурье - обычная классическая функция частоты. Пусть наряду с этим сигнал не удовлетворяет условию абсолютной интегрируемости и в обычном классическом смысле преобразование Фурье не существует. Однако можно расширить понятие спектральной плотности, допустив, что является обобщенной функцией в том смысле, который был установлен в § 1.2. Для этого в соответствии с обобщенной формулой Рэлея достаточно положить, что - функционал, который, действуя на известную функцию , дает следующий результат:

Приемы вычисления спектров неинтегрируемых сигналов целесообразно рассмотреть на конкретных примерах.

Спектральная плотность постоянного во времени сигнала. Простейший неинтегрируемый сигнал - это постоянная величина и . Предположим, что - произвольный вещественный абсолютно интегрируемый сигнал с известной спектральной плотностью

Раскрывая формулу (2.43), имеем

Но, как легко заметить,

Отсюда на основании фильтрующего свойства дельтафункции приходим к выводу, что равенство (2.43) возможно лишь при условии, что

Физический смысл полученного результата нагляден - неизменный во времени сигнал имеет спектральную составляющую только на нулевой частоте.

Спектральная плотность комплексного экспоненциального сигнала.

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

Отсюда искомая спектральная плотность S (со), выражается таким образом:

Отметим следующее:

1. Спектральная плотность комплексного экспоненциального сигнала равна нулю всюду, кроме точки где она имеет дельта-особенность.

2. Спектр данного сигнала несимметричен относительно точки и сосредоточивается в области либо положительных, либо отрицательных частот.

Спектральная плотность гармонических колебаний. Пусть По формуле Эйлера

Найденный выше спектр комплексного экспоненциального сигнала, а также свойство линейности преобразования Фурье позволяют сразу записать выражение спектральной плотности косинусоидального сигнала:

Читатель может легко проверитьсамостоятельно, что для синусоидального сигнала справедливо соотношение

Следует заметить, что выражение (2.46) представляет собой четную, а выражение (2.47) - нечетную функцию частоты.

Спектральная плотность произвольного периодического сигнала.

Ранее периодические сигналы исследовались методами теории рядов Фурье. Теперь можно расширить представления об их спектральных свойствах, описав периодические сигналы с помощью преобразования Фурье.

Периодический сигнал, заданный своим рядом Фурье в комплексной форме. На основании формулы (2.45), принимая во внимание свойство линейности преобразования Фурье, сразу получаем выражение спектральной плотности такого сигнала:

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

Спектральная плотность функции включения.

Вычислим спектральную плотность функции включения , которую для простоты определим во всех точках, кроме точки t = 0 [ср. с (1.2)]:

Заметим прежде всего, что функция включения получается путем предельного перехода из экспоненциального видеоимпульса:

Поэтому можно попытаться получить спектральную плотность функции включения, выполнив предельный переход при а- О в формуле спектральной плотности экспоненциального колебания:

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

Прежде всего выделим в спектральной плотности экспоненциального сигнала вещественную и мнимую части:

Можно убедиться в том, что

Действительно, предельное значение этой дроби при любых обращается в нуль, и в то же ремя

независимо от величины а, откуда и следует сделанное утверждение.

Итак, получено взаимно однозначное соответствие функции включения и ее спектральной плотности:

Дельта-особенность при свидетельствует о том, что функция включения имеет постоянную составляющую, равную 1/2.

Спектральная плотность радиоимпульса.

Как известно, радиоимпульс задается в виде произведения некоторого видеоимпульса играющего роль огибающей, и неинтегрируемого гармонического колебания: .

Чтобы найти спектральную плотность радиоимпульса, будем полагать известной функцию - спектр его огибающей. Спектр косинусоидального сигнала с произвольной начальной фазой получается путем элементарного обобщения формулы (2.46):

Спектр радиоимпульса есть свертка

Приняв во внимание фильтрующее свойство дельтафункции, получаем важный результат:

Рис. 2.8 иллюстрирует трансформацию спектра видеоимпульса при умножении его на высокочастотный гармонический сигнал.

Рис. 2.8. Частотные зависимости модуля спектральной плотности: а - видеоимпульса; б - радиоимпульса

Видно, что переход от видеоимпульса к радиоимпульсу при спектральном подходе означает перенос спектра видеоимпульса в область высоких частот - вместо единственного максимума спектральной плотности при наблюдаются два максимума при абсолютные значения максимумов сокращаются вдвое.

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

Пример 2.3. Спектральная плотность прямоугольного радиоимпульса.

Для простоты положим начальную фазу нулевой и запишем математическую модель радиоимпульса в виде

Зная спектр соответствующего видеоимпульса [см. формулу (2.20)], на основании (2.50) находим искомый спектр:

На рис. 2.9 изображены результаты расчета спектральной плотности по формуле (2.51) для двух характерных случаев,

В первом случае (рис. 2.9,а) импульс огибающей содержит 10 периодов высокочастотного заполнения частота здесь достаточно высока для того, чтобы избежать «перекрытия». Во втором случае (рис. 2.9, б) радиоимпульс состоит всего лишь из одного периода заполнения Наложение составляющих, которые соответствуют областям положительных и отрицательных частот, приводит к характерной асимметрии лепестковой структуры графика спектральной плотности радиоимпульса.

Рис. 2.9. Графики спектральных плотностей радиоимпульса с прямоугольной огибающей: а - при ; б - при

При исследовании автоматических систем управления удобно пользоваться еще одной характеристикой стационарного случайного процесса, называемой спектральной плотностью. Во многих случаях, особенно при изучении преобразования стационарных случайных процессов линейными системами управления, спектральная плотность оказывается более удобной характеристикой, чем корреляционная функция. Спектральная плотность случайного процесса определяется как преобразование Фурье корреляционной функцией , т. е.

Если воспользоваться формулой Эйлера то (9.52) можно представить как

Так как нечетная функция то в последнем выражении второй интеграл равен нулю. Учитывая, что четная функция получаем

Так как то из (9.53) следует, что

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

Если спектральная плотность известна, то по формуле обратного преобразования Фурье можно найти соответствующую ей корреляционную функцию:

Используя (9.55) и (9.38), можно установить важную зависимость между дисперсией и спектральной плотностью случайного процесса:

Термин «спектральная плотность» обязан своим происхождением теории электрических колебаний. Физический смысл спектральной плотности можно пояснить следующим образом.

Пусть - напряжение, приложенное к омическому сопротивлению 1 Ом, тогда средняя мощность рассеиваемая на этом сопротивлении за время равна

Если увеличивать интервал наблюдения до бесконечных пределов и воспользоваться (9.30), (9.38) и (9.55) при то можно формулу для средней мощности записать так:

Равенство (9.57) показывает, что средняя мощность сигнала может быть представлена в виде бесконечной суммы бесконечно малых слагаемых , которая распространяется на все частоты от 0 до

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

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

Взаимная спектральная плотность двух стационарных случайных процессов определяется как преобразование Фурье от взаимной корреляционной функции т. е.

По взаимной спектральной плотности можно, применяя к (9.58) обратное преобразование Фурье, найти выражение для взаимной корреляционной функции:

Взаимная спектральная плотность является мерой статистической связи между двумя стационарными случайными процессами: Если процессы некоррелированы и имеют равные нулю средние значения, то взаимная спектральная плотность равна нулю, т. е.

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

рассмотрим некоторые свойства спектральных плотностей

1 Спектральная плотность чистого случайного процесса, или белого шума, постоянна во всем диапазоне частот (см. рис. 9.5, г):

Действительно, подставляя в (9.52) выражение (9.47) для корреляционной функции белого шума, получим

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

Происхождение термина «белый шум» объъясняется аналогией такого процесса с белым светом, имеющим одинаковые интенсивности всех компонент, и тем, что случайные процессы типа белого шума впервые были выделены при исследовании тепловых флуктуациоиных шумов в радиотехнических устройствах.

2. Спектральная плотность постоянного сигнала представляет собой -функцию, расположенную в начале координат (см. рис. 9.5, а), т. е.

Чтобы доказать это, допустим, что спектральная плотность имеет вид (9.62), и иандем по (9.55) соответствующую ей корреляционную функцию. Так как

то при получаем

Это (в соответствии со свойством 5 корреляционных функций) означает, что сигнал, соответствующий спектральной плотности, определяемой (9.62), является постоянным сигналом, равным

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

3. Спектральная плотность периодического сигнала представляет собой две -функции, расположенные симметрично относительно начала кординат при (см. рис. 9.5, д), т. е.

Чтобы доказать это, допустим, что спектральная плотность имеет вид (9.63), и найдем по (9.55) соответствующую ей корреляционную функцию:

Это (в соответствии со свойством 6 корреляционных функций) означает, что сигнал, соответствующий спектральной плотности определяемой (9.63), является периодическим сиг налом, равным

Тот факт, что спектральная плотность представляет собой две -функции, расположенные при означает, что вся мощность периодического сигнала сосредоточена на двух частотах: Если рассматривать спектральную плотность только в области положительных частот, то получим,

что вся мощность периодического сигнала будет сосредоточена на одной частоте .

4. Спектральная плотность временной функции, разлагаемой в ряд Фурье имеет на основании изложенного выше вид

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

которая полностью совпадает с корреляционной функцией, определяемой по (9.45).

Из рис. 9.5, б, в видно, что чем шире график спектральной плотности тем уже график соответствующей корреляционной функции и наоборот. Это соответствует физической сущности процесса: чем шире график спектральной плотности, т. е. чем более высокие частоты представлены в спектральной плотности, тем выше степень изменчивости случайного процесса и тем же графики корреляционной функции. Другими словами, связь между видом спектральной плотности и видом функции времени получается обратной по сравнению со связью между корреляционной функцией и видом функции времени. Это особенно ярко проявляется при рассмотрении постоянного сигнала и белого шума. В первом случае корреляционная функция имеет вид горизонтальной прямой, а спектральная плотность имеет вид -функции (см. рис. 9.5, а). Во втором случае (см. рис. 9.5, г) имеет место обратная картина.

6. Спектральная плотность случайного процесса, на кото рой наложены периодические составляющие, содержит непрерывную часть и отдельные -функции, соответствующие частотам периодических составляющих.

Отдельные пики на графике спектральной плотности указывают на то, что случайный процесс смешан со скрытыми периодическими составляющими, которые могут и не обнаруживаться при первом взгляде на отдельные записи процесса. Если, например, на случайный процесс наложен один периодический сигнал с частотой то график; сцектральной плотности имеет вид, показанный на рис. 9.10,

Иногда в рассмотрение вводят нормированную

спектральную плотность являющуюся изображением Фурье нормированной корреляционной функции (9.48):

Нормированная спектральная плотность имеет размерность времени.

В статистической радиотехнике и физике при изучении детерминированных сигналов и случайных процессов широко используется их спектральное представление в виде спектральной плотности, которая базируется на преобразовании Фурье .

Если процесс имеет конечную энергию и квадратично интегрируем (а это нестационарный процесс), то для одной реализации процесса можно определить преобразование Фурье как случайную комплексную функцию частоты:

X (f) = ∫ − ∞ ∞ x (t) e − i 2 π f t d t . {\displaystyle X(f)=\int \limits _{-\infty }^{\infty }x(t)e^{-i2\pi ft}dt.} (1)

Однако она оказывается почти бесполезной для описания ансамбля. Выходом из этой ситуации является отбрасывание некоторых параметров спектра, а именно спектра фаз, и построении функции, характеризующей распределение энергии процесса по оси частот. Тогда согласно теореме Парсеваля энергия

E x = ∫ − ∞ ∞ | x (t) | 2 d t = ∫ − ∞ ∞ | X (f) | 2 d f . {\displaystyle E_{x}=\int \limits _{-\infty }^{\infty }|x(t)|^{2}dt=\int \limits _{-\infty }^{\infty }|X(f)|^{2}df.} (2)

Функция S x (f) = | X (f) | 2 {\displaystyle S_{x}(f)=|X(f)|^{2}} характеризует, таким образом, распределение энергии реализации по оси частот и называется спектральной плотностью реализации. Усреднив эту функцию по всем реализациям можно получить спектральную плотность процесса.

Перейдем теперь к стационарному в широком смысле центрированному случайному процессу x (t) {\displaystyle x(t)} , реализации которого с вероятностью 1 имеют бесконечную энергию и, следовательно, не имеют преобразования Фурье. Спектральная плотность мощности такого процесса может быть найдена на основании теоремы Винера-Хинчина как преобразование Фурье от корреляционной функции:

S x (f) = ∫ − ∞ ∞ k x (τ) e − i 2 π f τ d τ . {\displaystyle S_{x}(f)=\int \limits _{-\infty }^{\infty }k_{x}(\tau)e^{-i2\pi f\tau }d\tau .} (3)

Если существует прямое преобразование, то существует и обратное преобразование Фурье , которое по известной определяет k x (τ) {\displaystyle k_{x}(\tau)} :

k x (τ) = ∫ − ∞ ∞ S x (f) e i 2 π f τ d f . {\displaystyle k_{x}(\tau)=\int \limits _{-\infty }^{\infty }S_{x}(f)e^{i2\pi f\tau }df.} (4)

Если полагать в формулах (3) и (4) соответственно f = 0 {\displaystyle f=0} и τ = 0 {\displaystyle \tau =0} , имеем

S x (0) = ∫ − ∞ ∞ k x (τ) d τ , {\displaystyle S_{x}(0)=\int \limits _{-\infty }^{\infty }k_{x}(\tau)d\tau ,} (5)
σ x 2 = k x (0) = ∫ − ∞ ∞ S x (f) d f . {\displaystyle \sigma _{x}^{2}=k_{x}(0)=\int \limits _{-\infty }^{\infty }S_{x}(f)df.} (6)

Формула (6) с учетом (2) показывает, что дисперсия определяет полную энергию стационарного случайного процесса, которая равна площади под кривой спектральной плотности. Размерную величину S x (f) d f {\displaystyle S_{x}(f)df} можно трактовать как долю энергии, сосредоточенную в малом интервале частот от f − d f / 2 {\displaystyle f-df/2} до f + d f / 2 {\displaystyle f+df/2} . Если понимать под x (t) {\displaystyle x(t)} случайный (флуктуационный) ток или напряжение, то величина S x (f) {\displaystyle S_{x}(f)} будет иметь размерность энергии [В 2 /Гц] = [В 2 с]. Поэтому S x (f) {\displaystyle S_{x}(f)} иногда называют энергетическим спектром . В литературе часто можно встретить другую интерпретацию: σ x 2 {\displaystyle \sigma _{x}^{2}} – рассматривается как средняя мощность, выделяемая током или напряжением на сопротивлении 1 Ом. При этом величину S x (f) {\displaystyle S_{x}(f)} называют спектром мощности случайного процесса.

Свойства спектральной плотности

  • Энергетический спектр стационарного процесса (вещественного или комплексного) – неотрицательная величина:
S x (f) ≥ 0 {\displaystyle S_{x}(f)\geq 0} . (7)
  • Энергетический спектр вещественного стационарного в широком смысле случайного процесса есть действительная и четная функция частоты:
S x (− f) = S x (f) {\displaystyle S_{x}(-f)=S_{x}(f)} . (8)

Величина, характеризующая распределение энергии по спектру сигнала и называемая энергетической спектральной плотностью, существует лишь для сигналов, У которых энергия за бесконечный интервал времени конечна и, следовательно, к ним применимо преобразование Фурье.

Для незатухающих во времени сигналов энергия бесконечно велика и интеграл (1.54) расходится. Задание спектра амплитуд невозможно. Однако средняя мощность Рср, определяемая соотношением

оказывается конечной. Поэтому применяется более широкое понятие "спектральная плотность мощности". Определим ее как производную средней мощности сигнала по частоте и обозначим Сk(щ):

Индексом k подчеркивается, что здесь мы рассматриваем спектральную плотность мощности как характеристику детерминированной функции u(t), описывающей реализацию сигнала.

Эта характеристика сигнала менее содержательна, чем спектральная плотность амплитуд, так как лишена фазовой информации [см. (1.38)]. Поэтому однозначно восстановить по ней исходную реализацию сигнала невозможно. Однако отсутствие фазовой информации позволяет применить это понятие к сигналам, у которых фаза не определена.

Для установления связи между спектральной плотностью Сk(щ) и спектром амплитуд воспользуемся сигналом u(t), существующим на ограниченном интервале времени (-T<. t

где - спектральная плотность мощности сигнала, ограниченного во времени.

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

Функция автокорреляции детерминированного сигнала

Теперь в частотной области имеется две характеристики: спектральная характеристика и спектральная плотность мощности. Спектральной характеристике, содержащей полную информацию о сигнале u(t), соответствует преобразование Фурье в виде временной функции. Выясним, чему соответствует во временной области спектральная плотность мощности, лишенная фазовой информации.

Следует предположить, что одной и той же спектральной плотности мощности соответствует множество временных функций, различающихся фазами. Советским ученым Л.Я. Хинчиным и американским ученым Н. Винером практически одновременно было найдено обратное преобразование Фурье от спектральной плотности мощности:


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

В статистической радиотехнике и физике при изучении детерминированных сигналов и случайных процессов широко используется их спектральное представление в виде спектральной плотности, которая базируется на преобразовании Фурье .

Если процесс x(t) имеет конечную энергию и квадратично интегрируем (а это нестационарный процесс), то для одной реализации процесса можно определить преобразование Фурье как случайную комплексную функцию частоты:

Функция S_x(f)=|X(f)|^2 характеризует, таким образом, распределение энергии реализации по оси частот и называется спектральной плотностью реализации. Усреднив эту функцию по всем реализациям можно получить спектральную плотность процесса.

Перейдем теперь к стационарному в широком смысле центрированному случайному процессу x(t), реализации которого с вероятностью 1 имеют бесконечную энергию и, следовательно, не имеют преобразования Фурье. Спектральная плотность мощности такого процесса может быть найдена на основании теоремы Винера-Хинчина как преобразование Фурье от корреляционной функции:

Если полагать в формулах (3) и (4) соответственно f=0 и \tau=0, имеем

5
6

Формула (6) с учетом (2) показывает, что дисперсия определяет полную энергию стационарного случайного процесса, которая равна площади под кривой спектральной плотности. Размерную величину S_x(f)df можно трактовать как долю энергии, сосредоточенную в малом интервале частот от f-df/2 до f+df/2. Если понимать под x(t) случайный (флуктуационный) ток или напряжение, то величина S_x(f) будет иметь размерность энергии [В 2 /Гц] = [В 2 с]. Поэтому S_x(f) иногда называют энергетическим спектром . В литературе часто можно встретить другую интерпретацию: \sigma_x^2 – рассматривается как средняя мощность, выделяемая током или напряжением на сопротивлении 1 Ом. При этом величину S_x(f) называют спектром мощности случайного процесса.

Свойства спектральной плотности

  • Энергетический спектр стационарного процесса (вещественного или комплексного) – неотрицательная величина:
  • Корреляционная функция k_x(\tau) и энергетический спектр S_x(f) стационарного в широком смысле случайного процесса обладают всеми свойствами, характерными для пары взаимных преобразований Фурье . В частности, чем «шире» спектр S_x(f) тем «уже» корреляционная функция k_x(\tau), и наоборот. Этот результат количественно выражается в виде принципа или соотношения неопределенности.

См. также

Напишите отзыв о статье "Спектральная плотность"

Литература

  1. Зюко, А. Г. Теория передачи сигналов / А. Г. Зюко [и др.]. - М .: Связь, 1980. - 288 с.
  2. Тихонов, В. И. Статистический анализ и синтез радиотехнических устройств и систем / В. И. Тихонов, В. Н. Харисов. - М .: Радио и связь, 2004. - 608 с. - ISBN 5-256-01701-2 .
  3. Тихонов, В. И. Статистическая теория радиотехнических устройств / В. И. Тихонов, Ю. Н. Бакаев. - М .: Академия им. проф. Н. Е. Жуковского, 1978. - 420 с.

Отрывок, характеризующий Спектральная плотность

«Ну и пускай такой то обокрал государство и царя, а государство и царь воздают ему почести; а она вчера улыбнулась мне и просила приехать, и я люблю ее, и никто никогда не узнает этого», – думал он.
Пьер все так же ездил в общество, так же много пил и вел ту же праздную и рассеянную жизнь, потому что, кроме тех часов, которые он проводил у Ростовых, надо было проводить и остальное время, и привычки и знакомства, сделанные им в Москве, непреодолимо влекли его к той жизни, которая захватила его. Но в последнее время, когда с театра войны приходили все более и более тревожные слухи и когда здоровье Наташи стало поправляться и она перестала возбуждать в нем прежнее чувство бережливой жалости, им стало овладевать более и более непонятное для него беспокойство. Он чувствовал, что то положение, в котором он находился, не могло продолжаться долго, что наступает катастрофа, долженствующая изменить всю его жизнь, и с нетерпением отыскивал во всем признаки этой приближающейся катастрофы. Пьеру было открыто одним из братьев масонов следующее, выведенное из Апокалипсиса Иоанна Богослова, пророчество относительно Наполеона.
В Апокалипсисе, главе тринадцатой, стихе восемнадцатом сказано: «Зде мудрость есть; иже имать ум да почтет число зверино: число бо человеческо есть и число его шестьсот шестьдесят шесть».
И той же главы в стихе пятом: «И даны быта ему уста глаголюща велика и хульна; и дана бысть ему область творити месяц четыре – десять два».
Французские буквы, подобно еврейскому число изображению, по которому первыми десятью буквами означаются единицы, а прочими десятки, имеют следующее значение:
a b c d e f g h i k.. l..m..n..o..p..q..r..s..t.. u…v w.. x.. y.. z
1 2 3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160
Написав по этой азбуке цифрами слова L"empereur Napoleon [император Наполеон], выходит, что сумма этих чисел равна 666 ти и что поэтому Наполеон есть тот зверь, о котором предсказано в Апокалипсисе. Кроме того, написав по этой же азбуке слова quarante deux [сорок два], то есть предел, который был положен зверю глаголати велика и хульна, сумма этих чисел, изображающих quarante deux, опять равна 666 ти, из чего выходит, что предел власти Наполеона наступил в 1812 м году, в котором французскому императору минуло 42 года. Предсказание это очень поразило Пьера, и он часто задавал себе вопрос о том, что именно положит предел власти зверя, то есть Наполеона, и, на основании тех же изображений слов цифрами и вычислениями, старался найти ответ на занимавший его вопрос. Пьер написал в ответе на этот вопрос: L"empereur Alexandre? La nation Russe? [Император Александр? Русский народ?] Он счел буквы, но сумма цифр выходила гораздо больше или меньше 666 ти. Один раз, занимаясь этими вычислениями, он написал свое имя – Comte Pierre Besouhoff; сумма цифр тоже далеко не вышла. Он, изменив орфографию, поставив z вместо s, прибавил de, прибавил article le и все не получал желаемого результата. Тогда ему пришло в голову, что ежели бы ответ на искомый вопрос и заключался в его имени, то в ответе непременно была бы названа его национальность. Он написал Le Russe Besuhoff и, сочтя цифры, получил 671. Только 5 было лишних; 5 означает «е», то самое «е», которое было откинуто в article перед словом L"empereur. Откинув точно так же, хотя и неправильно, «е», Пьер получил искомый ответ; L"Russe Besuhof, равное 666 ти. Открытие это взволновало его. Как, какой связью был он соединен с тем великим событием, которое было предсказано в Апокалипсисе, он не знал; но он ни на минуту не усумнился в этой связи. Его любовь к Ростовой, антихрист, нашествие Наполеона, комета, 666, l"empereur Napoleon и l"Russe Besuhof – все это вместе должно было созреть, разразиться и вывести его из того заколдованного, ничтожного мира московских привычек, в которых, он чувствовал себя плененным, и привести его к великому подвигу и великому счастию.
Пьер накануне того воскресенья, в которое читали молитву, обещал Ростовым привезти им от графа Растопчина, с которым он был хорошо знаком, и воззвание к России, и последние известия из армии. Поутру, заехав к графу Растопчину, Пьер у него застал только что приехавшего курьера из армии.