Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Радиотехнические системы. Лабораторный практикум. Учебно-методическое пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

где 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]