Используйте рабочий режим сброса mlseeq эквалайзер. Демодулируйте сигнал и проверьте частоту битовых ошибок.

Укажите порядок модуляции, глубину отслеживания корректора, количество выборок на символ и длину сообщения.

M = 2;
tblen =  10; 
nsamp = 2;
msgLen = 1000;

Создайте эталонное созвездие.

const = pammod([0:M-1],M);

Создание сообщения со случайными данными. Модулируйте и увеличивайте частоту сигнала.

msgData = randi([0 M-1],msgLen,1); 
msgSym = pammod(msgData,M);
msgSymUp = upsample(msgSym,nsamp); 

Фильтрация данных по каналу искажения и добавление гауссова шума к сигналу.

chanest = [0.986; 0.845; 0.237; 0.12345+0.31i]; 
msgFilt = filter(chanest,1,msgSymUp); 
msgRx = awgn(msgFilt,5,'measured');

Выровнять и затем демодулировать сигнал для восстановления сообщения. Чтобы инициализировать эквалайзер, предоставьте оценку канала, опорную совокупность, глубину отслеживания эквалайзера, количество выборок на символ и установите режим работы на сброс. Проверьте частоту битовых ошибок сообщения. Результаты могут отличаться, поскольку в этом примере используются случайные числа.

eqSym = mlseeq(msgRx,chanest,const,tblen,'rst',nsamp);
eqMsg = pamdemod(eqSym,M);

[nerrs ber] = biterr(msgData, eqMsg)

nerrs = 1

ber = 1.0000e-03

Восстановите сообщение, содержащее преамбулу, выровняйте сигнал и проверьте частоту ошибок символов.

Укажите порядок модуляции, глубину отслеживания корректора, количество выборок на символ, преамбулу и длину сообщения.

M = 4; 
tblen = 16;
nsamp = 1;
preamble = [3;1];
msgLen = 500;

Создайте эталонное созвездие.

const = pskmod(0:3,4);

Создание сообщения с использованием случайных данных и добавление преамбулы к сообщению. Модулировать случайные данные.

msgData = randi([0 M-1],msgLen,1);
msgData = [preamble; msgData];
msgSym = pskmod(msgData,M);

Фильтрация данных по каналу искажения и добавление гауссова шума к сигналу.

chcoeffs = [0.623; 0.489+0.234i; 0.398i; 0.21];
chanest = chcoeffs;
msgFilt = filter(chcoeffs,1,msgSym);
msgRx = awgn(msgFilt,9,'measured');

Выровнять принятый сигнал. Для конфигурирования эквалайзера предоставьте оценку канала, опорную совокупность, глубину отслеживания эквалайзера, режим работы, количество выборок на символ и преамбулу. Те же самые символы преамбулы появляются в начале вектора сообщения и в синтаксисе для mlseeq. Поскольку система не использует postamble, пустой вектор указывается в качестве последнего входного аргумента в этом mlseeq синтаксис.

Проверьте частоту символьных ошибок выровненного сигнала. Результаты выполнения варьируются из-за использования случайных чисел.

eqSym = mlseeq(msgRx,chanest,const,tblen,'rst',nsamp,preamble,[]);
[nsymerrs,ser] = symerr(msgSym,eqSym)

nsymerrs = 8

ser = 0.0159

Использовать непрерывный режим работы mlseeq эквалайзер. Демодуляция принятых сигнальных пакетов и проверка статистики ошибок символов.

Укажите порядок модуляции, глубину отслеживания корректора, количество выборок на символ, длину сообщения и количество обрабатываемых пакетов.

M = 4;
tblen =  10; 
nsamp = 1;
msgLen = 1000; 
numPkts = 25;

Создайте эталонное созвездие.

const = pskmod(0:M-1,M);

Установите начальные входные параметры для метрики, состояний и входов эквалайзера в пустые векторы. Эти начальные назначения представляют параметры для первого переданного пакета.

eq_metric = [];
eq_states = [];
eq_inputs = [];

Назначение переменных для статистики ошибок символов.

ttlSymbErrs = 0;
aggrPktSER = 0;

Отправка и получение нескольких пакетов сообщений в цикле моделирования. Между передачей пакета и приемом каждый пакет фильтруется через канал искажения и добавляется гауссов шум.

for jj = 1:numPkts

Создание сообщения со случайными данными. Модулируйте сигнал.

    msgData = randi([0 M-1],msgLen,1); 
    msgMod = pskmod(msgData,M); 

Фильтрация данных по каналу искажения и добавление гауссова шума к сигналу.

    chanest = [.986; .845; .237; .12345+.31i]; 
    msgFilt = filter(chanest,1,msgMod);
    msgRx = awgn(msgFilt,10,'measured');

Выравнивание полученных символов. Для конфигурирования эквалайзера предоставьте оценку канала, опорную совокупность, глубину отслеживания эквалайзера, режим работы, количество выборок на символ и информацию инициализации эквалайзера. Для эквалайзера задан непрерывный режим работы. В непрерывном режиме работы информация инициализации эквалайзера (метрика, состояния и входы) возвращается и используется в качестве входов в следующей итерации for цикл.

    [eqSym,eq_metric,eq_states,eq_inputs] = ...
        mlseeq(msgRx,chanest,const,tblen,'cont',nsamp, ...
        eq_metric,eq_states,eq_inputs);

Сохраните статистику ошибок символов. Обновите статистику ошибок символов с помощью агрегированных результатов. Отображение общего количества ошибок. Результаты могут отличаться, поскольку в этом примере используются случайные числа.

    [nsymerrs,ser] = symerr(msgMod(1:end-tblen),eqSym(tblen+1:end));
    ttlSymbErrs = ttlSymbErrs + nsymerrs;
    aggrPktSER = aggrPktSER + ser;
end
printTtlErr = 'A total of %d symbol errors over the %d packets received.\n';
fprintf(printTtlErr,ttlSymbErrs,numPkts);

A total of 155 symbol errors over the 25 packets received.

Отображение совокупной частоты ошибок символов.

printAggrSER = 'The aggregate symbol error rate was %6.5d.\n';
fprintf(printAggrSER,aggrPktSER/numPkts);

The aggregate symbol error rate was 6.26263e-03.