Сгенерируйте последовательность данных.

data = 2:2:12

data = 1×6

     2     4     6     8    10    12

Сожмите последовательность данных при помощи компрессора mu-закона. Установите значение для mu к 255. Сжатая последовательность данных теперь располагается между 8,1 и 12.

compressed = compand(data,255,max(data),'mu/compressor')

compressed = 1×6

    8.1644    9.6394   10.5084   11.1268   11.6071   12.0000

Расширьте сжатую последовательность данных при помощи расширителя mu-закона. Расширенная последовательность данных почти идентична исходной последовательности данных.

expanded = compand(compressed,255,max(data),'mu/expander')

expanded = 1×6

    2.0000    4.0000    6.0000    8.0000   10.0000   12.0000

Вычислите разность между исходной последовательностью данных и расширенной последовательностью.

diffvalue = expanded - data

diffvalue = 1×6
10-14 ×

   -0.0444    0.1776    0.0888    0.1776    0.1776   -0.3553

Сгенерируйте последовательность данных.

data = 1:5;

Сожмите последовательность данных при помощи компрессора A-закона. Установите значение для к 87,6. Сжатая последовательность данных теперь располагается между 3,5 и 5.

compressed = compand(data,87.6,max(data),'A/compressor')

compressed = 1×5

    3.5296    4.1629    4.5333    4.7961    5.0000

Расширьте сжатую последовательность данных при помощи расширителя A-закона. Расширенная последовательность данных почти идентична исходной последовательности данных.

expanded = compand(compressed,87.6,max(data),'A/expander')

expanded = 1×5

    1.0000    2.0000    3.0000    4.0000    5.0000

Вычислите разность между исходной последовательностью данных и расширенной последовательностью.

diffvalue = expanded - data

diffvalue = 1×5
10-14 ×

         0         0    0.1332    0.0888    0.0888

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

Создайте экспоненциальный сигнал и вычислите его максимальное значение.

sig = exp(-4:0.1:4);
V = max(sig);

Квантуйте сигнал при помощи интервалов равной длины. Установите раздел и значения книги шифров, приняв 6-битное квантование. Вычислите среднеквадратическое искажение.

partition = 0:2^6 - 1;
codebook = 0:2^6;
[~,qsig,distortion] = quantiz(sig,partition,codebook);

Сожмите сигнал при помощи compand функция, сконфигурированная, чтобы применить метод mu-закона. Примените квантование и расширьте квантованный сигнал. Вычислите среднеквадратическое искажение сигнала companded.

mu = 255; % mu-law parameter
csig_compressed = compand(sig,mu,V,'mu/compressor');
[~,quants] = quantiz(csig_compressed,partition,codebook);
csig_expanded = compand(quants,mu,max(quants),'mu/expander');
distortion2 = sum((csig_expanded - sig).^2)/length(sig);

Сравните среднеквадратическое искажение для квантования по сравнению с объединенным компандированием и квантования. Искажение для companded и квантованного сигнала является порядком величины ниже, чем искажение квантованного сигнала. Интервалы равной длины хорошо подходят для логарифма экспоненциального сигнала, но не хорошо подходят для самого экспоненциального сигнала.

[distortion, distortion2]

ans = 1×2

    0.5348    0.0397

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

plot([sig' qsig' csig_expanded']);
title('Comparison Between Original, Quantized, and Expanded Signals');
xlabel('Interval');
ylabel('Apmlitude');
legend('Original','Quantized','Expanded','location','nw');
axis([0 70 0 20])