Радиотехнические системы. Лабораторный практикум. Учебно-методическое пособие
.pdf
где rik – безразмерные коэффициенты, определяющие значимость (це-
ну, риск) принятия решения в пользу k-го сигнала, если передавался i-й сигнал.
Если известно вероятностное описание помех, действующих в канале передачи информации, можно записать условные ПРВ для соответ-
ствующих событий: w U Ai , i 1... m , где U U (t), 0 t c – наблюдаемый на входе приемника сигнал, который может принадлежать одной из возможных областей линейного пространства k , k 1... m .
Таким образом, легко найти вероятность перепутывания сигналов:
pik p Ak* | Ai p U k | Ai w U Ai dU .
k
Исходя из этого, байесовский риск запишется в виде
m m |
|
|
|
|
|
|
r r ik pi w U |
|
Ai dU |
||||
|
||||||
i 1k 1 |
k |
|||||
m |
m |
|
|
|
|
|
|
r ik piw U |
|
Ai dU min. |
|||
|
||||||
k 1 k |
i 1 |
|
|
|
|
k |
m
Сумма r ik piw U Ai dU показывает условный риск принять ре- i 1
шение Ak* при наблюдении на входе приемника сигнала U (t) . Решение
принимается в пользу той гипотезы, для которой этот риск минимален. Таким образом, оптимальный приемник различения можно построить в виде многоканальной системы. Количество каналов соответствует количеству используемых для передачи информации сигналов. В каждом канале происходит вычисление условного риска, затем производится их сравнение, и решение принимается в пользу гипотезы с ми-
нимальным условным риском.
Частным случаем является различение равновероятных сигналов
S1(t) и S2 (t) с одинаковыми энергиями ЭS1 ЭS2 Э , наблюдаемых на фоне белого шума. При матрице рисков, имеющей вид
11
0, |
i k, |
|
|
r ik |
i k, |
1, |
|
|
|
принятие решения заключается в выборе сигнала, для которого значение корреляционного интеграла будет максимальным [2]:
A1*
z1(U ) z2 (U ) ,
где
A0*
c |
|
zi (U ) |
U (t)Si (t) dt , i 1, 2 . |
0 |
|
При этом решающее правило можно переписать в виде
A1*
z1(U ) z2 (U ) 0 .
A0*
Для анализа ошибок различения необходимо рассмотреть ПРВ значений сигнала (рис. 5), полученных после обработки. Очевидно [2], что при гауссовых ПРВ вероятность перепутывания сигналов при опти-
мальном пороговом уровне z0 0 будет определяться по выражению
|
0 |
|
|
pош p12 p21 w A2 |
d w A1 |
d . |
|
|
|
0 |
|
|
w( ) |
|
|
w A2 |
|
w A1 |
|
2 |
0 |
1 |
|
Рис. 5. Пример формирования условных ПРВ при различении сигналов с одинаковыми энергиями на фоне белого шума
Найдя ПРВ w A2 , получим, что при различении двух сигналов
в когерентной системе вероятность ошибки можно оценить по форму-
ле [2]
12
pош 1 Ф q |
0,5 (1 ) , |
T
где Э 1 S1(t)S2 (t) dt – коэффициент корреляции сигналов, исполь-
зуемых для0 передачи «0» и «1». Для некогерентной системы точные формулы для вероятностей ошибок различения достаточно громоздки [2], поэтому в лабораторной работе предлагается в качестве теоретической границы оценить вероятность ошибки для ортогональных в усиленном смысле сигналов, которая может быть вычислена по формуле [2]
q2
pош 0,5e 4 .
Типичным примером применения приемника различения сигналов является цифровая система передачи информации с частотной манипуляцией (СПИ с ЧМн). Сигнал с частотной манипуляцией (рис. 6) можно представить как сумму двух сигналов с амплитудной манипу-
ляцией с частотами |
f0 и f1 : |
|
|
|
|
|
|
|||
|
a cos(2 f |
0 |
t ), |
0 t |
с |
, если передается «0»; |
||||
S(t) |
|
|
|
|
|
|
|
|
||
|
|
|
|
0 t , еслипередается «1». |
||||||
|
a cos(2 f t ), |
|||||||||
|
|
1 |
|
|
с |
|
|
|
||
S(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
1 |
0 |
|
|
1 |
1 |
0 |
|
1 |
0 |
1 |
|
Рис. 6. Сигнал с частотной манипуляцией |
|
||||||||
Коэффициент корреляции для прямоугольных радиоимпульсов, из которых состоит указанный сигнал, определяется частотами f0 и f1 ,
атакже длительностью тактовой посылки с :
sin 2 ( f0 f1) с . 2 ( f0 f1) с
13
Минимальное значение коэффициента корреляции 0,212 (рис. 7) в системе с частотной манипуляцией достигается при
2 f с 1,5 .
2f с
2 |
|
|
2 |
–0,212
Рис. 7. Коэффициент корреляции двух прямоугольных радиоимпульсов
Структура СПИ с ЧМн приведена на рис. 8 и состоит из следующих частей:
|
|
|
|
|
s2(t) |
|
s3(t) |
|
s5(t) |
|
|
|
|
|
|||
|
u(t) |
|
|
|
СФ1 |
АД1 |
s7(t) |
y(t) |
|||||||||
|
|
s1(t) |
|
|
|
|
|
|
|||||||||
ИС |
|
ЧМ |
|
К |
|
|
|
s4(t) |
|
|
s6(t) |
– |
|
|
ПУ |
|
ПС |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
СФ2 |
|
АД2 |
|
|
|
|
|
|
||
|
|
Г |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 8. Структурная схема СПИ с частотной манипуляцией
ИС – источник случайного бинарного сообщения, описывающийся вероятностью появления нуля или единицы;
Г – генератор высокочастотного гармонического сигнала; ЧМ – модулятор, осуществляющий изменение частоты гармониче-
ского сигнала в соответствии с модулирующим сообщением; К – канал передачи информации, в котором сигнал подвергается
воздействию помехи; СФ – согласованный фильтр, выполняющий вычисление корреля-
ционного интеграла. Так как частотно-манипулированный сигнал можно представить суммой двух амплитудно-манипулированных сиг-
14
налов с частотами 0 и 1 соответственно, то в приемнике имеется
два канала, каждый из которых настроен на свою частоту; АД – амплитудный детектор;
ПУ – пороговое устройство, включающее в себя блок стробирования и синхронизации, необходимых для того, чтобы сигнал на пороговое устройство проходил только в момент окончания тактовой посылки; ПС – получатель сообщения, который является оконечным устрой-
ством, осуществляющим обработку информации.
Описание работы в моделирующей программе
Моделирующая программа реализована в пакете «Mathematica», который является комплексом для проведения инженерных и научных вычислений. Программирование в «Mathematica» осуществляется на Си-подобном языке. Рабочая область программы поделена на ячейки, каждая из которых ограничена квадратной скобкой справа. Запуск ячейки на выполнение производится нажатием клавиш Shift+Enter. Основными объектами являются: текст программы; результат вычислений; графики, рисунки, диаграммы и пр.
Система передачи информации моделируется с помощью набора подпрограмм (см. таблицу), реализующих определенные преобразования сигнала.
Подпрограммы, реализующие функции основных блоков системы передачи информации
Синтаксис |
Описание и параметры выходного сигнала |
|
datasource[ ] |
Генератор бинарного равновероятного сообщения: |
|
|
длина сообщения N 60 бит, |
|
|
длительность информационной посылки с 1 мс |
|
modulatorAM[u, f] |
Генератор амплитудно-манипулированного сигна- |
|
|
ла: |
сигнала |
|
u – модулирующее сообщение, амплитуда |
|
|
a 1 В, f – несущая частота (1…40 кГц) |
|
modulatorFM[u, f0, f1] |
Генератор частотно-манипулированного сигнала: |
|
|
u – модулирующее сообщение, амплитуда |
сигнала |
|
a 1 В, f0, f1 – несущие частоты (1…40 кГц) |
|
noise[s1, q] |
Канал передачи информации: |
|
|
s1 – входной сигнал, |
|
|
q – отношение сигнал/шум |
|
15
|
Окончание табл. |
|
|
Синтаксис |
Описание и параметры выходного сигнала |
optimfilter[s2, f] |
Согласованный фильтр: |
|
s2 – входной сигнал, |
|
f – несущая частота |
detector[s3] |
Амплитудный детектор: s3 – входной сигнал |
threshold[s4, z0] |
Пороговое устройство: s4 – входной сигнал, |
|
z0 – порог |
signalview[s, t1, t2] |
Оператор просмотра сигнала s на промежутке времени |
|
от момента t1 до момента t2 |
dataview[u, n1, n2] |
Оператор просмотра бинарного сообщения u от бита |
|
n1 до бита n2 |
В пакете «Mathematica» для предотвращения вывода на экран результата вычисления строка должна заканчиваться точкой с запятой.
Для оценки вероятности ошибки необходимо осуществить статистический эксперимент, состоящий из K наблюдений. Удобным инструментом для этого является программный цикл, синтаксис которого приведен ниже.
N1=0;
For[j=1,j<=K,j++,
Тело цикла;
U=datasource[];
…
V=threshold[s4,z0]; N1=N1+Sum[Abs[V[[i]]-U[[i]]],{i,60}];]
OSH=N[N1/(60*K)]
Домашнее задание
1.Определить необходимый объем эксперимента для измерения вероятностиошибки впределахот0,1 до0,5 спогрешностью не более 5 %.
2.Определить оптимальную величину порога для СПИ с АМн при равновероятном источнике.
3.Рассчитать теоретическую зависимость вероятности ошибки от отношения сигнал/шум для СПИ с АМн и построить график.
4.Для СПИ с ЧМн, задаваясь частотой 0 согласно своему вари-
анту, определить частоту 1 из условия 0 и min.
16
5. Рассчитать теоретическую зависимость вероятности ошибки от отношения сигнал/шум для СПИ с ЧМн для 0 и min и постро-
ить график.
Лабораторное задание
1. Исследовать когерентную СПИ с АМн.
1.1. Пронаблюдать и зарисовать сигнал во всех точках системы передачи информации при отношении сигнал/шум q2 100 .
1.2.Измерить зависимость вероятности ошибки от величины отношения сигнал-шум при оптимальном пороге, пороге больше и меньше оптимального.
1.3.Исследовать величину ошибки при оптимальном пороге и частотной расстройке согласованного фильтра.
2. Исследоватьнекогерентную СПИ сАМн по алгоритму пп. 1.1–1.3. 3. Сравнить полученные характеристики для когерентной и некоге-
рентной систем и сделать выводы о причинах их различия. 4. Исследовать СПИ с ЧМн.
4.1.Пронаблюдать и зарисовать сигнал во всех точках системы пе-
редачи информации при отношении сигнал/шум q2 100 .
4.2. Исследовать зависимость вероятности ошибки от величины отношения сигнал/шум при частотах, соответствующих минимальному коэффициенту корреляции и соответствующих коэффициенту корре-
ляции 0 .
Контрольные вопросы
1.Какую функцию выполняет каждый блок структурной схемы системы передачи информации с амплитудной манипуляцией?
2.Для чего нужна модуляция в радиотехнических системах? В каких системах можно не использовать модулятор в передатчике? Каким образом сигнал с частотной манипуляцией можно представить в виде суммы двух амплитудно-манипулированных сигналов с разными частотами?
3.Для чего нужен согласованный фильтр (СФ) в радиотехнических системах? Какие параметры сигнала необходимо знать, чтобы построить СФ? Какое условие должно выполняться для обеспечения каузаль-
17
ности СФ? Каков механизм увеличения отношения сигнал/шум на выходе согласованного фильтра?
4.Чем отличаются когерентный и некогерентный приемы? В каком случае нужен амплитудный детектор? В чем состоит необходимость наличия амплитудного детектора в системе с частотной манипуляцией?
5.Чем отличаются различные критерии выбора оптимального порога для порогового устройства?
6.От чего зависит вероятность ошибки в системе передачи информации? Как можно ее уменьшить?
7.Характеристикой какого блока структурной схемы является априорная вероятность?
8.В чем состоит основное различие обнаружения и различения сигналов?
9.Для чего нужны два канала обработки сигнала в приемнике системы передачи информации с частотной манипуляцией? Каким образом происходит разделение этих каналов? Для чего нужно вычитающее устройство?
10.Каким образом формируются ПРВ сигналов для системы с частотной манипуляцией? Каким образом выбирается оптимальный порог?
18
Лабораторная работа № 2
ОШИБКИ ИЗМЕРЕНИЯ ВРЕМЕННОГО ПОЛОЖЕНИЯ ИМПУЛЬСНОГО СИГНАЛА
Цель работы: изучение механизмов образования нормальных и аномальных ошибок измерения временного положения радиоимпульсного сигнала и способов обеспечения требуемой точности измерения.
Введение
В лабораторной работе, выполняемой методом компьютерного моделирования, изучаются ошибки измерения временного положения импульсного радиосигнала s(t) , поступающего на вход измерительно-
го приемника в смеси с белым шумом uш(t) :
u(t) s(t t0 ) uш(t), 0 t T .
Входное воздействие u(t) наблюдается на интервале времени T ,
обычно значительно большем, чем длительность импульсного сигнала s(t) . Оптимальный приемник (рис. 1), вырабатывающий оценку
максимального правдоподобия, состоит из фильтра СФ, согласованного с сигналом, амплитудного детектора АД и устройства, определяющего координату абсолютного максимума напряжения на выходе на интервале Т [2].
u(t) |
|
|
|
|
|
t0 |
СФ |
|
АД |
|
max |
||
|
|
|
|
|||
|
|
|
|
|
|
|
Рис. 1. Оптимальный приемник
19
Высокочастотный импульсный сигнал s(t) без внутриимпульсной модуляции (простой сигнал) можно представить в виде
s(t) S(t)cos(2 f0t ) ,
где – начальная фаза, а S(t) – импульсная огибающая, ширина
спектра которой П обычно на несколько порядков меньше несущей частоты: П f0 .
Моделирование фильтрации и детектирования, как правило, выполняют на основе метода комплексной огибающей [3], поскольку
имитация сигнала на высокой частоте f0 (несущей или промежуточной) на временном участке T потребовала бы использования порядка 2 f0T отсчетов. При имитации же процессов на основе метода огибающей число отсчетов уменьшается до 2ПT . Выигрыш в объеме вычислений определяется при этом отношением f0 / П и достигает де-
сятков, сотен или даже тысяч раз.
В работе имитируется сигнал с «колокольной» (гауссовской) огибающей (рис. 2):
|
S(t) S |
0 |
exp (t / )2 . |
|
||
|
|
|
и |
|
||
Параметр и имеет смысл длительности «колокольного» импульса, |
||||||
определяемой по уровню (рис. 2) |
|
|||||
|
|
|
|
S(t и / 2) |
exp( / 4) 0,46 . |
|
|
|
|
|
|
||
|
|
|
|
S0 |
|
|
|
Нетрудно убедиться [3], что нормирован- |
|||||
|
ный амплитудный спектр S( f ) данного сиг- |
|||||
|
нала описывается выражением, по форме по- |
|||||
|
добным выражению самого сигнала: |
|||||
Рис. 2. Гауссовская оги- |
|
|
S( f ) exp π f (1/ τи) 1 2 . |
|||
|
|
|
|
|
||
бающая радиосигнала |
Здесь теперь параметр П 1/ |
имеет смысл |
||||
|
|
|
|
|
и |
|
ширины спектра огибающей. Очевидно, что на «граничных» частотах |
||||||
f П 2 1/(2 и) , |
как и для самого сигнала во временной области, |
|||||
20
