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

Lab_code

.pdf
Скачиваний:
90
Добавлен:
11.04.2015
Размер:
1.48 Mб
Скачать

Вероятность ошибки в канале равна p

указана в таблице 3.1.

 

 

Таблица 3.1 Варианты вероятностей ошибок в канале

 

 

 

Вариант

0

1

2

3

4

5

6

7

8

9

Р

0,02

0,018

0,016

0,014

0,012

0,001

0,008

0,006

0,004

0,002

5. Методом статистических испытаний получить и построить по экспериментальным данным зависимость выигрыша от кодирования при изменении отношения сигнал/шум в таких пределах, чтобы вероятность ошибки изменялась в пределах от 0,1 до 0,001.

Выполнить статистические испытания СПИ с заданным корректирующим кодом не менее, чем для пяти значений отношения сигнал/шум.

Выигрыш от кодирования можно определить по формуле

kис f (H p )

pв

,

pвы

 

 

где pв и pвы – вероятности ошибок в кодовых символах на входе и выходе декодера (в таблицах результатов испытаний АРМ этим вероятностям соответствуют эквивалентные вероятности на входе и выходе УЗО).

6. Построить (или вывести на принтер ЭВМ) таблицу экспериментальных результатов для одной точки статистических испытаний. Для этой же точки построить (или вывести на принтер ЭВМ) распределение кратностей ошибок на входе и выходе декодера.

5. Содержание отчета

Отчет должен содержать:

1.результаты измерений в виде таблиц, графиков с соответствующими заголовками и пояснениями,

2.оценку погрешности измерений,

3.краткие выводы.

21

Приложение 3.1

ЗНАЧЕНИЯ ФУНКЦИЙ

 

 

 

1

 

e

x 2

 

1

 

 

1

 

 

 

t 2

 

 

w(x )

 

 

2 ; V (x )

[1 (x )]

 

e

2 dt

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2 x

 

 

 

x

w(x)

 

 

 

 

 

V(x)

 

 

 

 

x

 

w(x)

 

 

 

 

V(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,00

0,39894

 

 

 

 

0,50000

 

 

 

 

2,75

0,009094

 

0,002980

0,10

0,39695

 

 

 

 

0,46017

 

 

 

 

2,80

0,007915

 

0,002555

0,20

0,39104

 

 

 

 

0,42074

 

 

 

 

2,85

0,006873

 

0,002186

0,30

0,38139

 

 

 

 

0,38209

 

 

 

 

2,90

0,005953

 

0,001866

0,40

0,36827

 

 

 

 

0,34458

 

 

 

 

2,95

0,005143

 

0,001589

0,50

0,35207

 

 

 

 

0,30854

 

 

 

 

3,00

0,004432

 

0,001350

0,60

0,33322

 

 

 

 

0,27425

 

 

 

 

3,05

0,003810

 

0,001144

0,70

0,31225

 

 

 

 

0,24196

 

 

 

 

3,10

0,003267

 

0,000968

0,80

0,28969

 

 

 

 

0,21186

 

 

 

 

3,15

0,002794

 

0,000816

0,90

0,26609

 

 

 

 

0,18406

 

 

 

 

3,20

0,002384

 

0,000687

1,00

0,24197

 

 

 

 

0,15866

 

 

 

 

3,25

0,002029

 

0,000577

1,10

0,21785

 

 

 

 

0,13567

 

 

 

 

3,30

0,001723

 

0,000483

1,20

0,19419

 

 

 

 

0,11507

 

 

 

 

3,35

0,001459

 

0,000404

1,30

0,17137

 

 

 

 

0,09680

 

 

 

 

3,40

0,001232

 

0,000337

1,40

0,14973

 

 

 

 

0,08076

 

 

 

 

3,45

0,001038

 

0,000280

1,50

0,12952

 

 

 

 

0,06681

 

 

 

 

3,50

0,000873

 

0,000233

1,60

0,11092

 

 

 

 

0,05480

 

 

 

 

3,55

0,000732

 

0,000193

1,70

0,09405

 

 

 

 

0,04457

 

 

 

 

3,60

0,000612

 

0,000159

1,80

0,07895

 

 

 

 

0,03593

 

 

 

 

3,65

0,000510

 

0,000131

1,90

0,06562

 

 

 

 

0,02872

 

 

 

 

3,70

0,000425

 

0,000108

2,00

0,05399

 

 

 

 

0,02275

 

 

 

 

3,75

0,000353

 

0,000088

2,05

0,04879

 

 

 

 

0,02018

 

 

 

 

3,80

0,000292

 

0,000072

2,10

0,04398

 

 

 

 

0,01786

 

 

 

 

3,85

0,000241

 

0,000059

2,15

0,03955

 

 

 

 

0,01578

 

 

 

 

3,90

0,000199

 

0,000048

2,20

0,03547

 

 

 

 

0,01390

 

 

 

 

3,95

0,000163

 

0,000039

2,25

0,03174

 

 

 

 

0,01222

 

 

 

 

4,00

0,0001338

 

0,000032

2,30

0,02833

 

 

 

 

0,01072

 

 

 

 

4,10

0,00008926

 

0,0000207

2,35

0,02522

 

 

 

 

0,00939

 

 

 

 

4,20

0,00005894

 

0,0000133

2,40

0,02239

 

 

 

 

0,00820

 

 

 

 

4,30

0,00003854

 

0,0000085

2,45

0,01984

 

 

 

 

0,00714

 

 

 

 

4,40

0,00002494

 

0,0000054

2,50

0,01753

 

 

 

 

0,00621

 

 

 

 

4,50

0,00001598

 

0,0000034

2,50

0,017528

 

 

 

 

0,006210

 

 

 

 

4,60

0,00001014

 

0,0000022

2,55

0,015449

 

 

 

 

0,005386

 

 

 

 

4,70

0,00000637

 

0,0000

2,60

0,013583

 

 

 

 

0,004661

 

 

 

 

4,80

0,00000396

 

0,0000

2,65

0,011912

 

 

 

 

0,004025

 

 

 

 

4,90

0,00000244

 

0,00000

2,70

0,010421

 

 

 

 

0,003467

 

 

 

 

5,00

0,00000151

 

0,0000003

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

Лабораторная работа №4 ИССЛЕДОВАНИЕ СВЕРТОЧНЫХ КОДОВ И ДЕКОДЕРА ВИТЕРБИ

В СРЕДЕ MATLAB

1. Цель работы

Ознакомление с методами построения структурных схем в пакете расширения Simulink. Экспериментальное исследование исправляющей способности сверточных кодов с декодированием по Витерби.

2. Предварительная подготовка 2.1 Ознакомиться с описанием работы и изучить по указанной ниже литературе следующие вопросы:

-методика построения структурных схем в пакете расширения Simulink

среды Matlab

-сверточные коды, параметры и основные характеристики

-методы кодирования и декодирования сверточных кодов, энергетический выигрыш от кодирования.

-Декодер Витерби, выбор параметров, вероятность ошибки декодирования.

3.Исследование сверточных кодов в Simulink

Для исследования сверточных кодов запускаем файл модели vitharddec.mdl, производим выбор параметров сверточного кода и декодера Витерби в соответствии с вариантом. Структурная схема модели показана на рисунке 3.1

Схема также может быть построена в рабочем окне Симулинк используя библиотеки элементов.

Блок Data Source генерирует псевдослучайную последовательность с распределением Бернулли. Настройка блока выполняется при двойном щелчке правой кнопкой мыши.

Значения параметров:

Probability of a zero – вероятность появления нулевых значений на выходе блока

Initial seed – начальное значение для генератора случайных чисел Sample time – интервал дискретизации

Frame-based outputs – настройка выхода блока в виде фреймов Samples per frame – количество отсчетов во фрейме

23

Рисунок 4.1 Структурная схема модели исследования сверточных кодов и декодера Витерби

Рисунок 4.3 Блок Data Source

Блок Convolutional Encoder, рисунок 4.4, формирует кодированную последовательность двоичного сверточного кода в соответствии с выбранной схемой кодера и порождающими полиномами.

24

Рисунок 4.4 Блок Convolutional Encoder

Параметры блока Convolutional Encoder

Для определения сверточного кода применяется его решетчатая структура. Параметр Trellis structure задает сверточный код через его длину кодового ограничения, генераторные полиномы, и, при необходимости, через обратные связи. В поле блока Trellis structure задается команда poly2trellis, которая преобразует генераторные полиномы в решетчатую форму представления кодера. Например, для задания кода с кодовым ограничением 7, и генераторными полиномами в восьмеричной форме 171 и 133 применяется команда poly2trellis(7, [171 133]). Первоначально регистры кодера находятся в нулевом состоянии. В непрерывном режиме Continuous mode блок сохраняет состояние кодера в конце каждого фрейма для использования в последующем фрейме.

Выход кодера сверточного кода подается на блок двоичной фазовой модуляции BPSK Modulator Baseband.

25

Рисунок 4.5 Блок фазовой модуляции BPSK Modulator Baseband

Phase offset (rad) – фаза нулевой точки в сигнальном созвездии, Output data type – тип данных выходного сигнала.

Далее модулированный сигнал подается на вход канала с белым гауссовским шумом AWGN Channel

Рисунок 4.6 Блока канала с гауссовым шумом AWGN Channel

26

Блока канала с гауссовым шумом AWGN Channel добавляет к входному сигналу шум с нормальным распределением значений.

Параметр Initial seed инициализирует генератор шума. Дисперсию шума можно регулировать в нескольких режимах.

Режим Signal to noise ratio (Eb/No)

Eb/No – отношение энергии на один бит к спектральной плотности мощности шума, в децибелах,

Number of bits per symbol – количество бит на символ,

Input signal power – мощность входного сигнала, Symbol period – период символов.

Режим Signal to noise ratio (Es/No)

Es/No – отношение энергии сигнала к спектральной плотности мощности, в децибелах,

Input signal power – мощность входного сигнала, Symbol period – период символов.

Режим Signal to noise ratio (SNR)

SNR – отношение мощности сигнала к мощности шума Input signal power – мощность входного сигнала.

Изменение длительности периода сигнала влияет на дисперсию шума, добавляемого к отсчетам сигнала.

2

 

Ps Ts

 

n

 

 

Es / N0

 

 

 

T 10 10

 

 

 

0

 

 

Ps - мощность сигнала,

Ts - Длительность символа, T0 - интервал дискретизации.

Для комплексного входного сигнала связь между Eb/No, Es/No, и SNR выражается следующим образом:

Es/N0 = (Ts /T0 ) · SNR Es/N0 = Eb/N0 + 10log10(k), дБ

К – количество информационных бит на один символ. Для реальных сигналов

Es/N0 = 0.5(Ts /T0 ) · SNR.

Выход канала подается на вход демодулятора двоичной ФМ

BPSK Demodulator Baseband.

27

Рисунок 4.7 Демодулятор двоичной ФМ BPSK Demodulator Baseband.

Параметры блока

Phase offset (rad) – начальная фаза нулевой точки,

Samples per symbol – количество отсчетов на один период входного сигнала.

Выход демодулятора после преобразования типов сигнала (если это требуется) подается на вход блока декодера Витерби - Viterbi Decoder.

Параметры блока декодера Витерби

Trellis structure – решетчатая структура сверточного кода, применяется то же значение, что и в кодере сверточного кода.

Punctured code – устанавливается режим выкалывания.

Enable erasures input port – разрешение на применение стирания в канале.

Decision type – тип выхода демодулятора, входа декодера. Возможно неквантованные, жесткие и мягкие решения.

Number of soft decision bits – разрядность представления квантованных значений (при мягких решениях)

Traceback depth – глубина ветвей кодовой решетки, применяемой при нахождении выжившего пути

28

Рисунок 4.8 Блок декодера Витерби - Viterbi Decoder

Выход декодера подается на блок вычисления вероятности ошибки Error Rate Calculation. В данном блоке происходит подсчет количества ошибок и вычисляется вероятность ошибки декодирования. На один вход блока подается неискаженный каналом сигнал, с выхода источника дискретного сигнала, на второй вход – сигнал с выхода декодера. В блоке происходит подсчет несовпадающих символов, с учетом задержки на декодирование в декодере Витерби.

Параметры блока Error Rate Calculation

Receive delay – количество отсчетов сигнала, на которое выходной сигнал задержан относительно переданных данных.

Output data – передача результата в рабочее пространство.

Target number of errors – требуемое количество ошибок, необходимое для достижения заданной точности.

29

Рисунок 4.9 Блок вычисления вероятности ошибки Error Rate

Calculation

4.Порядок выполнения работы

4.1Включить ЭВМ. Вызвать и запустить программу Matlab

4.2Исследование сверточных кодов. Открыть файл модели vitharddec.mdl, настроить параметры кодера, декодера. С помощью BerTool получить зависимость вероятности ошибки на выходе от отношения сигнал/шум в канале. С помощью BerTool получить на одном графике результаты декодирования с применением кода и без него. Выполнить исследования вероятности ошибки декодирования для нескольких сверточных кодов, определяемых следующими параметрами

poly2trellis(9, [561 753]) poly2trellis(7, [171 133]) poly2trellis(3,[6 7])

30

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