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

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

КЫТИН Евгений Алексеевич

УДК 621.396.2

005009802

ОПРЕДЕЛЕНИЕ СПЕКТРАЛЬНЫХ ХАРАКТЕРИСТИК ПОЛОСОВЫХ ФИЛЬТРОВ ДЛЯ АНАЛИЗА КАЧЕСТВА НЕОДНОРОДНОГО КАНАЛА СВЯЗИ

Специальности:

05.13.01 - Системный анализ, управление и обработка информации (в науке и технике)

05.12.13 - Системы, сети и устройства телекоммуникаций

1 В ФЕВ 20Ї2

Автореферат диссертации на соискание ученой степени кандидата технических наук

Ижевск-2012

005009802

Работа выполнена на кафедре «Конструирование радиоэлектронной аппаратуры» в ФГБОУ ВПО «Ижевский государственный технический университет» (ИжГТУ, г. Ижевск)

Научный руководитель:

заслуженный изобретатель РФ,

доктор технических наук, профессор Лялин В.Е.

(ФГБОУ ВПО «ИжГТУ»)

Официальные оппоненты:

доктор технических наук, профессор Костарев С.Н.

(ФГБОУ ВПО «Пермский национальный политехнический университет»)

доктор технических наук, профессор НистюкА.И.

(ФГБОУ ВПО «ИжГТУ»)

Ведущая организация:

ФГБОУ ВПО «Вятский государственный университет»

Защита состоится 01 марта 2012 г. в 16:00 часов на заседании диссертационного совета Д 212.065.06 в ИжГТУ по адресу: 426069, г. Ижевск, ул. Студенческая, 7.

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

С диссертацией можно ознакомиться в библиотеке университета. С авторефератом можно ознакомиться на официальном сайте Министерства образования и науки РФ http://mon.gov.ru.

Автореферат разослан 31 января 2012 г. •

Ученый секретарь диссертационного совета, У ✓

кандидат технических наук, доцент

(/ / В.Н. Сяктерев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Одной из основных задач при повышении качества работы информационной радиотехнической системы является определение по-меховой обстановки в реальных рабочих условиях. Для этого необходимо определить минимальное число параллельных фильтров, которые позволили бы проанализировать зашумленность как в поддиапазонах, так и по всей полосе канала связи (КС). При этом количество фильтров такого анализатора должно удовлетворять заданной точности оценивания зашумленности канала. Нужна методика определения необходимого числа полосовых фильтров для анализатора качества канала связи при заданной точности измерений в канале связи.

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

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

Для решения вышеуказанных задач необходим спектральный анализ, который является мощным средством исследования сигналов. Методы спектрального анализа в настоящее время широко применяются в самых разных областях науки и практики. Для обеспечения анализа и статистической обработки проведено огромное количество работ видными учеными разных стран. Большой вклад в решение указанной проблемы внесли ученые Бокс Дж., Дженкинс Г., Бриллинджер Д., Вате В., Макс Ж., Робинсон Э.А., Отнес Р., Эноксон Л., Ста-

тулявичус В.В., Бенткус Р.Ю., Кэдзоу Дж.А., Поспелов Г.С., Поспелов Д.А., Яглом А.Н., Журбленко И.Г., Ибрагимов И.А., Линник Ю.В,

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

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

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

Известно, что параметры ленточного носителя, такие как жесткость и вязкость, в зависимости от длины ленты в бобине изменяются по случайному закону. Из этого следует, что при транспортировании ленты в каждый момент времени формируется вектор мгновенных собственных частот (СЧ), который отличается от векторов СЧ, сформированных в предыдущие моменты времени. Отклонения СЧ от номинальных, определяемых детерминированными параметрами МТЛ, синтезированными по существующим критериям и решающим функциям, могут привести к тому, что одно или несколько СЧ, располагаясь вблизи границ резонансно-опасной зоны (РОЗ), может зайти в запретные области, повышая уровень колебаний звеньев. С целью упреждения этого факта необходимо откорректировать значения границ РОЗ в соответствии с оценками дисперсии, закона распределения и асимптотики разложения оценки отклонений СЧ.

В связи с вышеизложенным тема настоящего диссертационного исследования весьма актуальна.

Объектом исследования являются: анализ и имитационное моделирование состояния неоднородного КС.

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

Цель работы состоит в разработке и научном обосновании математиче-

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

Для достижения поставленной цели требуется решить следующие задачи:

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

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

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

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

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

- определить условия инвариантности МТЛ стримеров по отношению к возмущающим факторам для поддержания заданных точностных характеристик информационной части записываемого сигнала при условии рассмотрения стримера как аналога или составной части канала передачи информации;

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

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

Методы исследования. В работе применялись теоретические и экспериментальные методы исследования.

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

Теоретические исследования базируются на использовании методов ста-

б

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

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

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

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

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

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

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

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

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

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

цесса, удовлетворяющего разностному уравнению АРСУ;

- анализ условий инвариантности МТЛ стримеров по отношению к паразитным колебаниям для поддержания заданных точностных характеристик информационной части регистрируемого сигнала;

- установление инвариантности к помехам по допустимой вероятности ошибки чтения-записи, при условии рассмотрения МТЛ стримеров как аналога или составной части канала передачи информации;

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

- рассмотрение колебаний (сигналов), зависящих от фактора времени в стримерах как инвариантные множества, поскольку МТЛ является динамической системой;

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

- методика коррекции границ РОЗ с целью ликвидации возможности захода СЧ в РОЗ в течение всего времени транспортирования носителя, построенная на основе изучения статистических характеристик отклонений СЧ МТЛ.

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

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

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

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

- в качестве одной из возможных математических моделей представления сигнала в созданном пакете программ "СПЕКТР" выбрана модель стационарного случайного процесса, удовлетворяющего разностному уравнению АРСУ. Для данной модели решаются две задачи: оценка параметров среднего, диспер-

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

- предложен метод проекционных оценок, заключающийся в аппроксимации неизвестной плотности распределения некоторым отрезком ее ряда Фурье по подходящей системе функций; в пакете "СПЕКТР" алгоритм определения весовой последовательности, при которой проекционная оценка проста в реализации, требует небольших ресурсов компьютера;

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

- показано, что описание конкретных видов МТЛ периодическими функциями возможно при использовании полученных классическими методами АЧХ и ФЧХ. Компоненты таких периодических функций, осуществимых МТЛ, являются парами преобразования Гильберта, позволяющими осуществлять переход от АЧХ к ФЧХ и обратно; в том числе не исключены случаи, когда передаточная функция может быть недостаточно гладкой для интегрирования по Риману, поэтому проводить операции над ней в общем виде можно только при помощи интеграла Лебега;

- получены формулы для среднего оценок отклонений СЧ, являющегося критерием риска захода СЧ в РОЗ; определены точные и асимптотические выражения для дисперсии оценки отклонений СЧ, установлен ряд ее оценок сверху через реологические константы и матрицу ковариации изменений параметров ленты, а также интервалы их корреляции;

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

- на основе изучения статистических характеристик отклонений СЧ МТЛ предложена методика коррекции границ РОЗ с целью ликвидации возможности захода СЧ в РОЗ в течение всего времени транспортирования носителя.

Практическая ценность. Разработана методика формирования ширины полосы пропускания, учитывающая, что форма спектральной плотности сигналов на выходе каждого из фильтров близка к прямоугольной, а АЧХ имитатора канала не имеет в частотном диапазоне возбуждения резонансов и антирезонансов, то есть спектральная плотность моделирующих сигналов представляет собой кусочно-постоянную функцию, а ее варьируемыми переменными являются уровни спектральной плотности составляющих процессов, сформированных к-ым фильтром, и координаты границ частотных диапазонов £-го и (&+1)-го фильтров.

Пакет "СПЕКТР" отличается своим функциональным назначением: определение спектральных характеристик временных рядов. Он позволяет проводить спектральный анализ временных рядов сравнительно систематично, пользователю на выбор представлены практически все методы спектрального анализа. В пакет "СПЕКТР" заложен и ряд методов, таких как, оценки спектральной плотности с адаптивными ковариационными окнами, ряд новых окон данных, спектральных и ковариационных окон.

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

Указано, что при транспортировании ленты в каждый момент времени формируется вектор мгновенных собственных частот (СЧ), а отклонения СЧ от номинальных, определяемых детерминированными параметрами МТЛ, могут привести к тому, что одно или несколько СЧ, располагаясь вблизи границ РОЗ, может зайти в запретные области, повышая уровень колебаний звеньев (помех в КС). С целью упреждения этого факта предложен алгоритм корректировки значений границ РОЗ в соответствии с оценками дисперсии, закона распределения и асимптотики разложения оценки отклонений СЧ. Установлены достаточно точные и удобные для применения экспоненциальные оценки сверху для вероятности, что отклонения СЧ превысят заданный уровень. Изучена асимптотика распределения отклонений СЧ при неограниченно возрастающем времени транспортирования ленточного носителя, предложены формулы, корректирующие границы РОЗ.

Реализация работы в производственных условиях. Методики анализа и имитации частотных характеристик неоднородного канала связи, повышения динамической точности функционирования стримеров, как составной части канала передачи информации, а также созданные алгоритмы и комплексы программ для статистического анализа сигналов в КС переданы в ОАО «Ижевский радиозавод» для практического использования.

Апробация работы. Основные положения и результаты диссертации докладывались и обсуждались на российских и международных научно-технических конференциях и конгрессах: 7-й Международной конференции «Цифровая обработка сигналов и ее применение» (Москва, 2005); VI и VII Международной научнотехнической конференции «Проблемы техники и технологии телекоммуникаций» (Самара, 2006); Российской научно-технической конференции «Приборостроение в XXI веке. Интеграция науки, образования и производства» (Ижевск, 2006); Международном симпозиуме «Надежность и качество» (Пенза, 2007—2009); 35-й Международной конференции «Информационные технологии в науке, образовании, телекоммуникации и бизнесе» (Украина, Крым, Ялта-Гурзуф, 2008, 2009); V Международ-

ной научно-практической конференции «Перспективы развития информационных технологий» (Новосибирск, 2011).

Публикации. Результаты работы отражены в 14 научных трудах в региональных журналах, сборниках научных трудов и материалов конференций. Автор имеет 3 научных труда в изданиях, выпускаемых в РФ и рекомендуемых ВАКом для публикации основных результатов диссертаций.

Структура и объем работы. Диссертация содержит введение, 4 главы и заключение, изложенные на 179 с. машинописного текста. В работу включены 37 рис., 7 табл., список литературы из 127 наименований и приложение, в котором представлен акт об использовании результатов работы.

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

В первой главе проведен обзор принципов анализа и имитации работы КС. Рассмотрены методики определения спектральных характеристик полосовых фильтров, математическое и программное обеспечение для спектрального анализа сигналов, методы оценки случайных возмущений в МТЛ стримеров, как составной части канала передачи информации, и методики снижения их уровня. Исследованы основные принципы контроля состояния КС по информационным признакам сигнала. ; . ..■■■•

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

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

Предложенная методика основана на учете физических параметров фильтров и степени их согла- ' ~ "

СОДЕРЖАНИЕ РАБОТЫ

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

сованности.

Сутью ее является

-►

(0

Рис. 1. Квантование частотной оси <а переменным шагом 25^

ния числа фильтров п в диапазоне Дсо необходимо было найти алгоритм квантования частотной оси со переменным шагом 28, (рис. 1). Степень разделения спектральной плотности для каждой пары смежных фильтров характеризуется отношением ординаты стыка У0ц к ординате максимума)^*, т.е. У0-У^ <с, (1) что справедливо в случае независимости Ко, Ум и с от положения характеристики фильтра на частотной оси <а. При учете вклада в резонансную характеристику /-го фильтра характеристик только /-1 и /+1 при постоянном затухании сі фильтров выражение для Ум можно, в частности, для /=2 представить:

УМ2=(і~1 +[(1-0)2М'2)2 +Ю2(й['2^2]'05+[(1-0)2(Вз2)г +Ю2юз2^2Г°5- (2)

Проведены некоторые преобразования для упрощения выражения (2). При этом отмечено, что при Зі«г52~53 справедливы выражения <а0, =<о, -8,, со02 = ю2 - 5,, й)03 = <у3 - Л,, и введены обозначения а = 8, • ш^}, Ь = 8, • со^ • (3)

Тогда выражения для Ум и У0 принимают вид соответственно

Гм = <Г' +(1-а)2(4аГ' + (1 + За)2[4а(1 + 2а)]"', (4)

у (1~а) 10 ~

а

-а\2 (Iі

і * ?

1-а; а а

-0.5

, 0 + °)

2 + яУ сі2

-0.5

(5)

Учитывая, что Ъ < 1, то с12(\-а)2 « 1 и выражение (5) можно упростить:

У0 =(1-а)2И2-а)Г’ +(1+а)2[а(2 + а)Г1. (6)

Принимая во внимание соотношения (4) и (6) из (1), получим приближенную формулу для определения параметра а: а - й'(1-0,5с)(с-2й()" . А если с < 0,5, с1 "* > 20, то а »<1 • с'1. Тогда из (3) следует, что шаг квантования частотной оси определяется из выражения 28, « 2©,„ • с/ • с", / = 1,и (7)

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

П

рительно заметив, что со0(л+о =шо1 +28! + ... + 28„; Дсо = со0(л+1) - со01 = 2^8,.

При учете,(7) Дсо = 2^] со0/-с/-сч. (8)

/=1

Очевидно, что со0, = со0((_!) + 25,_, = о)0(,_1)(1 +2с/ • с"1). (9)

При учете (9) выражение (8) принимает вид:

Дсо = со0, (1 +2с?-сГ1 У -со01. (Ю)

Из (10) определим число фильтров: п = ^(Дсо-сйо' +1) 1в-1 (1 + 2с1-с~1).

Если граничные частоты со01 и со0(„-ц) рассматриваемого частотного диапазона Дсо совпадают с резонансными частотами фильтров Ю| и со„, то выражение для определения п приобретает следующий вид:

п = 2 +1§-1 (1 + 2с/ • с-1) ■ {^[(Дсо • Юо’) +1 - (Л ■ с~')] -1§[1 + (с/ • с-1)]}. Оценка ошибки ДС(со) воспроизведения спектральной плотности С(а>) получена следующим образом. Верхней границей ошибки ДС(со) может слу-

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

г "

переданного набором п фильтров: є = Ум ■ Аш - У^со, со є Дш, (11)

где У, = [(1 - со2а> 2

У +ш2ш/"2^2]’05. Для определения є огибающую резонанс-

ных характеристик п фильтров по диапазону Дм аппроксимируем ломаной

(рис. 2). Тогда Е = Х(ГМ -Г0)-28,. 12 = (Ум -У0)Дсо/2.

(12)

/=і

Принимая во внимание (4) и (5), выражение (12) запишем в виде: є = [2(1 + 2 а)(а - сІ) + <і(1 + а)3] ■ [Ш{\ + 2а)]'1 • Дш

(13)

Оценка ошибки воспроизведения спектральной плотности получена на основании соотношения:

АС?(со) =| Д5(ш) |2 (яГ)-1.

Причем Д5(со) можно определить из (11), тогда Д5(ш)=в-Д(о“1 =

:УМ~ ДШ'

гш.

Учитывая (13), получим:

Рис. 2. Аппроксимация огибающей набора резонансных характеристик

ДС(со) = {[2(1 + 2 а)(а - </) + </( 1 + + 2а)]'1 }2(пТ)

-\

(14)

Как следует из (14), ДС’(ш) является функцией добротности, значения отношения Уо-Уц/] на частоте стыка резонансных характеристик и длительности реализации Т, т.е. ДС(ш) = /(с!,с,т).

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

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

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

(15)

спектральных плотностей 5^) и Бх(у), моделируемых случайных колебаний

N

и воспроизводимых в канале сигналов *(/)= ^х^(/). (16)

*=1

Предполагаем, что форма спектральной плотности сигналов на выходе каждого из N фильтров близка к прямоугольной, а также что амплитудночастотная характеристика имитатора не имеет в частотном диапазоне возбуждения [утт, утах] резонансов и антирезонансов, либо они скорректированы. При этих условиях спектральная плотность Ях(у) моделирующих колебаний (16) представляет собой кусочно-постоянную функцию (рис. 3) а ее варьируемыми переменными являются уровни 2* (А: = 1,Л^) спектральной плотности составля-

ляющих процессов х*(/), сформированных &-ым фильтром, и координаты у* (к = 1,— 1) границ частотных диапазонов к-го и (£+1)-го фильтров. Таким образом, погрешность моделирования (15) принимает вид:

е(у1.-"»Ч'лгч.б1--.6») = (Ч'* - УоГ1 х

> V

х£ Ж

(17)

Рис. 3. Спектральная плотность моделирующих колебаний

где \|/0 и ц>н не варьируются, а зафиксированы у0 =Ушш, Ч'а’ =Ущах- Необходимые условия экстремума функции (17): де/8ц1к= О (А = 1,// — 1); де/д()к=0 (к -1, Лг), - приводят к системе уравнений, связывающей варьируемые переменные Цк и 0$

(18)

(бы - 0* )[&+, + б* - 2$, (V*)] = О, (Л = 1, ЛГ -1),

(19)

Условие (18) удовлетворяется при равенстве нулю первого или второго сомножителя. Причем равенство ординат Q|C и соседних горизонтальных участков функции Бх(у) фактически объединяет их в одно целое и тем самым

максимизирует функцию (17), так как задача сводится к случаю с числом фильтров, на единицу меньшим заданного. Следовательно, для отыскания минимума е необходимо приравнять нулю второй сомножитель. Выполнение этого условия означает, что каждый вертикальный участок спектральной плотности ^(у) (кроме крайних левого при V = и правого при V = утах) должен пересекать график моделируемой спектральной плотности 52(у) и делиться точкой пересечения пополам. Подставляя в (18) £?*, выраженные через у* из уравнений

(19), приходим к системе нелинейных уравнений для нахождения у* (к = 1,ДО-1), минимизирующих критерий моделирования (17):

После того, как у* найдены, с помощью (19) определяются £?*.

В простейшем случае, когда N=2, задача сводится к решению одного трансцендентного уравнения, для которого с гарантией можно найти все корни, что обеспечивает отыскание глобального минимума погрешности (17). Совершенно иначе обстоит дело в случаях трех и более ступеней (N>3) в спектральной плотности ^(у), поскольку приходится решать не одно, а систему нелинейных трансцендентных уравнений (20). Известные численные методы не обеспечивают отыскания глобального минимума (17). Однако, учитывая особенность системы (20), заключающуюся в том, что в к-с уравнение входят только к-я, (Ы)-я и (к+1 )-я неизвестные, удалось разработать алгоритм решения системы (20), сводящий ее к одному нелинейному уравнению.

Предлагаемый алгоритм построен следующим образом. В качестве аргумента получаемого нелинейного уравнения выбирается у,. Остальные неизвестные ц/2, уз>---,Ч,лм последовательно выражаются через \]/|. Для этого каждое из неизвестных ум выражается через два предыдущие у* и у*., решением

При к = 2,3,...,(N-1) величина Л* определяется на предыдущем шаге при решении предыдущего уравнения вида (21). Так как (Л/-1)-ое уравнение вида (21) является последним уравнением в системе (20), найденное из него \\ifj должно равняться известной частоте утах в спектральной плотности моделируемого случайного процесса. Это условие и дает нелинейное уравнение отно-

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

М*.<-

- V -

15, (у>Лг + —-— У (уУV = 2^0,),

(20)

уравнения

(21)

сительно ц/р

Нч>і)= ^{ул--,[-¥2(ч/і)-]}-Утш =0,

(22)

линейной интерполяцией, можно найти аналитическое решение уравнения (21)

(Ш l)h

относительно у*+1. Действительно, левая часть уравнения (21) геометрически интерпретируется площадью заштрихованной на рис. 4 фигуры с основанием длиной (у*+1 - у*), боковыми сторонами ^(у*), 5г(^и) и ограниченной кусочно-линейной функцией 5Ду). Суммируя площади трапеций, подставляя результат в уравнение (21) и, разрешая его относительно получаем

Рис. 4. Кусочно-линейная функция Sz(v)

Vi+1 = Rk +S*+ /Г1 (vmjn +иЛХ5І+, -Snz)±{(«* +SD2 +2r'(vmin +nh) ■(Л* + OC.-S‘)~-h-\S!n+i -S„z){[vmin +(m + l)A-v*]-

*0.5

(23)

. + + Л 1(5; + 5;+1) + 2[х\,кЯк -(утш + «Л)5„х]}}°

у=т + 1 .

где Л - шаг интерполяции; тип- целые части чисел А '('К* -утт) и й"'(ч,*+1 -Ут,п)> соответственно,5/приу = 0,1,2,.,.,/Г1 (угмх - у|пт ) -заданные

значения 52(у) в узлах интерполяции Утш+у/г.

Формула (23) правомерна, если п > т + 2. В случае п = т + 1 в ней необходимо исключить слагаемое Л Если же п = т, то \|/,м определя-

,;=т+1

ется так = \»т1„ +т/» + А(5^+,-5*)"'[4.й*-Я*~5ДуА)]. (24)

Применение формул (23), (24) существенно сокращает затраты машинного времени, поскольку отпадает необходимость применения численного метода для уравнения (21), решаемого (ЛИ) раз при каждом вычислении функции У(ц/|), которое, в свою очередь, многократно повторяется при численном отыскании корней уравнения (22).

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

Сигналы идентифицированы временным рядом:

хк=х(кА), k = l,2...N, .

где А - шаг дискретизации аргумента t, NA<T.

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

x(t) = ц + ^Rj cos {2л у t + <Pj)+Z(t), -оо</ <оо,

■Н

(25)

(26)

где /и - постоянная составляющая процесса х(1); т - число скрытых периодичностей; V, и ср] - соответственно амплитуда, частота (в периодах за единицу времени 0 и фазау-ой периодичности, у =1,2, ..., т; г(() - стационарный случайный процесс со средним Е - 0 и дисперсией аг = ау2 = Е2г({).

Модули пакета позволяют по реализации (25) оценить ц , 7?,, у,, (р] и сг2 . Для определения частот V,, то в пакете предусмотрены два случая: случай известных частот у, и случай известных интервалов [ау, /?,], таких что а,Щ<Д, у =1,2, ..., т- Реализацию (26) перепишем в виде:

усоз (о^ + В^та)/к)+ г,.,

где 2к-2(кК), юу = 2я-ууДД= 1,2,И, А j= RJ соз ф1нВ]= ср). Амплитуды /? и фазы вычисляются по значениям А и В согласно формулам Л = 7Л2+Я2 ; ф = аОД (-В / А), А > 0; агс!ё (-5 / Л)-я; Л < 0, £>0; апЛ§(—2?//4)+я; Л < 0, Б>0; -я-/2, А = 0 В<0', л 12, неопределенно, А =0, В=0.

В случае известных частот считается, что частоты vJ (а следовательно и а]) известны. В таком случае параметры ц, А/, В], ..., А„„ Вт можно оценивать, минимизируя сумму квадратов:

2

,(27)

*-1 У"1 ' ' ' ’ '

В качестве оценки сг принимается статистика:

б (и, А\,ВХ,.... Ат,Вт )= £ хк - ц - £ (а/ соб аук + В] бш со ]к)

{1 1 —е(м,Д....я,А)| •

Однако вычислительная реализация формулы (27) при больших т является сравнительно сложной, поэтому в пакет "СПЕКТР" заложен другой алгоритм, так называемый метод циклического спуска. При его применении параметры Ар В] (/= 1,2,..., т) и /и оцениваются последовательно, при этом применяется несколько циклов уточнения оценок.

В случае известных интервалов частот исходными являются начальные значения амплитуд и фаз Яр и ср р , число Ыс циклов уточнения оценок (/= 1,2, аг а также интервалы [а^Р]\, удовлетворяющие условиям: 0<а,, рт< 1/Д , при у = 1, 2, ..., т-\ и Ъ<рга,<\1{ИА) при

у-1, 2,т. Для выбора интервалов {ап /?,]: вычисляются значения периодограммы для некоторого окна данных; находятся аргументы X/, х2, х„, ра-

зумного числа максимумов периодограммы; подбираются а, и так, чтобы 2жа]<х^2лР1 и РГа]<\1(ИА). Одновременно определяется и число т.

Во время работы алгоритма частоты могут выйти из исходных интервалов [а>, /3>]. Однако всегда

I 1 1 1

2 12 3 Ис

Поэтому при выборе параметра Ис и интервалов [а у, ру] желательно обеспечить нужные промежутки между соседними интервалами [ау, Ру], между нулем и а/, а также между рт и 1/Д.

В качестве одной из возможных математических моделей ряда (25) в пакете "СПЕКТР" (рис.5) выбрана модель стационарного случайного процесса, удовлетворяющего разностному уравнению АРСУ:

оГ, -//)+ ]>>,(*,., -/|) = г, + , (28)

где х = х ((А), Г = ...,-1,0,+1, ... ; ц - среднее процесса; 2, - случайный шум со средним Е 2,- 0 и дисперсией аг2= Е 2,2; Ц, Ь/ - неизвестные коэффициенты; р и (/-порядка соответственно авторегрессии (АР) и скользящего усреднения (СУ), р,<7>0.

Рис. 5. Возможные пути соединения модулей при спектральном анализе в рамках модели АРСУ

В пакете "СПЕКТР" для данной модели решаются две задачи: оценить параметры:

<гг={Е[х(*)-м¥)П>Р> Ч, А = (а,,а2,...,ар), 5 = (б„62,...,бДсг2; по этим параметрам вычислить оценки нормированной спектральной плотности и корреляционной функции.

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

В работе использован метод проекционных оценок, идея которого состоит в аппроксимации неизвестной плотности распределения некоторым отрезком ее ряда Фурье по подходящей системе функций. Тогда каждый коэффициент Фурье представляет собой линейный функционал от плотности и его можно оценить по наблюдениям X/,..., х^ ■ Когда плотность распределения/(х) аппроксимируется отрезком тригонометрического ряда Фурье и сосредоточена

в интервале (я; я), проекционной оценкой плотности распределения называется статистика вида:

1

Л*)=т- Х^с.е-

2.Л *=-оо

(29)

~ /V

где Ск - М - оценка коэффициента Фурье плотности распределения.

1~ I

Ск=Е/е"[*', к = 0,±1,...,а{сок} - некоторая весовая последовательность.

В пакете "СПЕКТР" используется следующий алгоритм определения весовой последовательности {&>*} в проекционной оценке (29). Обозначим А,={1}, А2={2,3}, А3={4,5,6},... ;Д5 = [5(5-1)/2+1, 5(5+ 1)/2]п7У

=5 Ск | . Положим Ж„=1и при | к\ е Л$,

К

[1-(м5)"‘, еслиN0> 1 + 51/4 и 5(5-1)<^,

[О, в противном случае.

(30)

Проекционная оценка с такой весовой последовательно-

стью проста в реализации и требует небольших ресурсов компьютера. На рис. 6 приведены

графики теоретической плостно-сти распределения, гистограммы и проекционной оценки плотно-_ , „ , „ сти распределения (29) с весовой

Рис. 6. График теоретической плотности распределения, г г ' '

гистограммы и проекционной оценки (29) плотности распре- Последовательностью (30) при

деления с весовой последовательностью (30) при N = 1000 N ~ 1000.

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

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

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

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

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

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

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

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

Р =/>(?,£). _ (31)

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

Р(Л>4) = Р(П>°) = Р^Л) Для всех £ е £ . (32)

Иначе говоря, р = йтг£ . (33)

В работе также рассмотрен МТЛ как тракт с постоянными коэффициентами. Тогда достигнуть абсолютной инвариантности можно в случае функциональной и линейной связи между собой случайных коэффициентов разложения паразитного колебания. Условия (32) и (33) требуют, чтобы число независимых переменных в разложении реализаций паразитных колебании было меньше раз-

мерности пространства этих колебаний. Достаточным условием относительной инвариантности к £ будет

= 0) + е(ч,е), (34)

где е - разность между левой частью (34) и первым слагаемым правой части. Известно, что помехи, вызванные рядом паразитных колебаний носителя в МТЛ, обладают достаточно четко очерченной структурой с ограниченным числом неопределенных параметров. Математической моделью такой помехи является квазидетерминированный случайный процесс с реализациями в виде функций ряда случайных параметров. Тогда, путем решения задачи минимизации коэффициентов разложения помехи и с учетом условий (32) и (34), может быть синтезирован МТЛ, инвариантный к оговоренной выше аддитивной ква-зидетерминированной помехе. '

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

В главе получены оценки случайных возмущений в МТЛ и методики снижения их уровня. Разработана математическая модель оценки отклонений СЧ МТЛ. Получены формулы для среднего оценок отклонений СЧ, являющегося критерием риска захода СЧ в РОЗ; определены точные и асимптотические выражения для дисперсии оценки отклонений СЧ. Предложены сравнительно точные и удобные для построения доверительных интервалов экспоненциальные оценки сверху для вероятности, что оценки отклонений СЧ превысят заданный уровень. На основе изучения статистических характеристик отклонений СЧ МТЛ предложена методика коррекции границ РОЗ с целью ликвидации возможности захода СЧ в РОЗ в течение всего времени транспортирования носителя.

При разработке математической модели оценки отклонений СЧ МТЛ учтено, что при малых изменениях параметров МТЛ отклонение СЧ р, от номинального значения р1 определяется первыми членами ряда Тейлора:

АР,(0 = Р, ~Р, = Е(Ф</5м(«/ ~ «/) = ЕСФ/Щ),,,*щ Аи/ •

/»1 ы

где и, - номинальное значение/-го параметра; Ар,(/)=Дд (/, и(з),0 <.у < ?)=[А/7/з]/=—

- случайное слагаемое, зависящее от времени.

Процесс отклонения г-го СЧ от номинального является реализацией триерного случайного процесса:

АР/(0 = /п(0 + ^(0, - °о </..< +аз, где т(/) = [ту(/)]у=^ - некоторая неслучайная функция, максимальная частота которой много меньше минимальной СЧ МТЛ; Х(() =[АГ^(/)] ^ - г-мерный

среднеквадратически непрерывный стационарный в узком смысле действительный случайный процесс со средним ЕХ(() и спектральной плотностью с низкочастотным диапазоном.

В качестве величины отклонения / -го СЧ примем случайную величину:

0* =^1{[(Ф, /5ц)„1=,7, Ли,(0]2 + [(Ф,/£Ы,г=,72Л“2(')]г +...+КФ,/диг),^Пг Аиг(1)]2}сЬ,

'о 1 '

представляющую собой среднеквадратическое отклонение /-го СЧ от номинального. Здесь Т - время транспортирования ленты.

Реально изменение параметров МТЛ в процессе транспортирования ленты приводит к отклонению не только одной СЧ, а одновременно вызывает смещение всего спектра СЧ. Поэтому естественно в качестве оценки отклонений частотного спектра принять следующую случайную величину:

Qr = \[£Р&) + £р2(?) + ..Л£р„(!)№■ (35)

* О

Дисперсия оценки отклонений СЧ 0()-р рассмотрена в работе в случае т(г) = 0. Предполагая, что ЕХ*(0<<х>, у = 1 ,г. Тогда с помощью непосредст-

1 тт г 2 венного подсчета получаем, что 0<^] = — // (*)}<*<* > (36)

Т 0 07Л=1

Если процесс гауссовский, то соу{Лг2(/),А'*(0)} = 2Яд(/), и формула (36) принимает более простой вид:

](Н'17,-1 >!*<') 12Л- (37)

I _г

Применение формулы (37) осложняется тем, что надо полностью знать ковариационную функцию /?(/)• Однако, подставляя в (37) более или менее точные оценки этих функций, можем получить и более или менее точные оценки дисперсии

Если процесс Х(1) гауссовский и имеет конечное время перемешивания, то из (37) получаем:

1 ],к=1

где ст2, - дисперсия компонент у-ой и к-он вектора Х(1); - оценка ин-

тервала корреляции, ],к = \,г.

В работе получены экспоненциальные оценки вероятности распределения оценки отклонений СЧ динамической модели МТЛ, которые особенно точны при больших х и Т (здесь х - абсцисса графика вероятности распределения). Это их свойство связано с тем, что при Т -> оо распределение оценки отклонений становится нормальным.

В случае, когда изменение параметров Х(1) подчиняется гауссовскому закону с необязательно конченым т(интервал корреляции), имеют место доста-

В работе указано всего 8 экспоненциальных оценок вероятности Р{<2^ >:о2 +х} через разные характеристики оценки отклонений СЧ X(I): 4 оценки с функцией С и 4 более простого вида. Заметим, что оценки с функцией б существенно точнее, поэтому в случае применения вычислительной техники надо выбирать их. С другой же стороны, оценки без функции б гораздо нагляднее. Фигурирующими в оценках вероятностными характеристиками оценки отклонений СЧ Х(/) являются: в оценках (38) а и йг,в (39) а2 и т. В случае применения установленных неравенств эти характеристики надо знать заранее. При этом они могут быть и завышенными, в таком случае приводимые оценки останутся справедливыми, но менее точными. Характеристики а, Вг, т, о2, и у являются весьма простыми и наглядными, и более или менее завышенные их оценки всегда можно определить экспериментально или просто указать заранее, исходя из конкретных соображений о МТЛ.

Разработана методика учета случайного закона изменения СЧ при синтезе МТЛ по частотным спектрам для управления коррекцией границ РОЗ.

Учитывая вышеизложенное, для устранения заходов СЧ в РОЗ границы последних следует раздвинуть и сделать равными:

них, то с вероятностью 0,95 можно считать, что СЧ в процессе транспортирования ленты не попадут в РОЗ. В качестве дополнительного критерия синтеза МТЛ по частотным спектрам можно выбрать функционал, минимизирующий случайную величину <2Т (35). Для этого необходимо экспериментальным путем изменять жесткость и вязкость носителя и выбирать ту марку носителя, которая обладает меньшим разбросом параметров.

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

р*(Г)<ю*(0<сх*+1(0. к = \,т-\.

Однако в процессе транспортирования ленточного носителя СЧ пред-

ты следующие обозначения

Тогда для всех х > 0 и Т > 0 справедливы оценки:

(39)

(38)

а* = ау - 3 шах д/; Р* = + 3 шах д/.

Если взять интервалы (а*,р*), ] = \,т и решить задачу отстройки для

ставляют собой в общем случае нестационарный случайный процесс. За время Т движения ленты СЧ в каждый момент времени претерпевают отклонения относительно средних значений со*(/). Эти отклонения могут привести к нарушения неравенств (40).

Оценим вероятность захода СЧ в РОЗ:

/у = 1 — Рфк (0 < Ю* (0 < а*+, ). А = 1, и, о < / < Г}.

Введем некоторые обозначения:

$*(0=<М0-®*(0,- и <41>

где (5* (0 - детерминированная функция изменения А-ой СЧ при движении ленты; со* (/) - процесс изменения А-ой СЧ при движении ленты.

Обозначим дисперсию процесса изменения к -ой СЧ:

/ВД=°*(0- _ (42)

Для удобства в дальнейшем будем оперировать безразмерной функцией.

4*(0/а*(0 = Л*(0 (43)

в предположении, ЧТО 1>Г|* (0 = м\ (/).

С учетом (42) и (43) неравенства (40) трансформируются в условия: *к(0<Цк(?)<ик(1), где и*(0=[а*+1(0-®*(0]/ст*(0; У*(/) = [со|1(/)-Р*(/)]/а4(0-а вероятность (41) определяется выражением:

рт = 1 - />*(/) < лДО < 0. А = 1,и, 0 < / < Г} •

Принимая во внимание теорему для абсолютно непрерывной функции

г

со*(0, для которой |м*(/)с* <00 и для всех 0 </<Г М4(0>0>

. О

имеем оценку сверху:

Рт * I {(0)]+т (0)] + {[<7* (0+Рк (№ | , (44)

*»1 [ О }

где qk{t) = ф[и*(0]{ц*(0ф["*(0/Ц*(01 - и'к(0/ц*(0]}«

рк (0=Фь (/)]{ц* (ОфК (0/ц* (0] - ^(О^К (0/ц* (/)]};

функции ф и V)/ определяются по следующим известным формулам:

ОО

ф(лг) = ехр(-х2/2)/л/2тг; Ч/(х) = 1-Ф(д:),гдеФ(х)= /ф(г)йЬ.

X

Рассмотрим частные случаи.

1. Пусть ю*(0 и а*(0 - медленно меняются. Тогда и[(/) г 0 и у'к(1) = О и д*(0 = фМ0К(')/2тг = Ц*('К"*(')/2 12*;

рМ=Фь№к(0/2п = 11к«)е-'^),2/2п. (45)

2. Стационарный случай

со*(0эш*; ст*(0 = ст4; Ц*(0 = И*; Щ^) = ик\ у4(/) = у*; оценка вероятности захода СЧ в РОЗ принимает следующий вид:

Рт ^Ет«*] + ^] + 7и*(в-*г/2 +е-"‘/2)(2я)-1} (46)

ы

Оценки вероятности (44)-(46) позволяют судить о степени динамической точности МТЛ и могут являться критериями спектрального синтеза стохастических динамических систем.

ЗАКЛЮЧЕНИЕ

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

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

3. Решена задача отыскания глобального минимума погрешности моделирования в случае трех и более ступеней спектральной плотности путем разработки алгоритма решения системы нелинейных трансцендентных уравнений, учитывая ее особенность, которая заключается в том, что в к-е уравнение входят только к-я, (£-1)-я и (&+1)-я неизвестные, что сводит ее решение к одному нелинейному уравнению.

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

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

6. При анализе сигналов последние рассматриваются как суперпозиция совокупности детерминированных сигналов и случайного шума. Разработанный пакет программ "СПЕКТР" включает алгоритм выделения гармонических составляющих, который предоставляет пользователю самостоятельно выбирать их число.

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

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

8. В качестве одной из возможных математических моделей представления сигнала в созданном пакете программ "СПЕКТР" выбрана модель стационарного случайного процесса, удовлетворяющего разностному уравнению АРСУ. Для данной модели решены две задачи: оценка параметров среднего, дисперсии сигнала и порядков АР и СУ и по этим параметрам вычисление оценок спектральной плотности и корреляционной функции.

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

10. В пакет "СПЕКТР" включены пять модулей, позволяющих оценивать одномерное распределение временного ряда, каждый из которых, соответственно, определяет выборочную функцию распределения сигнала, вычисляет гистограмму, реализует алгоритм адаптивной проекционной оценки плотности распределения с весовой последовательностью, строит гауссовскую функцию распределения или (по желанию пользователя) гауссовскую плотность распределения, формирует графики указанных функций.

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

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

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

14. Даны оценки вероятности захода СЧ в РОЗ в случае, когда инерцион-

ные и упруго-вязкие свойства ленточных носителей изменяются в зависимости от их длины по случайному закону. Оценки получены для стационарной и нестационарной динамической модели MTJI. На основе изучения статистических характеристик отклонений СЧ MTJI предложена методика коррекции границ РОЗ с целью ликвидации возможности захода СЧ в РОЗ в течении всего времени транспортирования носителя.

НАУЧНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Кытин Е.А. Оценка параметров деполяризованных сигналов // Вестник Московской Академии рынка труда и информационных технологий, 2005. - № 9 (21). - С. 106-108.

2. Кытин Е.А., Лялин В.Е. Синтез полосовых фильтров и фазокорректирующих цепей радиотелефонной системы на основе распознавания взаимного расположения спектра собственных частот И частот передающих станций // Научно-практический журнал «Интеллектуальные системы в производстве» № 2(16), 2010, С.142-154.

3. Лялин В.Е., Кытин Е.А. Критерии интеллектуальных систем синтеза приемо-передающих радиотелефонных устройств // Научнопрактический журнал «Интеллектуальные системы в производстве» № 2(16), 2010, С. 154-166.

4. Кытт Е.А. Преобразование электромагнитных волн элементами линейных трактов телекоммуникационных систем // Надежность и качество. Труды международного симпозиума: В 2-х томах / Под ред. Н.К. Юркова. - Пенза: Изд-во Пенз. гос. ун-та, 2007. - Т. 1.-С. 513-515.

5. Кытин Е.А. Оценка матрицы когерентности и вектора средних значений поляризованного сигнала // Информационные технологии в науке, образовании, телекоммуникации и бизнесе. Материалы 35-й междунар. конф. - Украина, Крым, Ялта-Гурзуф: Прилож. к журн. «Открытое образование», 2008. -С. 424-426.

6. Кытин Е.А. Оценка инварианта матрицы когерентности поляризованного сигнала // Информационные технологии в науке, образовании, телекоммуникации и бизнесе. Материалы 35-й междунар. конф. - Украина, Крым, Ялта-Гурзуф: Прилож. к журн. «Открытое образование», 2008. - С. 427-428.

7. Кытин Е.А. Адаптация поляризационных характеристик линейного тракта телекоммуникационной системы к состоянию канала связи // Надежность и качество. Труды международного симпозиума: В 2-х томах / Под ред. Н.К. Юркова. - Пенза: Изд-во Пенз. гос. ун-та, 2008.-Т. 2. - С. 184-185.

8. Кытин Е.А. Адаптация поляризационных характеристик приемного и передающего трактов к статистическим характеристикам канала связи // Надежность и качество. Труды международного симпозиума: В 2-х томах / Под ред. Н.К. Юркова. - Пенза: Изд-во Пенз. гос. ун-та, 2008. - Т. 2. - С. 185-186.

9. Кытин Е.А., Марков М.М. Анализ информационных признаков сигналов с амплитудной модуляцией для анализа качества канала связи // Надежность и качество. Труды международного симпозиума: В 2-х томах / Под ред. Н.К. Юркова. - Пенза: Изд-во Пенз. гос. ун-та, 2009. - Т. 1. - С. 253-256.

10. Кытин Е.А. Исследование информационных признаков сигнала для контроля качества работы канала связи // Информационные технологии в науке, образовании, телекоммуникации и бизнесе. Материалы 36-й междунар. конф. -Украина, Крым, Ялта-Гурзуф: Прилож. к жури. .«Открытое образование», 2009. -С. 341-342.

11. Кытин Е.А. Информационные признаки сигналов с частотной модуляцией для контроля качества работы канала связи // Информационные технологии в науке, образовании, телекоммуникации и бизнесе. Материалы 36-й междунар. конф. - Украина, Крым, Ялта-Гурзуф: Прилож. к журн. «Открытое образование», 2009. - С. 343-344.

12. Кытин Е.А. Оценка влияния статических и динамических ошибок анализатора качества канала связи // Научно-практический журнал «Отраслевые аспекты технических наук» - № 4. - Москва, 2011. - С. 2-5.

13. Кытин Е.А., Лялин В.Е. Методика выбора числа полосовых фильтров для анализатора качества канала связи // Научно-практический журнал «Приволжский научный вестник». - № 3. - Ижевск, 2011. - С. 11-17.

14. Кытин Е.А., Лялин В.Е. Спектральный анализ сигналов в рамках модели авторегрессии и скользящего усреднения // Научно-практический журнал «Приволжский научный вестник». - № 4. - Ижевск, 2011. - С. 19-24.

Подписано в печать 30. 01.12. Уел. печ. л. 1,0. Заказ № 17. Тираж 100 экз.

Издательство Ижевского государственного технического университета Отпечатано в типографии Издательства ИжГТУ. 426069, Ижевск, Студенческая, 7

61 12-5/3329

ФГБОУ ВПО «ИЖЕВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

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

КЫТИН Евгений Алексеевич

УДК 621.396.2

ОПРЕДЕЛЕНИЕ СПЕКТРАЛЬНЫХ ХАРАКТЕРИСТИК ПОЛОСОВЫХ ФИЛЬТРОВ ДЛЯ АНАЛИЗА КАЧЕСТВА НЕОДНОРОДНОГО КАНАЛА СВЯЗИ

Специальности:

05.13.01 - Системный анализ, управление и обработка информации

(в науке и технике); 05.12.13 — Системы, сети и устройства телекоммуникаций

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель:

заслуженный изобретатель РФ, доктор технических наук, профессор Лялин В.Е.

Ижевск-2012

ОГЛАВЛЕНИЕ

стр.

Введение............................................................

6

Глава 1. Теоретические основы проектирования неоднородных

каналов связи............................................................................... * ^

1.1. Организация телекоммуникационной системы резервного копирования данных ....................................................................

18

1.2. Программные продукты управления резервным копированием....

1.3. Ленточные библиотеки........................................................................................................................................25

1.4. Ленточные накопители........................................................................................................................................27

1.4.1. Серпантинная запись..................................................................................................................28

1.4.2. Линейная запись....................................................................................................................................31

1.5. Критерии отстройки спектра частот собственных колебаний относительно возмущающих частот............................................ 35

1.5.1. Консервативные динамические системы........................ 35

1.5.2. Диссипативные динамические системы.......................... 49

1.6. Частотно-избирательные устройства.......................................... ^

1.6.1. Виды АЧХфильтров, передаточные функции

и чувствительность........................................................... 58

1.6.2. Многопетлевые активные фильтры.......................................... 64

1.6.3. Реализация нулей передаточной функции фильтра.............. 72

1.6.4. Достоинства многопетлевых фильтров.............................. 72

1.7. Технические средства и методическое обеспечение анализа

канала связи..............................................

......................................... 73

1.8. Выводы, постановка цели и задач исследования.......................

Глава 2. Определение частотных характеристик полосовых фильтров, предназначенных для имитации и анализа сигналов в канале связи............................................................. "

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

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

2.3. Полученные результаты и выводы..............................................

Глава 3. Программный комплекс "Спектр" как инструмент для , 04 анализа канала связи..................................................................

3.1. Описание пакета "Спектр"........................................................... 105

3.1.1. Функциональные возможности пакета....................................................105

3.1.2. Структура пакета..........................................................................................................................106

3.1.3. Условия применения пакета "Спектр"........................................................107

3.2. Предварительная обработка сигналов..................................................................................107

3.3. Непараметрический спектральный анализ..................................................................108

3.3.1. Базовые статистики......................................................... ПО

3.3.2. Оценки типа сглаженной периодограммы...................... 111

3.3.3. Оценки, построенные с помощью

ковариационных окон......................................................... ^ ^

3.3.4. Оценки с адаптивными ковариационными окнами........ 116

3.3.5. Оценки с использованием окон данных............................ 117

3.3.6. Оценки типа осреднения периодограмм.......................... 119

3.3.7. Состав группы модулей для непараметрического спектрального анализа...................................................... ^ ^

3.4. Выделение скрытых периодичностей......................................... 121

3.4.1. Случай известных частот................................................ 122

3.4.2. Случай известных интервалов частот........................... 124

3.4.3. Состав группы модулей для выделения скрытых периодичностей..................................................................

3.5. Спектральный анализ в рамках модели авторегрессии -скользящего усреднения............................................................... 126

3.5.1. Метода Юла-Уокера......................................................... 126

3.5.2. Метод наименьших квадратов......................................... 128

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

3.6. Анализ одномерного распределения........................................... 131

3.6.1. Проекционные оценки........................................................ 132

3.6.2. Состав группы модулей для анализа одномерного распределения.....................................................................

3.7. Выбор методов спектрального анализа и построение

программ пользователя................................................................. 135

3.7.1. Пример программы пользователя.................................... 136

3.8. Полученные результаты и выводы.............................................. 137

Глава 4. Оценки случайных возмущений и методы снижения их

уровня в динамических системах стримеров.......................... 139

4.1. Анализ инвариантности воздействия механических возмущений при транспортировании ленточного носителя в стримерных устройствах на их параметрическую надежность...................... 139

4.2. Методика оценки отклонений частот собственных колебаний в неоднородном канале связи.......................................................

4.3. Дисперсия оценок отклонений частот

собственных колебаний................................................................ 150

4.4. Оценки плотности распределения оценок отклонений частот собственных колебаний................................................................ 153

4.5. Асимптотическое поведение оценок плотности распределения оценок отклонений частот собственных колебаний........... 155

4.6. Учет случайного закона изменения частот собственных колебаний в неоднородном канале связи............................................ 159

4.7. Оценка вероятности попадания частот собственных колебаний в резонансно-опасные зоны.................................................. 160

4.8. Полученные результаты и выводы................................................ 162

Заключение....................................................................................................... 164

Литература........................................................................................................ 168

ВВЕДЕНИЕ

Актуальность темы. Одной из основных задач при повышении качества работы информационной радиотехнической системы является определение по-меховой обстановки в реальных рабочих условиях. Для этого необходимо определить минимальное число параллельных фильтров, которые позволили бы проанализировать зашумленность как в поддиапазонах, так и по всей полосе канала связи (КС). При этом количество фильтров такого анализатора должно удовлетворять заданной точности оценивания зашумленности канала. Нужна методика определения необходимого числа полосовых фильтров для анализатора качества канала связи при заданной точности измерений в канале связи.

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

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

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

Для решения вышеуказанных задач необходим спектральный анализ, который является мощным средством исследования сигналов. Методы спектрального анализа в настоящее время широко применяются в самых разных областях науки и практики. Для обеспечения анализа и статистической обработки проведено огромное количество работ видными учеными разных стран. Большой вклад в решение указанной проблемы внесли ученые Бокс Дж., Дженкинс Г., Бриллинджер Д., Вате В., Макс Ж., Робинсон Э.А., Отнес Р., Эноксон Л., Ста-тулявичус В.В., Бенткус Р.Ю., КэдзоуДж.А., Поспелов Г.С., Поспелов Д.А., Яглом А.Н., Журбленко И.Г., Ибрагимов И.А., Линник Ю.В.

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

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

участки переприема.

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

Известно, что параметры ленточного носителя, такие как жесткость и вязкость, в зависимости от длины ленты в бобине изменяются по случайному закону. Из этого следует, что при транспортировании ленты в каждый момент времени формируется вектор мгновенных собственных частот (СЧ), который отличается от векторов СЧ, сформированных в предыдущие моменты времени. Отклонения СЧ от номинальных, определяемых детерминированными параметрами МТЛ, синтезированными по существующим критериям и решающим функциям, могут привести к тому, что одно или несколько СЧ, располагаясь вблизи границ резонансно-опасной зоны (РОЗ), может зайти в запретные области, повышая уровень колебаний звеньев. С целью упреждения этого факта необходимо откорректировать значения границ РОЗ в соответствии с оценками дисперсии, закона распределения и асимптотики распределения оценки отклонений СЧ.

В связи с вышеизложенным тема настоящего диссертационного исследования весьма актуальна.

Объектом исследования являются: анализ и имитационное моделирование состояния неоднородного КС.

Предметом исследования является разработка математических и инст-

рументальных средств для определения спектральных характеристик полосовых фильтров неоднородного КС.

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

Для достижения поставленной цели требуется решить следующие задачи:

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

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

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

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

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

- определить условия инвариантности МТЛ стримеров по отношению к возмущающим факторам для поддержания заданных точностных характеристик информационной части записываемого сигнала при условии рассмотрения стримера как аналога или составной части канала передачи информации;

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

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

Методы исследования. В работе применялись теоретические и экспериментальные методы исследования.

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

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

Экспериментальные исследования проводились путем имитационного мо-

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

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

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

Достоверность вычислительного эксперимен�