WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     || 2 | 3 |
-- [ Страница 1 ] --

Ярославский государственный университет им. П.Г. Демидова

На правах рукописи

Якимов Иван Михайлович ИССЛЕДОВАНИЕ И ОПТИМИЗАЦИЯ ЦЕПИ ПОСЛЕДОВАТЕЛЬНО СИНХРОНИЗИРУЕМЫХ ГЕНЕРАТОРОВ В УСЛОВИЯХ

КОМБИНИРОВАННЫХ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЙ Специальность: 05.12.04 – Радиотехника, в том числе системы и устройства радионавигации, радиолокации и телевидения Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель доктор технических наук Л.Н. Казаков Ярославль 2004 г.

-2ОГЛАВЛЕНИЕ стр. Введение _ 6 Глава 1. Математическое описание цепи последовательно синхронизируемых генераторов в условиях комбинированных случайных воздействий 1.1. Постановка задачи _ 1.2. Математические модели случайных воздействий _ 1.3. Математическая модель звена цепи на основе дискретной СФС в условиях комбинированных случайных воздействий _ 1.3.1. Модель звена в форме стохастических разностных уравнений 1.3.2. Модель звена в форме расширенного векторного уравнения Колмогорова–Чепмена 1.4. Модель цепи последовательно синхронизируемых генераторов в форме системы расширенных векторных уравнений Колмогорова–Чепмена и уравнений преобразования координат 1.5. Линеаризованная модель цепи последовательно синхронизируемых генераторов _ 1.5.1. Линеаризованная дискретная модель 1.5.2. Линеаризованная аналоговая модель 1.6. Выводы Глава 2. Анализ и оптимизация статистических характеристик звена цепи на основе дискретной СФС 2-го порядка в условиях комбинированных случайных воздействий 2.1. Постановка задачи _ 2.2. Оптимизация статистических характеристик в линейном приближении _ 2.3. Анализ статистических характеристик с помощью численного решения расширенного векторного уравнения Колмогорова– Чепмена _ 16 16 18 20 21 30 35 35 41 46 46 -32.3.1. Характеристики фазовой ошибки и фазовых флуктуаций выходного сигнала звена на основе бесфильтровой СФС _ 2.3.2. Анализ фазовой ошибки звена 2-го порядка 2.3.3. Характеристики фазовых флуктуаций на выходе звена 2-го порядка 2.4. Анализ статистических характеристик в линейном приближении _ 2.4.1. Применение аналоговой модели 2.4.2. Применение дискретной модели 2.5. Выводы Глава 3. Анализ и оптимизация статистических характеристик цепи последовательно синхронизируемых генераторов 3.1. Постановка задачи _ 3.2. Анализ и оптимизация цепи последовательно синхронизируемых генераторов с помощью аппарата марковских процессов 3.3. Статистические характеристики цепи последовательно синхронизируемых генераторов в линейном приближении для случая идентичных аналоговых звеньев _ 3.3.1. Случай бесфильтровых колец синхронизации _ 3.3.2. Случай колец синхронизации 2-го порядка _ 3.4. Статистические характеристики цепи последовательно синхронизируемых генераторов в линейном приближении для индивидуальной настройки аналоговых звеньев 3.4.1. Случай бесфильтровых колец синхронизации _ 3.4.2. Случай колец синхронизации 2-го порядка _ 3.5. Статистические характеристики цепи последовательно синхронизируемых генераторов в линейном приближении для случая идентичных дискретных звеньев 3.5.1. Случай бесфильтровых колец синхронизации _ 3.5.2. Случай колец синхронизации 2-го порядка _ 54 63 67 71 72 84 95 105 106 126 126 136 136 -43.6. Статистические характеристики цепи последовательно синхронизируемых генераторов в линейном приближении для индивидуальной настройки дискретных звеньев _ 3.6.1. Случай бесфильтровых колец синхронизации _ 3.6.2. Случай колец синхронизации 2-го порядка _ 3.7. Выводы Глава 4. Разработка и исследование имитационной модели цепи последовательно синхронизируемых генераторов на основе дискретных СФС 4.1. Постановка задачи _ 4.2. Описание структурной схемы имитационной модели методики проведения исследований 4.3. Разработка статистических характеристик сигналов на выходе цепи произвольной длины _ 4.4. Исследование имитационной модели для различных случайных воздействий 4.4.1. Случай комбинированного случайного воздействия со спектральной плотностью 1-го порядка _ 4.4.2. Случай комбинированного случайного воздействия со спектральной плотностью 2-го и 3-го порядков 4.4.3. Исследование телекоммуникационных характеристик качества сигнала на выходе системы _ 4.5. Выводы и сравнительный анализ результатов, полученных различными методами _ Заключение Библиографический список _ 144 145 150 154 154 163 168 168 172 177 183 185 Приложение 1. Текст программы анализа характеристик качества сигнала на выходе цепи последовательно синхронизируемых генераторов в линейном приближении -5Приложение 2. Текст программ обработки результатов имитационного моделирования цепи последовательно синхронизируемых генераторов -6ВВЕДЕНИЕ Актуальность работы Задача взаимодействия двух и более генераторов для области радиотехники является достаточно традиционной [46] и связана, как правило, с обеспечением синхронных режимов нескольких колебательных процессов. Примером является задача о стабилизации частоты, в которой менее стабильный генератор синхронизируется от более стабильного [1–6, 11–14]. Сюда же можно отнести современные системы частотного синтеза, строящиеся по принципу последовательно–параллельной синхронизации генераторов на кратных частотах [6, 9, 10, 34, 35, 43–45]. В подобных структурах число связанных генераторов может быть достаточно большим. В последние годы с развитием синхронных цифровых систем передачи информации цепочки последовательно соединенных генераторов приобрели новые области применения. В некоторых случаях они используются для синхронизации работы разнесенных в пространстве радиотехнических систем. Так, например, в сети передачи данных синхронной цифровой иерархии (СЦИ) используется цепь синхронизации, включающая в себя десятки последовательно синхронизируемых генераторов [90, 102]. В связи с широким распространением и важностью цифровых сетей передачи данных, исследования, посвященные улучшению качества их работы, представляют большой научный и практический интерес. В том числе важным является вопрос об улучшении качества сигнала синхронизации. Для решения данного вопроса необходимы исследования работы цепочки последовательно соединенных дискретных систем фазовой синхронизации. В основном существующие алгоритмы улучшения качества сигнала синхронизации используют специальные метки качества сигнала, передаваемые по специальным служебным каналам. На их основе производится выбор наилучшего из имеющихся сигналов синхронизации. Подобные алгоритмы на сегодняшний день разработаны достаточно хорошо и позволяют создать надежную сеть синхронизации [107–109, 110–113]. Но у подобного подхода существует и ряд существенных недостатков. Вопервых, он не учитывает реального качества сигналов синхронизации. Вся информация о том, насколько этот сигнал хорош, заключена в метке. Вовторых, этот метод не позволяет сказать, как можно улучшить качество сигнала -7синхронизации, как настроить связующие звенья, чтобы на выходе цепочки был сигнал наилучшего качества. Решение этих задач весьма важно для проектирования и управления цепью синхронизации. В то же время для анализа поведения цепи необходимо более детально изучить работу отдельных звеньев, в роли которых выступают дискретные системы фазовой синхронизации [102]. Системы фазовой синхронизации с элементами дискретизации интенсивно исследуются в последние годы [54, 67– 69, 81]. Такие системы обладают рядом преимуществ по сравнению с чисто аналоговыми устройствами, такими как повышенная помехоустойчивость, простота реализации, малая потребляемая мощность. Указанные достоинства совместно с неуклонным ростом рабочих частот цифровой схемотехники позволяют строить высокоэффективные системы обработки информации, включающие в себя системы фазовой синхронизации [7–10, 15, 27, 31]. Большое количество областей применения данных устройств определяет интерес к исследованиям данного класса систем [27, 32–36, 43–45, 47–49]. Значительный интерес вызывают исследования, посвященные поведению этих устройств в условиях помеховых воздействий [19]. Данный факт обусловлен постоянным ростом требований на качество обработки информации со стороны пользователей. С другой стороны помеховые воздействия всегда имеют место. Зачастую ими нельзя пренебречь, и они определяют качество работы системы в целом. Поэтому анализ реакции на действие помех достаточно важен для практики [1–5, 12–14, 41, 50, 51, 53]. Основная часть работ в данной области посвящена исследованию влияния наиболее простых и распространенных помех в виде аддитивного широкополосного гауссовского шума [16, 17, 22–25, 88]. Однако не всегда этот вид помех является определяющим. В связи с этим последнее время ряд авторов проводит исследования, посвященные анализу поведения систем фазовой синхронизации при наличии на входе помех различного вида. Так в ряде работ изучается влияние детерминированных и случайных сосредоточенных по частоте помех на работу указанных устройств, в том числе и помех, по структуре повторяющих полезный сигнал [15–18, 52, 56, 57, 60, 61, 63–65]. Учет данного типа помех позволяет ответить на многие вопросы функционирования современных систем передачи информации в условиях сложной электромагнитной обстановки. Под данный тип помех подходят помехи других -8станций, работающих на близкой частоте, что характерно для сотовых систем, условий многолучевого распространения сигналов и т.д. В то же время малоизученным или совсем неизученным остается поведение СФС при многих типах внешних воздействий, которые при некоторых условиях являются определяющими для качества работы систем фазовой синхронизации. В частности в случае анализа цепи последовательно синхронизируемых генераторов к таким воздействиям относятся фазовые шумы, присутствующие в различных сигналах. Их источники могут быть самыми различными [103, 104]. В некоторых случаях данные шумы являются основными, оказывающими влияние на функционирование системы. В ряде работ [90, 97, 102] показано, что на качество работы цифровых сетей передачи информации оказывают существенное влияние именно фазовые флуктуации сигналов. В связи с этим изучение данного вида воздействий является весьма актуальным. Однако круг работ, посвященных данной проблематике, весьма ограничен. К ним следует отнести работы [100, 102] а так же работы автора диссертации. Для анализа поведения и статистических характеристик систем при наличии сложных воздействия необходимо обладать определенным математическим аппаратом. Исследование работы систем фазовой синхронизации при наличии случайных воздействий ведется уже достаточно долго. За это время был разработан ряд подходов и методов к анализу данных устройств. Однако на этом пути существуют серьезные сложности. Данный факт связан с тем, что система фазовой синхронизации является сугубо нелинейной системой. Аналитические методики исследования таких устройств в большинстве случаев носят приближенный характер. В качестве примера таких подходов можно указать различные методы линеаризации и усреднения [20, 44]. На сегодняшний день одним из самых прогрессивных механизмов, позволяющих точно исследовать динамику нелинейных систем, является аппарат марковских процессов. Данный аппарат позволяет получить многие важные характеристики стохастических систем, такие как плотность распределения вероятности координат, среднее время достижения синхронизма, среднее время до срыва синхронизма. Ввиду явных достоинств этого метода его применению к анализу различных систем посвящено достаточно много работ. К их числу относятся труды Тихонова В.И., Миронова М.А., Казакова В.А., Стратоновича Р.Л. [11, 21, 28]. Применением данного -9метода к системам фазовой синхронизации, в том числе дискретным, занимались Шахтарин Б.И., Витерби А., Разевиг В.Д., Казаков Л.Н. и другие авторы [17, 25, 26, 56–58]. Ими наработаны основные методики и подходы, позволяющие применить аппарат марковских процессов к системам синхронизации, достаточно подробно проведен анализ поведения систем фазовой синхронизации в условиях аддитивных широкополосных шумов. В работах Башмакова М.В. [56–59, 61–65] рассмотрены статистические характеристики СФС при наличии детерминированных или случайных узкополосных аддитивных помех. В то же время очень мало работ посвящено анализу систем фазовой синхронизации в условиях присутствия фазовых шумов. Данный вопрос только начинает изучаться. Кроме того, аппарат марковских процессов в том виде, в котором он применяется сегодня, позволяет получить характеристики фазовой ошибки. При анализе же цепочки последовательно соединенных дискретных систем фазовой синхронизации интерес представляют характеристики фазы выходного сигнала перестраиваемого генератора. На сегодняшний день подходы к решению данной проблемы не известны. Отдельным вопросом является проблема оптимизации параметров систем и устройств при определенных входных воздействиях. В современных условиях постоянного роста требований на качество обработки сигналов эта проблема зачастую выходит на первое место. В то же время существует ряд проблем, связанных с тем, что системы фазовой синхронизации являются нелинейными устройствами. В случае линейных систем с целью построения оптимальных устройств широко используются винеровская и калмановская теории фильтрации. Существуют попытки применить данную теорию и к синтезу оптимальных систем фазовой синхронизации [115–117], но данные работы носят единичный характер. На основании вышесказанного тема диссертации, посвященная исследованию статистических характеристик цепи последовательно синхронизированных с помощью дискретных систем фазовой синхронизации генераторов и отдельных ее звеньев в условиях комбинированных случайных аддитивных и фазовых воздействий и оптимизации цепи, обеспечивающей качественный выходной сигнал, является актуальной.

- 10 Цели и задачи диссертации Целью диссертационной работы является моделирование, исследование и оптимизация цепи последовательно синхронизируемых генераторов на основе дискретных систем фазовой синхронизации в условиях комбинированных случайных аддитивных и фазовых воздействий. Для достижения поставленной цели в диссертации решаются следующие основные задачи: 1. Разработка математической модели цепи последовательно синхронизируемых генераторов, а так же отдельных ее звеньев в условиях комбинированных случайных воздействий. 2. Разработка методики анализа статистических характеристик сигналов в системе в условиях комбинированных случайных воздействий. 3. Исследование и параметрическая оптимизация отдельных звеньев цепи, представляющих собой системы фазовой синхронизации, с учетом специфики флуктуационных воздействий. 4. Исследование и параметрическая оптимизация цепи последовательно синхронизируемых генераторов с учетом специфики флуктуационных воздействий. 5. Построение имитационной модели цепи последовательно синхронизируемых генераторов с использованием компьютерных пакетов динамического моделирования. 6. Разработка методики и исследование имитационной модели цепи последовательно синхронизируемых генераторов и сравнение результатов с результатами анализа математической модели. Общая методика исследований Разрабатываемые в диссертации методы исследования цепочки последовательно соединенных дискретных систем фазовой синхронизации основаны на общих положениях качественных методов теории дискретных систем, теории нелинейных разностных уравнений, на аппарате теории вероятности и в частности марковских процессов и цепей, на прикладной теории математической статистики, на статистической теории радиотехнических систем и устройств.

- 11 Для решения поставленных задач используются также компьютерное моделирование, численное решение нелинейных стохастических разностных уравнений. Разработанные методы и алгоритмы анализа статистических характеристик дискретных СФС и цепочки из дискретных СФС ориентированы на использование персональных компьютеров. Научная новизна результатов 1. Построены математические модели дискретных систем фазовой синхронизации при наличии комбинированных случайных воздействий в форме стохастических уравнений и векторных уравнений Колмогорова– Чепмена. 2. Получены математические модели цепи последовательно синхронизируемых генераторов на основе дискретных СФС в виде комбинации векторных уравнений Колмогорова–Чепмена и уравнений перехода к плотности распределения вероятности выходных координат. 3. Разработана методика анализа фазовых флуктуаций на выходе дискретной системы фазовой синхронизации, основанная на расширенной марковской модели. 4. С помощью аппарата марковских процессов проведено исследование статистических характеристик фазовой ошибки и фазовых флуктуаций на выходе дискретной СФС 2-го порядка в условиях близких к белым частотных шумов входного сигнала и сигнала перестраиваемого генератора и белого аддитивного канального шума. Получены зависимости дисперсии фазовых ошибок и выходных фазовых флуктуаций от различных параметров воздействий и системы. 5. Проведено исследование и оптимизация статистических характеристик фазовых флуктуаций сигнала на выходе цепи последовательно синхронизируемых генераторов на основе систем фазовой синхронизации. Рассмотрено несколько алгоритмов оптимизации параметров звеньев и проведено их сравнение. Практическая ценность 1. В диссертации предложены методики исследования, позволяющие определить основные статистические характеристики систем фазовой - 12 синхронизации, а так же цепи последовательно соединенных СФС в условиях комбинированных флуктуационных воздействий. На основе методик разработаны алгоритмы для расчета статистических характеристик отдельных звеньев и цепи в целом, в том числе телекоммуникационных характеристик качества TDEV, TIE. 2. Разработанные программы позволяют оптимизировать параметры систем фазовой синхронизации с целью обеспечения заданных статистических свойств сигналов в условиях комбинированных воздействий и параметры цепи последовательно соединенных СФС. 3. Полученные в диссертации результаты позволили сформулировать предложения по повышению эффективности работы цепи последовательно синхронизируемых генераторов, функционирующей в условиях сложных флуктуационных воздействий на систему. 4. Предложенные и развитые в диссертации методики и разработанные на их основе алгоритмы и программы можно использовать в научно– исследовательских и опытно–конструкторских работах для анализа статистических свойств систем фазовой синхронизации и синтеза СФС различного назначения в условиях комбинированных случайных воздействий. Методики и алгоритмы могут быть использованы при разработке и исследовании цепей последовательно синхронизируемых генераторов различного назначения. Часть материалов, включая разработанное программное обеспечение, используется в учебном процессе ЯрГУ г. Ярославль, МГТУ им. Н.Э. Баумана г. Москва. Положения, выносимые на защиту 1. Математическая модель дискретной СФС в форме векторного уравнения Колмогорова–Чепмена для случая комбинированного случайного воздействия, представляющего собой аддитивный широкополосный гауссовский шум, фазовые флуктуации входного сигнала и фазовые флуктуации сигнала перестраиваемого генератора. 2. Математическая модель цепи последовательно синхронизируемых с помощью дискретных СФС генераторов в виде комбинации векторных уравнений Колмогорова–Чепмена отдельных звеньев и уравнений перехода к плотности распределения вероятности выходных координат.

- 13 3. Методика анализа фазовых флуктуаций на выходе дискретной системы фазовой синхронизации, основанная на применении расширенной марковской модели. 4. Результаты исследования и оптимизации однокольцевой дискретной СФС 2-го порядка в условиях комбинированных случайных воздействий с полиномиальной спектральной плотностью. 5. Результаты исследования и оптимизации цепи последовательно синхронизируемых генераторов одного и различных типов для трех вариантов оптимизации: гомогенной, позвенной, смешанной. 6. Имитационная модель цепи последовательно синхронизируемых генераторов разного типа на базе дискретных СФС, выполненная в среде динамического моделирования Simulink пакета Matlab. 7. Методика и результаты исследования имитационной модели цепи последовательно синхронизируемых генераторов разного типа. Объем и структура диссертации Диссертация состоит из введения, четырех глав, заключения, списка литературы и приложений. Во введении обоснована актуальность темы и ее практическая значимость, сформулированы цели и задачи исследования, дан критический анализ работ в области исследования цепей последовательно соединенных генераторов, а так же различных классов систем фазовой синхронизации. В первой главе предложены линейные и нелинейные математические модели цепи последовательно синхронизируемых генераторов на основе систем фазовой синхронизации и модели звеньев этой цепи в условиях комбинированных случайных воздействий на систему. Нелинейные модели представлены в виде разностных стохастических уравнений и уравнений Колмогорова–Чепмена. Предложена методика исследования характеристик фазовых флуктуаций сигнала на выходе цепи СФС на основе аппарата марковских процессов. Линейные модели получены в виде коэффициентов передачи и квадратов их АЧХ для различных источников флуктуационных воздействий в системе. Вторая глава посвящена анализу и оптимизации систем фазовой синхронизации в условиях комбинированного случайного воздействия.

- 14 Последнее представляет собой аддитивную смесь широкополосного шума и псевдогармонического сигнала, фаза которого флуктуирует по сложному закону. Кроме того, учитывается, что сигнал перестраиваемого генератора СФС так же имеет свои флуктуации фазы. В главе рассмотрена аналитическая оптимизация параметров аналоговых и дискретных систем фазовой синхронизации для различных типов входных воздействий. В качестве критерия оптимальности использовались минимумы дисперсий фазовой ошибки и фазовых флуктуаций на выходе СФС. С помощью аппарата марковских процессов исследованы характеристики фазовой ошибки и фазовых флуктуаций выходного сигнала дискретных систем фазовой синхронизации для случая белых частотных шумов сигнала на входе и сигнала ПГ. Получены зависимости оптимального значения параметров СФС от характеристик входных воздействий. В главе рассмотрено поведение различных линеаризованных моделей систем при наличии сложных флуктуационных воздействий. Исследована зависимость дисперсии фазовых флуктуаций на выходе СФС от параметров входных воздействий, а так же от параметров самих систем. Рассмотрены вопросы оптимизации параметров СФС с целью минимизации дисперсии фазового шума выходного сигнала. Проведено сравнение поведения аналоговых и дискретных СФС при наличии комбинированного входного воздействия. Глава 3 посвящена анализу поведения цепи последовательно синхронизируемых генераторов на основе систем фазовой синхронизации. Рассмотрена зависимость дисперсии фазовых флуктуаций на выходе цепи от ее длины, от параметров звеньев и шумовых воздействий. Получены ограничения на параметры звеньев, которые можно использовать в цепи. Изучено поведение системы, состоящей из аналоговых и дискретных СФС. Проведено исследование различных алгоритмов оптимизации параметров звеньев цепочки, а так же сравнение качества работы этих методов. Рассмотрены особенности поведения гомогенной и негомогенной цепочек последовательно синхронизируемых генераторов. Четвертая глава посвящена разработке и анализу имитационной модели цепи последовательно синхронизируемых генераторов, а так же сравнению результатов исследования имитационной и математической моделей. В рамках пакетов динамического моделирования SystemView и Simulink построены модели системы. Изучено и представлено несколько методов формирования - 15 случайных процессов, близких по своим свойствам к фликкер–шумам. Для формирования во временной области использовалась система нелинейных дифференциальных уравнений Ланжевена, для формирования в спектральной области использовался формирующий фильтр нижних частот 8-го порядка. На их основе созданы имитационные модели генераторов флуктуационных сигналов с полиномиальной спектральной плотностью. Разработана методика исследования статистических характеристик сигналов на выходе системы по их временным реализациям, полученным в результате моделирования. Проведен анализ дисперсии фазовых флуктуаций выходного сигнала отдельных звеньев и всей цепи в целом. Путем перехода к исследованию дисперсии частотных флуктуаций выходного сигнала решена проблема анализа статистических характеристик при неэргодичности процесса флуктуаций фазы. Получены зависимости телекоммуникационных параметров качества работы системы TDEV и TIE от времени наблюдения и показано влияние на эти характеристики основных источников шумов в цепи. Проведено сравнение результатов, полученных с помощью имитационной модели, с результатами анализа математических моделей. В заключении подведены итоги диссертации и показаны направления дальнейшего развития идей, предложенных в работе.

- 16 Глава 1. Математическое описание цепи последовательно синхронизируемых генераторов в условиях комбинированных случайных воздействий 1.1. Постановка задачи Цепь последовательно синхронизируемых генераторов является важной частью многих синхронных цифровых сетей передачи информации. Наиболее ярким примером таких сетей является система передачи синхронной цифровой иерархии (СЦИ). Роль цепи последовательно синхронизируемых генераторов в данной системе весьма важна, она обеспечивает систему сигналом синхронизации, необходимым для надежного функционирования [90]. В то же время на сегодняшний день цепь последовательно синхронизируемых генераторов изучена недостаточно. Частично эту ситуацию можно объяснить отсутствием адекватных моделей, описывающих ее функционирование при произвольных воздействиях. Данному вопросу посвящено ограниченное число работ, рассматривающих частные задачи [92, 100, 102]. Еще меньше освещены проблемы выбора параметров отдельных звеньев, обеспечивающих наилучшее качество работы системы в целом. Одной из задач, которые необходимо решить для ответа на поставленные вопросы, является изучение поведения цепи и отдельных ее звеньев в условиях комбинированных случайных воздействий. Как известно, последовательно соединенные генераторы можно синхронизировать различными способами. Однако в современных телекоммуникационных системах широкое распространение получил способ, когда соединение генераторов происходит средствами систем фазовой синхронизации (СФС) [102, 106]. Изучение поведения и характеристик дискретных систем фазовой синхронизации ведется уже достаточно давно. В рамках данных исследований были построены различные модели данных систем. Широкое распространение получили модели в виде систем разностных уравнений [10, 27, 32–34], позволяющих описать широкий класс устройств, к которому относятся как ряд цифровых СФС, так и импульсные системы, имеющие в себе аналоговые звенья. Возможность описания с помощью данных моделей такого обширного круга устройств обуславливает популярность этого подхода.

- 17 Рост интереса к статистической динамике систем фазовой синхронизации привел к появлению моделей в виде стохастических разностных уравнений. На ранних этапах исследований в основном учитывалось влияние аддитивного широкополосного шума, подмешиваемого к сигналу в канале передачи. В ряде работ показано, что такой шум может быть эффективно пересчитан на выход фазового детектора [37, 39, 40]. В дальнейшем шумовое воздействие, оказывающее влияние на систему, усложнялось под влиянием требований к учету различных факторов. Так в некоторых работах поднимается вопрос моделирования систем фазовой синхронизации при наличии сосредоточенных по частоте помех [52, 56, 57, 60, 61]. В то же время при моделировании цепочки дискретных СФС следует учитывать ряд специфических помех. В частности существенную роль играют флуктуации фазы сигналов присутствующих в системе генераторов. Данные воздействия учитывались лишь в малом числе работ [86]. Кроме того, сами по себе данные воздействия имеют достаточно сложную структуру, что требует построения математических моделей их самих. С развитием применения аппарата марковских процессов к анализу систем фазовой синхронизации распространение получили марковские модели систем. Несмотря на то, что указанный аппарат является весьма удобным, его применение возможно лишь в том случае, если удается описать исследуемую систему таким образом, чтобы процесс изменения ее состояния являлся марковским. Это удается не всегда. Применению метода марковских процессов посвящено достаточно много работ, в том числе и работ автора диссертации [25, 28, 37, 57, 89]. Однако, в силу некоторых трудностей при построении марковских моделей, для некоторых входных воздействий данный аппарат пока не применялся. В частности фазовое воздействие на систему учитывается в малом количестве работ. Кроме того, в основном аппарат марковских процессов применялся для анализа статистических характеристик фазовой ошибки в дискретной СФС. При исследовании же цепочки последовательно соединенных систем фазовой синхронизации интерес представляют статистические характеристики флуктуаций фазы перестраиваемого генератора. Это требует внесения определенных изменений в методику применения аппарата марковских процессов. Как уже говорилось, возможности применимости аппарата процессов Маркова ограничены. Не для любых входных воздействий можно представить - 18 систему фазовой синхронизации в виде марковской модели. В связи с этим требуются иные формы представления моделей дискретных СФС. Кроме всего прочего, данные модели должны быть удобны для использования их при анализе цепочки. Из курса радиотехнических систем известно [79, 80], что цепочечные структуры удобно описывать, если известны коэффициенты передачи отдельных звеньев. В силу того, что СФС является нелинейной системой, в общем случае ее нельзя описать с помощью коэффициента передачи. Однако при некоторых допущениях это возможно. Так в ряде работ [77, 82, 85] автором диссертации предложена и апробирована именно линейная модель цепи. В связи с этим задачей данной главы является построение математической модели цепи последовательно синхронизируемых генераторов, связанных с использованием дискретных систем фазовой синхронизации. Необходимо получить модели звеньев цепи в форме стохастических уравнений и марковских моделей при наличии комплексного случайного воздействия, а так же в форме коэффициентов передачи и квадратов их АЧХ, необходимых при анализе прохождения шумов через систему. Под комплексным случайным воздействием понимается следующее. Входной сигнал представляет собой сумму широкополосного случайного процесса и информационного сигнала, чья фаза флуктуирует по некоторому закону. Кроме того, в данной главе в математическом описании систем учитываются фазовые флуктуации перестраиваемого генератора СФС. 1.2. Математические модели случайных воздействий Пусть входной сигнал системы фазовой синхронизации представим в виде:

s вх (t ) = A H ( t + + (t )) + n(t ), (1.2.1) где A, и – постоянные амплитуда, частота и начальная фаза сигнала, n(t) – широкополосный случайный процесс, H(·) – периодическая функция своего аргумента с периодом 2, определяющая форму входного сигнала, (t) представляет собой флуктуации фазы входного сигнала.

- 19 Пусть сигнал перестраиваемого синхронизации представим в виде: генератора системы фазовой s пг (t ) = U G ( пг 0 t + пг (t ) + пг (t )), где U, пг0 – амплитуда и собственная частота перестраиваемого генератора, G(·) – периодическая функция своего аргумента с периодом 2, определяющая форму сигнала перестраиваемого генератора, пг(t) – добавка к фазе ПГ, обусловленная каналом управления, пг(t) – собственный фазовый шум перестраиваемого генератора. Функции H и G могут иметь различный вид в зависимости от применения систем. Наиболее распространены следующие виды сигналов: sin, cos, rect. С ними могут применяться различные типы фазовых детекторов, однако, в данной главе рассматривается общий случай, и вид фазового детектора не конкретизируется. Относительно флуктуационных процессов n(t), (t) и пг(t) в работе используются следующие предположения. Случайный процесс n(t) представляет собой широкополосный шум с постоянной спектральной плотностью: Sn() = c0. (1.2.2) Случайные процессы (t) и пг(t) представляют собой собственные флуктуации фазы сигналов на входе СФС и перестраиваемого генератора. В качестве модели фазовых шумов генераторов, используемых в цепочке, применяется полиномиальная модель. В работах [103, 104], посвященных анализу статистических характеристик систем фазовой синхронизации, показано, что фазовые шумы перестраиваемых генераторов описываются моделью Лисона. При этом их спектральная плотность имеет вид:

2 2 2 K n kT f f ( + + + 1), S ( ) = Ps 3 - 20 где k – постоянная Больцмана, Т – температура в Кельвинах, Kn – коэффициент шума, PS – мощность сигнала, – частота, на которой уровень фликкер–шума достигает уровня белого шума, f – половина полосы пропускания цепи обратной связи по уровню половинной мощности. Данная модель учитывает виды шумов, перечисленные в таблице 1.1. Таблица 1.1 Слагаемое, содержащее отвечает за белый фазовый шум фликкер фазовый шум белый частотный шум фликкер частотный шум 0 –1 –2 – Из модели Лисона видно, что спектральная плотность фазовых флуктуаций определяется тремя независимыми параметрами. Ниже в работе используется более общая модель спектральной плотности фазовых флуктуаций сигнала генератора [102]:

a a a S ( ) = a0 + 1 + 2 + 3.

(1.2.3) В отличие от модели Лисона, эта модель позволяет отдельно задавать уровень каждой составляющей спектральной плотности. 1.3. Математическая модель звена цепи на основе дискретной СФС в условиях комбинированных случайных воздействий Описание дискретной СФС в виде стохастических разностных уравнений основано на ее функциональной схеме, приведенной на рис. 1.1. Особенностью схемы в отличии от известных [16, 37] является учет фазовых флуктуаций входного сигнала и сигнала перестраиваемого генератора. На схеме приняты следующие обозначения: [k] – фазовая ошибка, представляющая собой разность полной фазы входного сигнала и полной фазы сигнала перестраиваемого генератора, n[k] – пересчитанный на выход фазового - 21 детектора аддитивный широкополосный шум, [k] – случайная составляющая фазы входного сигнала, пг[k] – случайная составляющая фазы перестраиваемого генератора, F(·) – 2–периодическая функция, описывающая дискриминационную характеристику фазового детектора, K(z) – передаточная функция ФНЧ. Величина вх[k] представляет собой детерминированную составляющую фазы входного сигнала:

вх [k ] = T k +, где T – период дискретизации.

вх[k] [k] + T / (z–1) + – [k] E·F( [k]) n[k] + K(z) + S пг[k] пг Рис. 1.1. Функциональная схема дискретной системы фазовой синхронизации при наличии флуктуаций фазы входного сигнала и сигнала перестраиваемого генератора, а так же широкополосного аддитивного воздействия. 1.3.1. Модель звена в форме стохастических разностных уравнений В терминах фазовой ошибки стохастическое разностное уравнение данной системы в Z–области имеет вид [37]:

( z 1) ( z ) = ( z 1) вх ( z ) + ( z 1) ( z ) ( z 1) пг ( z ) T у K ( z ) F ( z ) T S K ( z ) n( z ) T пг 0 ( z ), (1.3.1) где у = E·S. Переход в область дискретного времени с помощью обратного Z– преобразования позволяет получить из уравнения (1.3.1) стохастическое - 22 разностное уравнение системы фазовой синхронизации, записанное в терминах фазовой ошибки. Для определения флуктуационной составляющей фазы выходного сигнала используется следующий подход. Полная фаза выходного сигнала имеет вид:

пг(z) = вх пол(z) – (z).

где вх пол(z) = вх(z) + (z) – полная фаза входного сигнала. В соответствии с этим запишем выражение для флуктуационной составляющей фазы выходного сигнала:

пг(z) = (z) – (z).

(1.3.2) Выражение (1.3.2) послужит в дальнейшем основой для расчета ПРВ фазовых и частотных флуктуаций выходного сигнала звена при заданной совместной ПРВ фазовых флуктуаций входного сигнала и фазовой ошибки. Рассмотрим конкретные случаи уравнения (1.3.1). 1. Пусть система является бесфильтровой: K(z) = 1. В этом случае уравнение (1.3.1) принимает вид:

( z 1) ( z ) = ( z 1) вх ( z ) + ( z 1) ( z ) ( z 1) пг ( z ) T у F ( z ) T S n( z ) T пг 0 ( z ).

(1.3.3) Проделав уравнение:

обратное Z–преобразование [73], получим следующее [k + 1] [k ] = вх [k + 1] вх [k ] + [k + 1] [k ] пг [k + 1] + пг [k ]. T у F ([k ]) T S n[k ] T пг (1.3.4) Используя выражение для детерминированной составляющей фазы входного сигнала вх, получим из (1.3.4):

[k + 1] [k ] = T н + [k + 1] [k ] пг [k + 1] + пг [k ] T у F ([k ]) T S n[k ], (1.3.5) - 23 где н = – пг0 – начальная частотная расстройка. Необходимо отметить следующую особенность уравнения (1.3.1). Фазовые шумы и пг сигналов в системе вне зависимости от типа используемого фильтра входят в уравнение в виде (z–1)·(z). В ряде случаев бывает удобно ввести в рассмотрение новый случайный процесс (z) = (z–1)·(z). По физическому смыслу этот процесс представляет собой частотные флуктуации. Во временной области соотношение между и имеет вид:

[k] = [k+1] – [k].

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

( z 1) ( z ) = ( z 1) вх ( z ) + ( z ) пг ( z ) T у K ( z ) F ( z ) T S K ( z ) n( z ) T пг 0 ( z ).

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

( z 1) ( z ) = ( z ).

(1.3.7) Стохастическое разностное уравнение бесфильтровой дискретной системы фазовой синхронизации (1.3.5) в данном случае записывается в виде:

[k + 1] [k ] = T н + [k ] пг [k ] T у F ([k ]) T S n[k ]. (1.3.8) [k + 1] = [k ] + [k ] Система уравнений (1.3.8) в случае слабокоррелированных случайных процессов [k], пг[k] и n[k] (значения процессов в соседние моменты времени независимы) описывает двумерную марковскую случайную - 24 последовательность. Саму систему (1.3.8) с учетом новой координаты в дальнейшем будем рассматривать в качестве расширенной марковской модели. Как видно, такой подход позволяет несколько упростить запись уравнения. Кроме того, он позволяет выразить текущий отсчет сигнала фазовой ошибки только через значения сигналов в предыдущие моменты времени. Данный факт зачастую является немаловажным, особенно при построении марковских моделей систем, что будет показано позднее. 2. Рассмотрим уравнение (1.3.1) с K(z), имеющим следующий вид:

a z + a0 K ( z) = 1, zd где 0 d 1. Подставив выражение для передаточной функции в (1.3.6), получим:

( z d ) ( z 1) ( z ) = ( z d ) ( z 1) вх ( z ) + ( z d ) ( z ) ( z d ) пг ( z ) T у (a1 z + a0 ) F ( z ) T S (a1 z + a0 ) n( z ) ( z d ) T пг 0 ( z ).

(1.3.9) После перехода во временную область имеем стохастическое уравнение 2го порядка следующего вида:

[k + 2] (1 + d ) [k + 1] + d [k ] = (1 d ) T н + ([k + 1] пг [k + 1]) d ([k ] пг [k ]). T y a1 F ([k + 1]) T y a0 F ([k ]) T S a1 n[k + 1] T S a0 n[k ] (1.3.10) С учетом уравнения для входного фазового шума в итоге получаем следующую систему стохастических уравнений:

[k + 2] (1 + d ) [k + 1] + d [k ] = (1 d ) T н + ([k + 1] пг [k + 1]) d ([k ] пг [k ]). T y a1 F ([k + 1]) T y a0 F ([k ]) T S a1 n[k + 1] T S a0 n[k ] [k + 1] = [k ] + [k ] (1.3.11) - 25 Случайный векторный процесс с компонентами ([k], [k]) не является марковским. Для перехода к марковской дополнительную координату. модели необходимо ввести 1.3.2. Модель звена в форме расширенного векторного уравнения Колмогорова–Чепмена Целью данного раздела является построение марковских моделей рассматриваемых дискретных систем фазовой синхронизации в форме уравнений Колмогорова–Чепмена. Для этого необходимо получить выражения для условных ПРВ на основе соответствующих систем стохастических разностных уравнений, и преобразовать их с учетом существования инвариантных движений в системе, обусловленных периодической нелинейностью [11, 14, 37]. Кроме того, необходимо разработать методику, позволяющую определить плотность распределения вероятности флуктуаций фазы на выходе СФС по результатам решения уравнения Колмогорова– Чепмена. Рассмотрим бесфильтровую дискретную систему фазовой синхронизации. Переписав (1.3.8) в виде:

[k + 1] = [k ] + T н + [k ] пг [k ] T у F ([k ]) T S n[k ], [k + 1] = [k ] + [k ] (1.3.12) получим марковскую модель системы в виде стохастического отображения. Динамика марковской модели описывается с помощью уравнения Колмогорова–Чепмена, векторная форма которого представлена следующим уравнением [28]:

Wk +1 (Q ) = q k (Q | Z ) Wk ( Z ) d Z, (1.3.13) где Z – вектор состояния системы на k-ом шаге, Q – вектор состояния системы на k+1-ом шаге, Wk – плотность вероятности координат системы на k-ом шаге, qk – условная плотность распределения вероятности на k-ом шаге, – область - 26 интегрирования, определяемая интервалами изменения координат вектора состояния. Определим условную ПРВ qk. Для этого необходимы плотности распределения вероятностей случайных процессов, пг и n. Будем считать, что они являются независимыми гауссовскими процессами с нулевыми математическими ожиданиями и СКО, пг и n соответственно. В этом случае плотности распределения вероятностей будут иметь вид:

1 2 1 2 пг 1 2 n p ( ) = p( пг ) = p (n) = exp( 2 ) ).

exp( пг 2 2 пг (1.3.14) exp( n 2 2 n ) Для дальнейших преобразований введем новый флуктуационный процесс:

[k ] = пг [k ] + T S n[k ].

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

2 2 2 = пг + T 2 S 2 n.

Следует особо отметить, что параметры процесса зависят не только от параметров порождающих его процессов, но и от параметров самой системы T и S. С учетом замены уравнение (1.3.12) примет вид:

[k + 1] = [k ] + T н + [k ] [k ] T у F ([k ]). [k + 1] = [k ] + [k ] (1.3.15) - 27 Введя обозначение = –, получим выражение для двумерной плотности распределения вероятности p(,):

1 2 p (, ) = exp( 2 ) exp( ( ) 2 2 ).

С учетом полученного выражения, условная ПРВ будет иметь вид:

rr q k (Q | Z ) = 1 2 exp( (Q2 Z 2 ) 2 2 ) exp( [Q2 Z 2 Q1 + Z1 + T н T y F ( Z1 )]2 2 ), (1.3.16) где Z1 = [k], Z2 = [k], Q1 = [k+1], Q2 = [k+1]. С учетом инвариантных движений в системе согласно (1.3.13) получим следующее уравнение [37]:

r rr r r Wk +1 (Q) = q' k (Q | Z ) Wk ( Z ) dZ.

(1.3.17) Область интегрирования учитывает изменение разности фаз в пределах от – до. Условная ПРВ qk связана с (1.3.5) следующим образом:

r q 'k (Q1, Q2 | Z ) = r qk (Q1 + 2l, Q2 | Z ).

l = По интегрирование ведется в бесконечных пределах. Необходимо отметить следующий важный момент. В полученном выражении присутствует операция бесконечного суммирования. Но замена пределов суммирования на конечные не приводит к качественным изменениям получаемого результата. Пределы суммирования в надо выбирать, исходя из оценки вероятности перескока изображающей точки на несколько периодов за одну итерацию. Так как условная ПРВ перехода (1.3.5) представляет собой экспоненциальную зависимость, вероятность такого события резко уменьшается с ростом дистанции, поэтому полученный ряд очень быстро сходится. Как показывают - 28 расчеты, для типичных параметров входного воздействия и параметров СФС при суммировании можно ограничиться 7–11 слагаемыми без заметной потери точности результата. Результатом вычислений по формуле (1.3.17) является совместная плотность распределения вероятности фазовой ошибки и флуктуаций фазы входного сигнала на любом шаге. Для нахождения плотности распределения вероятности (ПРВ) флуктуаций фазы на выходе системы воспользуемся формулой (1.3.2). Используя известную формулу, связывающую ПРВ разности двух случайных величин с совместной ПРВ этих величин, получим:

pk ( 'пг ) = Wk ( 'пг, )d.

(1.3.18) Рассмотрим дискретную систему фазовой синхронизации 2-го порядка. Ее уравнение имеет вид (1.3.11). Как было указано, данная система не является марковской моделью. Для перехода к марковской модели введем новую переменную:

x[k + 1] = d ([k + 1] [k ]) d ([k ] пг [k ]) T у a0 F ([k ]) TSa0 nk + (1 d )T н.

Система уравнений (1.3.11) преобразуется к виду:

[k + 1] = [k ] + [k ] пг [k ] T у a1 F ([k ]) TSa1n[k ] + x[k ] x[k + 1] = dx[k ] T у (da1 + a0 ) F ([k ]) TS (da1 + a0 )n[k ] + (1 d )T н. [k + 1] = [k ] + [k ] (1.3.19) Система уравнений (1.3.19) в случае слабокоррелированных процессов, пг и n описывает 3-хмерную марковскую случайную последовательность. Саму систему (1.3.19) с учетом новой координаты будем рассматривать в качестве расширенной марковской модели 3-го порядка. Для получения условной ПРВ введем новые флуктуационные процессы:

- 29 1[k ] = [k ] пг [k ] T S a1 n[k ] 2 [k ] = T S (d a1 + a0 ) n[k ]. 3 [ k ] = [ k ] Можно показать, что совместная плотность распределения вероятности этих новых величин имеет вид:

1 (2 ) T S | d a1 + a0 | n пг p( 1, 2, 3 ) = exp( 2 ) exp( 2 2 [T S (d a1 + a0 )]2 n ( 3 1 + ) exp(. a1 2 ) 2 d a1 + a0 ) 2 пг 2 (1.3.20) Согласно (1.3.2) введем обозначения: Z1 = [k];

Z2 = x[k];

Z3 = [k];

Q1 = [k+1];

Q2 = x[k+1];

Q3 = [k+1].

Тогда процессы 1, 2 и 3 можно выразить через эти величины с учетом (1.3.19) следующим образом:

1 = Q1 Z1 + T у a1 F ( Z1 ) Z 2 2 = Q2 d Z 2 + T у (d a1 + a0 ) F ( Z1 ) (1 d ) T н. 3 = Q3 Z Подставив эти выражения в (1.3.20), получим условную плотность распределения вероятности:

q k (Q | Z ) = q k (Q1, Q2, Q3 | Z1, Z 2, Z 3 ) == p ( 1 (Q, Z ), 2 (Q, Z ), 3 (Q, Z )).

(1.3.21) В итоге уравнение Колмогорова–Чепмена для дискретной системы 2–го порядка при наличии указанных воздействий будет иметь следующий вид:

- Wk +1 (Q1, Q2, Q3 ) = [ [ q' k (Q1, Q2, Q3 | Z1, Z 2, Z 3 ) Wk ( Z1, Z 2, Z 3 ) dZ1 ] dZ 2 ] dZ где q' k (Q1, Q2, Q3 | Z1, Z 2, Z 3 ) = l = qk (Q1 + 2l, Q2, Q3 | Z1, Z 2, Z 3 ).

Полученное уравнение будем называть расширенным векторным уравнением Колмогорова–Чепмена. Выражение для нахождения плотности распределения вероятности фазовых флуктуаций выходного сигнала в данном случае имеет вид:

p k ( ' пг ) = [ Wk ( ' пг, x, ) dx] d.

(1.3.22) 1.4. Модель цепи последовательно синхронизируемых генераторов в форме системы расширенных векторных уравнений Колмогорова– Чепмена и уравнений преобразования координат Цепь последовательно синхронизируемых генераторов представляет собой цепочку последовательно соединенных дискретных систем фазовой синхронизации, показанную на рис. 1.2.

ОГ ФД1 ФНЧ1 ГУН1 ФДN ФНЧN ГУНN Рис. 1.2. Структурная схема цепи последовательно синхронизируемых генераторов. Здесь ОГ – опорный генератор, ФДi – фазовый детектор, ФНЧi – фильтр нижних частот, ГУНi – генератор, управляемый напряжением (перестраиваемый генератор). На основе описания отдельных ячеек цепи последовательно синхронизируемых генераторов, представляющих собой дискретные СФС, - 31 построена общая модель этой системы. Она представляет собой комбинацию векторных уравнений Колмогорова–Чепмена следующего вида:

r rr r r Wk +1 (Q1 ) = q '1 (Q1 | Z1 ) Wk ( Z1 ) dZ1 r p k ( ' пг1 ) = L[Wk (Q1 )] L r r r r r W (Q ) = q ' (Q | Z ) W ( Z ) dZ k +1 N N N N N k N r p k ( ' пгN ) = L[Wk (Q N )] (1.4.1) Здесь Zi и Qi – переменные состояния i–го звена цепи. q'i – условная ПРВ компонент вектора Qi, p( ' пгi ) – ПРВ фазовых флуктуаций сигнала на выходе i– го звена. Оператор L представляет собой преобразование многомерной совместной ПРВ переменных состояния звена вида (1.3.18) или (1.3.22). Система уравнений (1.4.1) имеет ряд особенностей. Во-первых, она содержит векторные уравнения Колмогорова–Чепмена для каждого звена цепи, а так же уравнения для вычисления ПРВ фазовых флуктуаций выходного сигнала звена. Последние необходимы, поскольку выходной сигнал одной ячейки цепи является входным сигналом для следующей ячейки. Во-вторых, r условная ПРВ q'i для i–ого звена зависит от p( ' пгi 1 ). Это означает, что для решения уравнений, отвечающих за i–ую ячейку, необходимо решить уравнения из (1.4.1), отвечающие за все предыдущие ячейки. В свою очередь, это несколько усложняет процедуру поиска решения, поскольку заранее не известен вид условной ПРВ q'i при i > 1. В-третьих, согласно (1.3.19) в уравнения i-ой ячейки входят не фазовые, а частотные флуктуации, связанные с фазовыми следующим соотношением:

i [k ] = ' пгi 1 [k + 1] 'пгi 1 [k ].

(1.4.2) Согласно (1.4.2) для определения ПРВ i[k] необходима двумерная совместная плотность распределения вероятности координат ' пгi 1 [k+1] и ' пгi 1 [k]. Для ее нахождения предлагается следующий подход. Запишем - 32 совместную ПРВ координат системы в два соседних момента времени через условную ПРВ q' (Qk +1 | Qk ) :

r r r r r W (Qk +1, Qk ) = q' (Qk +1 | Qk ) W (Qk ), r r (1.4.3) где Wk (Q) = W (Qk ) – ПРВ координат системы на k-ом шаге. Интегрируя (1.4.3), можно перейти к совместной ПРВ фазовых флуктуаций выходного сигнала на соседних шагах ' пг [k+1] и ' пг [k] – W ( 'пгk +1, 'пгk ). Поскольку Qk содержит в себе k и k, то используя (1.3.2) и формулы преобразования ПРВ при переходе к новым координатам, можно получить искомую функцию W ( 'пгk +1, 'пгk ), из которой получается плотность r r r распределения вероятности частотных флуктуаций выходного сигнала:

p ( k ) = W ( 'пгk + k, 'пгk ) d 'пгk (1.4.4) В итоге уточненная модель цепи будет иметь вид:

r rr r r Wk +1 (Q1 ) = q '1 (Q1 | Z1 ) Wk ( Z1 ) dZ1 r p k ( ' пг1 ) = L[Wk (Q1 )] r p k ( 2 ) = M [Wk (Q1 )] L=L r r r r r W (Q ) = q ' (Q | Z ) W ( Z ) dZ k +1 N N N N k N N r p k ( ' пгN ) = L[Wk (Q N )] r p ( k N +1 ) = M [Wk (Q N )] (1.4.5) Здесь оператор M представляет собой описанное выше расширенное преобразование, приводящее к расчету плотности распределения вероятности частотных флуктуаций сигнала на выходе соответствующей ячейки. Модель (1.4.5) является более полной, по сравнению с моделью (1.4.1). Таким образом, процесс анализа поведения цепи последовательно синхронизируемых генераторов сводится к следующему. Для текущего звена - 33 определяются параметры входных воздействий. С помощью решения векторного уравнения Колмогорова–Чепмена находятся плотности распределения вероятности координат состояния системы. Затем производится преобразование ПРВ данных координат в ПРВ фазовых флуктуаций выходного сигнала и определяются статистические характеристики этого процесса. При данных входных воздействиях он не является стационарным, его дисперсия увеличивается со временем. Определяется ПРВ частотных флуктуаций выходного сигнала звена с помощью описанных преобразований. Исходя из этой ПРВ, находится условная плотность распределения вероятности координат системы для следующего звена. Затем аналогичная процедура выполняется для следующей ячейки цепи. Блок–схема данного алгоритма показана на рис. 1.3. Как следует из приведенных выше рассуждений, важнейшим фактором для анализа поведения цепи последовательно синхронизируемых с помощью дискретных СФС генераторов является исследование поведения одной ячейки этой цепи, находящейся под воздействием аддитивного белого шума и слабо коррелированных частотных флуктуаций входного сигнала и сигнала свободного перестраиваемого генератора. Существенно упростить процедуру анализа цепи может предположение, что частотные флуктуации выходных сигналов ячеек являются гауссовскими. В этом случае общий вид условных ПРВ является известным для всех ячеек. Неизвестными являются только отдельные их параметры. Такое допущение позволяет, в частности, достаточно эффективно проводить анализ цепи с помощью вычислительной техники [70, 71]. Полученная марковская модель цепи в форме комбинации векторных уравнений Колмогорова–Чепмена для отдельных ячеек и уравнений преобразования координат позволяет анализировать поведение цепи для произвольных режимов отдельных ячеек. В частности, речь может идти о переходных процессах, для которых является характерным изменение координат в широких пределах. Для анализа синхронных и квазисинхронных режимов, особенно в условиях ограниченных воздействий, весьма полезными могут оказаться линейные модели. При этом они требуют значительно меньших вычислительных затрат по сравнению с марковскими моделями. Ниже рассмотрен ряд таких моделей.

- 34 Задать значения параметров флуктуационных воздействий: ог, пгi, ni. Рассчитать ПРВ координат состояния текущего звена с помощью уравнения КЧ (1.3.17) на следующем шаге Рассчитать ПРВ фазовых пг и частотных флуктуаций на выходе звена с помощью (1.4.3)–(1.4.4) Определить стат. характеристики фазовых и частотных флуктуаций на выходе звена Определить условную ПРВ q для следующего звена Вернуться к первому звену Нет Достигнут конец цепи? Да Нужно вычислять стат. характеристики на следующем шаге?

Нет Закончить обработку Рис. 1.3. Блок–схема алгоритма анализа цепи последовательно соединенных генераторов с помощью аппарата марковских процессов.

- 35 1.5. Линеаризованная модель синхронизируемых генераторов цепи последовательно Линеаризованная модель рассмотрена автором диссертации в работе [82]. В данном разделе представлено более полное описание цепи, являющееся развитием результатов приведенной работы. На рис. 1.4 представлена схема исследуемой цепи, на которой изображены источники внешних воздействий.

n1 ОГ ог + ФД1 ФНЧ1 ГУН1 пг1 nN + ФДN ФНЧN ГУНN пгN Рис. 1.4. Сема модели цепи последовательно синхронизируемых генераторов с рассматриваемыми источниками случайных воздействий. К таковым относятся фазовые шумы перестраиваемых генераторов пгi, фазовый шум опорного генератора ог, и канальные шумы ni. Отсюда следует, что для описания функционирования одного звена цепи необходимо знание коэффициентов передачи на выход перестраиваемого генератора входных фазовых шумов, канальных шумов и шумов самого перестраиваемого генератора. С помощью этих коэффициентов передачи можно определить спектральную плотность фазовых флуктуаций на выходе любого звена системы. Используя данную спектральную плотность можно вычислить важные характеристики качества работы цепи, такие как дисперсию фазовых флуктуаций выходного сигнала. Кроме того, с помощью этой спектральной плотности можно вычислять и некоторые телекоммуникационные характеристики качества, такие как TDEV, TIE, дисперсию Алана, и др.

1.5.1. Линеаризованная дискретная модель На рис. 1.5 представлена линеаризованная модель дискретной системы фазовой синхронизации. Здесь E – крутизна характеристики фазового детектора - 36 в рабочей точке, S – крутизна перестраиваемого генератора, K(z) – передаточная функция фильтра в цепи управления, пг0 – собственная частота перестраиваемого генератора. Параметр T является периодом дискретизации в системе.

A ФД ФНЧ – вх[k] ГУН С B E K(z) T/(z–1) пг[k] S пг Рис. 1.5. Линейная модель дискретной СФС. Данная схема использует так называемые «фазовые» сигналы. На ней отмечены основные источники флуктуаций в системе. пг – фазовые флуктуации сигнала перестраиваемого генератора СФС, вх – фаза входного сигнала звена. В состав вх входят как детерминированная составляющая, так и флуктуации фазы сигнала на выходе предыдущего звена. Кроме того, в ее состав входит и шум, вносимый каналом передачи. Обычно аддитивный широкополосный шум канала пересчитывается на выход фазового детектора [17]. Однако использование линейной модели системы позволяет просто пересчитать его на вход СФС (точка А). Использование такого подхода позволяет свести описание звена цепи к двум коэффициентам передачи [55]: из точки A в точку C (описывает прохождение через систему шумов канала передачи и фазовых флуктуаций входного сигнала) и из точки B в точку C (описывает прохождение через систему фазовых шумов перестраиваемого генератора):

ESTK ( z ) K AC ( z ) = z 1, ESTK ( z ) 1+ z K BC ( z ) = 1 ESTK ( z ) 1+.

z (1.5.1) - 37 Введем обозначение = E·S·T. В этом случае для бесфильтрового звена имеем следующие выражения для этих коэффициентов передачи:

K AC ( z ) = 1 D, zD K BC ( z ) = z 1, zD где D = 1 –. Для определения спектральной плотности шума на выходе системы по спектральной плотности шума на ее входе необходимы квадраты модулей коэффициентов передачи. В данном случае они имеют вид:

1 2D + D 2 1 2 D cos( ) + D L AC ( ) =, LBC ( ) = 2 [1 cos( )] 1 2 D cos( ) + D.

(1.5.2) Здесь = ·T – цифровая частота. Пусть В этом случае:

K AC ( z ) = K BC ( z ) = K ( z) = m + 1. zd m z + (1 m d ) z + ( m 1 d ) z + (1 m d ) + d z 2 (1 + d ) z + d z 2 + ( m 1 d ) z + (1 m d ) + d,.

Для удобства введем следующие обозначения: A = ·m;

Тогда K AC ( z ) = K BC ( z ) = B = ·(1–m·d).

A z + B z 2 + ( A 1 d) z + B + d z 2 (1 + d ) z + d z 2 + (A 1 d ) z + B + d,.

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

B < 1 d B > A. B > A 2 (1 + d ) Функции LAC и LBC в данном случае будут иметь вид:

A 2 + B 2 + 2 A B cos( ) 1 + C 2 + D 2 + 2 C ( D + 1) cos( ) + 2 D cos(2 ) 1 + C 2 + D 2 + 2 C ( D + 1) cos( ) + 2 D cos( 2 ) L AC ( ) = LBC ( ) =,, (1.5.3) 1 + (1 + d ) 2 + d 2 2 (1 + d ) 2 cos( ) + 2 d cos( 2 ) где C = A – 1 – d, D = B + d. Используя функции LAC и LBC, запишем выражения для спектральной плотности выходных фазовых флуктуаций.

S вых ( ) = L AC ( ) S вх ( ) + LBC ( ) Sпг ( ) (1.5.4) где Sвх() – спектральная плотность фазовых флуктуаций сигнала на входе системы, Sпг() – спектральная плотность фазовых флуктуаций сигнала перестраиваемого генератора. Используя соотношение (1.5.4) можно построить линейную модель цепи последовательно синхронизируемых генераторов в виде:

- 39 S вых1 ( ) = L AC1 ( ) [ Sог ( ) + S n1 ( )] + LBC1 ( ) Sпг1 ( ) S вых 2 ( ) = L AC 2 ( ) [ S вых1 ( ) + S n 2 ( )] + LBC 2 ( ) Sпг 2 ( ) L, (1.5.5) S выхi ( ) = L ACi ( ) [ S выхi 1 ( ) + S ni ( )] + LBCi ( ) Sпгi ( ) L S выхN ( ) = L ACN ( ) [ S выхN 1 ( ) + S nN ( )] + LBCN ( ) SпгN ( ) где Sог() – спектральная плотность фазовых шумов опорного генератора, Sогi() – спектральная плотность фазовых шумов i-го перестраиваемого генератора, Sni() – спектральная плотность фазовых шумов, вносимых i–ым каналом передачи. Использование специальных допущений позволяет значительно упростить описание системы. Предположим, что все звенья и каналы передачи одинаковы: Sn1() = Sn2() = … = Sn(), Sпг1() = Sпг2() = … = Sпг(), LAC1() = LAC2() = … = LAC(), LBC1() = LBC2() = … = LBC(). (1.5.6) В таком случае вклад в спектральную плотность фазовых флуктуаций сигнала на выходе цепочки от опорного генератора равен S1 ( ) = Lог ( ) S ог ( ), Lог ( ) = L AC ( ) N.

[ ] (1.5.7) где N – длина цепи. Вклад от i–го перестраиваемого генератора:

S 2, i ( ) = LBC ( ) [ L AC ( )] N i S пг ( ) Просуммировав вклады спектральных плотностей от перестраиваемых генераторов, получим, что их общий вклад составит всех - 40 S 2 ( ) = Lпг ( ) Sпг ( ), N Lпг ( ) = LBC ( ) L AC ( ) N i. i = [ ] (1.5.8) Вклад, вносимый спектральной плотностью аддитивного шума на входе i– го звена, имеет вид:

S 3, i ( ) = [ L AC ( )] N i S n ( ).

Тогда общий вклад в спектральную плотность фазовых шумов на выходе цепи от аддитивных шумов составит S 3 ( ) = Lс ( ) S n ( ), N Lc ( ) = [ L AC ( )]i. i = (1.5.9) Просуммировав вклады от всех источников шумов, получим, что итоговая спектральная плотность фазовых шумов на выходе цепи длиной N имеет вид:

S выз ( ) = S1 ( ) + S 2 ( ) + S 3 ( ) = Lог ( ) Sог ( ) + Lпг ( ) Sпг ( ) + LC ( ) S n ( ) (1.5.10) Входящие в (1.5.10) функции Lог(), Lпг(), LC() выражаются в соответствии с формулами (1.5.7)–(1.5.9). Их можно назвать квадратами модулей эквивалентных коэффициентов передачи соответствующих шумов на выход системы. В данном случае все многочисленные источники флуктуационных воздействий сведены лишь к трем, отвечающим за канальный шум, фазовый шум опорного генератора и фазовый шум генератора цепи. Функции Lог(), Lпг() и LC() описывают передачу этих шумов на выход. Выражение (1.5.10) справедливо только в предположении (1.5.6). В ином случае следует использовать более общее выражение (1.5.5).

- 41 1.5.2. Линеаризованная аналоговая модель Наличие тригонометрических функций в выражениях для квадратов модулей коэффициентов передачи (1.5.2) и (1.5.3) значительно усложняет аналитическое исследование цепи. В то же время если предположить, что случайные воздействия имеют спектральную плотность, в основном сосредоточенную в области низких частот (это относится к фазовым шумам), то можно применить для анализа аналоговую модель звеньев. Из формул (1.2.2) и (1.2.3) видно, что основная мощность фазовых шумов перестраиваемых генераторов сосредоточена именно в низкочастотной области. Что касается аддитивного белого шума, то при анализе процессов на больших временах его роль будет несущественна. Поэтому применение аналоговой модели в этом случае оправдано. На рис. 1.6 изображена схема линейной модели аналоговой СФС. Здесь K(p) – передаточная функция фильтра в цепи управления. Остальные обозначения соответствуют обозначениям на рис. 1.5. Для данной модели ES 1 K ( p) K ( p) p K AC ( p ) = =, 1 + ES 1 K ( p ) p + K ( p) p K BC ( p) = 1 p =, 1 + ES 1 K ( p) p + K ( p) p (1.5.11) где = ES – параметр системы, характеризующий ее полосу удержания. От самой полосы удержания он отличается на множитель равный крутизне характеристики фазового детектора в рабочей точке.

- 42 A ФД ФНЧ – вх(t) ГУН С B E K(p) 1/p пг(t) S пг Рис. 1.6. Линейная модель аналоговой СФС. В случае звеньев на основе бесфильтровых СФС формулы (1.5.11) принимают вид:

K AC ( p) =, p+ K BC ( p) = p. p+ Квадраты их модулей:

2, 2 + L AC ( ) = LBC ( ) = 2 2 + (1.5.12) Данная система является устойчивой при >0. В случае, когда в основе звена цепи лежит СФС с RC–фильтром в цепи управления, имеющим коэффициент передачи вида K ( p) = 1, имеем 1+ T p K AC ( p ) = T p2 + p +, K BC ( p) = p + T p2 T p2 + p +.

Квадраты их модулей:

2 ( T 2 ) 2 + L AC ( ) =, L BC ( ) = 2 + T 2 ( T 2 ) 2 +.

(1.5.13) - 43 В случае, когда в основе звена цепи лежит СФС с пропорционально интегрирующим фильтром в цепи управления, имеющим коэффициент передачи вида K ( p) = 1+ m T p, имеем 1+ T p p + T p2 T p2 + p + + m T p K AC ( p) = + m T p T p2 + p + + m T p, K BC ( p ) =.

Квадраты их АЧХ представимы как:

2 + ( m T ) 2 ( T 2 ) 2 + ( + m T ) L AC ( ) = LBC ( ) =,. (1.5.14) 2 + T 2 ( T 2 ) 2 + ( + m T ) При положительных параметрах T, m и система устойчива. Формулы, позволяющие получить спектральную плотность фазовых флуктуаций сигнала на выходе любого звена цепи, аналогичны выражениям (1.5.4)–(1.5.10) с той разницей, что вместо цифровой частоты используется обычная циклическая частота. Кроме того, функции LAC и LBC выбираются из (1.5.12)–(1.5.14).

1.6. Выводы В данной главе рассмотрен ряд математических моделей цепи последовательно синхронизируемых генераторов. Введены модели входных воздействий, а так же модели звеньев цепи, являющихся дискретными системами фазовой синхронизации. Получены математические модели звеньев в форме стохастических разностных уравнений. Произведено сведение указанных уравнений к марковскому виду. Так же получены выражения для условных ПРВ, используемых в расширенных уравнениях Колмогорова– Чепмена, и записаны сами эти уравнения. Кроме того, построена - 44 линеаризованная модель цепи, позволяющая при определенных условиях выполнить исследование с меньшими затратами. Основные результаты главы можно свести к следующим: 1. Разработана математическая модель цепи последовательно синхронизируемых генераторов в форме комбинации расширенных уравнений Колмогорова–Чепмена и уравнений преобразования координат. Для бесфильтровых звеньев уравнение Колмогорова–Чепмена имеет второй порядок, для звеньев с фильтром 1-го порядка в цепи управления – 3-ий порядок. 2. Для получения характеристик фазовых флуктуаций на выходе звена предложен подход, основанный на увеличении размерности марковской модели. Выражение для ПРВ фазового шума на выходе вычисляется из многомерной ПРВ, являющейся результатом решения расширенного уравнения Колмогорова–Чепмена. 3. Для расчета ПРВ частотных флуктуаций, необходимых для марковской модели, предложен подход, основанный на использовании совместной ПРВ координат состояния системы на двух последовательных шагах работы. 4. Математическую модель цепи последовательно синхронизируемых генераторов в форме комбинации расширенных уравнений Колмогорова– Чепмена и уравнений преобразования координат системы в выходные координаты целесообразно использовать для цепей ограниченной длины. Для исследования цепей произвольной длины требуются значительные вычислительные ресурсы. Некоторым ограничением к использованию является также тот факт, что структура уравнений заранее неизвестна и определяется решением на предыдущем шаге. Имеется в виду вид условной ПРВ. В то же время в случае выполнения ряда условий аппарат процессов Маркова дает возможность точно исследовать поведение цепи. 5. Использование линеаризованной модели цепи последовательно синхронизируемых генераторов в условиях ограниченных воздействий позволяет получить результаты при существенно меньших вычислительных затратах. Кроме того, она дает возможность вычислять спектральные плотности фазовых флуктуаций в различных точках цепи. Знание данных спектральных плотностей позволяет вычислять важные радиотехнические и телекоммуникационные характеристики качества работы системы. Линейная - 45 модель позволяет описать поведение каждого звена с помощью только лишь двух коэффициентов передачи. 6. Аналитическое исследование дискретной линеаризованной модели цепи может быть весьма затруднено наличием тригонометрических функций циклической частоты в квадратах модулей коэффициентов передачи звеньев. В этом случае полезным может быть рассмотрение аналоговой модели цепи. Такая модель так же является адекватной в том случае, если основная мощность случайных воздействий сосредоточена в области низких частот.

- 46 Глава 2. Анализ и оптимизация статистических характеристик звена цепи на основе дискретной СФС 2-го порядка в условиях комбинированных случайных воздействий 2.1. Постановка задачи Основой цепи последовательно синхронизируемых генераторов являются ее звенья, представляющие собой дискретные системы фазовой синхронизации. Следует ожидать, что поведение всей цепи при сложных входных воздействиях на систему тесным образом связано с поведением отдельных ее звеньев. Этой тематике посвящена достаточно много работ, написанных как российскими, так и зарубежными авторами [27, 31–36]. Например, большое внимание исследователей привлекли вопросы динамики дискретных СФС в условиях шумовых воздействий [37–40, 43, 44]. Разработаны специальные методы исследования статистической динамики. К ним относится мощный аппарат теории марковских процессов, позволяющий исследовать динамику нелинейных систем, каковыми являются дискретные СФС [17, 20, 57, 58]. В то же время существует определенная специфика при исследовании систем фазовой синхронизации, являющихся звеньями цепи последовательно синхронизируемых генераторов. В первую очередь она связана с характером флуктуационных воздействий на систему. Анализу поведения дискретных СФС в условиях шумов посвящены работы [16, 17, 22–26, 38]. Однако большинство из них направлено на изучение влияния простейших флуктуационных воздействий. Широко рассмотрены вопросы, касающиеся воздействия широкополосного аддитивного шума. Значительно меньше работ посвящено более сложным воздействиям на систему [52, 56, 60]. Практически без внимания на сегодняшний день оставлен класс воздействий в виде фазовых флуктуаций сигналов в СФС. Между тем именно такие шумы способны играть решающую роль в формировании качества работы цепи последовательно синхронизируемых генераторов [106]. Следующей особенностью, которую необходимо учитывать при исследовании дискретных СФС в цепочке, является предмет исследования. Подавляющее большинство работ по анализу характеристик систем фазовой синхронизации посвящено изучению фазовой ошибки в системе. Уравнения, описывающие СФС, так же написаны в терминах фазовой ошибки [34]. Однако - 47 при исследовании цепи последовательно синхронизируемых генераторов интерес представляют характеристики фазы выходного сигнала, поскольку именно они определяет качество работы системы. В последнее время большой интерес вызывают вопросы оптимизации и синтеза оптимальных систем. Возросшие требования на качество обслуживания заставляют разработчиков добиваться наилучшего в определенном смысле функционирования от устройств. Это касается и систем фазовой синхронизации [42, 86]. Широкий круг применения данных устройств еще больше подогревает интерес к вопросам оптимизации их работы. В случае использования дискретных СФС в цепи последовательно синхронизируемых генераторов вопросы выбора наилучших параметров выдвигаются на первое место. Интуитивно понятно, что чем лучше сигнал на выходе одного звена, тем качественнее он будет на выходе всей цепи. Однако, как выяснил в своих работах автор диссертации [75], это выполняется далеко не всегда. В связи с этим вопросы оптимизации параметров звеньев в условиях комбинированных случайных воздействий выходит на первый план. Целью данной главы является анализ и оптимизация характеристик фазовой ошибки и фазы выходного сигнала одного звена цепи последовательно синхронизируемых генераторов на основе дискретных СФС 2-го порядка в условиях комбинированных флуктуационных воздействий. В качестве исследуемой характеристики рассматривается дисперсии фазовой ошибки и фазовых флуктуаций выходного сигнала. Критерием оптимальности является минимум дисперсий. 2.2. Оптимизация приближении статистических характеристик в линейном Рассмотрим вначале вопрос оптимизации линейной дискретной модели звена цепи на основе дискретной СФС Автором диссертации опубликовано несколько работ [76, 78, 86], в которых он рассмотрел различные подходы к оптимизации дискретных СФС при наличии сложных флуктуационных воздействий на систему. На рис. 2.1 приведена схема линейной модели бесфильтровой дискретной системы фазовой синхронизации с источниками флуктуационных воздействий.

- 48 вх[k] A E + C пг[k] 1 / (z–1) S – [k] D n[k] + B Рис. 2.1. Схема линейной модели бесфильтровой дискретной системы фазовой синхронизации с источниками флуктуационных воздействий. Данная схема не содержит источников детерминированных сигналов (таких как пг0), поскольку они не вносят вклад во флуктуационные составляющие фазовой ошибки и фазы выходного сигнала. Жирными буквами обозначены основные точки схемы, в которых вводятся входные воздействия или снимаются выходные сигналы. Поскольку данная схема является линейной, то ее можно описывать с помощью передаточных функций. В таблице 2.1 приведены передаточные функции от входных воздействий на каждый из выходов. Таблица 2.1. На фазовую ошибку (точка D) K AD ( z ) = z 1, z + S 1 S, K BD ( z ) = z + S 1 ( z 1). K CD ( z ) = z + S На фазу выходного сигнала (точка E) K AE ( z ) = S, z + S 1 S, K BE ( z ) = z + S 1 z 1. K CE ( z ) = z + S Спектральная плотность флуктуаций на выходе вычисляется по формуле:

S X ( z ) = K AX ( z ) Sвх ( z ) + K BX ( z ) S n ( z ) + K СX ( z ) Sпг ( z ) 2 2 (2.2.1) где X = D для фазовой ошибки и X = E для фазы выходного сигнала. Sвх, Sпг, Sn – спектральные плотности входных воздействий. Дисперсия сигнала в соответствующей точке вычисляется по формуле:

- 49 2 X = 1 2j | z | = S X ( z) z dz (2.2.2) В дальнейшем будем считать, что шум n[k], представляющий собой канальный шум, пересчитанный внутрь кольца, является широкополосным:

2 S n ( z) = n (2.2.3) а) Оптимизация дисперсии фазовой ошибки при белых фазовых шумах генераторов В данном случае спектральные плотности фазовых шумов:

2 Sвх ( z ) = вх 2 Sпг ( z ) = пг (2.2.4) Тогда из (2.2.1) и (2.2.2) следует, что дисперсия фазовой ошибки выражается формулой:

2 = S 2 2 2 2 n + вх + пг 2S 2S [ ] (2.2.5) Производная (2.2.5) по S всегда положительна:

2 d 2 2 2 2 вх + пг + n dS = [ ] (2 S ) (2.2.6) Это означает, что в данном случае оптимальным будет нулевое значение S, соответствующее границе области устойчивости. б) Оптимизация дисперсии фазовой ошибки при белых частотных шумах генераторов В этом случае спектральные плотности фазовых шумов:

- 50 Sвх ( z ) = 2 вх ( z 1)( z 1 1) Sпг ( z ) = 2 пг ( z 1)( z 1 1) (2.2.7) Подставив (2.2.7) в (2.2.1) и (2.2.2), получим:

2 = S 1 2 2 2 n + вх + пг 2S (2 S ) S [ ] (2.2.8) В данном случае 2 d 2 2 2 2 2 2 n S 2 + ( вх + пг ) S ( вх + пг ) dS = [ ] [ S (2 S )] (2.2.9) Введем обозначение:

2 2 вх + пг = 2 n (2.2.10) В этом случае оптимальное значение параметра S внутри области устойчивости имеет значение:

2 + 4 S opt = (2.2.11) в) Оптимизация дисперсии фазовых флуктуаций выходного сигнала при белых фазовых шумах генераторов Спектральные плотности фазовых шумов в этом случае выражаются формулами (2.2.4). Здесь 2 = 2 S 2 2 2 вх + n + пг 2S 2S [ ] (2.2.12) - 51 Как и в случае анализа фазовой ошибки при тех же видах шумов имеем, что производная (2.2.12) по S всюду положительна:

2 d 2 2 2 2 вх + пг + n dS = [ ] (2 S ) (2.2.13) То есть оптимальным так же будет являться нулевое значение S. г) Оптимизация дисперсии фазовых флуктуаций выходного сигнала при белых частотных шумах генераторов В этом случае провести процедуру оптимизации не удается, поскольку при белых частотных шумах генераторов фазовый шум выходного сигнала является нестационарным. Более сложной задачей является оптимизация параметров фазового шума выходного сигнала системы фазовой синхронизации, если входные воздействия имеют более сложный характер. Для этого предлагается следующий подход. В качестве критерия оптимальности будем использовать минимум дисперсии этого фазового шума. В случае использования линейной модели бесфильтрового аналогового звена спектральная плотность фазовых флуктуаций на его выходе выражается соотношением:

2 2 S ( ) + S ( ) 2 + 2 вх 2 + 2 пг S ( ) = (2.2.14) Здесь Sвх – спектральная плотность фазовых флуктуаций на входе звена, включающая вх и n, Sпг – спектральная плотность собственных фазовых флуктуаций перестраиваемого генератора звена. Отметим, что эти величины не зависят от. Искомый критерий оптимальности связан с S() соотношением:

2 = S ( ) d.

l h (2.2.15) - 52 Здесь h – верхняя частота, определяемая полосой пропускания системы и частотой дискретизации. l – нижняя частота, определяемая временем наблюдения за процессом. S() – спектральная плотность флуктуаций в точке, где вычисляется дисперсия. Наличие нижней граничной частоты является принципиальным, что будет показано ниже. Чтобы определить, при каком значении полосы удержания дисперсия фазового шума на выходе будет минимальной, подставим (2.2.14) в (2.2.15) и продифференцируем по. Получено выражение:

d 2 = d h l ( 2 + ) [S вх ( ) S пг ( )] d (2.2.16) Обозначим для краткости G() = Sвх() – Sпг(). Проанализируем данное выражение. Функция 2 (2.2.17) ( + ) является всюду положительной при >0. Поэтому если G() сохраняет свой знак во всей области интегрирования, то (2.2.16) нигде не обращается в 0. Т. е. дисперсия фазового шума не имеет экстремального значения. При G()<0 (Sвх() < Sпг()) производная всюду отрицательна. Значит, с ростом значение дисперсии непрерывно уменьшается. Это соответствует случаю, когда входной сигнал является более качественным, чем собственный сигнал перестраиваемого генератора. Тогда увеличение полосы удержания приводит к большему ослаблению сигнала ПГ и лучшему пропусканию входного сигнала. При = на выход будет подаваться непосредственно входной сигнал. При G()>0 (Sвх() > Sпг()) производная всюду положительна. Т. е. с ростом значение дисперсии увеличивается. В данном случае собственный сигнал перестраиваемого генератора более высококачественный, чем входной сигнал. Уменьшая, мы ослабляем входной сигнал. При =0 на выходе будет чистый сигнал собственного генератора ячейки.

- 53 Более интересен вариант, когда в интервале интегрирования [l;

h] имеются области, в которых G()>0, и области, где G()<0. В этом случае возможен вариант, когда (2.2.16) обратится в ноль при значениях полосы удержания, отличных от 0 и. Существует качественное различие между аналоговой и дискретной моделями. В случае если фазовые шумы перестраиваемых генераторов являются белыми, формула (2.2.13) показывает, что при любом соотношении мощностей внешних воздействий оптимальным будет нулевое значение параметра S, определяющего полосу удержания СФС. В случае же аналоговой модели формула (2.2.16) дает такой же результат только при Sвх() > Sпг(). Физически данное различие можно объяснить следующим образом. Цепь обратной связи аналоговой системы не обладает запаздыванием. Поэтому сигнал управления может быть коррелирован с текущими значениями шумов. Таким образом, может происходить их гашение. В дискретной системе всегда присутствует запаздывание сигнала управления. При белом фазовом шуме ПГ сигнал управления способен только увеличить мощность фазовых флуктуаций на выходе. Именно поэтому нулевое значение S, при котором данный сигнал равен 0, является оптимальным. С точки зрения передаточных функций системы такое отличие объясняется следующим образом. С ростом полосы пропускания входной фазовый шум лучше проходит на выход как в случае аналоговой, так и дискретной СФС. С фазовым же шумом перестраиваемого генератора дело обстоит иным образом. В обоих случаях с ростом усиливается подавление низкочастотных составляющих этого сигнала. Однако в высокочастотной области наблюдается существенная разница. У дискретной системы происходит усиление ВЧ– составляющих шума, чего нет у аналоговой системы. Таким образом в интервале учитываемых частот от l до h площадь под кривой АЧХ у дискретной системы увеличивается, а у аналоговой – уменьшается. Этим и объясняется различная реакция на белые фазовые шумы перестраиваемого генератора. В то же время реакция на флуктуации, чья спектральная плотность спадает достаточно быстро, одинакова у обеих систем. Это вызвано тем, что обе системы ведут себя одинаково в низкочастотной области. Поведение же в высокочастотной области в данном случае не важно, поскольку мощность флуктуаций там пренебрежимо мала.

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

2(cos( ) 1) d 2 (1 D 2 ) S вх ( ) + 2( D cos( )) S пг ( ) d = 22 dD l [1 2 D cos( ) + D ] [ ] где D = 1 – ·T, – цифровая частота. Sвх() – спектральная плотность фазовых флуктуаций входного сигнала, Sпг() – спектральная плотность фазовых флуктуаций сигнала перестраиваемого генератора. l в данном случае – нижняя цифровая частота, определяемая временем наблюдения за процессом. 2.3. Анализ статистических характеристик с помощью численного решения расширенного векторного уравнения Колмогорова–Чепмена В этом разделе рассмотрены результаты анализа статистических характеристик дискретных систем фазовой синхронизации с помощью аппарата марковских процессов. Часть данного материала была опубликована автором диссертации в работе [74]. Математической основой исследования, проводимого здесь, являются соотношения, приведенные в разделе 1.3. Особенностью входных воздействий на звено является то, что фазовые флуктуации входного сигнала и сигнала перестраиваемого генератора дискретной СФС представлены белыми частотными шумами с дисперсиями 2вх и 2пг соответственно. Это означает, что в модели спектральной плотности фазовых шумов (1.2.3) отличен от 0 только коэффициент a2. Канальный шум, пересчитанный внутрь кольца, так же является белым с дисперсией 2n. 2.3.1. Характеристики фазовой ошибки и фазовых выходного сигнала звена на основе бесфильтровой СФС флуктуаций Математическая модель бесфильтровой дискретной СФС в форме уравнений Колмогорова–Чепмена имеет вид (1.3.17). Для определения характеристик фазовых флуктуаций выходного сигнала требуется использовать - 55 совместную плотность распределения вероятности координат и вх. Типичный пример этой функции показан на рис. 2.2.

W(, вх) вх 0 – – Рис. 2.2. Совместная ПРВ фазовой ошибки и входного фазового шума при вх = 0.1, пг = 0.7, n = 0.1, T = 1, S = 0.84.

В работе рассмотрена зависимость дисперсии фазовых флуктуаций на выходе звена от параметра S (рис. 1.1). Данный параметр отвечает за величину полосы удержания в системе. В исследовании предполагалось, что T = 1 (система цифровая без учета эффектов квантования). Кроме того, считалось, что начальная расстройка равна 0. Последнее допущение основано на том, что в цепочке последовательно синхронизируемых генераторов их собственные частоты обычно равны или близки. На рис. 2.3 представлена график данной зависимости для различных значений дисперсии белого частотного шума.

- 56 S Рис. 2.3. Зависимость дисперсии фазовых флуктуаций на выходе дискретной СФС при вх = 0.1, n = 0.1 для 1) пг = 0.1, 2) пг = 0.4, 3) пг = 0.7.

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

- 57 1 2 S Рис. 2.4. Зависимость дисперсии фазовых флуктуаций на выходе дискретной СФС при пг = 0.1, n = 0.1 для 1) вх = 0.1, 1 2 2) вх = 0.4, 3) вх = 0.7.

S Рис. 2.5. Зависимость дисперсии фазовых флуктуаций на выходе дискретной СФС при вх = 0.1, пг = 0.1 для 1) n = 0.1, 2) n = 0.4, 3) n = 0.7.

- 58 Очень важным является факт, что фазовые флуктуации сигнала на выходе дискретной СФС не являются стационарными. Флуктуации фазы входного сигнала вх связаны с флуктуациями его частоты формулой вх[k+1] = вх[k] + вх[k].

В случае если вх – белый гауссов шум, то вх в пределе имеет бесконечную дисперсию. Аналогичными свойствами обладает и фазовый шум выходного сигнала. На рис. 2.6 показано, как изменяется ПРВ фазовой ошибки и фазового шума выходного сигнала со временем. Как видно, ПРВ фазовой ошибки практически не меняется, она достигла стационарного состояние. В то же время ПРВ ’пг постоянно расширяется.

W() а) - 59 W(’пг) ’пг б) Рис. 2.6. ПРВ а) фазовой ошибки и б) флуктуаций фазы выходного сигнала на 8-ом шаге (1) и на 16-ом шаге (2) при вх = 0.1, пг = 0.4, n = 0.1, T = 1, S = 1.

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

- 60 1 2 S Рис. 2.7. Зависимость дисперсии фазовых флуктуаций выходного сигнала от S для вх = 0.1, пг = 0.4, n = 0.1, T = 1, S = 1 на 1) 8-ом шаге;

1 2 2) 16-ом шаге;

3) 24-ом шаге.

S Рис. 2.8. Зависимость дисперсии фазовой ошибки от параметра S при вх = 0.1, n = 0.1 для 1) пг = 0.1, 2) пг = 0.4, 3) пг = 0.7.

- 61 1 2 S Рис. 2.9. Зависимость дисперсии фазовой ошибки от параметра S при вх = 0.1, пг = 0.1 для 1) n = 0.1, 1 2 2) n = 0.4, 3) n = 0.7.

S Рис. 2.10. Зависимость дисперсии фазовой ошибки от параметра S при n = 0.1, пг = 0.1 для 1) вх = 0.1, 2) вх = 0.4, 3) вх = 0.7.

- 62 На рис. 2.9 показана зависимость дисперсии фазовой ошибки от S для различных мощностей канального шума. С ростом n оптимальное значение полосы удержания уменьшается. Графики на рис. 2.8 и 2.10 показывают данные зависимости при различных мощностях фазового шума ПГ и фазового шума на входе соответственно. Графики идентичны, поскольку влияния данных шумов на фазовую ошибку эквивалентны. Данный факт следует из уравнения (1.3.1) в которое образы этих фазовых шумов в Z–области входят идентично. Было проведено сравнение оптимального значения параметра S, рассчитываемого по формуле (2.2.11), с оптимальным значением, полученным с помощью метода марковских процессов. Результаты представлены на рис. 2.11. Как видно из графиков, с увеличением мощности флуктуационных воздействий реальное оптимальное значение S становится больше расчетного. Данный результат объяснятся влиянием нелинейной характеристики фазового детектора. В работе была получена более общая формула (2.3.1), позволяющая учесть влияние данной нелинейности.

S opt = (L ) 2 + 4 (L ) (2.3.1) Здесь L – коэффициент усиления фазового детектора. В случае если мощности флуктуационных воздействий на систему малы, он равен производной характеристики ФД в рабочей точке. В данном исследовании использовалась синусоидальная характеристика, поэтому L=1. С ростом мощности воздействий на систему L падает, что приводит к росту Sopt.

- 63 Sopt Sopt S S а) Sopt б) S в) Рис. 2.11. Зависимость дисперсии фазовой ошибки от параметра S при n = 0.1, вх = 0.1 для а) пг = 0.1, б) пг = 0.2, в) пг = 0.3.

2.3.2. Анализ фазовой ошибки звена 2-го порядка Был проведен анализ статистических характеристик сигналов дискретной системы фазовой синхронизации с фильтром 2-го порядка [78]. В качестве фильтра был использован цифровой ПИФ. В данном разделе рассмотрено поведение дисперсии фазовой ошибки в таком звене. На рис. 2.12 показана зависимость данного параметра от S для различных мощностей белого частотного шума перестраиваемого генератора. Как видно, качественно графики не отличаются от аналогичных кривых для бесфильтровой дискретной - 64 СФС. Так же происходит увеличение оптимального значения полосы удержания при росте пг. На рис. 2.13 представлены аналогичные зависимости для разных мощностей белого частотного шума входного сигнала. Графики на рис. 2.12 и 2.13 не отличаются друг от друга, что говорит о том, что влияния фазовых шумов входного сигнала и сигнала перестраиваемого генератора на фазовую ошибку эквивалентны. В пользу этого факта говорит и выражение для образа фазовой ошибки в Z–области (1.3.1), в которое данные шумы входят симметрично.

1 2 S Рис. 2.12. Зависимость дисперсии фазовой ошибки от параметра S при d = 1, m = 2, вх = 0.1, n = 0.1 для 1) пг = 0.1, 2) пг = 0.4, 3) пг = 0.7.

- 65 1 2 S Рис. 2.13. Зависимость дисперсии фазовой ошибки от параметра S при d = 1, m = 2, пг = 0.1, n = 0.1 для 1) вх = 0.1, 2) вх = 0.4, 3) вх = 0.7.

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

- 66 1 2 S Рис. 2.14. Зависимость дисперсии фазовой ошибки от параметра S при d = 1, m = 2, вх = 0.1, пг = 0.1 для 1) n = 0.1, 2) n = 0.4, 3) n = 0.7.

Было проведено исследование влияния параметров фильтров на характеристики фазовой ошибки. На рис. 2.15 показана зависимость дисперсии фазовой ошибки от параметра S при различных значениях m. Как видно, с увеличением значения этого параметра фильтра происходит уменьшение оптимальной величины полосы пропускания. С изменением m происходит изменение области значений S, в которой система устойчива. Максимально допустимое значение:

S max = 2 (1 + d ). (1 + d ) m Как видно, с ростом m оно уменьшается. Этим объясняется тот факт, что кривые для различных значений этого параметра перекрывают не весь диапазон изменения S.

- 67 3 2 S Рис. 2.15. Зависимость дисперсии фазовой ошибки от параметра S при d = 1, вх = 0.1, пг = 0.7, n = 0.1 для 1) m = 2, 2) m = 4, 3) m = 6.

Еще одной немаловажной особенностью является то, что с увеличением m падает величина дисперсии фазовой ошибки для оптимального значения полосы удержания. Дело в том, что с ростом параметра m увеличивается подавление низкочастотных составляющих фазовых шумов входного сигнала и сигнала перестраиваемого генератора. Поскольку основная их мощность сосредоточена именно в этой области, происходит уменьшение дисперсии фазовой ошибки. 2.3.3. Характеристики фазовых флуктуаций на выходе звена 2-го порядка Так же был произведен анализ статистических характеристик фазовых флуктуаций на выходе звена 2-го порядка. На рис. 2.16–2.18 показаны зависимости дисперсии этого процесса от S для различных мощностей входных воздействий.

- 68 1 2 S Рис. 2.16. Зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S при d = 1, m = 2, вх = 0.1, n = 0.1 для 1) пг = 0.1, 1 2 2) пг = 0.4, 3) пг = 0.7.

S Рис. 2.17. Зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S при d = 1, m = 2, пг = 0.1, n = 0.1 для 1) вх = 0.1, 2) вх = 0.4, 3) вх = 0.7.

- 69 1 2 S Рис. 2.18. Зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S при d = 1, m = 2, вх = 0.1, пг = 0.1 для 1) n = 0.1, Характер зависимостей 2) n = 0.4, аналогичен 3) n = 0.7. графикам, полученным для бесфильтрового звена. На рис. 2.19 продемонстрированы аналогичные кривые, снятые для разного числа шагов системы. Видно, что фазовые флуктуации выходного сигнала дискретной СФС являются нестационарным процессом. С течением времени их дисперсия увеличивается. В то же время положение оптимального значения S практически не изменяется.

- 70 1 2 S Рис. 2.19. Зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S при d = 1, m = 2, n = 0.1, вх = 0.1, пг = 0.1 и времени анализа 1) 16 шагов, 2) 24 шага, 3) 32 шага.

На рис. 2.20 представлена зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S для различных m. С ростом m так же наблюдается уменьшение оптимального значение полосы удержания, как и для случая фазовой ошибки. Рост данного параметра фильтра вызывает усиление ростом m. подавления низкочастотных составляющих фазового шума перестраиваемого генератора. Этим объясняется падение уровня дисперсии с - 71 3 2 S Рис. 2.20. Зависимость дисперсии фазовых флуктуаций выходного сигнала дискретной СФС от параметра S при d = 1, вх = 0.1, пг = 0.4, n = 0.1 для 1) m = 2, 2) m = 4, 3) m = 6.

2.4. Анализ статистических характеристик в линейном приближении Как уже говорилось ранее, метод марковских процессов позволяет проанализировать поведение системы только для некоторых типов входных воздействий. В данном случае приходится ограничиваться широкополосными частотными шумами входного сигнала и сигнала перестраиваемого генератора СФС, а так же широкополосным канальным шумом. В общей формуле (1.2.3) марковская модель позволяет рассмотреть только случай, когда a20, а все другие коэффициенты очень малы по сравнению с a2. Для анализа поведения звена цепи в более общей ситуации приходится использовать линеаризованные модели звеньев, описанные в разделе 1.5. Автором диссертации опубликовано несколько работ [83, 86, 87], в которых им рассматривается поведение СФС при наличии сложных флуктуационных воздействий. Материал этого раздела основан на указанных работах и является их продолжением.

- 72 2.4.1. Применение аналоговой модели В начале было рассмотрено поведение линеаризованных аналоговых звеньев при наличии комплексного воздействия на них. Аналоговые звенья являются более простыми для анализа. Рассмотрение поведения дискретных звеньев велось путем сравнения с аналоговыми моделями и вычленения отличий. а) Анализ бесфильтровой модели Воспользовавшись моделями входных воздействий из раздела 1.2, были получены следующие результаты для работы бесфильтровой системы. На рис. 2.21 показана зависимость нормированной дисперсии ошибки на выходе одного звена от для случая белых фазовых шумов генераторов.

, рад/с, рад/с а) б) Рис. 2.21. Зависимость нормированной дисперсии ошибки на выходе одного звена от при а) a0 = 100 рад с;

b0 = 1 рад с;

б) a0 = 1 рад с;

b0 = 100 рад с.

Результаты показывают полное соответствие с выводами, сделанными в разделе 2.2. В случае если входной сигнал ячейки «хуже» собственного, оптимальным является нулевое значение, и наоборот. На рис. 2.22 представлена зависимость нормированной дисперсии ошибки на выходе одного звена от, для случая фазовых флуктуаций 1-го порядка.

- 73 Считается, что на вход системы поступает сумма фазовых шумов канала и предыдущего звена, которую можно описывать моделью (1.2.3). Так же присутствует собственный шум генератора звена. Здесь для фазового шума перестраиваемого генератора: b0 = 1 рад с, b1 = 30000 рад2. Как видно, с ростом белого фазового шума на входе значение оптимальной полосы удержания уменьшается, что обусловлено необходимостью лучше фильтровать входной сигнал. Аналогичное влияние на оптимальное значение оказывает изменение мощности фликкер–фазового шума. Увеличение его мощности во входном сигнале приводит к уменьшению положения минимума на кривой.

, рад/с, рад/с а) б), рад/с, рад/с в) г) Рис. 2.22. Зависимость нормированной дисперсии фазового шума на выходе звена от полосы удержания для а) a0 = 100 рад·с;

a1 = 20000 рад2;

б) a0 = 340 рад·с;

a1 = 20000 рад2;

в) a0 = 200 рад·с;

a1 = 10000 рад2;

г) a0 = 200 рад·с;

a1 = 25000 рад2.

- 74, рад/с, рад/с а) б), рад/с, рад/с в) г) Рис. 2.23. Зависимость нормированной дисперсии фазового шума на выходе звена от полосы удержания для а) a3 = 0 рад4·/ с2;

b2 = 29000 рад3 / c;

b3 = 0 рад4 / c2;

б) a3 = 0 рад4·/ с2;

b2 = 11000 рад3 / c;

b3 = 0 рад4 / c2;

в) a3 = 20000 рад4·/ с2;

b2 = 25000 рад3 / c;

b3 = 23300 рад4 / c2;

г) a3 = 20000 рад4·/ с2;

b2 = 25000 рад3 / c;

b3 = 24350 рад4 / c2. В работе так же проведены исследования зависимости значения дисперсии фазового шума на выходе звена от полосы удержания для шумов 2го и 3-го порядков. На рис. 2.23 показано, какой сложный характер может носить эта зависимость. Как видно из данного рисунка, в зависимости от соотношения параметров входного воздействия и шума собственного генератора звена мы можем получить самые различные характеры - 75 зависимостей. На кривых 2() могут иметься как максимумы, так и минимумы, что еще раз подтверждает необходимость проведения процедуры оптимизации параметров звена. Здесь a0 = 140 рад·с;

a1 = 20000 рад2;

b0 = 1 рад·с;

b1 = 30000 рад2;

a2 = 50000 рад3 / с. Серьезным вопросом при проведении оптимизации является зависимость оптимального значения полосы удержания от времени анализа. Существование данной зависимости объясняется рис. 2.24. На нем изображен вариант графика функции G() для области низких частот. Как известно, минимальная фиксируемая в исследовании частота обратно пропорциональна времени анализа. В случае если нижний предел интегрирования в (2.2.2) равен 1l, во всей области интегрирования функция G() положительна. Это означает, что минимальная дисперсия фазовой ошибки будет при нулевой полосе удержания. При увеличении времени анализа (уменьшении нижнего предела интегрирования) появляется область, в которой данная функция имеет уже другой знак. В этом случае может существовать отличное от 0 значение оптимальной полосы удержания.

2 l 1 l Рис. 2.24 График функции G(). На рис. 2.25 показаны графики, отражающие зависимость этого оптимального значения от времени анализа. В данном случае в представлении спектральных плотностей шума мы ограничились 1–ой степенью частоты. Следует отметить, что дальнейшее увеличение времени анализа незначительно - 76 влияет на характер зависимости производной фазового шума на выходе звена от. Это связано с тем, что функция 2 ( + ), играющая роль окна в выражении (2.2.3), быстро стремится к 0 при уменьшении частоты. Влияние может оказать только составляющая спектральной плотности шумов, стоящая в (1.2.3) при a3, но она обычно мала [102].

, рад/с, рад/с а) б) Рис. 2.25. Зависимость нормированной дисперсии фазового шума от полосы удержания СФС для различных времен анализа при a0 = 140 рад·с;

a1 = 20000 рад2;

b0 = 1 рад·с;

b1 = 30000 рад2: а) = 1.57·10–2 с;

б) = 3.79·10–2 с. б) Анализ модели с RC–цепочкой Анализ поведения моделей систем фазовой синхронизации с фильтрами в цепи управления будем проводить путем сравнения квадратов их АЧХ с квадратами АЧХ бесфильтровой системы. Используя выражения (1.5.12) и (1.5.13), было исследовано поведение данных систем по отношению к шумам различных частотных диапазонов. На рис. 2.26 показаны графики LAС() для бесфильтровой системы и системы с RC–цепочкой при различных постоянных времени T. Видно, что с увеличением Т АЧХ становится немонотонной. Функция LAС() похожа на квадрат АЧХ параллельного колебательного контура. В данном случае аналогом резонансной частоты является p = T. С - 77 ростом Т величина пика АЧХ увеличивается, однако уменьшается полоса пропускания и возрастает скорость спада на больших частотах. В результате фильтр лучше пропускает составляющие шума на частотах, близких к p, но лучше подавляет высокочастотный шум.

1 LAC() 2а 2б 2в, рад/с Рис. 2.26. Графики LAС() для бесфильтровой СФС (1) и СФС с RC–цепочкой (2) для = 30 рад/с для а) Т = 10 / ;

б) Т = 5 / ;

в) Т = 1 /.

Поскольку LAС() представляет собой квадрат АЧХ ФНЧ, то можно, как и в случае бесфильтровой системы, рассмотреть величину шумовой полосы. На рис. 2.27 показана зависимость шумовой полосы от параметра T. Так как при T = 0 система становится бесфильтровой, использование RC–цепочки в качестве фильтра в цепи управления увеличивает шумовую полосу СФС. Анализ функции LBC(), отвечающей за прохождение на выход собственных фазовых шумов перестраиваемых генераторов в цепи, дал следующие результаты. На рис. 2.28 показаны графики этой функции для бесфильтровой СФС и для системы с RC–цепочкой при различных постоянных времени. С увеличением постоянной времени Т частота среза фильтра уменьшается. В то же время усиливается пропускание на выход СФС низкочастотных составляющих фазового шума перестраиваемого генератора. Подавление высокочастотных составляющих данного шума при этом несколько - 78 увеличивается. Однако, при любом T, значение LBС() для системы с RC– цепочкой больше значения этой функции для бесфильтровой системы.

Пш, рад/с T, с Рис. 2.27. Зависимость шумовой полосы системы, описываемой LAC(), от постоянной времени фильтра T.

LBC() 2а 2б 2в, рад/с Рис. 2.28. Графики LBС() для бесфильтровой СФС (1) и СФС с RC–цепочкой (2) для = 30 рад/с для а) Т = 0.1 / ;

б) Т = 0.3 / ;

в) Т = 0.5 /.

- 79 В результате сделаны следующие выводы: 1. Одно звено (одна СФС) лучше пропускает на выход низкочастотные составляющие фазового шума на входе, но больше подавляет высокочастотные составляющие этого процесса по сравнению с бесфильтровой системой. 2. Одно звено при любом выборе параметра фильтра Т пропускает больше фазовых шумов перестраиваемого генератора на выход, чем бесфильтровая система с той же полосой удержания. Данные выводы подтверждаются полученными результатами компьютерного анализа. На рис. 2.29 показаны зависимости нормированной дисперсии фазовых флуктуаций на выходе СФС с RC–фильтром от постоянной времени фильтра. Шумовые воздействия в данном случае имеют следующие параметры: a0 = 149 рад·с;

b0 = 1 рад·c;

b1 = 30000 рад2.

T, с T, с а) б) Рис. 2.29. Зависимость нормированной дисперсии фазового шума на выходе одного звена от постоянной времени RC–фильтра для а) a1 = 20000 рад2;

б) a1 = 10000 рад2. Полоса удержания была выбрана оптимальной для бесфильтровой системы при таких же входных воздействиях. Как видно, нулевое значение постоянной времени Т является наилучшим, что соответствует бесфильтровой системе.

- 80 в) Анализ модели с фильтром с пропорциональным каналом Было проанализировано поведение звена с фильтром в цепи управления системы фазовой синхронизации, коэффициент передачи которого имеет вид:

K ( p) = 1+ m T p. В случае если 01. Для создания такого фильтра можно использовать RC– цепочку, в комбинации с усилителями. На рис. 2.30 представлены зависимости LAC() для бесфильтровой СФС и СФС с ПИФ при различном выборе параметров. Как видно, СФС с данным типом фильтра обладает большим подавлением низкочастотных составляющих фазового шума на входе (в точке А). В то же время она больше пропускает высокочастотные составляющие этого шума. Увеличение параметров T и m приводит к усилению этого эффекта. Однако рост этих параметров приводит к увеличению шумовой полосы системы. Соответствующие графики показаны на рис. 2.31. Рост величины шумовой полосы объясняется увеличением значения модуля коэффициента пропускания в области высоких частот. При расчете Пш учитывается весь бесконечный диапазон частот, в то время как в реальной системе верхняя частота ограничена определенными факторами (частотой дискретизации и др.). На рис. 2.32 приведены графики для функций LВC(), отвечающих за передачу на выход фазовых шумов перестраиваемых генераторов. Из рисунков видно, что при данном выборе параметров фильтра обеспечивается большее подавление фазовых шумов перестраиваемого генератора звена.

- 81 1 L A C ( ) L A C ( ) 2I 2II 2III 2I 2II 1 2III, рад/с, рад/с а) б) Рис. 2.30. Графики LAC() бесфильтровой СФС (1) и СФС с ПИФ (2) при = 30 рад/с для а) m = 2, Т = (I – 0.5/, II – 0.3/, III – 0.1/);

б) T = 0.3/, m = (I – 6, II – 4, III – 2).

Пш, рад/с Пш, рад/с T, с а) б) Рис. 2.31. Графики зависимости шумовой полосы СФС с ПИФ от параметров фильтра при = 30 рад/с для а) m = 2;

б) T = 0.3/.

m - 82 1 L В C ( ) L В C ( ) 2I 2II 2III 2I 2II 2III, рад/с, рад/с а) б) Рис. 2.32. Графики LВC() бесфильтровой СФС (1) и СФС с ПИФ (2) при = 30 рад/с для а) m = 2, Т = (I – 0.5/, II – 0.3/, III – 0.1/);

б) T = 0.3/, m = (I – 6, II – 4, III – 2). В результате исследования системы с фильтром с пропорциональным каналом в цепи управления сделаны следующие выводы: 1. Одно звено (одна СФС) больше подавляет низкочастотные составляющие фазового шума на входе, но пропускает больше высокочастотных составляющих этого процесса по сравнению с бесфильтровой системой. 2. Одно звено при любом выборе параметров фильтра Т и m лучше подавляет фазовые шумы перестраиваемого генератора, чем бесфильтровая система с той же полосой удержания. Результаты компьютерного анализа [85] представлены на следующих рисунках. На рис. 2.33 показана зависимость нормированной дисперсии фазовой ошибки от параметра T2 = m·T для различных значений Т. Как видно, оптимальное значение Т2 увеличивается с ростом Т, поскольку параметр m должен быть больше 1. С другой стороны, по-видимому, существует оптимальное значение параметра T, при котором выигрыш по сравнению с бесфильтровой системой будет максимален.

- 83 T2, с T2, с а) б) Рис. 2.33. Зависимость нормированной дисперсии фазового шума на выходе одного звена от параметра Т2 при а) T = 0.3/;

б) T = 0.5/.

Как видно из рис. 2.34 значение параметра m тоже изменяется. Причем чем меньше T, тем больше становится оптимальное значение m. Это вызвано необходимостью увеличить подавление низкочастотных составляющих шумов, поскольку уменьшение параметра Т усиливает их пропускание.

m m а) б) Рис. 2.34. Зависимость нормированной дисперсии фазового шума на выходе одного звена от параметра m при а) T = 0.1/;

б) T = 0.5/.

- 84 2.4.2. Применение дискретной модели Рассмотрим поведение линеаризованных дискретных звеньев при наличии комплексного воздействия на них. а) Анализ бесфильтровой модели В данном случае система описывается передаточными функциями для входных флуктуационных воздействий, имеющими вид (1.5.2). Анализ поведения данного звена будем проводить из сравнения с аналоговой бесфильтровой моделью. На рис. 2.35 проведено сравнение функций LAC и LBC для аналоговой и дискретной систем.

= 0.9 = 0.5 = 0. = 0. = 0. = 0. L BC () L АC () а) б) Рис. 2.35. Графики функций LAC() и LBC() для аналогового (толстые линии) и дискретного (тонкие линии) бесфильтровых звеньев. Здесь параметр = E·S·T. Как видно из рисунка, чем ближе к 0 параметр, тем меньше различие между графиками. В то же время на любой частоте графики для дискретной системы проходят выше графиков аналоговой СФС. Это означает, что дискретная система пропускает больше шума, чем эквивалентная ей аналоговая система в той же полосе частот. При этом существуют некоторые отличия в зависимости данных функций от, по сравнению с аналоговым случаем. Увеличение данного параметра приводит к увеличению пропускания входных шумов (в точке А). Однако сильно растет - 85 пропускание высокочастотных составляющих этого шума. Аналогичный эффект наблюдается и для функции LBC(). С ростом усиливается подавление низкочастотных составляющих собственных фазовых шумов перестраиваемого генератора. Но имеется усиление высокочастотных составляющих данных шумов, чего не наблюдалось в аналоговом случае. Анализ показывает, что 2 LBC ( ) =. В результате максимальное усиление спектральной плотности фазовых шумов перестраиваемого генератора может составлять 4 раза. На основе анализа дискретной бесфильтровой системы сделаны следующие выводы: 1. При стремлении параметра к 0 (случай, когда частота дискретизации существенно больше полосы удержания СФС) поведение дискретной модели становится аналогичным поведению аналоговой бесфильтровой модели. 2. При любом значении параметра дискретная система фазовой синхронизации пропускает больше шума, чем эквивалентная ей аналоговая система. В качестве подтверждения данных выводов на рис. 2.36 приведены графики зависимости нормированной дисперсии фазовой ошибки от величины параметра. Как видно из представленных кривых, значение дисперсии фазовых флуктуаций на выходе дискретной ячейки больше, чем на выходе аналоговой. Причем для больших значений параметра данный эффект сильнее. Здесь параметры флуктуационных воздействий: a0 = 140 рад·с;

b0 = 1 рад·с;

b1 = 30000 рад2. В то же время качественный характер зависимостей нормированной дисперсии фазовых флуктуаций на выходе звена от его параметров остается прежним. Так на рис. 2.37 показаны графики зависимости этой величины от для шумовых воздействий 2-го порядка. Здесь a0 = 140 рад·с;

a1 = 20000 рад2;

b0 = 1 рад·с;

b1 = 30000 рад2;

a2 = 50000 рад3 / с. Аналогичные графики для - 86 аналоговых бесфильтровых звеньев показаны на рис. 2.23. Как видно из сравнения, характер кривых не изменился.

, рад/с, рад/с а) б) Рис. 2.36. Зависимость нормированной дисперсии фазового шума на выходе звена от для аналоговой (1) и дискретной (2) ячеек при а) a1 = 20000 рад2;

б) a1 = 10000 рад2;

, рад/с, рад/с а) б) Рис. 2.37. Зависимость нормированной дисперсии фазового шума на выходе звена от полосы удержания для б) b2 = 11000 рад3 / c. а) b2 = 29000 рад3 / c;

Однако результаты могут иметь и качественно иной характер. На рис. 2.38 показаны графики зависимости нормированной дисперсии фазового шума на - 87 выходе звена от полосы удержания при наличии только белых фазовых шумов генераторов.

, рад/с, рад/с а) б) Рис. 2.38. Зависимость нормированной дисперсии фазовых флуктуаций на выходе одного звена от при а) a0 = 1 рад с;

b0 = 100 рад с;

б) a0 = 100 рад с;

b0 = 1 рад с.

, рад/с, рад/с а) б) Рис. 2.39. Зависимость нормированной дисперсии фазовых флуктуаций на выходе одного звена от при а) b1 = 10000 рад2;

б) b1 = 5000 рад2.

На рис. 2.39 показана зависимость дисперсии фазовых флуктуаций от при наличии только фазовых шумов перестраиваемого генератора, которые - 88 имеют 1-ый порядок. Здесь b0 = 100 рад с. Объяснение такому поведению системы дано в конце раздела 2.2 и связано с характером зависимости АЧХ системы от величины полосы удержания. б) Анализ модели с фильтром с пропорциональным каналом В работе рассмотрено поведение звена с ПИ–фильтром в цепи управления. Коэффициент передачи такого фильтра в Z–области имеет вид:

K ( z) = m + 1. zd В данном случае математическая модель данного звена имеет вид, представленный выражением (1.5.3). При фиксированных и d изменение параметра m может привести к качественно различным ситуациям, показанным на рис. 2.40. При больших значениях m LAC() для системы с ПИФ на всех частотах больше аналогичной функции системы без ПИФ. Это означает, что такая дискретная СФС пропускает больше входных шумов, чем бесфильтровая дискретная СФС. В то же время пропускается меньше низкочастотных составляющих шумов перестраиваемого генератора. Высокочастотные составляющие этих шумов пропускаются лучше. Уменьшение параметра m приводит к тому, что в области высоких частот ( близка к ) система с ПИФ начинает подавлять фазовый шум входного сигнала лучше, чем в бесфильтровая система. При этом увеличивается пропускание низкочастотных составляющих шумов перестраиваемого генератора, а высокочастотные составляющие этих шумов подавляются лучше, чем в бесфильтровом случае. Наконец, при отрицательных m возможен вариант, когда функция LAC() для системы с ПИФ на всех частотах меньше этой функции для бесфильтровой дискретной СФС. Но лучше пропускаются низкочастотные составляющие шумов перестраиваемого генератора.

- 89 1 2 L B C ( ) ) L А C ( ) ) 1 а) 1 L BC () L АC () 2 1 б) 1 L BC () L АC () 2 1 в) Рис. 2.40. Графики функций LAC() и LBC() для бесфильтровой дискретной СФС (1) и дискретной СФС с ПИФ (2) при = 0.1, d = 0.1 для а) m = 4;

б) m = 1;

в) m = –0.8.

- 90 Исследования показали, что существуют оптимальные значения параметра m. На рис. 2.41 показаны графики зависимости нормированной дисперсии фазовых флуктуаций на выходе одного звена от величины m. Здесь параметры шумовых воздействий: a0 = 140 рад·с;

a2 = a3 = 0;

b0 = 1 рад·с;

b2 = b3 = 0. Величины в обоих случаях выбирались оптимальными для бесфильтровых систем при данных внешних воздействиях. Параметр фильтра d равен 0.

m m а) б) m в) Рис 2.41. Зависимость нормированной дисперсии фазового шума на выходе звена от m при 2 б) a1 = 10000 рад2;

b1 = 30000 рад2;

а) a1 = 20000 рад ;

b1 = 30000 рад2;

в) a1 = 20000 рад2;

b1 = 25000 рад2. Как видно из рисунка, изменение параметров входных воздействий практически не влияет на оптимальное значение m. Это обусловлено тем, что значения различны для каждого из экспериментов. Изменение само по себе - 91 компенсирует изменения во входных воздействиях. В то же время при фиксированном значении полосы удержания наблюдается изменение оптимального значения m (рис. 2.42).

2 m m а) б) m в) Рис. 2.42. Зависимость нормированной дисперсии фазового шума на выходе звена от m при = 43.2 рад/с для а) a1 = 20000 рад ;

b1 = 30000 рад2;

б) a1 = 10000 рад2;

b1 = 30000 рад2;

в) a1 = 20000 рад2;

b1 = 25000 рад2. Увеличение фликкер–составляющей фазовых флуктуаций входного сигнала ведет к уменьшению m, поскольку это позволяет лучше отфильтровать низкочастотные составляющие этого шума. Аналогичный эффект вызывает уменьшение фликкер–составляющей фазовых флуктуаций сигнала генератора ячейки, поскольку уменьшение m вызывает большее пропускание низкочастотных составляющих этого сигнала.

- 92 Существенным отличием от случая аналогового ПИ–фильтра является следующее. Изменением параметров аналогового ПИ–фильтра в цепи управления ячейки можно было добиться одновременного усиления подавления низкочастотных составляющих фазовых шумов сигнала на входе ячейки и сигнала перестраиваемого генератора. В случае дискретного ПИ– фильтра уменьшение m подавляет низкочастотные составляющие фазовых шумов входного сигнала, но увеличивает пропускание аналогичных составляющих фазовых шумов перестраиваемого генератора. В связи с этим в проведенных экспериментах цепочка из звеньев с ПИФ показала худшее качество работы по сравнению с цепочкой бесфильтровых звеньев. 2.5. Выводы На основании результатов, полученных в данной главе, сделаны следующие выводы: 1. В случае если все фазовые флуктуации в системе являются белыми шумами, оптимальным значением полосы удержания в дискретной СФС будет нулевое значение. В аналоговой СФС оптимальным будет либо нулевое, либо бесконечное значение. 2. В случае воздействия на СФС белых частотных шумов входного сигнала и/или сигнала перестраиваемого генератора может иметься оптимальное значение параметров системы, при котором дисперсия фазовой ошибки будет минимальна. 3. В случае воздействия на СФС белых частотных шумов входного сигнала и сигнала перестраиваемого генератора может иметься оптимальное значение, при котором будет наблюдаться минимум дисперсии фазовых флуктуаций на выходе системы. Аналитически получить связь между параметрами шумов и оптимальным значением полосы удержания не удается, поскольку фазовые флуктуации на выходе системы являются нестационарным процессом. Их дисперсия со временем увеличивается. Однако, несмотря на это, значение, при котором наблюдается минимум этой дисперсии в каждый момент времени, практически не изменяется. 4. Воздействия фазовых флуктуаций входного сигнала и сигнала перестраиваемого генератора на фазовую ошибку эквивалентны. То есть их - 93 замена друг на друга не приводит к изменению характеристик фазовой ошибки. Для дисперсии фазовых флуктуаций на выходе СФС это не так. 5. Рост мощности канальных шумов приводит к уменьшению оптимального значения полосы удержания, как для фазовой ошибки, так и для фазовых флуктуаций выходного сигнала. 6. Рост мощности фазовых флуктуаций сигнала перестраиваемого генератора приводит к увеличению оптимального значения полосы удержания, как для фазовой ошибки, так и для фазовых флуктуаций выходного сигнала. 7. Рост мощности фазовых флуктуаций входного сигнала СФС приводит к увеличению оптимального значения полосы удержания для фазовой ошибки и к ее уменьшению для фазовых флуктуаций выходного сигнала. 8. При наличии комплексных флуктуационных воздействий на систему, спектральные плотности которых выражаются формулой (1.2.3), возможны сложные зависимости дисперсии фазовых флуктуаций на выходе СФС от ее параметров. Кривые этой зависимости могут иметь несколько экстремумов (как минимумов, так и максимумов), что подтверждает необходимость проведения процедуры оптимизации. 9. Большое влияние на значение оптимальных параметров оказывает время анализа. При малых временах анализа оптимальным оказывается нулевое значение полосы удержания. С ростом времени наблюдения эта величина может принять конечное, отличное от нуля значение. 10. Аналоговая СФС с RC–цепочкой показывает худшие результаты по сравнению с бесфильтровой системой в случае рассматриваемого воздействия. Это вызвано тем, что данная система пропускает больше фазовых шумов перестраиваемого генератора и больше низкочастотных составляющих фазовых шумов входного сигнала и канала передачи. Поскольку именно здесь сосредоточена основная мощность шумов, то данная система проигрывает бесфильтровой СФС. 11. Аналоговая СФС с ПИ–фильтром может иметь выигрыш перед бесфильтровой СФС при условии, что m>1. Такая система лучше подавляет основные источники шумов. 12. Дискретная СФС всегда пропускает больше шумов, чем эквивалентная ей в смысле параметров аналоговая система. В то же время если частота дискретизации существенно превышает полосу удержания, то поведение дискретной СФС мало отличается от поведения аналоговой СФС и - 94 характеристики качества их работы близки. С ростом качество работы дискретной СФС по сравнению с аналоговой СФС падает. 13. Введение дискретного ПИ–фильтра в дискретной СФС позволяет либо лучше фильтровать низкочастотные составляющие фазовых шумов перестраиваемого генератора, либо низкочастотные составляющие шумов канала передачи и фазовых шумов входного сигнала. Но одновременно это сделать не удается. Уменьшение влияния одних шумов приводит к увеличению влияния других. Поэтому введение такого фильтра не позволяет получить превосходства над бесфильтровой системой, какое давало введение ПИФ в аналоговом случае.

- 95 Глава 3. Анализ и оптимизация статистических характеристик цепи последовательно синхронизируемых генераторов 3.1. Постановка задачи В главе исследуется ряд проблем, связанных с исследованием и параметрической оптимизацией цепи последовательно синхронизируемых с помощью дискретных СФС генераторов. По аналогии со 2-ой главой, посвященной изучению одного звена в условиях комбинированных случайных воздействий, будет рассмотрен общий случай нелинейной цепи без ограничения на уровень входных воздействий и соответственно диапазон изменения координат системы. Для этого исследования будет применена математическая модель цепи в виде комбинации расширенных векторных уравнений Колмогорова-Чепмена и уравнений преобразования координат. Подобная модель особенно эффективна для цепей, состоящих из небольшого числа генераторов. Для ограниченных входных воздействий и близких к синхронному режимов цепи, допускающих применение линейных моделей, выполнены исследования для произвольной длины цепи. Заметим, что задача исследования и оптимизации цепей произвольной длины имеет важное значение для специалистов в области синхронных цифровых сетей передачи данных. Можно утверждать, что на сегодняшний день тщательно изучен вопрос о допустимых уровнях фазовых флуктуаций сигналов синхронизации [91, 93–102], выработаны рекомендации, указывающие допустимое количество последовательно соединяемых элементов цепи, маски на допустимые уровни шумов [97, 106], рассмотрен ряд других проблем надежности функционирования цепи [99, 114]. В то же время имеющиеся на сегодняшний день алгоритмы управления сетью синхронизации основаны на метках качества сигналов, которые передаются вместе с информационным потоком. Эти метки только весьма условно передают информацию об уровне качества сигналов синхронизации. Вместе с тем, очевидно, что более эффективно можно было бы управлять сетью на основе знаний о реальном качестве этих сигналов. Кроме того, полученные нормы на фазовые дрожания, на число последовательно соединяемых генераторов были получены эмпирическим путем. Недостаточно исследований посвящено процессам, происходящим внутри цепи последовательно соединенных генераторов и определяющим качество ее работы [92, 100, 102].

- 96 Отсутствие модели цепи тактовой синхронизации не позволяет выявить основные закономерности работы этой системы, влияние параметров ее ячеек на общее качество сигнала на выходе. Дополнительную сложность в задачу вносят специфические воздействия на систему. При анализе цепи последовательно синхронизируемых генераторов нельзя ограничиваться случаем широкополосных флуктуационных воздействий. В связи с вышесказанным целью данной главы является анализ модели цепи последовательно соединенных генераторов на основе дискретных систем фазовой синхронизации с учетом комбинированных флуктуационных воздействий. Исследуемой величиной является дисперсия фазовых флуктуаций на выходе системы. Рассмотрены также различные алгоритмы оптимизации цепи. В качестве критерия оптимальности выбран минимум указанной дисперсии. Часть результатов исследований, представленных в главе автор диссертации опубликовал в работах [75, 77, 85, 87]. В них, в частности, достаточно подробно рассматриваются вопросы исследования и оптимизации цепи в условиях ограниченных воздействий. Ниже в систематизированном виде представлены результаты исследований для произвольных воздействий. 3.2. Анализ и оптимизация цепи последовательно синхронизируемых генераторов с помощью аппарата марковских процессов В разделе выполнен анализ цепи, состоящей из последовательно синхронизируемых генераторов, на основе методики, описанной в первой главе. Математической основой для исследования является система уравнений (1.4.5), позволяющая получить ПРВ фазовых и частотных флуктуаций сигнала на выходе любого звена цепи. Предполагается, что частотные флуктуации на выходе любого звена являются гауссовскими. Это допущение позволяет записать уравнения, содержащиеся в (1.4.5) в явном виде с точностью до некоторых параметров (дисперсии частотных флуктуаций). На рис. 3.1–3.5 показаны ПРВ фазовых и частотных флуктуаций на выходе первого звена цепи, полученные на основе решения (1.4.5) для различных воздействий и параметров звена. При расчетах предполагалось, что спектральная плотность входных фазовых флуктуаций имеет квадратичный вид, а спектральная плотность частотных флуктуаций соответствует белому шуму.

- 97 На рис 3.1 приведено семейство ПРВ частотных флуктуаций выходного сигнала бесфильтрового звена для различных усилений S и параметров частотных и аддитивного воздействий (расширенное уравнение Колмогорова– Чепмена в (1.4.5) имеет 2-ой порядок). С хорошей степенью достоверности можно утверждать, что ПРВ имеет вид, близкий к гауссовскому. На рис 3.2–3.5 приведены семейства ПРВ фазовых флуктуаций на выходе звена 2-го порядка (расширенное уравнение Колмогорова–Чепмена в (1.4.5) имеет 3-ий порядок). Кривые получены для различных усилений, параметров фильтра в цепи управления, значений СКО входных частотных и аддитивного воздействий. Аналогично бесфильтровому звену следует отметить гауссовский характер всех кривых ПРВ.

p() S=0.63 0.95 1. p() S=0.63 0.95 1. а) б) Рис. 3.1. ПРВ частотных флуктуаций выходного сигнала дискретной СФС при вх = 0.1, n = 0.1 для а) пг = 0.4;

б) пг = 0.7.

- 98 p(пг) S=0.05 0.40 0. p(пг) S=0.13 0.93 1. пг пг а) б) Рис. 3.2. ПРВ фазовых флуктуаций на выходе звена 2-го порядка при d = 1, вх = 0.1, n = 0.1, пг = 0.1 для а) m = 4;

б) m = 2.

пг=0. 0.4 0. p(пг) p(пг) пг=0. 0.4 0. пг пг а) б) Рис. 3.3. ПРВ фазовых флуктуаций на выходе звена 2-го порядка при d = 1, вх = 0.1, n = 0.1 для а) m = 4, S = 0,29;

б) m = 2, S = 0.66.

- 99 p(пг) m=6 4 p(пг) m=6 4 пг пг а) б) Рис. 3.4. ПРВ фазовых флуктуаций на выходе звена 2-го порядка при d = 1, вх = 0.1, n = 0.1, S = 0.13 для б) пг = 0.7. а) пг = 0.4;

p(пг) d=0.2 0.6 p(пг) d=0.2 0.6 пг пг а) б) Рис. 3.5. ПРВ фазовых флуктуаций на выходе звена 2-го порядка при вх = 0.1, n = 0.1, пг = 0.4, S = 0.3 для а) m = 4;

б) m = 6. На рис. 3.6 приведены зависимости дисперсий фазовых и частотных флуктуаций сигнала на выходе звена от усиления в кольце.

- 100 пг = 0. 0. пг = 0. 0. S S а) б) Рис. 3.6. Зависимость дисперсии фазовых а) и частотных б) флуктуаций сигнала на выходе дискретной СФС от параметра S при вх = 0.1, n = 0.1,. Видно, что с ростом усиления в кольце мощность частотных флуктуаций на выходе увеличивается. Данный факт объясняется следующим образом. Зависимость дисперсии частотных флуктуаций сигнала на выходе СФС от параметров системы при наличии белого частотного шума на входе такая же, как зависимость дисперсии фазовых флуктуаций выходного сигнала от этих параметров при наличии входного белого фазового шума. В главе 2 подробно рассмотрен и объяснен характер последней зависимости. Показано, что с ростом S увеличивается дисперсия фазовых флуктуаций выходного сигнала. По аналогичным причинам происходит и увеличение дисперсии частотных флуктуаций. На рис. 3.7 приведены результаты оптимизации цепи последовательно синхронизируемых генераторов, выполненной с помощью гомогенного алгоритма. Цифрами обозначены длины цепочки. Результаты получены на основе марковской модели (1.4.5) и блок–схемы алгоритма оптимизации, приведенной на рис. 3.8. Предполагается, что все ячейки имеют одинаковые параметры, параметры аддитивных воздействий и собственных частотных и фазовых шумов также одинаковы. Гомогенный алгоритм характеризуется тем, что всем звеньям цепи дается одинаковое значение настраиваемого параметра. Таким образом производится настройка сразу всех звеньев.

- 101 2 а) S S б) Рис. 3.7. Зависимость дисперсии фазовых флуктуаций выходного сигнала цепи от параметра S звеньев при гомогенной оптимизации а) вх = 0.07, n = 0.07, пг = 0.28;

б) вх = 0.1, n = 0.1, пг = 0.4.

- 102 Задать значения параметров флуктуационных воздействий: ог, пгi, ni. Установить значение полосы удержания у всех ячеек цепи равным у.

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

Изменить у.

Нужно проводить расчет для другого у? Нет Из всех использованных у выбрать то, для которого дисперсия фазовых флуктуаций на выходе цепи 2 минимальна.

Да Установить значение полосы удержания у всех звеньев цепи равной полученной величине у опт.

Закончить обработку Рис. 3.8. Блок–схема гомогенного алгоритма оптимизации цепи последовательно соединенных генераторов. Из рисунков видно, что оптимальное значение полосы удержания быстро падает с ростом длины цепочки. Это объясняется эффектом накопления шумов в системе. На входе следующего звена дисперсия частотных флуктуаций больше, чем на входе предыдущего. Как уже говорилось, с ростом S значение - 103 2 увеличивается. В результате при больших значениях полосы удержания рост дисперсии частотных флуктуаций происходит быстрее, что видно на обоих графиках. Однако при больших мощностях шумов начинает существенно сказываться влияние нелинейности, которое приводит к ограничению мощности фазовых флуктуаций выходного сигнала. Этот эффект ясно прослеживается на графиках для цепи из 3-х элементов. При малых же S имеет место другой эффект. Поскольку при узкой полосе удержания входной фазовый шум отфильтровывается хорошо, то эффект его накопления сказывается минимально. Кроме того, очень быстро система приходит к некоторому балансу между накоплением шумов и их подавлением. Особенно хорошо этот эффект выражен на рис. 3.7а. На рис. 3.9 приведены результаты оптимизации цепи последовательно синхронизируемых генераторов, выполненной с помощью позвенного алгоритма. Цифрами на графике обозначены длины цепи.

S Рис. 3.9. Зависимость дисперсии фазовых флуктуаций выходного сигнала цепи от параметра S звеньев при позвенной оптимизации вх = 0.1, n = 0.1, пг = 0.4.

- 104 Задать значения параметров флуктуационных воздействий: ог, пгi, ni. Установить значение полосы удержания у текущей ячейки цепи равным у.

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

Изменить у.

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

Да Нет Установить значение полосы удержания у текущей ячейки цепи равным полученной величине у опт.

Pages:     || 2 | 3 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.