
- •1. Описание сигнала и помехи
- •2.Проверка статистических гипотез
- •2.1.Определения
- •2.2. Критерий Байеса и правило принятия решений
- •2.3. Критерий максимума апостериорной вероятности и правило принятия решений
- •2.4. Критерий и правило максимума отношения правдоподобия
- •2.5. Критерий Неймана-Пирсона и правило принятия решений
- •2.6. Минимаксный критерий
- •3. Расчет вероятностей ошибок
- •4. Функционал отношения правдоподобия
- •5. Вероятность ошибки и вероятность правильного принятия решения
- •6. Методика эксперимента
5. Вероятность ошибки и вероятность правильного принятия решения
Определим вероятность ошибки
и вероятность правильного обнаружения
как функция порога
.
Как известно, ошибка первого рода
возникает тогда, когда принимается
решение о наличии сигнала
в принятой реализации в то время каксигнал отсутствует.
Вероятность этого события равна
. (5.1)
Для упрощения
расчетов произведем нормирование
случайной величины V , сдвинув V на
величину условного математического
ожидания при отсутствии сигнала
и пронормировав полученное выражение
относительно условного среднеквадратического
отклонения s
v0 при
отсутствии сигнала:
(5.2)
где
,
.
(5.3)
Вероятность правильного обнаружения
как функция порога С*есть
вероятность превышения случайной
величиной V порога С*, но
при наличии сигнала в принятой реализации:
(5.4)
где
,
.(5.5)
В интегральной форме вероятность ошибки aи вероятность правильного решения D = 1-bимеют вид
, ( 5.6)
, (5.7)
где
и
вычисляются при помощи формул (4.11),
(5.3), (5.5).
Формулы (5.6) и (5.7) позволяют построить
рабочую характеристику приемника при
использовании непрерывных наблюдений
на интервале (0,
).
6. Методика эксперимента
1. Для того, чтобы снять рабочую
характеристику, необходимо измерить
вероятность правильного обнаружения
D=1-bи вероятность
ошибки первого родаa.
Эти вероятности измеряются методом
статистических испытаний. Для измерения
D=1-bв течение времени
периодически воспроизводится передача
“1” по каналу связи при заданном
отношении сигнал/шум и фиксированном
значении порога. (Переключатель П1в положении “Сигн.+Шум”, а переключатель
П2в положениях “Е1“, “Е2“, “Е3“).
В качестве оценки вероятности правильного
обнаружения D=1-bберется отношение числа импульсов на
выходе приемника Niк общему числу
сигналов N*, поступивших на вход
приемника за время измерения
,
т.е. D*= Ni/ N*.
Для измерения aв
течение времени
периодически воспроизводится передача
“0” по каналу связи с определенным
уровнем шума. (Переключатель П1в
положении “Сигн.+Шум”, а переключатель
П2в положении “0”). Вероятность
ложной тревоги вычисляется как отношение
числа ошибочных решений Niк общему
числу испытаний за время
:a*= Ni/ N*.
При N*®¥оценкиa*и D*стремятся к своим истинным значениямaи D. Следовательно,
результаты измерений тем точнее, чем
больше время измерения. Отношение Ni/N*измеряется частотомером. Для этого
импульсы с выхода приемника Niподаются на вход А, а с выхода N*-
на вход Б частотомера. Переключатель
частотомера “Род работы” устанавливается
в положение “Отношение частот A/B”.
Время измерения
устанавливается максимальным при помощи
переключателя “Время счета”. При этом
на цифровом индикаторе частотомера
будет высвечиваться значениеa*или D*в зависимости от положения
переключателя П1.
Рекомендуется регистрировать значения aи D при одном и том же значении порога С*. Для этого при фиксированном значении порога С*переключатель П1поставить в положение “Сигн.+Шум” и манипулировать переключателем П2.
Для того, чтобы вычислить значение порога триггера Шмидта приведем правило решения (4.16) к виду
Выходное напряжение интегратора pавно
Здесь К(CT) множитель, учитывающий коэффициент передачи перемножителя и интегратора. Поэтому величина порога СТтриггера Шмидта связана с порогом С выражением
(6.1)
Значение К(C T) – указано на установке.
2. В пункте 8 ЗАДАНИЯ необходимо построить правила d1,d2проверки гипотез о состоянии источника. Согласно правилуdjинтервал значений порогов (С Т min , C T max ) разбивается точкой СT j на область (С Т min ,СT j ) - приема гипотезы Н0 , и область (СT j , C T max ) - приема гипотезы Н1 . Примем СT1 = С TБ , где С TБ - порог Байеса, найденный в пункте 7, а порог C Т2 определим по формуле
C Т2 = С Т min + k0 ( C T max - С Т min ),
где 0 < k0 < 1 и задается преподавателем.
3. При расчетах вероятностей a, D=1-bи порога С можно воспользоваться как экспериментальными данными , так и построенными графиками. Погрешность представленияa, D и С зависит от качества интерполяции и качества построения графиков. Вероятностиaи D по графику определяются следующим образом:
- для известного С в произвольном месте
чертежа под углом qк оси абсцисс проводится прямая,
удовлетворяющая уравнению
,
- построенная прямая переносится параллельно самой себе до касания с выбранной рабочей характеристикой,
- в точке касания определяются вероятности aи D.
З А Д А Н И Е
1. Ознакомиться с описанием экспериментальной установки, методикой эксперимента, получить от преподавателя значения априорной вероятности P0состояния источника, матрицу потерь П, отношение сигнал/шум и коэффициент k0 .
2. Поставить переключатель П1в положение “Сигн.”. Какому состоянию канала соответствует это положение ?
Поставить переключатель П2в положение “Е1“.
Пронаблюдать и зарисовать осциллограммы напряжений в контрольных точках КТ1, КТ2, КТ3. Объяснить изменение формы сигнала в контрольной точке КТ3при изменении порога.
Измерить амплитуду и длительность
радиосигнала в контрольной точке КТ1во всех положениях переключателя П2.
По результатам измерений вычислить
энергию сигнала
.
3. Поставить переключатель П1в положение “Сигн.+Шум”.
Пронаблюдать форму напряжений в контрольных точках КТ1, КТ2, КТ3 при всех положениях переключателя П2. Что изменяется и почему ?
Для каждого положения переключателя П2вычислить отношение сигнал/шум. В положении переключателя П2“0” пронаблюдать форму напряжений в контрольных точках. Объяснить изменение формы сигналов в точке КТ3при изменении порога.
4. Снять зависимость вероятности правильного обнаружения D и вероятности ошибки aот величины порога для всех значений энергии сигнала, вычисленных в пункте 2.
5. По полученным данным построить рабочие характеристики приемника на плоскости (a, D) для всех значений энергии сигнала, вычисленных в пункте 2, пересчитав предварительно пороги СТтриггера Шмидта в пороги С, пользуясь соотношением (6.1).
6. Вычислить вероятности aи D по формулам (5.6) и (5.7) для заданного преподавателем отношения сигнал/шум .
Построить теоретическую рабочую характеристику на том же графике, что и экспериментальная рабочая характеристика и сравнить теоретическую кривую с экспериментальной.
7. Критерий Байеса. Для заданной преподавателем матрицы потерь П, априорной вероятности состояния источника P0и отношения сигнал/шум построить график среднего риска как функцию порога С по формуле (2.7), приведя его к виду
Определить по графику порог СБ,
при котором функция риска достигает
минимума, и соответствующие ему
вероятностиaи D.
Определить порог триггера Шмидта С
TБ, соответствующий порогу СБ.
Для найденных aи D проверить формулу (3.16).
8. Минимаксный критерий. Задать два правила d1 ,d2 проверки гипотез о состоянии источника. Выбрать наилучшее правило.
9. Критерий максимума апостериорной вероятности. По заданной априорной вероятности P0состояния источника найти порог СМАВи вероятности (a(di) , D(di)) по экспериментальным кривым как функцию отношения сигнал/шум. Построить график D(di) и объяснить его поведение.
10. Критерий максимума правдоподобия. Найти порог СМП. Для каждого отношения сигнал/шум по экспериментальным кривым найти вероятности (a(d i), D(d i )). Построить график D(di) и объяснить его поведение.
11. Критерий Неймана-Пирсона. Положим вероятность ошибки первого рода равна вероятности a, найденной в п. 8. Определить порог СН-Пи вероятность правильного обнаружения D по экспериментальной кривой и теоретически для отношения сигнал/шум, заданного в п.1. Сравнить полученные результаты.
12. Построить таблицу для заданного в п.1. отношения сигнал/шум
Критерий |
Априорные данные |
Порог |
a, D |
Байеса |
|
|
|
minmax |
|
|
|
МАВ |
|
|
|
МП |
|
|
|
Неймана-Пирсона |
|
|
|
Сравнить полученные результаты. При каких условиях критерий Байеса можно заменить критериями minmax, МАВ, МП.