Метод статистической стабилизации частот независимо функционирующих генераторов

Содержание

Основная часть

выводы

Библиографический список


Основная часть

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

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

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

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

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

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

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

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

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


,

где  – функция Хевисайда.

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

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

, (1)

в котором корреляционная матрица  данной векторной случайной величины имеет вид

 (2)

Элементами вектора  являются отклонения частот в каждом из  генераторов от номинальных значений.

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

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

 (3)

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

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

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

Рассмотрим систему из  генераторов, соединенных как показано на рис. 1. Формируемые каждым генератором сигналы подаются не только на функциональные элементы, определяемые назначением устройств и систем телекоммуникаций, но и дополнительно поступают на вход измерительного устройства И.

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

Рис. 1 - Блок-схема из k+1 генераторов


Для каждого из данной совокупности генераторов одновременно, в течение одного и того же интервала измерений длительностью  с использованием измерителя Иопределяется полная фаза колебаний, формируемых этими генераторами. Значения фаз колебаний генераторов определяются соотношениями

. (4)

В соотношении (4)  – частота k-го генератора на интервале измерений длительностью .

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

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


,(5)

где ,…, - составляющие отклонения фаз измеренных колебаний от номинального значения для k-го генератора, обусловленные только нестабильностью временного интервала;  – составляющие отклонения фаз измеренных колебаний от номинального значения для k-го генератора, обусловленные только нестабильностью собственных частот генераторов;  – составляющие отклонений фаз колебаний генераторов, обусловленные нестабильностью как собственных частот генераторов, так и нестабильностью временного интервала одновременно. Однако из-за существенной малости значений последними составляющими можно пренебречь по сравнению с остальными слагаемыми.

С учетом сделанного замечания перепишем выражения (5) в виде

,(6)

где  – номинальное значение фазы k-го генератора;  – отклонение фазы k-го генератора от номинального значения вследствие собственной нестабильности частоты;  – отклонение фазы k-го генератора от номинального значения вследствие нестабильности длительности временного интервала.

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

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

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

.(7)

При этом значения фазы , как отмечалось ранее, являются измеряемыми, а величины  – априорно известными .

Подстановка выражения (4) в (7) позволяет получить следующие равенства


.(8)

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

 - отклонение длительности временного интервала от номинального значения;  - собственные отклонения частот от номинальных значений формируемых колебаний каждого из  генераторов.

Каждая из этих величин является неизвестной, но их сумма, равная , определяется с помощью измерительного устройства И.

Выразим из (8) величины отклонений фаз генераторов, обусловленных только их собственными нестабильностями  в виде:

.(9)

Подстановка уравнений из системы (7) в уравнения из (9) позволяет получить следующую систему уравнений:


. (10)

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

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

Перепишем с учетом сделанного допущения выражения из системы (10) в виде:

.(11)

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

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

.(12)

Запишем формулу для плотности нормального закона распределения величины :

,(13)

где  – относительная нестабильность k-го генератора .

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

. (14)

После подставки выражения (8) в (14) и получим следующее выражение относительно переменной :


.(15)

После логарифмирования данной функции получаем:

,(16)

или с учетом свойства логарифмической функции

. (17)

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

. (18)

В качестве оценки нестабильности временного интервала  берется такое значение параметра , при котором производная (17) обращается в нуль. Приравняем данное выражение к нулю, что позволяет получить:

. (19)

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


. (20)

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

.

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

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

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

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

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

.(21)

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

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


.(22)

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

.(23)

При записи данного равенства было учтено свойство, что

.

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

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

.


Выполняя несложные, но громоздкие преобразования, получаем:

.(24)

С учетом того, что отклонения  и  каждого из  генераторов являются некоррелированными, получаем:

.(25)

Выражение (25) позволяет записать следующее равенство:

.

С учетом этого представление для дисперсии  имеет вид

.(26)

Преобразуем данное выражение при условии, что:

.

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


.(27)

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

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

.(28)

.(29)

Таким образом, увеличение числа генераторов приводит к уменьшению дисперсии оценки нестабильности временного интервала.

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

стабилизация частота генератор

.(30)


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

Анализ рассмотренных случаев показывает, что наиболее высокая стабильность обеспечивается при идентичных значениях частот и относительных нестабильностей (30) в системе генераторов, которая составляет величину, равную .

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

Воспользуемся соотношением (21), на основании которого представим выражение

. (31)

С учетом свойств математического ожидания (8–15) преобразуем зависимость (30) следующим образом

. (32)

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


, (33)

Используем соотношение , которое легко преобразуется к виду . С учетом последней формулы и ранее сделанных предположений о статистических характеристиках отклонений частот генераторов от номинальных значений в различных генераторах зависимость (32) принимает вид

, . (34)

Дисперсии получаемых оценок определяются зависимостью

, .(35)

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

.(36)

Таким образом, соотношения (23), (24) и (34), (36) определяют статистические характеристики оценок отклонения длительности временного интервала измерений и частот генераторов от номинальных значений.

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

Некомпенсированное отклонение частоты k-го генератора  определяется как:

, (37)

или с учетом полученной оценки :

. (38)

В силу того, что измеренное отклонение числа импульсов k-го генератора обусловлено как собственной нестабильностью , так и нестабильностью генератора, задающего временной интервал измерений , запишем данное выражение в следующем виде:

.

Раскроем величины  и , что позволит получить равенство:

.(39)

После несложных преобразований окончательно запишем:


. (40)

Исследуем статистические характеристики нескомпенсированного отклонения частоты k-го генератора (математическое ожидание и дисперсию).

Перепишем выражение (40), раскрыв значение оценки отклонения частотыk-го генератора с использованием соотношения (21):

.(41)

Запишем выражение для нахождения математического ожидания некомпенсированного значения отклонения частоты k-го генератора:

, (42)

где учтено, что

и, следовательно, .

Запишем выражение для нахождения дисперсии нескомпенсированного значения отклонения частоты k-го генератора в виде:


.(43)

После выполнения преобразований получаем

. (44)

Соотношение (44) получено с использованием выражения

.

Такое значение (44) дисперсия нескомпенсированного отклонения частоты k-го генератора принимает лишь в том случае, когда относительные нестабильности  и частоты  идентичны.

Таким образом, из выражения (44) следует, что дисперсия нескомпенсированного значения частоты k-го генератора в  раз меньше дисперсии отклонения частоты k-го генератора при идентичных параметрах используемых генераторов.

Рассмотрим случай, при котором все генераторы работают на различных частотах и имеют одинаковые относительные нестабильности . В силу того, что математическое ожидание нескомпенсированного отклонения частоты k-го генератора, определяемое выражением (43), равно нулю, определим дисперсию нескомпенсированного отклонения частоты k-го генератора для данного случая. Перепишем выражение (43) с учетом заданных исходных данных:

.

После преобразования данного выражения получим:

 (45)

Предположим, что в выражении (45) , и, следовательно, решение данной задачи может быть сведено к частному случаю, рассмотренному выше:

.

Рассмотрим общий случай, при котором все генераторы работают на различных частотах и имеют отличные друг от друга значения относительных нестабильностей . В силу того, что математическое ожидание значения нескомпенсированного отклонения частоты k-го генератора равно нулю, определим дисперсию нескомпенсированного отклонения частоты k-го генератора для данного случая. Перепишем выражение (43) с учетом заданных исходных данных:


.

После проведенных преобразований запишем окончательно:

. (46)

Как несложно убедиться, в частном случае идентичных генераторов  и  решение рассматриваемой задачи сводится к случаю, рассмотренному ранее:

.

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

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

Рассмотрим систему из одиннадцати независимых генераторов , работающих на одинаковых частотах с равными относительными нестабильностями. Для численных расчетов использовались следующие исходные данные: номинальные частоты генераторов Гц, относительная нестабильность . Номинальная длительность временного интервала измерений выбиралась равной . Здесь и далее для анализа основных закономерностей рассматривались отклонения как длительности временного интервала, так и частот генераторов от номинальных значений, значения оценок отклонения длительности временного интервала и частот генераторов, а также значения некомпенсированных отклонений частот генераторов от номинальных значений. Данные результаты приведены в табл. 1.

Таблица 1 - Отклонения частот генераторов с одинаковыми значениями нестабильностей

№ генератора

Отклонение длительности интервала измерений

от номиналь-ного значения

Оценка отклонения длительности интервала измерений

от номиналь-ного значения

Отклонение частоты

генератора

от номиналь-ного значения

Оценка отклонения

частоты генератора

от номиналь-ного

значения

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

частоты генератора

от номиналь-ного значения

1–3,0617·10–9–3,3566·10–9–20,2894
Актуально: