ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧЕРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
«БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ
ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ»
( Н И У
« Б е л Г У » )
ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК
КАФЕДРА МАТЕМАТИЧЕСКОГО И ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ
ИНФОРМАЦИОННЫХ СИСТЕМ
РАЗРАБОТКА ТЕХНОЛОГИИ ФОРМИРОВАНИЯ СИГНАЛОВ,
УСТОЙЧИВЫХ К ВОЗДЕЙСТВИЮ КОРОТКОДВИЖУЩИХ ИЛИ
СОСРЕДОТОЧЕННЫХ ПО СПЕКТРУ ПОМЕХ
Выпускная квалификационная работа
обучающегося по направлению подготовки 02.04.01 Математика и
компьютерные науки, группы 07001631
Киселёва Юрия Игоревича
Научный руководитель
к.т.н., доцент
Муромцев В.В.
Рецензент
д.т.н., доцент
Жиляков Е.Г.
БЕЛГОРОД 2018
2
ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ........................................................... Error! Bookmark not defined.
ГЛАВА 1 ПОНЯТИЕ ШУМА И ПОМЕХИ СИГНАЛА.
СОВРЕМЕННЫЕ МЕТОДЫ СИНХРОНИЗАЦИИ СИГНАЛОВ ..... Error!
Bookmark not defined.
1.1 Виды шумов и помех сигнала ..................... Error! Bookmark not defined.
1.2 Кратковременные помехи ............................................................................. 9
1.3 Типы синхронизации ................................... Error! Bookmark not defined.
1.4 Роль синхронизации в современных сетях Error! Bookmark not defined.
1.5 Тактовая синхронизация ............................. Error! Bookmark not defined.
1.6 Символьная синхронизация ........................................................................ 20
ГЛАВА 2 ФОРМИРОВАНИЕ СИГНАЛОВ, УСТОЙЧИВЫХ К
ВОЗДЕЙСТВИЯМ КОРОТКОЖИВУЩИХ ИЛИ
СОСРЕДОТОЧЕННЫХ ПО СПЕКТРУ ПОМЕХ ...... Error! Bookmark not
defined.
2.1 О классе сигналов, допускающих восстановление одних фрагментов по
остальным ........................................................... Error! Bookmark not defined.
ГЛАВА 3 РАЗРАБОТКА ТЕХНОЛОГИИ ФОРМИРОВАНИЯ
СИГНАЛОВ НА БАЗИСЕ СОБСТВЕННЫХ ВЕКТОРОВ
СУБПОЛОСНОЙ МАТРИЦЫ УСТОЙЧИВЫХ К ВОЗДЕЙСТВИЮ
КОРОТКОЖИВУЩИХ ИЛИ СОСРЕДОТОЧЕННЫХ ПО СПЕКТРУ
ПОМЕХ ................................................................................................................. 37
3.1 Исследование разработанной технологии ................................................. 37
ЗАКЛЮЧЕНИЕ ................................................................................................... 63
3
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ ....................................... 65
ВВЕДЕНИЕ
C увеличением потребления цифрового контента каждый день в мире
увеличивается
потребность
в
качественном
обеспечении
передачи
информационных сигналов. В связи с этим данная система должна обладать
определенными критериями, одним из которых является корректная работа
всей системы в целом.
Под сигналом понимается некоторая функция времени, в параметрах
которой
закодирована
информация,
предназначенная
определенным
пользователям (коммуникационные сигналы). Сигналы нашли широкое
применение в телекоммуникациях, телемеханике, дистанционном управлении
объектами и т.д.
Наблюдающаяся тенденция возрастания интенсивности
информационного обмена между различными пользователями привела к
необходимости решения проблемы искажений одних сигналов другими,
которые в таком контексте называются помехами.
Достаточно часто искажениям подвергаются либо отдельные временные
фрагменты сигналов либо фрагменты их спектров. В первом случае речь идет
о
воздействиях
короткоживущих
помех,
тогда
как
во
втором-
о
сосредоточенных по спектру Фурье помехах. Короткоживущие по сравнению
с длительностью сигнала помехи достаточно часто возникают при замираниях
в каналах их передачи, например при многолучевом распространении
радиоволн
или
переотражениях
в
условиях
городской
постройки.
Возникновение узкополосных помех также может быть обусловлено
различными причинами, среди которых основное место занимают источники
электромагнитного излучения. Важно то, что искажения из- за воздействий
других сигналов возникают непредсказуемо, хотя их можно фиксировать в
процессе передачи/приема собственных сигналов. Таким образом, борьба с
этими помехами приобретает все большую актуальность
4
Актуальностью выпускной квалификационной работы выступает спрос
в постоянной и качественной работе беспроводной сети. Для этого необходима
согласованная работа всей беспроводной системы связи, что обеспечивается
благодаря борьбе с помехами.
Целью работы является разработка технологии формирования сигналов,
устойчивых к воздействию короткодвижущих или сосредоточенных по
спектру помех. Для достижения поставленной цели необходимо выполнить
следующие задачи:
1.
Провести
анализ
современных
методов
формирования
информационного сигнала устойчивого к воздействию короткодвижущих или
сосредоточенных по спектру помех;
2.
Изучить
информацию
о
классе
сигналов,
допускающих
восстановление одних фрагментов по остальным;
3. Исследовать возможность использования собственных векторов
субполосных ядер в поставленной задаче;
4.
Разработать
программный
код,
реализующих
технологию
формирования сигналов, устойчивых к воздействию короткодвижущих или
сосредоточенных по спектру помех на основе базиса собственных векторов
субполосных ядер;
5. Провести анализ полученных результатов реализации и выполненной
работы, сделать выводы;
Предметом исследования является процесс передачи информационного
сигнала в современных системах связи.
Объектом исследования являются методы формирования сигналов
устойчивых к воздействию помех.
В главе 1 рассмотрено понятие шума и помехи сигнала, последствия их
воздействия на информационный сигнал, методы борьбы с ними, а также
современные методы синхронизации сигналов в системах связи.
5
В главе 2 исследованы классы сигналов, допускающих восстановление
одних фрагментов по остальным на основе базиса собственных векторов
субполосной матрицы.
Глава
3
рассматривает
программную
реализацию
технологии
формирования сигналов, устойчивых к воздействию короткодвижущих или
сосредоточенных по спектру помех на основе базиса собственных векторов
субполосных ядер, а также вывод по полученным результатам программной
реализации
6
ГЛАВА I
ПОНЯТИЕ ШУМА И ПОМЕХИ СИГНАЛА. СОВРЕМЕННЫЕ
МЕТОДЫ СИНХРОНИЗАЦИИ СИГНАЛОВ
1.1 Виды шумов и помех сигнала
Шум и помехи одинаково влияют на прием удаленных сигналов,
наслаждаться высококачественной музыкой, просматривать фотографии и
видео. Шум формируется внутри радиоэлектронного устройства (телевизор,
радио, сотовый телефон или камера). Шум отличается от шума тем, что они
входят в электронное устройство снаружи. На рисунке 1.1 показаны основные
схемы, на которые воздействуют помехи.
Рис. 1.1 Основные электрические цепи, на которые воздействует
электрические помехи
Конечно, основным приемником помех является вход усилителя.
Именно отсюда усиливается сигнал. И если помеха направлена на вход
усилителя, он будет усилен, и в будущем будет очень сложно отделить
полезный сигнал.
7
Не менее опасно с точки зрения. В этом случае помехи могут быть
сформированы в источниках вторичного источника питания. В этом смысле,
как приемник помех, могут служить схемы с общими проводами. Когда токи
протекают через них из разных блоков, образуется падение напряжения,
которое может влиять на высокочувствительные компоненты устройства.
Основные источники помех:
• Магнитное поле;
• Электрическое поле;
• Радиоволны;
• Совместное падение напряжения на одном проводнике;
• Микрофонный эффект;
Магнитное поле действует в основном на индуктивные элементы, такие
как индукторы, дроссели или электродинамические микрофоны. При высоких
интенсивностях магнитного поля он будет вызывать заклинивающие токи на
соединительных проводах, в том числе на силовых цепях транзисторов и
микросхем.
Электрическое поле похоже на характер его воздействия на элементы
радиоэлектронного оборудования, такие как магнитное поле, но в отличие от
магнитного поля электрическое поле индуцирует напряжение или потенциал
помех.
Радиоволны представляют собой электромагнитное поле по своей
природе, поэтому они могут направлять на элементы радиоэлектронной
схемы, как токи, так и напряжения. Борьба с этим типом помех является
наиболее сложной задачей, поскольку высокочастотные токи могут проникать
глубоко в экранирующие материалы.
Когда ток протекает через омическое или индуктивное сопротивление,
на нем происходит падение напряжения. В радиоэлектронных устройствах
одновременно присутствуют мощные блоки или блоки с повышенным
уровнем помех и высокочувствительные блоки. Обычно они питаются от
8
одного и того же источника тока. Однако сопротивление проводника невелика,
но все же имеет сопротивление, при котором происходит падение напряжения.
Поскольку один и тот же проводник подключен к блоку с высокой
чувствительностью, помеха может проникать на выход этого устройства. Это
означает, что повышенные требования обычно устанавливаются на схему
силовых цепей.
Учитывая
огромное
разнообразие
помех
в
радиоэлектронном
устройстве, существует огромное количество способов борьбы с помехами.
Это и технические, и организационные методы борьбы с помехами.
Вмешательство может быть подавлено в месте их возникновения и
предотвратить их проникновение в электрические цепи радиоэлектронного
устройства. Обычно при рассмотрении вопроса об управлении помехами
рассматриваются технические методы снижения уровня помех на выходе
устройства.
Технические методы снижения шума:
• Экранирование;
• Заземление;
• Балансировка;
• Фильтрация;
• Разнесение и ориентация;
• Регулировка величины полного сопротивления схемы;
• Выбор кабеля;
• Подавление (в частотной или временной области);
Шум
отличается
от
помех
тем,
что
появляется
внутри
радиоэлектронного устройства. В качестве источника шума используются
радиоэлементы, из которых он собирается. Основными источниками шума
являются следующие элементы:
• Резисторы (тепловой шум);
• Транзисторы;
9
• Диоды;
1.2 Кратковременные помехи
В
каналах
связи,
помимо
флуктуационных
помех,
хорошо
аппроксимированных гауссовским шумом, часто встречается аддитивная
интерференция
различной
средневолновых
и
природы.
коротковолновых,
В
радиоканалах,
доминирующую
особенно
роль
играют
сосредоточенные и импульсные шумы, которые являются предметом этой
главы.
Концентрированная интерференция называется такой аддитивной
интерференцией, в которой основная часть мощности сосредоточена в
отдельных полосах частот, меньших или сравнимых с, где - длительность
сигнального элемента. Они чаще всего встречаются в радиоканалах в
результате
воздействия
принадлежащих
на
посторонним
принимающее
каналам
устройство
связи.
В
сигналов,
этом
случае
концентрированный шум также называется взаимной помехой.
Импульсивная
интерференция
называется
такой
аддитивной
интерференцией, которая отличается от нуля только на отдельные интервалы
времени, значительно меньшие, чем те, которые разделены гораздо более
длинными интервалами, без помех. Импульсный шум является регулярной
или случайной последовательностью мешающих импульсов. Источники
импульсного шума в радиоканалах очень разнообразны. К ним относятся
атмосферные выбросы, промышленные установки, в частности система
зажигания в двигателях внутреннего сгорания, медицинские и бытовые
электроприборы и т. Д.
Концентрированный и импульсный шум наблюдается также в других
каналах связи, например, кабельных, гидроакустических и т. Д.
Идеализированным предельным случаем концентрированного шума
является сумма монохроматического заклинивания со случайными (но не
10
изменяющимися во времени) амплитудами, частотами и фазами. Суммарная
ширина полосы, занимаемая такой интерференцией, имеет меру нуля и,
следовательно, не уменьшает пропускную способность канала. Из этого
следует, что в принципе должны быть методы получения сигнала, в котором
идеализированные концентрированные помехи могут быть полностью
подавлены, т. Е. Не вызывать ошибок.
Идеализированным предельным случаем импульсного шума является
последовательность дельта-функций, возникающих в случайные моменты
времени со случайными интенсивностями. Энергетический спектр такого
возмущения неограничен, но общее время его существования имеет меру нуля.
Следовательно, он также не уменьшает пропускную способность канала, и
должны быть методы приема, которые могут полностью подавлять
идеализированный импульсный шум.
Реальные компоненты концентрированного шума не являются точно
монохроматическими, так же как фактические импульсы не являются дельтафункциями. Поэтому полное подавление таких помех невозможно, но их
можно частично подавить. Мы объясним, что мы имеем в виду под этим.
Предположим, что на выходе канала имеется сигнал и некоторый
аддитивный негауссовский шум, который мы будем характеризовать средней
спектральной плотностью мощности. Предположим, что решающая схема
выбрана оптимальной для гауссовской интерференции. Вероятность ошибок в
этом случае может быть, в зависимости от характера помехи, больше или
меньше, чем в случае гауссовской интерференции с той же спектральной
плотностью. Однако во многих случаях, как будет показано ниже, он мало чем
отличается от вероятности ошибки для гауссовских помех. Если можно
выбрать сигналы и построить такую схему принятия решений, вероятность
того, что вероятность ошибочного приема сигнального элемента с конечной
длительностью для данной интерференции будет сколь угодно мала, то мы
говорим, что помехи полностью подавлены. Если, однако, вероятность
ошибки остается конечной, но существенно меньше, чем для гауссовской
11
интерференции с той же спектральной плотностью, интерференция частично
подавляется.
К сожалению, в настоящее время практически нет общей теории
оптимального приема негауссовского шума. Его конструкция связана с
большими трудностями, поскольку такое препятствие не может быть
полностью охарактеризовано первыми двумя моментами. Кроме того, даже
одномерное распределение вероятности негауссовского шума не является
инвариантным относительно линейных преобразований.
Выбор оптимальной системы сигналов и схемы решения для данной
структуры помех теоретически рассматривается только для некоторых
частных случаев с дополнительными упрощающими предположениями. Такие
особые случаи включают идеализированные концентрированные или
импульсные шумы. Огромные трудности возникают при изучении каналов, в
которых одновременно происходят концентрированные, импульсные и
флуктуационные помехи.
На основе интуитивных соображений разработаны многие методы
подавления концентрированного и импульсного шума. В настоящее время
теория этого вопроса в лучшем случае объясняет суть используемых методов
и способна прояснить некоторые из их деталей.
Следует отметить, что методы, позволяющие подавлять хорошо
сфокусированные
помехи,
обычно
ухудшают
условия
подавления
импульсных шумов и наоборот. Это будет показано ниже в ряде примеров.
Кроме того, оба метода приводят к тому, что схема принятия решений не
является оптимальной для флуктуационных помех.
Радиоканалы в диапазонах длинной, средней и короткой волн всегда
подвержены большому количеству концентрированных шумов. Это является
следствием условий распространения радиоволн этих диапазонов, что
приводит к тому, что любое излучение создает значительную напряженность
поля на значительных расстояниях от передатчика. Большинство из этих
нарушений
имеют
относительно
небольшую
интенсивность.
12
Складывающиеся, они образуют общий шумный фон, который по своим
характеристикам мало чем отличается от обычного белого шума. Для такого
типа помех все результаты, полученные в предыдущих главах при изучении
флуктуационных помех, действительны, но они не являются основным
предметом
этой
главы.
Здесь
нас
будет
интересовать
отдельный
концентрированный шум, который выделяется на общем фоновом шуме и
соизмерим по мощности с полезным сигналом. Такое вмешательство
встречается во всех диапазонах, а при проектировании систем и оборудования
радиосвязи их возможность всегда учитывается.
С точки зрения механизма влияния концентрированных помех на прием
сигналов их можно разделить на три типа:
1. помехи, спектр которых сосредоточен в полосе частот, полностью или
частично совпадающей с полосой частот, занимаемой сигналом;
2. помехи, спектр которых лежит вне полосы частот сигнала, их часто
называют помехами «по соседним каналам»;
3. помехи, которые на входе приемного устройства имеют спектр,
лежащий вне полосы частот сигнала, но в результате нелинейных
преобразований
в
приемном
устройстве
образуют
составляющие,
попадающие в ту же полосу частот, что и сигнал;
Уменьшение потенциала таких нелинейных эффектов помех от каналов
приема боковых частот является одной из основных задач при разработке
радиоприемников. Это подробно обсуждается во всех руководствах по
радиоприемникам, поэтому здесь мы рассмотрим его только в той мере, в
какой это связано с введением в канал приемника специальных нелинейных
элементов для защиты от других типов помех.
Быстрое
развитие
радиосвязи,
а
также
других
применений
радиоэлектроники привело к тому, что освоенные диапазоны радиочастот
были перегружены различными излучениями. В результате вмешательство во
многих
случаях
преобладает
над
всеми
остальными
факторами,
ограничивающими точность и реальную полосу пропускания радиоканалов.
13
Основные методы защиты от концентрированных помех, применяемые
с самого начала развития радиосвязи до настоящего времени, основаны на
использовании
пространственная
частотной
селективности.
избирательность,
Широко
обеспечиваемая
используется
направленными
передающими и приемными антеннами. Хотя еще в 1930-х годах было
показано, что частотная и пространственная избирательность отнюдь не
являются единственными методами выделения полезного сигнала от помех,
однако первые попытки реализовать другие возможности относятся только к
самому последнему времени.
Использование частотной избирательности для устранения взаимных
помех подразумевает наличие хорошо известной регулировки частот,
назначенных различным каналам связи. В «идеальном» случае казалось бы,
что можно полностью решить проблему взаимных помех, если каждый канал
связи снабжен определенной полосой частот, запрещая любое другое
излучение в этой полосе. Однако такое решение невозможно по ряду причин.
Основными
препятствиями
для
идеального
регулирования
являются
следующие
1. Существующие международные соглашения по регулированию
использования радиочастотных полос часто и по разным причинам
нарушаются;
2. Строгое регулирование, при котором каждому радиоканалу
назначается собственный диапазон частот, свободный от других выбросов,
невозможно, поскольку количество существующих каналов во многих
полосах частот намного больше, чем количество полос, которые можно
выделить в наиболее экономичном распределении , с учетом возможной
пространственной селективности. Кроме того, необходимое количество
каналов постоянно растет, и скорость передачи информации в этих каналах
также увеличивается, что обычно приводит к необходимости расширения
выделенных полос частот;
14
3. Полоса частот, в которой сосредоточена основная часть энергии
сигнала, существенно шире, чем условная полоса частот. Это необходимо
учитывать при распределении полос частот, обеспечивая наличие известного
запаса или «защитных» полос частот, которые не используются эффективно
для
передачи
информации.
Эти
защитные
полосы
должны
быть
дополнительно расширены из-за возможных неточностей и нестабильности
частот сигнала;
Во многих случаях помехи не генерируются первичным излучением
передатчика, а различными побочными излучениями, такими как гармоники
или комбинационные частоты, которые происходят в пути генерации сигнала.
Фактически, для любого радиочастотного распределения одни и те же
полосы распределяются по нескольким каналам. Это учитывает взаимное
расположение приемников и передатчиков, так что на входе каждого
приемника рассчитанные напряженности поля помех существенно меньше,
чем напряженность поля полезного сигнала. Однако условия распространения
радиоволн широко варьируются, особенно с радиосвязью ионосферы, и не
могут быть точно предсказаны. Поэтому никогда нельзя быть уверенным, что
уровень помех не будет увеличиваться до такой степени, когда точность
приема ниже приемлемой. Значительно более надежными в этом отношении
являются радиоканалы, использующие
распространение радиоволн
в
свободном пространстве или вдоль земли (например, микроволновые каналы)
и, конечно же, каналы проводной связи.
В тех случаях, когда помехи мешают работе канала, часто необходимо
изменить используемый диапазон частот, что часто приводит к помехам в
других рабочих каналах.
В связи с концентрированным шумом возникают различные проблемы,
которые важны для дальнейшего развития коммуникационных технологий.
Среди них следует отметить проблему снижения эффекта концентрации помех
в текущей ситуации в отношении регулирования радиочастот и используемых
в настоящее время сигналов. Эта проблема включает в себя вопросы
15
построения схемы принятия решений, применения приема разнообразия и т.
Д. Это также включает вопросы, связанные с устранением побочных
излучений и ложных каналов приема, а также вопросы применения
пространственной избирательности, которые не рассматриваются Вот.
Второй круг задач - изучение возможностей снижения взаимного
вмешательства в строительство новых систем связи.
Многие из этих задач еще не были удовлетворительно решены.
Существенным препятствием на пути развития теории в этой области является
отсутствие достаточных статистических данных по концентрированному
шуму.
Получение
этих
данных
осложняется
тем,
что
характер
концентрированного шума различен в разных диапазонах радиочастот, в
разное время суток и даже в разное время года. Очевидно, что распределение
концентрированных помех также зависит от наличия и степени соответствия
регулированию
использования
радиочастот.
Поэтому
мы
должны
использовать грубые приближенные представления, которые не всегда
приводят к достоверным результатам.
1.3 Типы синхронизации
Трудности, связанные с синхронизацией, начались сразу с появлением
цифровых методов передачи информации. Буквально каждая процедура
оцифровки, передачи и приема информации в виде двоичного сигнала или
закодированного двоичного сигнала требует согласованности частот
передачи и приема, в противном случае передаваемая информация будет
приниматься некорректно. Сложность синхронизации не ограничивается
только цифровой первичной сетью, но она важна при изучении сетей ISDN,
передачи данных (SPD), цифровой телефонии, сетей специального назначения
и других вторичных сетей. В технологии современной коммуникации
существуют три основные концепции синхронизации:
1. Частотная;
16
2. Фазовая;
3. Временная;
Наиболее важным типом синхронизации для первичной сети является
частотная синхронизация или синхронизация частоты (ES), которая
определяет согласованность генераторов цифровых устройств в сети по
частоте. В идеальных условиях в этом случае все сетевые генераторы
работают на одной и той же частоте, скорость передачи цифровой
информации с высокой степенью точности равна скорости приема, в
результате отсутствует потеря информации из-за проскальзывания в системы
связи, то есть ошибок, связанных с нарушениями синхронизации, нет.
Конечно, это основная цель эксплуатации - добиться наименьшего уровня
ошибок в сети. Поэтому частотная синхронизация является основным
интересом операторов связи.
Согласование фаз (PS) относится к соответствию фаз принимающих и
передающих сигналов. Эта синхронизация приобрела большое значение в
различных электронных устройствах. Достижение использует различные
компоненты (фазовращатели), фазовые детекторы и т. Д. В современной
практике систем связи FS часто используется в схемах приемника сигналов,
где достигается фазовая синхронизация между линейным сигналом и схемой
приемника.
Синхронизация времени (BPS) или синхронизация по времени
подразумевают, что все устройства в сети имеют единовременный набор
времени. Это время соответствует универсальному скоординированному
времени (UTC - скоординированное универсальное время). Чаще всего задача
синхронизации времени связана с различными вторичными сетями и
системами SPD (сети специального назначения, биллинговые системы и т. Д.).
Например, вы можете рассмотреть довольно известный метод защиты
электронных банковских сетей, когда каждому документу присваивается
точная метка времени. На принимающей стороне эта метка сравнивается с
текущим временем и анализируется основная возможность преобразования
17
электронного документа при его передаче по сети. Если задержка передачи
превышает установленный порог, документ считается недействительным. Эта
схема защиты является сильной, но необходимо, чтобы все рабочие станции в
сети синхронизировались во времени. Следует отметить, что BPS - это задача,
полностью независимая от частотной синхронизации. В современной
практике построения корпоративных сетей существует всего несколько
методов для комплексного решения обеих проблем.
1.4 Роль синхронизации в современных сетях
Основной
задачей
синхронизации
частот
является
достижение
совпадающих или множественных частот генераторов (тактовых частот) всех
цифровых устройств, поступающих в систему цифровой связи. Чтобы
добиться синхронизации в сети, вы должны:
1. установить единую тактовую частоту для всей системы связи, чтобы
система работала с одинаковой скоростью, обеспечить, чтобы все цифровые
устройства в сети работали синхронно;
2. Сетевая синхронизация должна поддерживаться в любое время
независимо от изменений в структуре сети, вызванных такими факторами, как
сбои в каналах и узлах, перестройка сети, расширение и т. Д .
3. Система синхронизации (CC) также должна быть достаточно
стабильной в отношении таких нарушений, как изменение частоты или фазы в
узле, изменение времени передачи на линии и т. Д.;
Из вышесказанного определяется интегральный смысл синхронизации для систем синхронизации нет абсолютно никакого деления на первичные и
вторичные сети, синхронизация должна охватывать все цифровые устройства
сети. В любом случае, современные методы проектирования и построения
телекоммуникационной системы определяют установленный приоритет при
проектировании систем синхронизации. Поскольку в телекоммуникационной
системе первичная сеть является ядром сети и создает каналы для вторичных
18
сетей, SS наиболее правильно построен в следующем порядке: сначала
создается первичная сеть CC, а затем вторичная сеть синхронизируется с
первичной сети. Таким образом, SS начинается с первичной сети.
На данный момент наиболее актуальными проблемами являются
синхронизация первичной сети. Современная первичная сеть ориентирована
на использование технологии SDH. В отличие от технологии PDH, где
предусматривался режим, плезиохронная («почти синхронная») работа
различных устройств систем передачи данных, технология SDH обеспечивает
синхронную работу всей сети и всех устройств, входящих в сеть.
Вторым фактором, влияющим на рост роли синхронизации в
современных системах цифровой связи, является развитие последних. В
общем случае для синхронизации нет первичной и вторичной сетей для
синхронизации, все цифровые устройства в телекоммуникационной системе
должны использоваться в системе синхронизации. В результате уровень
проблем, возникающих при синхронизации, определяется в зависимости от
количества цифровых устройств в сети. С небольшим количеством цифровых
устройств трудности с синхронизацией могут быть решены конфиденциально,
и CC не обрабатывается операторами как отдельная система внутри сети.
Однако, с увеличением числа цифровых устройств, проблемы синхронизации
больше не могут рассматриваться частично и должны рассматриваться
систематически. С этим связана определенная локальная революция в
подходе: появление определенной «критической массы» цифровых устройств
в сети приводит к необходимости рассматривать систему синхронизации как
отдельный компонент телекоммуникационной системы. При дальнейшем
увеличении числа цифровых устройств концепции построения систем
синхронизации начинают меняться.
1.5 Тактовая синхронизация
19
Широкая доступность устройств, обеспечивающих синхронную работу
приемных и передающих станций цифровых систем передачи, а также
высококачественную работу ретрансляторов станций и линий, называется
тактовыми синхронизаторами (TCB). Реализация TCB во многом зависит от
методов создания тактового сигнала и его использования для синхронизации
терминальных станций. Для этого вы можете использовать:
1.
Передача
синхронизирующего
сигнала
(или
периодической
последовательности импульсов с тактовой частотой) по специальной паре
сердечников кабеля или по отдельному каналу. Этот метод является
неэкономичным и, кроме того, связан со значительными трудностями,
связанными с точной коррекцией фазовых соотношений для получения того
же группового времени прохождения (GWP) для частотных компонентов,
несущих оперативную и синхронизирующую информацию;
2. Выбор тактовой частоты сигнала из группового сигнала ИКМ на
принимающей станции. Операция фильтрации или извлечения тактовой
частоты из группового сигнала основана на том факте, что в спектре
случайной однополярной последовательности импульсов, длительность
которых t короче, чем продолжительность их времени T, регулярная
(периодическая) составляющая, тактовая частота. Часть TCF, которая
обеспечивает фильтрацию тактовой частоты, называется распределителем
тактовой частоты (WTCH);
Сигнал групповой PCM подается на узкополосный фильтр тактовой
частоты (UPFTCH), который извлекает тактовую частоту из регулярного
компонента сигнала группы PCM. В простейшем случае UPFT представляет
собой параллельный резонансный контур. Под воздействием импульсов
сигнала ИКМ на выходе схемы формируются затухающие гармонические
колебания с тактовой частотой. Кроме того, эти гармонические колебания
усиливаются и ограничены по амплитуде усилителю-ограничителю (VO). На
выходе VO мы имеем гармонический сигнал, ограниченный по амплитуде (с
достаточно крутыми фронтами). После прохождения этого сигнала через
20
дифференцирующую
схему
(DC)
мы
получаем
биполярную
последовательность стробирующих импульсов, которая преобразуется в
последовательность однополярных импульсов, которые следуют за тактовой
частотой fm.
После
этого
эта
последовательность
отправляется
генератору
периодической последовательности прямоугольных импульсов (PFTP),
заданной амплитуды и следующей с тактовой частотой fm. FPPI, в простейшем
случае, представляет собой ожидающий мультивибратор, который генерирует
импульсы заданной формы, длительности и амплитуды, частота которых
определяется частотой возбуждающего сигнала.
Для настройки фаз тактовой частоты тактового сигнала и тактовой
частоты синхронизированного сигнала используется регулируемая линия
задержки (RLZ).
Рассматриваемая схема TCB использует так называемый резонансный
метод выделения тактовой частоты, преимуществами которой являются
простота реализации и улучшение экономических показателей системы, что
важно для систем передачи локальных или зональных первичных сетей.
Минусы резонансного типа UTS включают:
1. быстрая потеря тактовой частоты во время разрывов связи или при
появлении длинной серии нулей в принятом групповом сигнале ИКМ;
2. зависимость стабильности выведенной тактовой частоты от характера
кодовых комбинаций или структуры и скорости передачи линейного
цифрового сигнала.
1.6 Символьная синхронизация
При синхронизации символов, для оптимальной демодуляции, все
цифровые приемники должны быть синхронизированы с переходами
входящих цифровых символов.
21
Класс синхронизаторов, используемых в этом случае, называется
синхронизаторами без данных (NDA). Существует еще один класс
символьных синхронизаторов, которые используют известную информацию
об информационном потоке. Эта информация может быть извлечена из
решений обратной связи относительно полученных данных или из известной
последовательности, введенной в информационный поток.
В настоящее время методы, использующие данные (DA), становятся
более важными и доминируют в выборе модуляции, которые эффективно
используют полосу.
Рассматриваемые синхронизирующие символы можно разделить на две
основные группы. Первая группа состоит из открытых синхронизаторов. Эти
схемы обращают внимание на копию выходного сигнала передатчика
непосредственно из входящего информационного потока. Вторая группа закрытые синхронизаторы; они синхронизируют локальный генератор
тактовых импульсов с входящим сигналом путем сравнения локального и
поступающего сигналов. Закрытые синхронизаторы, как правило, более
точные, но в то же время более сложные и дорогие.
Открытые символьные синхронизаторы также иногда называются
нелинейными синхронизаторами на фильтрах; это имя говорит само за себя.
Синхронизаторы этого класса генерируют частотную составляющую со
скоростью символа, пропуская входящий узкополосный сигнал через
последовательность фильтров и нелинейное устройство. Работа этого
устройства аналогична восстановлению носителя в схеме поддержки с
подавленной несущей. В этом случае желаемый частотный компонент,
передаваемый со скоростью символа, выделяется полосовым фильтром, после
чего насыщающий усилитель с высоким коэффициентом насыщения дает
желаемую форму. В результате восстанавливается прямоугольный сигнал
генератора тактовых импульсов.
Основным недостатком открытых синхронизирующих символов
является наличие неизбежной ошибки обслуживания с ненулевым средним
22
значением. Эта ошибка может быть уменьшена при высоких соотношениях
сигнал-шум, но поскольку форма сигнала синхронизации напрямую зависит
от входящего сигнала, никогда не удастся исправить ошибку.
Замкнутые символьные синхронизаторы сравнивают входной сигнал с
локально
генерируемым
сигналом
и
последующей
синхронизацией
локального сигнала с переходами во входном сигнале. Фактически, эта
процедура ничем не отличается от процедуры, используемой в открытых
синхронизаторах.
Среди
наиболее
популярных
синхронных
синхронизаторов
с
замкнутым циклом вы можете выбрать синхронизатор раннего / позднего
затвора, который состоит в выполнении двух отдельных интеграций энергии
входного сигнала в течение двух разных интервалов интервала символов с
длительностью (T-d) секунд. интеграция (расширенная) начинается с
момента, указанного в начале периода передачи символа (условно - время 0),
и заканчивается на (Td) секунд. Вторая итерация (отставание) начинается с
момента задержки d секунд и заканчивается в конце символа передачи
(shareware - время T). Разница между абсолютными значениями, описанными
в выводах итератора y1 и y2, является мерой ошибки синхронизации символа
приемника, которая может быть возвращена для последующей коррекции
приема.
23
ГЛАВА II
ФОРМИРОВАНИЕ СИГНАЛОВ, УСТОЙЧИВЫХ К
ВОЗДЕЙСТВИЯМ КОРОТКОЖИВУЩИХ ИЛИ
СОСРЕДОТОЧЕННЫХ ПО СПЕКТРУ ПОМЕХ
2.1 О классе сигналов, допускающих восстановление одних
фрагментов по остальным
Пусть далее x(t ), 0 t T , вещественнозначная непрерывная функция
времени (сигнал) с ограниченной евклидовой нормой
T
|| x || = x 2 (t )dt .
2
(2.1)
0
Положим
x1 (t ) = x(t ), 0 t T1; x2 (t ) = x(t + T1 ), 0 t T2 = T − T1 . (2.2)
.
Задача заключается в разработке достаточно общего подхода к
построению класса сигналов, позволяющих вычислить один из определенных
в (2.2) фрагментов при наличии другого. Такие вычисления должны
осуществляться на основе некоторых функционалов, которые ниже
представлены в общем виде
x1 (t ) = F1 (t , x2 ( ), 0 T2 ), 0 t T1 .
(2.3)
x2 (t ) = F2 (t , x1 ( ), 0 T1 ), 0 t T2 .
(2.4)
24
Утверждение 1.
Сигналы, допускающих восстановление одних
фрагментов по остальным,. могут быть сформированы на основе
собственных функций симметричных операторов.
Доказательство. Без потери общности рассмотрим построение таких
сигналов и функционалов для восстановления начальных из их фрагментов
(2.3).
Пусть A(s,u), 0 s, u T - симметричное ядро
A(s,u) = A(u, s), 0 s, u T .
(2.5)
Некоторого интегрального оператора, удовлетворяющее условиям
равномерного разложения в ряд [3]
A(s,u) = k g k ( s) g k (u ) .
(2.6)
k =1
По ортонормальным собственным функциям
T
k g k ( s ) = A( s, u ) g k (u )du , 0 s T ; k = 1,2,... .
(2.7)
0
T
( g k , gi ) = g k (t ) gi (t )dt = 0, i k .
(2.8)
0
|| g k ||2 = (g k , g k ) = 1.
Для
дальнейшего
собственные
функции
изложения
(2.9)
существенным
симметричного
ядра
ортонормированный базис пространства L2 [3].
является
составляют
то,
что
полный
25
Положим для определенности
x(t ) cg1 (t ), 0 t T .
(2.10)
Здесь c - некоторая постоянная.
Тогда в соответствии с определением (2.7) нетрудно получить
следующее соотношение
T1
1 x1 (t ) − A11 (t , u ) x1 (u )du = y1 (t ) .
(2.11)
0
где
T2
y1 (t ) = A12 (t , u ) x2 (u )du .
(2.12)
A11 ( s, u ) = A( s, u ), 0 s, u T1 .
(2.13)
A12 ( s, u ) = A( s, u + T1 ), 0 s T1; 0 u T2 .
(2.13)
0
Ясно, что соотношение (2.11) представляет собой интегральное
уравнение Фредгольма второго рода [4] относительно фрагмента (2.3).
Возможность вычисления его решения определяется свойствами ядра
интеграла в (2.11) и величиной 1 .
Для завершения доказательства утверждения построим представление
для искомого фрагмента сигнала, которое при выполнении некоторых условий
позволяет осуществить соответствующие вычисления.
Очевидно, что ядро в (2.11) также является симметричным. Это
позволяет воспользоваться представлением вида (2.6)
26
A11(s,u) = 1k g1k ( s) g1k (u ), 0 s,u T1 .
(2.15)
k =1
T1
1k g1k ( s) = A11 ( s, u ) g1k (u )du, 0 s T1; k = 1,2,... .
(2.16)
0
T1
( g1k , g1i ) = g1k (t ) g1i (t )dt = 0, i k .
(2.17)
0
|| g1k ||2 = 1 .
(2.18)
Ясно, что рассматриваемые сигналы в рамках теории пространств L2
являются элементами линеалов [5]. Поэтому в виду полноты базиса
собственных функций ядра (2.13) справедливы представления
x1 (t ) = 1k g1k (t ) .
(2.19)
k =1
y1 (t ) = 1k g1k (t ) .
(2.20)
k =1
коэффициенты которых являются скалярными произведениями вида
1k = ( x1 , g1k ) ; 1k = ( y1 , g1k ) .
(2.21)
.
Подстановка представлений (2.15), (2.19) и (2.20) в уравнение (2.11) дает
тождество
27
k =1
[1k (1 − 1k ) − 1k ) g1k (t ) 0, 0 t T1
(2.22)
Необходимым и достаточным условием достижения которого является
выполнение равенств
1k (1 − 1k ) = 1k , k = 1,2,...
(2.23)
Предположим, что ни одно из собственных чисел ядра (2.13) не
совпадает с 1
1k 1 , k = 1,2...
(2.24)
Тогда с учетом равенств (2.23) представление(2.19) для искомого
фрагмента сигнала принимает вид
x1 (t ) = 1k g1k (t ) / (1 − 1k )
(2.25)
k =1
Таким
образом,
полученное
соотношение,
определяющее
при
выполнении условий (2.24) один из фрагментов сигнала вида (2.10) через
другой. Доказательство закончено.
Справедливо следующее.
Утверждение2.
Пусть
ядро
(2.6)
является
положительно
определенным [4] и его собственные числа упорядочены по убыванию
1 2 .... 0
(2.26)
причем максимальное из них 1 является единственным. Тогда для
собственных чисел ядра (2.15) выполняются неравенства
28
1k 0 ; 1 1k , k = 1,...
(2.27)
Причем знак равенства в последнем неравенстве соответствует
случаю совпадения собственных функций ядер (2.6) и (2.15).
Доказательство. Нетрудно показать, что при любом T1 T ядро (2.13)
является положительно определенным.
В самом деле, положим
z(t) = g1k (t ), 0 t T1; z (t ) 0, T1 t T
(2.28)
Тогда в виду положительной определенности исходного ядра нетрудно
получить неравенство
T
d (T1 ) =
0
T
0
T1 T2
A(t1 , t2 ) z (t1 ) z (t2 ) dt1dt2 =
0
A11 (t1 , t2 ) g1k (t1 ) g1k (t2 ) dt1dt2 0 (2.28.1)
0
В соответствии с соотношениями
(2.16) и (2.18) отсюда получаем
неравенства
d(T1 ) = 1k 0, k 1
что
(2.29)
является необходимым и достаточным условием положительной
определенности ядра (2.13). С другой стороны для любого нормированного
сигнала y(t), 0 t T ; || y ||= 1 выполняется вариационное неравенство
T
T
0
0
1 =
T
T
0
0
A(t1 , t2 ) g1 (t1 ) g1 (t2 )dt1dt2 = max
A(t1 , t2 ) y (t1 ) y (t2 )dt1dt2
(2.30)
29
Сопоставление (2.30) и (2.29) показывает, что неравенства (2.27)
выполняются. Это и завершает доказательство последнего утверждения.
Следствие. Если собственные функции ядер (2.6) и (2.15),
советующие
максимальным
собственным
числа
не
совпадают,
то
соотношение (2.25) является корректным математическим представлением
фрагмента сигнала вида (2.10) любой длительности, удовлетворяющей
неравенству T1 T .
Таким
образом,
представление
(2.25)
может
служить
теоретической основой для вычислений фрагментов сигналов вида (2.10) по
оставшемуся фрагменту. Ясно, что реализация таких вычислений требует
дискретизации области определения ядра.
В соответствие с этим ядра аппроксимируются с помощью следующих
матриц
В
A = {aik } = (t )2{ A(it , k t )}, i, k = 1,..,N
(2.31)
A11 = {aik }, i, k = 1,..., M N
(2.32)
A12 = (t )2{ A(it ,( M + n)t )}, i = 1,..., M ; n = 1,..., N − M
(2.33)
соответствии
с
этим
целесообразно
рассматривать
дискретизированные сигналы
xk = x(k t ), k = 1,.., N ; N t = T
(2.34)
x = ( x1,.., xM , xM +1,.., xN )' = ( x1' , x2' )' ; M t = T1
(2.35)
- где штрих означает транспонирование.
30
Выбор шага дискретизации осуществляется на основе анализа
конкретных ядер, с обеспечением, например, свойства положительной
определенности получаемых матриц. Важно, то, что матрицы (2.31)
получаются симметричными, наборы собственных векторов которых
образуют полные в смысле пространства N
-мерных
вещественных
векторов ортонормальные базисы[6,7].
Собственные числа матриц (2.31) и (2.32) ниже для простоты
обозначаются теми же символами, что и у аппроксимируемых ими ядер, если
это не вызывает затруднений в их интерпретации.
Если отсчеты сигнала выбираются на основе собственного вектора
матрицы (2.31), соответствующего максимальному собственному числу
x = cq1
(2.36)
то нетрудно показать, что уравнение (2.11) переходит в матричное
соотношение
(1I − A11 ) x1 = A12 x2
Если
соответствующие
максимальным
(2.37)
собственным
числам
собственные векторы матриц (2.31) и (2.32) не совпадают, то представление
решения уравнения (2.37) имеет вид
x1 = (1I − A11 )−1 A12 x2
(2.38)
Полученные соотношения можно обобщить на случай, когда среди
собственных чисел матрицы (2.31) имеются кратные. Не нарушаю общности,
предположим, что выполняются равенства
31
1 = 2 = ... = M =
(2.39)
При этом, как известно [6], каждому из этих собственных чисел будет
соответствовать ортогональные собственные векторы. Поэтому их любая
линейная комбинация вида
M
x = cm qm
(2.40)
m =1
будет также собственным вектором матрицы (2.31), соответствующим
собственному числу (2.39), так что соотношение (2.37) принимает вид
( I − A11 ) x1 = A12 x2
(2.41)
Входящие сюда переменные определены выше. В свою очередь, если
соответствующие максимальным собственным числам собственные векторы
матриц (2.31) и (2.32) не совпадают, то представление решения уравнения
(2.41) имеет вид
x1 = ( I − A11 )−1 A12 x2
(2.42)
Отметим, что по аналогии нетрудно также получить соотношения для
восстановления любого из фрагментов, занимающего сплошную подобласть
внутри области определения исходного сигнала. При этом очевидным образом
должны быть изменены определения ядра (2.13) и функции (2.12) и их
дискретизованные аналоги.
Следует только иметь в виду, что определить условия выполнения
неравенств типа (2.24) по аналогии с Утверждением 2 при наличии нескольких
одинаковых собственных чисел затруднительно.
32
Для иллюстрации предлагаемого подхода к формированию сигналов с
восстанавливаемыми
фрагментами
воспользуемся
вычислительными
экспериментами. В качестве элементов симметричных матриц дискретной
аппроксимации
рассматриваемой
задачи
восстановления
используем
следующие:
В
качестве
aik = exp(−4ik / N 2 );
(2.43)
aik = cos(4 (i − k ) / N );
(2.44)
aik = sin(4 (i − k ) / N ) / ( (i − k ))
(2.45)
aik = sin(32 (i − k ) / N ) / ( (i − k ))
(2.46)
важной
конкретизации
предлагаемого
подхода
к
формированию сигналов, позволяющих восстановить одни фрагменты по
остальным рассмотрим применение введенных в работе [8] субполосных ядер.
Целесообразность их использования обусловлена в частности тем, что они
построены на основе оптимизации процедур субполосного анализа и синтеза
сигналов и, как будет показано ниже, позволяют построить процедуры
восстановления одних фрагментов спектров Фурье по остальным.
Для сохранения цельности изложения представляется необходимым
привести здесь некоторые сведения из теории субполосных ядер.
Согласно [8] субполосное ядро определяется для частотного интервала
(полосы)
= [-2 , −1 ) [1 , 2 ), 2 1 0
на основе преобразования Фурье
(2.47)
33
A (t ) =
exp(− jzt )dz / 2
(2.48)
z
что дает
A (t − ) = 2 A0 (t − )cos(s (t − ), A (0) = A0 (0) = 2 D /
(2.49)
где
A0 (t − ) = sin( D(t − )) / (t − )
(2.50)
D = ( 2 − 1 ) / 2 ; s = ( 2 + 1 ) / 2 .
(2.51)
Отметим, что ядра вида (2.50) рассматривались ранее в связи с
исследованием
класса
функций
с
двойной
ортогональностью
(см.
например[9]).
Из определения (2.48) следует, что субполосное ядро (2.49) является
симметричным
и
положительно
определенным.
Поэтому
для
него
справедливы все выводы, изложенные в предыдущем разделе.
Рассмотрим
теперь
проблему
восстановления
пораженных
воздействиями узкополосных помех участков спектров сигнала вида (2.10).
На основе определений (2.7) и (2.48) нетрудно получить следующее
представление для собственных функций субполосного ядра
k g k (t ) =
Gk ( z )exp( jzt )dz / 2
(2.52)
z
через отрезок их преобразований (трансформанта, спектр) Фурье
34
T
Gk ( z ) = g k (t )exp( − jzt )dt
(2.53)
0
Справедливы следующие соотношения [8]
k =
| G ( z ) |2 dz / 2 1, k = 1,...
(2.54)
z
причем с высокой точностью выполняются равенства
k 0. k J
(2.55)
J = 2[TD / ] + 4
(2.56)
где
Справедливо следующее.
Утверждение 3. Фрагменты спектров собственных функций
ядра (2.49), соответствующих ненулевым собственным числам, вне
частотного интервала (2.47) можно вычислить по фрагментам
внутри него на основе следующего соотношения
k Gk ( y ) =
Gk ( z )(1 − exp(− jT ( y − z ))) / j 2 ( y − z )dz
(2.57)
z
Доказательство.
умножения
слева
и
Соотношение
справа
(2.56)
представления
получается
(2.52)
на
в
результате
exp ( − jyt )
и
интегрирования в соответствии с определением (2.53) по переменной t в
пределах 0 t T . Очевидно, что соотношение (2.56) определяет способ
35
аналитического продолжения спектров собственных функций из выбранного
интервала (2.47) на всю частотную ось.
Иными словами оказывается возможным вычислять фрагменты
спектров собственных векторов вне выбранного частотного интервала по
фрагменту внутри него, что и завершает доказательство Утверждения 3.
Более сложно решается задача восстановления фрагментов спектра
внутри выбранного частотного интервала (2.47).
Положим
Gk ( y) = exp( jTy / 2)Gk ( y)
(2.58)
Ty / 2 = u − T 1 / 2 = u − 1
(2.59)
Tz / 2 = v − T 1 / 2 = v − 1
(2.60)
Fk (u ) = Gk (2u /T + 1 ) , D = TD
(2.61)
Отметим, что представление (2.57) соответствует в (2.53) замене
(простому сдвигу на − T / 2 ) переменной интегрирования. При этом из (2.56)
получаем соотношение
k Fk (u ) =
D
Fk (v)sin(u − v) / ( (u − v))dv
(2.62)
−D
показывающее, что функции (2.60) являются собственными функциями
его подынтегрального ядра.
Как было показано в предыдущем разделе, возможность восстановления
одних фрагментов функций (2.60) по другим фрагментам определяется
свойствами собственных чисел ядра в уравнении (2.61).
36
Дискретизация области определения по обеим переменным с шагом
D = D / N
(2.63)
позволяет получить аппроксимирующую ядро матрицу
S = {sik }, i, k = 1,...,2 N
(2.64)
с элементами
sik = sin( D(i − k )) / (i − k )
(2.65)
С учетом (2.62) легко заметить, что представления (2.64), (2.45), (2.46)
идентичны. Отметим, что возможности этого определяются величиной
произведения вида D = TD .
37
ГЛАВА III
РАЗРАБОТКА ТЕХНОЛОГИИ ФОРМИРОВАНИЯ
СИГНАЛОВ НА БАЗИСЕ СОБСТВЕННЫХ ВЕКТОРОВ
СУБПОЛОСНОЙ МАТРИЦЫ УСТОЙЧИВЫХ К
ВОЗДЕЙСТВИЮ КОРОТКОЖИВУЩИХ ИЛИ
СОСРЕДОТОЧЕННЫХ ПО СПЕКТРУ ПОМЕХ
3.1 Исследование разработанной технологии
После формирования собственных векторов субполосной матрицы,
описанных
в
относительной
главе
2,
было
погрешности
проведено
восстановления
исследование
от
зависимости
количества
отчетов
неиспорченного фрагмента исходного вектора при добавление белого
гауссовского шума с различных соотношениях шум/сигнал:
N
q=
i =1
2
i
(3.1)
N
x
i =1
2
i
где - белый шум, x - неиспорченный фрагмент исходного вектора.
Относительная
погрешность
восстановления
рассчитывается
следующим образом:
N
r=
(x − x )
i
i =1
i
(3.2)
N
x
i =1
i
2
38
где
x - исходный собственный вектор,
x - восстановленный
собственный вектор.
В качестве исследуемых собственных векторов были выбраны первые 3
позиции собственных векторов субполосного ядра, описанного во главе 2.
В качестве соотношения шум/сигнал в исследовании были выбраны
следующие соотношения:
Таблица 3.1
Соотношения сигнал/шум, используемые в исследовании
Соотношение
шум/сигнал (q)
0
0.1
0.3
0.5
Номер собственного вектора – 1 ( K i = 1 )
В первую очередь
погрешности
было рассчитана зависимость относительной
восстановления
от
количества
отчетов
неиспорченного
фрагмента исходного вектора без добавления шума:
Таблица 3.2
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 1 , q = 0 )
M
511
R
M
2,74
E-14
501
1,89
E-14
491
1,21
E-14
481
8,00
E-15
471
5,73
E-15
461
4,16
E-15
451
3,28
E-15
441
2,55
E-15
431
2,12
E-15
421
1,81
E-15
411
1,65
E-15
401
1,53
E-15
391
1,52
E-15
R
381
1,32
E-15
371
1,25
E-15
361
1,27
E-15
351
1,81
E-15
341
2,04
E-15
331
2,84
E-15
321
2,97
E-15
311
3,39
E-15
301
3,49
E-15
291
4,66
E-15
281
5,62
E-15
271
7,23
E-15
261
7,89
E-15
M
251
R
8,90
E-15
241
1,16
E-14
231
1,44
E-14
221
1,71
E-14
211
2,13
E-14
201
2,51
E-14
191
3,03
E-14
181
3,82
E-14
171
4,81
E-14
161
6,02
E-14
151
7,76
E-14
141
9,86
E-14
131
1,23
E-13
M
121
111
101
91
81
71
61
51
41
31
21
11
1
39
R
1,60
E-13
На
2,05
E-13
2,65
E-13
рис.3.1
3,45
E-13
4,53
E-13
6,01
E-13
представлен
8,10
E-13
1,11
E-12
график
1,58
E-12
2,39
E-12
3,99
E-12
относительной
8,56
E-12
1,05
E-10
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 1 , q = 0 ):
Рис.3.1 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 1 , q = 0 )
На рис.3.2 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 1 , q = 0 ):
40
Рис.3.2 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 1 , q = 0 )
Как видно по графику 3.1 и по результатам из таблицы 3.2 погрешность
восстановления сигнала минимальна вплоть до использования минимального
количества отчетов неиспорченного фрагмента исходного вектора.
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.1
Таблица 3.3
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 1 , q = 0.1 )
M
511
R
M
0,38
88
501
0,03
24
491
0,16
27
481
0,06
78
471
0,02
68
461
0,02
70
451
0,01
69
441
0,08
42
431
0,06
55
421
0,06
33
411
0,03
12
401
0,01
48
391
0,01
64
R
381
0,06
67
371
0,05
04
361
0,02
00
351
0,02
93
341
0,03
54
331
0,03
30
321
0,02
32
311
0,01
92
301
0,08
31
291
0,00
94
281
0,03
12
271
0,02
64
261
0,02
28
M
251
R
0,00
32
241
0,02
07
231
0,00
38
221
0,02
35
211
0,00
33
201
0,02
70
191
0,01
39
181
0,04
87
171
0,02
34
161
0,01
30
151
0,01
69
141
0,01
83
131
0,05
51
41
M
121
R
0,00
32
111
0,02
07
На
101
0,00
38
рис.3.3
91
0,02
35
81
0,00
33
71
0,02
70
представлен
61
0,01
39
51
0,04
87
график
41
0,02
34
31
0,01
30
относительной
21
0,01
69
11
0,01
83
1
0,05
51
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 1 , q = 0.1 ):
Рис.3.3 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 1 , q = 0.1 )
На рис.3.4 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 1 , q = 0.1 ):
42
Рис.3.4 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 1 , q = 0.1 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.5
Таблица 3.4
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 1 , q = 0.5 )
M
511
R
M
1,05
87
501
1,20
80
491
0,35
55
481
0,70
63
471
0,61
35
461
0,07
38
451
0,11
07
441
0,21
37
431
0,33
43
421
0,04
12
411
0,24
65
401
0,07
56
391
0,07
61
R
381
0,06
26
371
0,06
03
361
0,17
52
351
0,09
73
341
0,08
27
331
0,18
03
321
0,02
55
311
0,07
83
301
0,10
04
291
0,02
09
281
0,00
86
271
0,40
56
261
0,13
57
M
251
R
0,08
21
241
0,13
21
231
0,03
67
221
0,19
31
211
0,09
95
201
0,22
55
191
0,03
95
181
0,13
74
171
0,15
23
161
0,30
03
151
0,20
50
141
0,12
77
131
0,19
10
M
121
R
0,14
69
111
0,01
55
101
0,03
21
91
0,14
59
81
0,27
32
71
0,14
87
61
0,23
42
51
0,25
23
41
0,04
76
31
0,01
61
21
0,07
93
11
0,13
07
1
0,25
96
43
На
рис.3.5
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума при М=511 ( K i = 1 , q = 0.5 ):
Рис.3.4 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 1 , q = 0.5 )
На рис.3.5 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 1 , q = 0.5 ):
44
Рис.3.5 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 1 , q = 0.5 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.8
Таблица 3.5
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 1 , q = 0.8 )
M
511
R
3,14
43
501
0,73
41
491
1,31
63
481
0,96
76
471
0,90
36
461
0,23
43
451
0,12
19
441
0,02
97
431
0,89
42
421
0,30
06
411
0,27
45
401
0,15
73
391
0,13
94
381
R 0,19
52
371
0,06
96
361
0,22
87
351
0,38
62
341
0,32
81
331
0,08
61
321
0,28
44
311
0,27
41
301
0,33
05
291
0,48
54
281
0,31
08
271
0,06
41
261
0,14
89
M
241
0,08
03
231
0,16
37
221
0,20
93
211
0,02
06
201
0,02
78
191
0,31
26
181
0,30
83
171
0,27
32
161
0,28
68
151
0,16
28
141
0,18
25
131
0,29
28
111
0,46
15
101
0,40
23
91
0,07
93
81
0,11
02
71
0,25
12
61
0,05
30
51
0,17
08
41
0,04
52
31
0,09
02
21
0,30
78
11
0,28
84
1
0,48
82
M
251
R 0,06
03
M
121
R
0,16
37
45
На
рис.3.6
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума при М=511 ( K i = 1 , q = 0.8 ):
Рис.3.6 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 1 , q = 0.8 )
На рис.3.7 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 1 , q = 0.8 ):
46
Рис.3.7 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 1 , q = 0.8 )
Номер собственного вектора – 2 ( K i = 2 )
В первую очередь была рассчитана зависимость относительной
погрешности
восстановления
от
количества
отчетов
неиспорченного
фрагмента исходного вектора без добавления шума:
Таблица 3.6
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 2 , q = 0 )
M
511
R
3,44
E-15
501
1,56
E-15
491
2,33
E-15
481
2,56
E-15
471
2,20
E-15
461
2,28
E-15
451
2,68
E-15
441
2,80
E-15
431
2,86
E-15
421
3,43
E-15
411
3,78
E-15
401
4,11
E-15
391
4,44
E-15
381
R 4,78
E-15
371
5,28
E-15
361
5,79
E-15
351
6,43
E-15
341
7,26
E-15
331
8,25
E-15
321
9,70
E-15
311
1,14
E-14
301
1,33
E-14
291
1,57
E-14
281
1,86
E-14
271
2,26
E-14
261
2,70
E-14
M
251
R
3,17
E-14
241
3,91
E-14
231
4,80
E-14
221
5,93
E-14
211
7,58
E-14
201
9,80
E-14
191
1,34
E-13
181
1,97
E-13
171
3,52
E-13
161
1,43
E-12
151
6,96
E-13
141
2,79
E-13
131
1,64
E-13
M
47
M
121
R
1,35
E-13
На
111
1,06
E-13
101
8,93
E-14
рис.3.8
91
7,21
E-14
81
6,59
E-14
71
6,19
E-14
представлен
61
6,75
E-14
51
7,90
E-14
график
41
9,86
E-14
31
1,31
E-13
21
1,93
E-13
относительной
11
3,39
E-13
1
7,62
E-13
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 2 , q = 0 ):
Рис.3.8 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 2 , q = 0 )
На рис.3.9 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 2 , q = 0 ):
48
Рис.3.9 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 2 , q = 0 )
Как видно по графику 3.8 и по результатам из таблицы 3.6 погрешность
восстановления сигнала минимальна вплоть до использования минимального
количества отчетов неиспорченного фрагмента исходного вектора.
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.1
Таблица 3.7
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 2 , q = 0.1 )
M
511
R
0,00
69
501
0,02
24
491
0,04
07
481
0,02
44
471
0,00
70
461
0,05
25
451
0,02
32
441
0,00
42
431
0,00
57
421
0,02
02
411
0,01
42
401
0,00
74
391
0,00
81
381
R 0,01
17
371
0,02
50
361
0,01
34
351
0,03
70
341
0,01
43
331
0,00
74
321
0,02
61
311
0,00
39
301
0,01
91
291
0,01
59
281
0,00
54
271
0,01
36
261
0,01
67
M
251
R
0,01
39
241
0,01
83
231
0,06
28
221
0,06
81
211
0,16
11
201
0,14
79
191
0,11
49
181
0,09
80
171
0,15
39
161
0,26
61
151
0,32
97
141
0,21
87
131
0,04
47
M
49
M
121
R
0,08
88
111
0,00
23
На
101
0,02
94
рис.3.8
91
0,02
97
81
0,06
84
71
0,03
22
представлен
61
0,02
95
51
0,01
97
график
41
0,01
84
31
0,00
88
21
0,01
07
относительной
11
0,04
07
1
0,01
56
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 2 , q = 0.1 ):
Рис.3.8 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 2 , q = 0.1 )
На рис.3.9 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 2 , q = 0.1 ):
50
Рис.3.9 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 2 , q = 0.1 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.5
Таблица 3.7
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 2 , q = 0.5 )
M
511
R
0,02
23
501
0,02
63
491
0,16
77
481
0,04
17
471
0,00
46
461
0,09
49
451
0,04
37
441
0,07
52
431
0,11
94
421
0,07
78
411
0,01
94
401
0,06
30
391
0,04
65
381
R 0,00
56
371
0,03
05
361
0,01
98
351
0,01
53
341
0,09
12
331
0,08
46
321
0,05
65
311
0,07
11
301
0,06
63
291
0,08
46
281
0,11
05
271
0,07
65
261
0,01
51
M
241
0,07
94
231
0,27
48
221
0,26
50
211
0,05
53
201
0,33
52
191
0,63
14
181
0,71
23
171
1,07
52
161
3,99
72
151
0,11
67
141
0,41
76
131
0,40
13
111
0,38
51
101
0,18
01
91
0,19
62
81
0,08
97
71
0,01
88
61
0,05
05
51
0,00
84
41
0,05
44
31
0,10
19
21
0,08
76
11
0,15
34
1
0,18
76
M
251
R 0,03
40
M
121
R
0,13
17
51
На
рис.3.10
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума при М=511 ( K i = 2 , q = 0.5 ):
Рис.3.10 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 2 , q = 0.5 )
На рис.3.11 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 2 , q = 0.5 ):
52
Рис.3.11 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 2 , q = 0.5 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.8
Таблица 3.8
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 2 , q = 0.8 )
M
511
R
0,52
85
501
0,03
00
491
0,46
88
481
0,07
64
471
0,04
31
461
0,12
41
451
0,14
77
441
0,02
32
431
0,04
86
421
0,14
93
411
0,08
39
401
0,15
01
391
0,06
65
381
R 0,05
98
371
0,11
62
361
0,03
21
351
0,19
25
341
0,08
85
331
0,10
89
321
0,16
09
311
0,22
35
301
0,08
47
291
0,10
77
281
0,08
57
271
0,03
27
261
0,27
27
M
241
0,08
18
231
0,01
33
221
0,57
61
211
0,57
16
201
1,17
49
191
0,68
29
181
1,24
58
171
0,34
93
161
9,22
27
151
2,70
42
141
0,97
15
131
0,37
09
111
0,19
54
101
0,07
96
91
0,15
67
81
0,05
36
71
0,06
75
61
0,17
84
51
0,19
06
41
0,02
27
31
0,05
09
21
0,17
51
11
0,23
36
1
0,32
77
M
251
R 0,36
83
M
121
R
0,77
62
53
На
рис.3.12
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума при М=511 ( K i = 2 , q = 0.8 ):
Рис.3.12 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 2 , q = 0.8 )
На рис.3.13 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 2 , q = 0.8 ):
54
Рис.3.13 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 2 , q = 0.8 )
Номер собственного вектора – 3 ( K i = 3 )
В первую очередь была рассчитана зависимость относительной
погрешности
восстановления
от
количества
отчетов
неиспорченного
фрагмента исходного вектора без добавления шума:
Таблица 3.9
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 3 , q = 0 )
M
511
R
3,56
E-15
501
2,15
E-15
491
2,09
E-15
481
2,02
E-15
471
2,10
E-15
461
2,21
E-15
451
2,17
E-15
441
2,26
E-15
431
2,57
E-15
421
2,70
E-15
411
2,68
E-15
401
2,94
E-15
391
3,20
E-15
381
R 3,49
E-15
371
3,76
E-15
361
4,04
E-15
351
4,50
E-15
341
4,80
E-15
331
5,21
E-15
321
5,46
E-15
311
6,03
E-15
301
7,15
E-15
291
2,98
E-14
281
8,02
E-15
271
7,12
E-15
261
6,38
E-15
M
251
R
7,97
E-15
241
6,43
E-15
231
6,83
E-15
221
6,74
E-15
211
5,78
E-15
201
6,87
E-15
191
5,44
E-15
181
7,85
E-15
171
9,17
E-15
161
1,01
E-14
151
1,88
E-14
141
9,72
E-14
131
7,18
E-14
M
55
M
121
R
4,42
E-14
На
111
2,49
E-14
101
2,95
E-14
рис.3.14
91
1,69
E-14
81
1,43
E-14
71
1,32
E-14
представлен
61
9,03
E-15
51
6,98
E-15
график
41
6,37
E-15
31
7,70
E-15
21
9,22
E-15
относительной
11
1,47
E-14
1
2,70
E-14
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 3 , q = 0 ):
Рис.3.14 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 3 , q = 0 )
На рис.3.15 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 3 , q = 0 ):
56
Рис.3.15 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 3 , q = 0 )
Как видно по графику 3.14 и по результатам из таблицы 3.9 погрешность
восстановления сигнала минимальна вплоть до использования минимального
количества отчетов неиспорченного фрагмента исходного вектора.
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.1
Таблица 3.10
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 3 , q = 0.1 )
M
511
R
0,00
59
501
0,00
05
491
0,00
12
481
0,01
47
471
0,00
33
461
0,01
29
451
0,00
17
441
0,00
48
431
0,00
03
421
0,00
95
411
0,00
16
401
0,02
55
391
0,00
59
381
R 0,00
29
371
0,00
77
361
0,02
37
351
0,01
57
341
0,01
37
331
0,02
50
321
0,00
40
311
0,06
76
301
0,11
39
291
0,44
74
281
0,04
97
271
0,02
59
261
0,04
53
M
251
R
0,05
31
241
0,01
19
231
0,01
93
221
0,00
63
211
0,01
00
201
0,03
98
191
0,01
23
181
0,01
25
171
0,05
13
161
0,01
54
151
0,04
76
141
0,06
50
131
0,29
81
M
57
M
121
R
0,01
11
111
0,03
00
На
101
0,01
12
рис.3.16
91
0,01
07
81
0,00
33
71
0,01
62
представлен
61
0,01
31
51
0,00
76
график
41
0,00
18
31
0,00
42
21
0,00
65
относительной
11
0,00
84
1
0,02
73
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 3 , q = 0.1 ):
Рис.3.16 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 3 , q = 0.1 )
На рис.3.17 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 3 , q = 0.1 ):
58
Рис.3.17 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 3 , q = 0.1 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.5
Таблица 3.11
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 3 , q = 0.5 )
M
511
R
0,02
45
501
0,04
68
491
0,01
33
481
0,03
25
471
0,07
25
461
0,02
83
451
0,00
74
441
0,01
91
431
0,00
86
421
0,03
68
411
0,02
72
401
0,02
45
391
0,04
68
381
R 0,06
11
371
0,10
56
361
0,17
04
351
0,24
19
341
0,23
72
331
0,16
06
321
0,21
14
311
4,05
12
301
0,12
56
291
0,10
60
281
0,45
41
271
0,06
11
261
0,10
56
M
241
0,03
90
231
0,05
88
221
0,06
15
211
0,09
31
201
0,23
32
191
0,02
43
181
0,16
00
171
0,28
27
161
0,90
95
151
0,66
16
141
0,06
42
131
0,03
90
111
0,02
49
101
0,08
42
91
0,10
62
81
0,02
58
71
0,02
18
61
0,11
98
51
0,03
45
41
0,07
75
31
0,04
03
21
0,07
26
11
0,22
90
1
0,02
49
M
251
R 0,06
42
M
121
R
0,22
90
59
На
рис.3.18
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума ( K i = 3 , q = 0.5 ):
Рис.3.18 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 3 , q = 0.5 )
На рис.3.19 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 3 , q = 0.5 ):
60
Рис.3.19 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 3 , q = 0.5 )
Далее к неиспорченному фрагменту исходного вектора было добавлен
белый шум с соотношением шум/сигнал – q = 0.8
Таблица 3.12
Зависимость относительной погрешности восстановления (M) от
количества отчетов неиспорченного фрагмента исходного вектора без
добавления шума (R) ( K i = 3 , q = 0.8 )
M
511
R
0,16
09
501
0,00
98
491
0,04
05
481
0,09
02
471
0,02
54
461
0,08
28
451
0,05
95
441
0,07
20
431
0,05
67
421
0,04
74
411
0,02
96
401
0,16
09
391
0,00
98
381
R 0,03
13
371
0,13
67
361
0,14
31
351
0,09
86
341
0,04
67
331
0,19
60
321
0,28
36
311
7,43
40
301
0,24
85
291
0,23
94
281
0,36
39
271
0,03
13
261
0,13
67
M
241
0,20
42
231
0,11
18
221
0,09
32
211
0,10
20
201
0,04
94
191
0,28
83
181
0,16
78
171
0,55
44
161
2,05
01
151
1,15
71
141
0,07
75
131
0,20
42
111
0,10
33
101
0,16
05
91
0,20
09
81
0,15
74
71
0,15
20
61
0,23
77
51
0,13
32
41
0,01
56
31
0,13
39
21
0,26
82
11
0,42
44
1
0,10
33
M
251
R 0,07
75
M
121
R
0,42
44
61
На
рис.3.20
представлен
график
относительной
погрешности
восстановления от количества отчетов неиспорченного фрагмента исходного
вектора без добавления шума при М=511 ( K i = 3 , q = 0.8 ):
Рис.3.20 График зависимости относительной погрешности
восстановления (R) от количества отчетов неиспорченного фрагмента
исходного вектора без добавления шума (M) ( K i = 3 , q = 0.8 )
На рис.3.21 представлен график исходного и восстановленного
собственного вектора без добавления шума при М=511 ( K i = 3 , q = 0.8 ):
62
Рис.3.21 График исходного и восстановленного собственного вектора
без добавления шума при М=511 ( K i = 3 , q = 0.8 )
По полученным результатам можно сделать следующие выводы:
1. Погрешность восстановления собственного вектора при отсутствии
добавления белого шума к отчетам неиспорченного фрагмента исходного
вектора дает полное восстановления вектора
2.
При
увеличении
соотношения
шум/сигнал
относительная
погрешность восстановления увеличивается
3. Чем более ортогонален собственный вектор субполосного ядра, тем
менее значимо на него влияет соотношение шум/сигнал, соответственно
использование данных векторов наиболее приемлемо в поставленной задаче
63
ЗАКЛЮЧЕНИЕ
В ходе выполнения выпускной работы была достигнута установленная
цель, а именно: была разработана технология формирования сигналов,
устойчивых к воздействию короткодвижущих или сосредоточенных по
спектру помех.
Для достижения цели были реализованы все поставленные задачи. Была
рассмотрена проблема формирования сигналов, одни из фрагментов которых
могут быть восстановлены по остальным. Показано, что класс таких сигналов
не пуст и может быть реализован на основе собственных функций
симметричных интегральных операторов. Получены интегральные уравнения
Фредгольма второго рода, которым удовлетворяют восстанавливаемые
фрагменты и установлены условия существования их решений.
По итогам реализации были получены следующие выводы:
1. Погрешность восстановления собственного вектора при отсутствии
добавления белого шума к отчетам неиспорченного фрагмента исходного
вектора дает полное восстановления вектора
2.
При
увеличении
соотношения
шум/сигнал
относительная
погрешность восстановления увеличивается
3. Чем более ортогонален собственный вектор субполосного ядра, тем
менее значимо на него влияет соотношение шум/сигнал, соответственно
использование данных векторов наиболее приемлемо в поставленной задаче
Рассмотрены вычислительные аспекты восстановления. Приведены
результаты вычислительных экспериментов, иллюстрирующие влияние на
доступные к восстановлениям длительности фрагментов свойств исходных
ядер (аппроксимирующих матриц).
В качестве важной конкретизации предлагаемого подхода рассмотрен
класс интегральных операторов на основе субполосных ядер [8]. Эти
интегральные операторы являются обобщением операторов, рассмотренных в
работе [9] Собственные функции этих ядер позволяют сформировать сигналы
64
не только с восстанавливаемыми фрагментами по другим фрагментам, но и
восстанавливать фрагменты их спектров Фурье. Получены соотношения,
позволяющие реализовать такие восстановления.
Актуальность полученных результатов обусловлена тем, что в
настоящее
время
существенно
возрастает
количество
источников
электромагнитных излучений, которые оказывают искажающее воздействие
друг на друга. В частности воздействиям электромагнитных излучений,
которые
генерируются
посторонними
источниками,
подвержены
управляющие и коммуникационные сигналы. Разработанный в статье подход
к формированию сигналов является вкладом в решение проблемы
электромагнитной совместимости при реализации взаимодействия на
расстоянии
65
СПИСОК ИСПОЛЬЗОВАННЫХ ЛИТЕРАТУРЫ
1.
Жиляков, Е.Г. Оптимальные канальные сигналы при цифровой
передаче с частотным уплотнением [Текст] / Е.Г. Жиляков, С.П. Белов, Д.В.
Урсол // Научные ведомости БелГУ Серия: Информатика, 2009. – № 7(62). ¬–
С. 166-172.
2.
Жиляков, Е.Г. Вариационные методы анализа и построения
функций по эмпирическим данным [Текст] / Е.Г. Жиляков. Белгород: Изд-во
БелГУ, 2007. – 160 с.
3.
Эммануил, С.А. Цифровая обработка сигналов: Практический
подход [Текст] / С. А. Эммануил, У. Барри. Москва: Вильямс, 2008. – 282 с.
4.
Тихвинский, В.В. Сети мобильной связи LTE. Технологии и
архитектура [Текст] / В.В. Тихвинский, С.А. Терентьев, А.А. Юрчук. Москва:
Эко-Трендз, 2013. - 200 с.
5.
Гмурман, В.Е. Теория вероятностей и математическая статистика
[Текст]: Учебное пособие для вузов / В.Е. Гмурман. Москва: Высшая школа,
2004. – 479 с.
6.
Берлин, А.В. Цифровые сотовые системы связи [Текст] / А.В.
Берлин. Москва: Эко-Трендз, 2012. – 133 с.
7.
Гордиенко, В. Н. Многоканальные телекоммуникационные
системы. Учебник для вузов [Текст] / В.Н. Гордиенко М.С Тверецкий. Москва:
Горячая линия-Телеком, 2007. – 88с.
8.
Слепов, Н.Н. Синхронные цифровые сети SDH [Текст] / Н.Н
Слепов. Москва: Эко-трендз, 2012. – 122 с.
9.
Соболев, Ю.В. Системы синхронизации в цифровых системах
связи [Текст] / Ю.В. Соболев. Москва: Эко-трендз, 2002. – 101 с.
10.
Курицын, С.А. Телекоммуникационные технологии и системы
[Текст] / С.А. Курицын. Москва: Academia, 2008. – 155 с.
11.
Берлин, А.В. Сотовые системы связи [Текст] / А.В. Берлин.
Москва: Эко-Трендз, 2011. - 121 с.
66
12.
Френкс, Л. Теория сигналов [Текст] / Френкс Л. М.: Сов. Радио,
1974. – 344 с.
13.
Малла, С. Вейвлеты в обработке сигналов [Текст] / Малла С. М.:
Мир, 2008. – 671 с.
14.
Смирнов, В.И. Курс высшей математики [Текст] / В.И. Смирнов.
М.:Наука, 1974. – 520 с.
15.
Васильева, А.Б. Интегральные уравнения [Текст] / А.Б. Васильева,
Н.А. Тихонов. М.: ФИЗМАТЛИТ, 2002. – 501 с.
16.
Ректорис, К. Вариационные методы в математической физике и
технике [Текст] / Ректорис К. М.: Мир, 1984
17.
Гантмахер, Ф.Р. Теория матриц [Текст] / Ф.Р. Гантмахер М.:
Наука, 1967. – 370 с.
18.
Хорн, Р. Матричный анализ [Текст] / Хорн Р., Джонсон Ч. М.: Мир,
1989. - 528 с.
19.
Жиляков, Е. Г. Оптимальные субполосные методы анализа и
синтеза сигналов конечной длительности [Текст] / Е.Г. Жиляков // Автомат. и
телемех, 2015. – № 4. – С. 51–66.
20.
Хургин, Я.И. Финитные функции в физике и технике [Текст] / Я.И.
Хургин, В.П. Яковлев. М.: Наука, 1971.
21. Rappaport, S. S. Spread-spectrum signal acquisition-Methods and
technology [Текст] / S.S. Rappaport, D.M. Grieco // IEEE Communications
Magazine. – 1984. – Т. 22. – С. 6-21.
22. Peterson, R. L. Introduction to spread-spectrum communications [Текст]
/ R. L. Peterson, R. E. Ziemer, D. E. Borth – New Jersey : Prentice hall, 1995. – Т.
995.
23.
Гмурман, В.Е. Теория вероятностей и математическая статистика
[Текст]: Учебное пособие для вузов / В.Е. Гмурман. Москва: Высшая школа,
2004. – 479 с.
67
24. Ozdemir E. A resampling–based markovian model for automated colon
cancer diagnosis [Текст] / E. Ozdemir, C. Sokmensuer, C. Gunduz-Demir // IEEE
Transactions on Biomedical Engineering. – 2012. – Vol. 59(1). P. 281–289.
25. Siegel, R. Cancer statistics [Текст] / R. Siegel, D. Naishadham, A. Jemal
// CA: A Cancer Journal for Clinicians. – 2013. – Vol. 63. – P. 11–30.
26. Жиляков, Е.Г., Веселых Н.К. Сжатие изображений на основе
субполосного анализа/синтеза [Текст] / Е.Г. Жиляков, Н.К. Веселых //
Научные ведомости Белгородского государственного университета. Серия:
Экономика. Информатика. 2014. №21-1 (192).
27. Жиляков, Е.Г. Вариационные алгоритмы анализа и обработки
изображений на основе частотных представлений: моногр. [Текст] / Е.Г.
Жиляков, А.А. Черноморец. – Белгород: Изд-во ГиК, 2009. – 146 c.
28.
Жиляков, Е.Г. Методы анализа и построения функций по
эмпирическим данным на основе частотных представлений [Текст] / Е.Г.
Жиляков. – Белгород, изд-во БелГУ, 2007. – 160 с.
29.
Жиляков, Е.Г. Метод определения точных значений долей энергии
изображений в заданных частотных интервалах [Текст] / Е.Г. Жиляков, А.А.
Черноморец, И.В. Лысенко // Вопросы радиоэлектроники. Сер. РЛТ. — 2007.
— Вып. 4. — С. 115-123.
30.
Жиляков,
Жиляков, Е.Г. Частотный анализ речевых сигналов [Текст] / Е.Г.
Е.И.
Прохоренко
//
Научные
ведомости
Белгородского
государственного университета. – Белгород, 2006. – №2 (3.118), выпуск 3. – С.
201-208. – (Серия: информатика и прикладная математика).
31.
Черноморец, А.А. Оценка количества значимых долей энергии
изображения, определяемых его проекциями на собственные векторы
субполосных матриц [Текст] / Е.Г. Жиляков, А.А. Черноморец //
Информационные системы и технологии. — № 3 (71). — 2012. — С. 51-58.
32.
Синхронизация по уровню сигнала в триггерах [Электронный
ресурс] / habrahabr.ru – Хабрахабр / URL: https://habrahabr.ru/post/115690/htm
(дата обращения 25.05.2018)
68
33.
Синхронизация в LTE [Электронный ресурс] / anisimoff.org / LTE,
WiMax, WiFi / URL: http://anisimoff.org/lte/lte_synch.html (дата обращения
25.05.2018)
34.
Синхронизация цифровых систем связи по сигналам спутниковых
радионавигационных систем [Электронный ресурс] / radio-systems.org/Радиотехнические
системы
/
URL:
systems.org/synchronization_of_communication_systems
http://radio-
(дата
обращения
25.05.2018)
35.
LTE [Электронный ресурс] / ru.wikipedia.org / Википедия –
свободная энциклопедия / URL: https://ru.wikipedia.org/wiki/LTE (дата
обращения 25.05.2018)
36.
LTE Advanced
[Электронный ресурс] /
ru.wikipedia.org /
Википедия – свободная энциклопедия / URL: https://ru.wikipedia.org/wiki/LTE
(дата обращения 25.05.2018)
37.
Корреляция [Электронный ресурс] / ru.wikipedia.org / Википедия –
свободная энциклопедия / URL: https://ru.wikipedia.org/wiki/ Корреляция (дата
обращения 25.05.2018)
38.
QPSK and OFDM with MATLAB System Objects [Электронный
ресурс] / mathworks.com / MathWork – Makers of Matlab and Simulink / URL:
http://www.mathworks.com/examples/matlab-communications/mw/commex74651154-qpsk-and-ofdm-with-matlab-system-objects
(дата
обращения
25.05.2018)
39.
MathWork
LTE System Toolbox [Электронный ресурс]/ mathworks.com /
–
Makers
of
Matlab
and
Simulink/
URL:
http://www.mathworks.com/products/lte-system/ (дата обращения 25.05.2018)
40.
Downlink
frame
timing
estimate
[Электронный
ресурс]
/
mathworks.com / MathWork – Makers of Matlab and Simulink/ URL:
http://www.mathworks.com/help/lte/ref/ltedlframeoffset.html (дата обращения
25.05.2018)
69
ПРИЛОЖЕНИЕ 1
clc
clear
close all
N=512;
Om=4*pi/N;
sig=2;
sigq=sig^2;
qArray = [0.1];
rSredArrayQ = zeros(1,length(qArray));
rSredArif = zeros(length(qArray),N-1);
figureIndex = 1;
for qIndex= 1:length(qArray)
q = qArray(qIndex);
%Восстановление фрагментов
A0 = zeros(N, N);
Q = zeros(N, N);
Q1 = zeros(N, N);
for i=1:N
for k=1:N
if i==k A0(i,k)=Om/pi;else
A0(i,k)=sin(Om*(i-k))/(pi*(i-k));end
%A0(i,k)=exp(-4*(i*k)/N);
%A0(i,k)=cos(4*pi*(i-k)/N);
end
end
[Q,L]=eig(A0);
for i=1:N
Ni=N+1-i;
Q1(:,i)=Q(:,Ni);
L1(i)=L(Ni,Ni);
end
Ki=1; %номер собственного вектора (собственного числа)
r = zeros(1,N-1);
%
M = 300
for M=1:N-1
B=zeros(M,M);
QQ1=zeros(M,M);
LL1=zeros(1,M);
Dv=zeros(M,N-M);
z2=zeros(1,M);
z1=zeros(1,M);
c=zeros(1,M);
for i=1:M
z1(i)=Q1(i,Ki);%+Q1(i,Ki+2); востанавливаемый фрагмент исходного вектора
for k=1:M %формирование матрицы для восстановления
if i==k
70
B(i,k)=L1(Ki)-A0(i,k);
else
B(i,k)=-A0(i,k);
end
end
% dd(M) = det(B);
z2 = zeros(1, M);
zzz = zeros(1, M);
v = zeros(1, M);
for k=1:M
v(k) = randn;
zzz(k)=Q1(k,Ki);%+Q1(k+M,Ki+2); неиспорченный фрагмент исходного вектора
end
for k=1:N-M
Dv(i,k)=A0(i,k+M); %блок исходной матрицы
end
%
z2(1:sh) = zzz(1:sh) +cc*v(1:sh);
end
vv = norm(v);
zz2 = norm(zzz);
cc = q*zz2./vv; %z2
z2 = zzz;
if M < 100 uu = M; else uu = 100; end
z2(1:uu) = zzz(1:uu) +cc*v(1:uu);
QQ = zeros(M, M);
LL = zeros(M, M);
QQ1 = zeros(M, M);
LL1 = zeros(M, M);
[QQ,LL]=eig(B);
for k=1:M
Mk=M+1-k;
QQ1(:,Mk)=QQ(:,Mk);
LL1(Mk)=LL(Mk,Mk);
end;
%псевдо решение СЛАУ восстановления
zv=Dv*z2';%правая часть системы слау
c1=QQ'*zv;
for k=1:M
c(k)=c1(k)/LL1(k);
end
z1v=inv(B)*zv;%QQ*c'; прямое решение СЛАУ, восстановленная часть вектора
z1v=z1v.*1*10^-6;
% dB(M)=det(B);
% yv=inv(B)*Dv*z2';
dy=norm(z1v-z1');
zz=norm(z1);
r(M)=dy/zz;%относительная погрешность восстановления
% среднеквадр погрешность
sum1=0;
71
sum2=0;
for ind=1:M
sum1=sum1+(z1(ind) - z1v(i))^2;
sum2=sum2+(z1(i))^2;
end
rSredArif(qIndex, M) = sqrt(sum1/sum2);
end
%конц
zx=B*z1';
zy=Dv*z2';
cc=QQ'*zx;
ccv=QQ'*zy;
ccz1=QQ'*z1';
for k=1:M
ccz1v(k)=ccz1(k)*LL1(k);
cccv(k)=(ccv(k))/LL1(k)-ccz1(k);
end
zz1v=QQ*cccv';
zzzv=QQ*ccz1;
[QQ,LL]=eig(B);
%Среднеарифметическая погрешность
rSredQ = sum(rSredArif(qIndex,:))/length(rSredArif(qIndex,:));
%figure(1)
f1 = figure(figureIndex);
%
subplot(2,2,1);
plot(z1)
legend('Исходный сигнал')
grid on
title(['Ki=' num2str(Ki) ', q=' num2str(q)])
%
subplot(2,2,2);
plot(z1v)
legend('Восстановленный сигнал')
grid on
title(['Ki=' num2str(Ki) ', q=' num2str(q)])
%
path = ['Ki' num2str(Ki)];
mkdir(path)
subplot(2,2,[3,4]);
hold on
plot(z1)
plot(z1v)
hold off
grid on
legend('Исходный сигнал','Восстановленный сигнал')
title(['Ki=' num2str(Ki) ', q=' num2str(q)])
saveas(f1, [path '\signals_q' num2str(q) '.png']);
% close(f1)
72
%figure(2)
f2 = figure(figureIndex + 1);
plot(r)
legend('Относительная погрешность восстановления')
title(['Ki=' num2str(Ki) ', q=' num2str(q)])
saveas(f2, [path '\r_q' num2str(q) '.png']);
% close(f1)
% close all
figureIndex = figureIndex + 2;
rSredArrayQ(1,qIndex) = rSredQ;
end
% %figure(3)
% f3 = figure(figureIndex+1);
% plot(qArray, rSredArrayQ)
% legend('Значение отношения средней среднеквадратической ошибки от шум/сигнал')
% xlabel('Шум/Сигнал')
% ylabel('Среднее среднеквадратической погрешности')
% title(['Ki=' num2str(Ki)])
% saveas(f3, [path '\rArifSred_q' num2str(q) '.png']);
% % close(f1)
% close all
Отзывы:
Авторизуйтесь, чтобы оставить отзыв