gold matlab 编解码_Matlab PCM编码解码

好久没发博文了，这两天正好在做通原的实验，正好学学怎么用Live Writer。

通常我们会把模拟信号抽样、量化，直到最后变换成二进制符号的集成过程称为脉冲编码调制(Pulse Code Modulation)，实现的主要包括三个步骤：抽样(sampling)、量化(quantization)和编码(coding)。

1、Sampling

Matlab仿真过程中选定模拟信号sin(4*t)，按抽样定理，抽样频率只需高于模拟信号最高频率的两倍即可，但是实际情况下会大于两倍的模拟最高频率。

2、Quantization

量化方法包括：均匀量化和非均匀量化。这里主要讨论非均匀量化，并模拟A率压缩，通过13折现近似实现。

代码如下：

%13折线量化

function y=zhe13(x)

x=x/max(x);

z=sign(x);

x=abs(x);

for i=1:length(x)

if ((x(i)>=0)&x(i)<1/64)

y(i)=16*x(i);

else if ((x(i)>=1/64)&x(i)<1/32)

y(i)=8*x(i)+1/8;

else if ((x(i)>=1/32)&x(i)<1/16)

y(i)=4*x(i)+2/8;

else if ((x(i)>=1/16)&x(i)<1/8)

y(i)=2*x(i)+3/8;

else if ((x(i)>=1/8)&x(i)<1/4) y(i)=x(i)+4/8;

else if ((x(i)>=1/4)&x(i)<1/2)

y(i)=(1/2)*x(i)+5/8;

else y(i)=(1/4)*x(i)+6/8;

end

end

end;

end;

end;

end;

end;

y=z.*y;

3、Coding

编码包括A率13折线编码和解码两个过程，具体原理可参考樊昌信所著《通信原理》第六版中‘模拟信号的数字传输’一章。

代码如下，包括A率13折线编码和A率13折线解码：

% A率13折线编码

function CodeofPCM = C_CodeofPCM(A_Crent,Amax)

%% 极性

Jixing = 1;

if A_Crent<0

Jixing=0;

end

Jixing = dec2bin(Jixing,1);

%% 段号

DuanHao = 0;

A_Crentmaxto1 = floor(abs(A_Crent)/Amax*2048);

LinghuaJG = 0; %分段操作

if A_Crentmaxto1 < 16

DuanHao = 0;

DuanluoQS = 0;

LinghuaJG = 1;

elseif A_Crentmaxto1 < 32

DuanHao = 1;

DuanluoQS = 16;

LinghuaJG = 1;

elseif A_Crentmaxto1 < 64

DuanHao = 2;

DuanluoQS = 32;

LinghuaJG = 2;

elseif A_Crentmaxto1 < 128

DuanHao = 3;

DuanluoQS = 64;

LinghuaJG = 4;

elseif A_Crentmaxto1 < 256

DuanHao = 4;

DuanluoQS = 128;

LinghuaJG = 8;

elseif A_Crentmaxto1 < 512

DuanHao = 5;

DuanluoQS = 256;

LinghuaJG = 16;

elseif A_Crentmaxto1 < 1024

DuanHao = 6;

DuanluoQS = 512;

LinghuaJG = 32;

else

DuanHao = 7;

DuanluoQS = 1024;

LinghuaJG = 64;

end

DuanHao = dec2bin(DuanHao,3);

%% 段内码

if A_Crentmaxto1 == 2048

A_Crentmaxto1 = A_Crentmaxto1-1;

end

DuanNeima1 = floor((A_Crentmaxto1-DuanluoQS)/LinghuaJG);

DuanNeima = dec2bin(DuanNeima1,4);

%% PCM码组 CodeofPCM

PCM = strcat(Jixing,DuanHao, DuanNeima);

CodeofPCM = str2num(PCM);

% A率13折线解码

function A_Crent=D_CodeofPCM(CodeofPCM)

CodeofPCM_int = CodeofPCM;

%% 极性判断

Jixing = 1;

if floor(CodeofPCM_int/10^7) == 0

Jixing = -1;

end

%% 段号

DuanHao = floor(mod(abs(CodeofPCM_int),10^7)/10^4);

DuanHao = num2str(DuanHao);

DuanHao = bin2dec(DuanHao);

switch DuanHao

case 0

DuanluoQS=0;

LinghuaJG=1;

case 1

DuanluoQS=16;

LinghuaJG=1;

case 2

DuanluoQS=32;

LinghuaJG=2;

case 3

DuanluoQS=64;

LinghuaJG=4;

case 4

DuanluoQS=128;

LinghuaJG=8;

case 5

DuanluoQS=256;

LinghuaJG=16;

case 6

DuanluoQS=512;

LinghuaJG=32;

case 7

DuanluoQS=1024;

LinghuaJG=64;

end

%%段内码提取

DuanNeima = mod(abs(CodeofPCM_int),10^4);

DuanNeima = num2str(DuanNeima);

DuanNeima = bin2dec(DuanNeima);

absofA_Crent = DuanluoQS + LinghuaJG*DuanNeima + 0.5*LinghuaJG;

A_Crent=absofA_Crent*Jixing/2048;

PCM具体实现：

1、选取模拟信号x = sin(4*t)，抽样频率为10Hz，得到：

2、量化，13折线量化：

3、编码：

编码得到：

Columns 1 through 8

11111110    11111011    11101110     1111001     1111101     1111111     1111111     1111111

Columns 9 through 16

1111101     1111010    10000000    11111010    11111101    11111111    11111111    11111111

Columns 17 through 21

11111101    11111001     1101110     1111011     1111110

解码后与编码比较：