Автокорреляция для функции времени f(t) представляет собой временно́е среднее от произведения f(t+τ) на f(t). Удобно вести комплексные функции времени, если даже в реальных случаях мы рассматриваем действительные функции. Тогда автокорреляция становится равной среднему произведению f(t+τ) на величину, сопряженную с f(t). Работаем ли мы с действительными или с комплексными функциями, спектр мощности функции f(t) равен преобразованию Фурье от ее автокорреляции.
Я уже говорил о непригодности чернильных записей для дальнейшей математической обработки. Прежде чем ожидать многого от идеи автокорреляции, необходимо было заменить чернильные записи какими-либо другими, более пригодными.
Одним из лучших способов фиксации малых флюктуирующих напряжений для дальнейшей обработки - применение магнитной ленты. Она позволяет сохранять флюктуирующее электрическое напряжение в виде постоянной записи, которую можно затем использовать когда угодно. Один из таких приборов был придуман около десяти лет тому назад в научно-исследовательской лаборатории электроники Массачусетсского технологического института под руководством проф. Уолтера А. Розенблита и д-ра Мэри А. Б. Бразье.
В этом приборе применяется запись на магнитную ленту с частотной модуляцией. Дело в том, что считывание всегда связано с некоторым стиранием магнитной ленты. При записи с амплитудной модуляцией стирание приводит к изменению хранимого сообщения, и при последовательных считываниях ленты мы по существу имеем дело с меняющимся сообщением.
При частотной модуляции также происходит некоторое стирание, но приборы, посредством которых мы читаем ленту, сравнительно нечувствительны к амплитуде и считывают только частоту. Пока лента не сотрется настолько, что станет совершенно неразборчива, частичное стирание ленты не искажает значительно сообщения, которое она хранит. Поэтому ленту можно [c.271] читать много раз почти с такой же точностью, как и при первом считывании.
Как следует из самого понятия автокорреляции, нам понадобится механизм, задерживающий считывание ленты на регулируемый интервал времени. Если отрывок записи длительности А пропустить через прибор с двумя последовательными считывающими головками, то образуются два одинаковых, но сдвинутых во времени сигнала. Временной сдвиг зависит от расстояния между считывающими головками и от скорости подачи ленты, и его можно менять по нашему желанию. Мы можем обозначить один сигнал через f(t), а другой - через f(t+τ), где τ - временной сдвиг. Произведение этих сигналов можно, например, получить при помощи квадратических детекторов и линейных смесителей, используя тождество
4ab = (a+b)-(a-b) (10.01)
Это произведение можно приближенно усреднить на интегрирующей реостатно-емкостной цепи, имеющей большую постоянную времени сравнительно с длительностью А нашей выборки. Полученное среднее [c.272] пропорционально значение автокорреляционной функции при задержке τ. Повторение процесса при различных τ даст некоторый ряд значений автокорреляции (или, вернее, выборочной автокорреляции за большое время включения А). На рис. 9 показан график одной реальной автокорреляции такого рода. Заметим, что здесь показана лишь половина кривой, так как автокорреляция для отрицательных времен совпадает с автокорреляцией для положительных времен, по крайней мере, в случае, когда мы отыскиваем автокорреляцию действительной кривой.
Рис. 9. Автокорреляция
Заметим, что подобные автокорреляционные кривые применялись уже много лет в оптике и что прибором, с помощью которого их получали, был интерферометр Майкельсона (рис. 10). Интерферометр Майкельсона посредством системы зеркал и линз разделяет световой луч на две части, которые посылаются по путям разной длины и затем вновь соединяются в один луч. Различные длины путей вызывают различные задержки во [c.273] времени, и результирующий луч будет равен сумме двух отражений входящего луча, которые можно опять обозначить через f(t) и f(t+τ). Если измерить чувствительным фотометром силу луча, то его показание будет пропорционально квадрату суммы f(t)+ f(t+τ) и, следовательно, должно содержать член, пропорциональный автокорреляции. Другими словами, яркость интерференционных полос даст нам автокорреляцию (с точностью до линейного преобразования).
Рис. 10. Интерферометр Майкельсона
Все это неявно содержалось в работе Майкельсона. Нетрудно видеть, что при выполнении преобразования Фурье над интерференционными полосами интерферометр дает нам энергетический спектр света и тем самым по существу является спектрометром. Более того, это самый точный из известных нам типов спектрометров.
Спектрометр такого типа получил должное признание лишь в последние годы. Мне говорили, что теперь он принят в качестве важного средства прецизионных измерений. Отсюда видно, что методы обработки автокорреляционных записей, которые я сейчас изложу, применимы также в спектроскопии и позволяют довести до предела ту информацию, которую может дать спектрометр.
Рассмотрим, как получить спектр мозговой электрической волны по автокорреляции. Пусть C(t) - автокорреляция функции f(t). Тогда C(t) можно записать в виде
(10.02)
Здесь F всегда является возрастающей или по меньшей мере неубывающей функцией от ω; мы будем называть ее интегральным спектром функции f. Вообще говоря, этот интегральный спектр состоит из трех аддитивных частей. Линейчатая часть спектра возрастает лишь на счетном множестве точек. После ее исключения останется непрерывный спектр, равный, в свою очередь, сумме двух частей: одна из них возрастает только на множестве меры нуль, а другая абсолютно непрерывна и является интегралом положительной интегрируемой функции.
Будем впредь полагать, что первые две части спектра: дискретная часть и непрерывная часть, возрастающая [c.274] на множестве меры нуль, - отсутствуют. В этом случае можно написать
(10.03)
где φ (ω) - спектральная плотность. Если φ (ω) принадлежит к классу Лебега L, то можно написать
(10.04)
Как видно по автокорреляционной кривой мозговых волн, преобладающая часть мощности спектра сосредоточена в окрестности частоты 10 гц. В таком случае φ (ω) будет иметь форму, подобную следующей диаграмме: