На правах рукописи
АМРО Аладдеин Мухаммед
РАЗРЕШЕНИЕ КЛАССОВ СИГНАЛОВ МЕТОДОМ
МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ В УСЛОВИЯХ
ВОЗДЕЙСТВИЯ ПОМЕХ
Специальность 05.12.04 - радиотехника, в том числе системы и
устройства радионавигации, радиолокации и телевидения.
Автореферат диссертации на соискание ученой степени кандидата физикоматематических наук
КАЗАНЬ 2003
Работа выполнена в Казанском государственном университете
Научный руководитель: доктор физико-математических наук, профессор Э. А. Ибатуллин
Официальные оппоненты: доктор технических наук, профессор В. А. Белавин кандидат технических наук, доцент Ю. Л. Комаров
Ведущая организация: Казанский научно-исследовательский институт радиоэлектроники Зашита состоится 01.07.2003 в 14 часов на заседании диссертационного совета Д212.081.18 при Казанском государственном университете по адресу:
Кремлевская, 18, физический факультет, ауд. 210.
С диссертацией можно ознакомиться в библиотеке Казанского государственного университета.
Автореферат разослан «30» мая 2003 г.
Ученный секретарь диссертационного совета, доктор физико-математических наук Карпов А. В.
{5^—
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Увеличение количества радиоэлектронных систем (РЭС), используемых в быту и в производстве, а также переполненность радиочастотного диапазона привели к осложнению электромагнитной обстановки, что вызвало в свою очередь ухудшение разрешаемости классов сигналов, — с одной стороны, а с другой стороны желание получить высокое быстродействие мешало использованию сложных алгоритмов разрешения классов сигналов. Но в настоящее время после существенного прогресса средств программирования и электроники появилась возможность использовать более сложные алгоритмы разрешения для повышения точности при допустимом быстродействии.
Определим вначале понятие класса сигналов. Под классом сигналов будем понимать совокупность сигналов от одного источника, когда существенные параметры сигнала меняются от сигнала к сигналу случайным образом.
Случайный характер изменения параметров сигнала может происходить вследствие аппаратурных нестабильностей и ошибок измерения. Например, в качестве случайных существенных параметров сигналов могут выступать время запаздывания сигнала или его несущая частота.
Класс сигналов представляет собой совокупность дискретных сигналов, в качестве которых могут использоваться простые радиоимпульсы, сигналы с внутриимпульсной модуляцией (ЛЧМ, ФМн, сигналы на базе составных последовательностей, сигналы на базе кодов Баркера и т. д.). Эти сигналы обычно используются в радиолокационных системах, в системах управления воздушным движением в аэропортах, в системах траекторных измерений, в системах активного запроса и ответа, в радиотехнических системах предупреждения столкновений самолетов.
В настоящее время актуальной является проблема оптимального разрешения классов сигналов, которое заключается в нахождении оценок параметров частных распределений, характеризующих классы, и в распределении сигналов по их классам. Процедуре разрешения классов сигналов предшествует разрешение сигналов. Говоря иначе, до процедуры разрешения классов сигналы должны быть предварительно обработаны, т е.
усилены, сжаты, продетектированы, декодированы и иметь метку (label) со значением того параметра, по которому будет производиться разрешение классов. Например, если в качестве такого параметра выступает время запаздывания, то должно быть значение временного положения сигнала относительно некоторого опорного момента времени. Если же в качестве параметра берется несущая частота, то с помощью частотно-избирательного устройства определяется значение несущей частоты или значение отклонения от номинала несущей частоты сигнала. Эти значения представляются в цифровом виде, так как рассматриваемый алгоритм разрешения классов сигналов, предполагает цифровую обработку.
При разработке алгоритма разрешения классов сигналов с учетом выше сказанного невозможно игнорировать проблему электромагнитной совместимости (ЭМС). Это крайне важно в данном исследовании, так как алгоритм имеет дело со случайными сигналами и должен корректно функционировать в условиях воздействия непреднамеренных помех.
В настоящее время в качестве критериев ЭМС часто используют энергетический критерий, при котором электромагнитное воздействие мешающего передатчика на исследуемый приемник определяется сравнением величины мощности помехи на входе (выходе) приемника с мощностью, при которой срабатывает исполнительное устройство на выходе приемника с вероятностью близкой к единице. Если величина мощности помехи равна или превышает определенную таким образом пороговую мощность, то считается, что условия ЭМС этих информационных систем (ИС) нарушены, и необходимо принимать меры по уменьшению мощности помех.
По нашему мнению, широко используемый энергетический критерий электромагнитной совместимости наряду с простотой не дает возможности провести количественную оценку ЙА^ОВб^нЪсиЙл^инймаелюй информации в условиях воздействия помех для большого класса информационных систем.
Действительно, довольно часто, даже в случае превышения мощности помехи порогового значения, она может не оказывать существенного влияния на работу ИС в целом. Примером может служить воздействие импульсной помехи большой скважности на ИС, использующую импульсный сигнал.
Поэтому в данной диссертации мы предлагаем определять пороговую мощность с учетом влияния помехи на эффективность РЭС в целом. Это позволяет избежать те случаи, когда критерии ЭМС бывают неоправданно завышены, что является актуальным в настоящее время.
Цель и задачи исследования. Цель настоящей работы заключается в разработке алгоритма разрешения классов сигналов на основе метода максимального правдоподобия, в синтезе помеховой обстановки, создаваемой непреднамеренными помехами, а также в анализе генератора класса сигналов.
Для достижения поставленной цели необходимо решить следующие задачи.
1. Определить понятие класса сигналов и дать его статистическое описание;
2. Определить необходимые и достаточные условия разрешения классов сигналов;
3. Синтезировать оптимальный алгоритм разрешения классов сигналов с помощью метода максимального правдоподобия;
4. Получить и решить систему уравнений максимального правдоподобия для определения средних значений в практически важном случае, когда частные плотности распределений моментов прихода сигналов подчиняются нормальному закону;
5. Синтезировать допустимую помеховую обстановку с учетом критерия эффективности информационной системы;
6. Разработать методику по определению пороговой мощности шумовой помехи для импульсной радиоэлектронной системы с учетом критерия эффективности этой системы;
7. Провести анализ генератора класса сигналов.
Научная новизна работы 1. Сформулирована и доказана теорема, определяющая необходимое и достаточное условия разрешения классов сигналов.
2. Получена система уравнений максимального правдоподобия для определения средних значений в практически важном случае, когда частные плотности распределений моментов прихода сигналов подчиняются нормальному закону. Система уравнений решена методом Ньютона.
Получены численные результаты, которые показывают работоспособность алгоритма.
3. Разработана новая методика синтеза допустимой помеховой обстановки с учетом критерия эффективности информационной системы.
4. Для конкретной РЭС разработана методика по определению пороговой мощности для шумовой помехи.
5. Проведен анализ генератора класса импульсных сигналов.
Публикации. По теме диссертации опубликовано десять печатных работ, в том числе 8 статей и 2 тезисов докладов.
Практическая ценность результатов работы 1. Оптимальный алгоритм разрешения классов сигналов, синтезированный с помощью метода максимального правдоподобия, позволяет создавать самонастраивающиеся устройства разрешения повышенной точности для информационных систем, которые к настоящему времени практически отсутствуют.
информационной системы позволяет качественно на новом уровне определять условия ЭМС информационных систем.
3. Разработка генераторов классов импульсных сигналов позволит создать соответствующие имитаторы для исследований и испытаний импульсных информационных систем.
Личный вклад. Автор принимал участие в синтезе алгоритма разрешения классов сигналов и в синтезе методики расчета пороговой мощности шумовой помехи на выходе радиоприемного устройства импульсной РЭС. Непосредственно им написана программа разрешения классов сигналов на алгоритмическом языке C++, получены численные результаты и исследованы возможности алгоритма разрешения классов сигналов, Автор также принимал участие в исследовании генератора класса сигналов.
Апробация работы. Основные результаты и положения диссертационной работы доложены и обсуждены: на международных научных сессиях, посвященных Дню Радио, г. Москва, 1999 + 2003 г.г.; на шестой российской научно-технической конференции «ЭМС-2000» г. Санкт-Петербург, на итоговых научных конференциях Казанского государственного университета 2000 + 2002 г.г.
Структура и объем диссертация Диссертация состоит из введения, четырех глав, заключения и приложений. Общий объем - 139 страниц, включая библиографию из 82 наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введения содержится обоснование актуальности разрабатываемой темы, изложена цель работы и задачи исследований, показана научная новизна, практическая ценность, апробация представленных к защите задач.
Первая глава диссертационной работы посвящена общему подходу к решению задач распознавания сигналов в условиях воздействия помех Определены основные понятия проблемы разрешения классов сигналов:
класс, сигнал, смесь сигналов, описана вероятностная модель сигналов, среди различных алгоритмов распознавания (разрешения) выделены следующие: метод моментов, метод минимума хи-квадрат, метод максимального правдоподобия.
Представлено общее состояние вопроса электромагнитной совместимости в настоящее время и имеющиеся методы анализа электромагнитной обстановки.
разрешения классов сигналов, сформулирована и доказана теорема, определяющая необходимое и достаточное условия разрешения классов сигналов. Эти условия сводятся к тому, чтобы семейство функций распределения, из которых состоит смесь, было линейно независимым множеством над полем действительных чисел.
При этом для статистического описания классов сигналов использовано понятие смеси распределений где k - количество классов стохастических сигналов; f(X/&t}- частная характеризуемого совокупностью параметров ©,;/?, -вес г'-го класса в смеси.
Для разрешения классов сигналов или для нахождения оценок pt и 0,, в выражении (1), когда известны виды плотностей вероятностей исходных распределений, а также их число в смеси был использован метод максимального правдоподобия.
При этом для логарифма функции правдоподобия пачки из п независимых сигналов будем иметь Была получена система уравнений максимального правдоподобия для практически важного случая, когда частные плотности распределений моментов прихода сигналов подчиняются нормальному закону. При этом распределения, а -дисперсия, одинаковая для всех частных распределений Предполагалось, что pi и известны, k=3. Требуется определить оценки для математических ожиданий at,i-\,k, которые пропорциональны дальностям до источников излучений. При этих условиях были получены следующие уравнения правдоподобия:
Таким образом, имеем систему нелинейных уравнений, для решения которой могут быть использованы только численные методы Одним из таких методов является итерационный метод Ньютона, используя который могут быть найдены оценки для неизвестных параметров.При этом для нашего случая была получена система уравнений где t - номер шага итераций, /,/e;(3M) = ^-/I(5/-1),« = U,y = U, В свою очередь, Aa'j^a'j-a™.
Теперь, решая совместно уравнения (4), находим приращения Да', Далее по найденным значениям приращений могут быть определены оценки для математических ожиданий на шаге t в соответствии с выражением (5) Итерационный процесс продолжается до тех пор, пока не будут одновременно выполнены неравенства для всехj = l,k,, где к- заданная точность вычислений Полученные оценки математических ожиданий и будут оценками максимального правдоподобия, вычисленными с точностью s Для проверки алгоритма разрешения классов сигналов был проведен вычислительный эксперимент При этом генерировались выборки с различным расположением математического ожидания классов сигналов a,,i = l,k относительно друг друга Сначала вычисления проводились, когда классы находились на достаточно большом расстоянии друг от друга, потом постепенно сближали их до пересечения. Результаты получены при распределений рг = 0,3, р2 = 0,3, р3 = 0,4; дисперсия а2 = 0,09; одинаковая для всех классов, размер выборки 15-3000; начальные значения в пределах (1-6) от истинных значений.
На рисунках 1, 2 представлен один из случаев, для которого проводились вычисления. На рисунке 1 показана гистограмма для трех классов сигналов в случае их пересечения, на рисунке 2 представлены результаты итерационного процесса, где истинные значения отображены пунктирными линями. Из рисунка видно, что приближение получаемых оценок к разыгрываемым значениям средних при перечисленных исходных данных происходит в районе 5-й итерации.
Эксперимент показывает, что правильность разрешения и номер итерации, на котором происходит приближение получаемых оценок к разыгрываемым значениям средних, зависит от нескольких факторов, а именно, размера выборки, степени пересечения классов, начальных значений и насколько они близки к истинным значениям. Важнейшим из этих факторов является последний, влияние которого становится более значительным в случае пересечения классов.
В случае, когда классы не пересекаются, разрешение происходит даже математического ожидания (на и более), а в случае пересечения классов разрешение происходит только при начальных значениях близких к математическому ожиданию (в пределах одной ).
Третья глава посвящена рассмотрению вопросов помехоустойчивости и ЭМС информационной системы, где может быть реализован алгоритм разрешения классов сигналов, описанный во второй главе. При этом использован новый подход для синтеза допустимой помеховой обстановки с учетом критерия эффективности информационной системы, потому что мы считаем, что широко используемый для определения возможностей работы различных систем в условиях воздействия непреднамеренных помех энергетический критерий наряду с простотой не дает возможности провести количественную оценку достоверности принимаемой информации в условиях воздействия помех для большого класса ИС. В предложенном подходе устанавливаются зависимости между показателями эффективности ИС и мощностью помехи известного вида, затем, задавшись допустимым снижением эффективности при воздействии помех, определяют пороговую мощность для данной помехи.
Таким образом, если установлена зависимость вида где Рт — мощность помехи на выходе приемного устройства, А,„ — средняя частота следования для хаотически-импульсной помехи (ХИЛ) на входе приемного устройства, рх — вероятность прохождения ХИЛ через приемное устройство, Рс,Ре — соответственно мощности сигнала и естественных шумов на выходе приемного устройства, Qn,Q0,Qp -соответственно операторы приема, обработки информации и принятия решения, то для заданного значения Эт из (9) находим Рл™, которое и будет являться пороговым значением мощности для данного вида помехи.
В случае использования векторного критерия эффективности для оценки ИС, когда z показателей качества зависят от воздействия помехи, имеем z выражений подобных (9) В этой ситуации из каждого выражения определяется значение мощности помехи и в качестве пороговой мощности выбирается ее минимальное значение При воздействии помех от нескольких источников случае скалярного эффективности при воздействии помех не хуже заданной, определяем требования к мощности помехи на выходе приемного устройства для каждого вида помехи, действующего на входе ИС Таким образом, для каждого источника помех будем иметь свое значение пороговой мощности помехи Э(2)>Э"Ч2)=>О2), При использовании векторного критерия для каждого из z показателей качества определяем требования к величине мощности помехи на выходе приемного устройства для каждого источника помехи Если число источников помех равно то в общем случае получаем zl значений мощностей помех, из которых необходимо выбрать минимальные значения для каждого источника помехи В результате будем иметь / пороговых мощностей Данный подход был использован при разработке методики синтеза электромагнитной обстановки для импульсной радиоэлектронной системы, расположенной на борту, и предназначенной для приема радиосигналов, передаваемых бинарным кодом. Исследуемый тракт состоит из антенны, циркулятора, приемника, порогового устройства (двоичный квантователь бинарный код с пассивным нулем, у которого синхрогруппа состоит из п позиций, а информационная группа — из т позиций, на которых могут подключается устройство принятия решения, которое необходимо для накопления дешифрированных бинарных сигналов с целью увеличения достоверности принимаемой информации. Для выполнения такой функции устройство принятия решения содержит обнаружитель начала пачки и накопитель дешифрированных бинарных сигналов В роли показателей качества векторного критерия, характеризующего эффективность такой РЭС, были приняты вероятности правильного и ложного решений, которые в соответствии с алгоритмом, заложенным в устройстве принятия решения, определяются следующими выражениями.
где Pm — полная вероятность искажения сигнала, Рпт — полная вероятность обнаружителя начала пачки, а, п — соответственно нижний и верхний пороги накопителя.
порогового устройства приемника от импульсных помех и выбросов шума подчиняется закону Пуассона будем иметь :
где N - среднее количество срабатываний порогового устройства приемника от импульсных помех и выбросов шума в единицу времени, icmf — интервал стробирования дешифратора, ^ — число нулевых позиций в синхрогруппе, - число единичных позиций в информационной группе, и/ -число позиций в синхрогруппе бинарного сигнала.
Выражения (14) и (15) показывают, что воздействие помех на вероятности Pnvvi /"„,„сказывается через среднее количество срабатываний порогового устройства N.
Дальнейшие преобразования выражения (14) с использованием бинома Ньютона позволили выразить в явной форме зависимость N от Рт :
Преобразования выражения (15) привели его к виду: