MATLAB——PCM编译码
- 一、实验原理
- - 1.掌握PCM编码原理和译码原理
  - 2. 练习使用Matlab编程实现PCM编码和译码
  - 3. 了解失真度的概念，能对译码结果进行失真度分析
- 二、实验原理
- 三、实验要求
- - 1、用Matlab产生一模拟信号，如：或者自己编写一信号，或者找某一语音信号……。
  - 2、在满足抽样定理要求下，对该信号进行抽样，用Matlab编程实现，画出抽样信号，并与原信号进行对比。
  - 3、对抽样信号进行PCM编码，用Matlab编程实现，画出编码波形（矩形脉冲、单极性不归零），并标出某个抽样值的8位编码结果。
  - 4、对经信道传输得到的PCM编码进行PCM解码，用Matlab编程实现，并画出解码波形。
  - 5、对解码结果进行失真度分析。
  - 6、（扩展练习）采用PCM编码，2PSK结合simulink实现语音信号（可使用Matlab自带的语音信号）的传输仿真。
- 四、实验内容
- - 4.1用Matlab产生一模拟信号
  - 4.2画出抽样信号，并与原信号进行对比
  - 4.3对抽样信号进行PCM编码
  - 4.4 PCM解码
  - 4.5失真度分析
  - 4.6 simulink实现语音信号的传输仿真
- 五、引用

MATLAB——PCM编译码

一、实验原理

1.掌握PCM编码原理和译码原理

2. 练习使用Matlab编程实现PCM编码和译码

3. 了解失真度的概念，能对译码结果进行失真度分析

二、实验原理

脉冲编码调制就是把一个时间连续，取值连续的模拟信号变换成时间离散，取值离散的数字信号后在信道中传输。脉冲编码调制就是对模拟信号先抽样，再对样值幅度量化，编码的过程。

抽样，就是对模拟信号进行周期性扫描，把时间上连续的信号变成时间上离散的信号，抽样必须遵循奈奎斯特抽样定理。该模拟信号经过抽样后还应当包含原信号中所有信息，也就是说能无失真的恢复原模拟信号。它的抽样速率的下限是由抽样定理确定的。

量化，就是把经过抽样得到的瞬时值将其幅度离散，即用一组规定的电平，把瞬时抽样值用最接近的电平值来表示,通常是用二进制表示。

量化误差，量化后的信号和抽样信号的差值。量化误差在接收端表现为噪声，称为量化噪声。量化级数越多误差越小，相应的二进制码位数越多，要求传输速率越高，频带越宽。为使量化噪声尽可能小而所需码位数又不太多，通常采用非均匀量化的方法进行量化。非均匀量化根据幅度的不同区间来确定量化间隔，幅度小的区间量化间隔取得小，幅度大的区间量化间隔取得大。一个模拟信号经过抽样量化后，得到已量化的脉冲幅度调制信号，它仅为有限个数值。

编码，就是用一组二进制码组来表示每一个有固定电平的量化值。然而，实际上量化是在编码过程中同时完成的，故编码过程也称为模/数变换，可记作A/D。
为解决均匀量化时小信号量化误差大，音质差的问题，在实际中采用不均匀选取量化间隔的非线性量化方法，即量化特性在小信号时分层密，量化间隔小，而在大信号时分层疏，量化间隔大。在实际中使用的是两种对数形式的压缩特性：A律和μ律。

三、实验要求

1、用Matlab产生一模拟信号，如：或者自己编写一信号，或者找某一语音信号……。

2、在满足抽样定理要求下，对该信号进行抽样，用Matlab编程实现，画出抽样信号，并与原信号进行对比。

3、对抽样信号进行PCM编码，用Matlab编程实现，画出编码波形（矩形脉冲、单极性不归零），并标出某个抽样值的8位编码结果。

4、对经信道传输得到的PCM编码进行PCM解码，用Matlab编程实现，并画出解码波形。

5、对解码结果进行失真度分析。

6、（扩展练习）采用PCM编码，2PSK结合simulink实现语音信号（可使用Matlab自带的语音信号）的传输仿真。

四、实验内容

4.1用Matlab产生一模拟信号

t=0:1/2000:0.1;
f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);plot(t,f);
xlabel('t');
title('f(t)的时域波形');

4.2画出抽样信号，并与原信号进行对比

clear;
clc;
T=0.0005;
t=-0.05:T:0.05;
fs=2000;
sdt=1/fs;
t1=-0.05:sdt:0.05;% 原始信号
f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
%抽样信号
fs=cos(100*pi*t1)+2*sin(200*pi*t1)+3*cos(300*pi*t1);figure;
subplot(2,1,1);plot(t,f);title('原始信号');grid on;
subplot(2,1,2);stem(t1,fs,'.');title('抽样信号');grid on;

4.3对抽样信号进行PCM编码

% 原始信号
f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
%抽样信号
fs=cos(100*pi*t1)+2*sin(200*pi*t1)+3*cos(300*pi*t1);
max = max(abs(fs));
% 原始信号figure;
subplot(2,1,1);plot(t,f);title('原始信号');grid on;
subplot(2,1,2);stem(t1,fs,'.');title('抽样信号');grid on;pcm_encode = PCMcoding(fs);
figure;
stairs(pcm_encode);%绘制信号的阶梯图
axis([0 328 -0.1 1.1]);
title('PCM 编码');axis([0 168 -0.1 1.1]);
set(gca,'xtick',[ 1 9  17 25 33 41 49 57 168]);
xticklabels({'1','9','17','25','33','41','49','57','168'})
grid on;function code=PCMcoding(S)z=sign(S);                                %判断S的正负MaxS=max(abs(S));                         %求S的最大值 S=abs(S/MaxS);                            %归一化Q=2048*S;                                 %量化code=zeros(length(S),8);                  %PCM编码存储矩阵%%% 段落码判断程序for i=1:length(S)if (Q(i)>128)&&(Q(i)<=2048)code(i,2)=1;            %在第五段与第八段之间，段位码第一位都为"1"endif (Q(i)>32)&&(Q(i)<=128)||(Q(i)>512)&&(Q(i)<=2048)code(i,3)=1;            %在第三四七八段内，段位码第二位为"1"endif (Q(i)>16)&&(Q(i)<=32)||(Q(i)>64)&&(Q(i)<=128)||(Q(i)>256)&&(Q(i)<=512)||(Q(i)>1024)&&(Q(i)<=2048)code(i,4)=1;            %在二四六八段内，段位码第三位为"1"endend%%% 段内码判断程序N=zeros(length(S));                              for i=1:length(S)N(i)=bin2dec(num2str(code(i,2:4)))+1;        %找到code位于第几段,bin2dec将二进制整数的文本表示转换为双精度值，enda=[0,16,32,64,128,256,512,1024];                 %13折线各段起始对应的量化单位数b=[1,1,2,4,8,16,32,64];                          %除以16，得到每段的最小量化间隔for i=1:length(S)  q=ceil((Q(i)-a(N(i)))/b(N(i)));              %求出在段内的位置,ceil将 X 的每个元素四舍五入到大于或等于该元素的最接近整数if q==0code(i,(5:8))=[0,0,0,0];                 %如果输入为零则输出"0"else k=dec2bin(q-1,4);                       %将字符数组或字符串转换为数值数组,编码段内码为二进制,dec2bin将十进制整数转换为其二进制表示字符向量code(i,5)=str2num(k(1));code(i,6)=str2num(k(2));code(i,7)=str2num(k(3));code(i,8)=str2num(k(4));end%符号位的判断if z(i)>0code(i,1)=1;elseif z(i)<0code(i,1)=0;end                                        endcode = reshape(code', 1, []);
end

我们选取观察的点是模拟信号中负值最大的点，同时也是抽样信号的第5个点。它的PCM编码结果应该是0111 1111，我们打开Workbench中pcm_encode，观察第33~40位。发现编码结果确实是0111 1111，符合编码逻辑，证明代码编写无误。

4.4 PCM解码

function s=PCMdecoding(encode, max)encode=(reshape(encode',8,length(encode)/8))';l=size(encode,1);a=[0,16,32,64,128,256,512,1024];b=[1 1 2 4 8 16 32 64];c=[0 1.5:15.5];for i=1:lx=encode(i,1);T=bin2dec(num2str(encode(i,(2:4))))+1;Y=bin2dec(num2str(encode(i,(5:8))));if Y==0k(i)=a(T)/2048;elsek(i)=(a(T)+b(T)*c(Y))/2048;endif x==0s(i)=-k(i);elses(i)=k(i);endends = s*max;
end

通过原始信号与PCM编码译码后的对比可以看到能够比较好地重建出原始信号。通过下图能够更好得出对比的效果。

4.5失真度分析

figure;
subplot(2,1,1);plot(t,f);title('原始信号');grid on;
subplot(2,1,2);stem(t1,fs,'.');title('抽样信号');grid on;pcm_encode = PCMcoding(fs);pcm_decode = PCMdecoding(pcm_encode, max);
figure;
plot(t, pcm_decode);hold on
title('PCM 译码');grid on;
plot(t,f);
title('原始信号');grid on;% 计算失真度
da=0;
for i=1:length(t)dc=(f(i)-pcm_decode(i))^2/length(t);da=da+dc;
end
fprintf('失真度是：%.6f\n',da);

失真度我们认为是前后信号之间的均方误差，经过计算得到为：
失真度是：0.010716

4.6 simulink实现语音信号的传输仿真

首先我设计了一个2psk的编码模块，可以将单极性基带二进制码转换为2PSK信号。从Scope模块中我们可以得到Bernoulli Binary模块生成的随机二进制基带信号生成的2PSK波形，验证模块的正确性。

接下来我将之前编写的PCM编码和解码模块和2PSK模块结合起来，模拟语音信号编码、传输和解码过程。我使用了MATLAB自带的语音包crisp，可以发出鸟鸣的声音。
接下来我将之前编写的PCM编码和解码模块和2PSK模块结合起来，模拟语音信号编码、传输和解码过程。语音信号的选择我使用了MATLAB自带的语音包crisp，可以发出鸟鸣的声音。

load chirp % matlab自带语音信号
y1=[y,y];
x=y1(1:20000); %取前20000个采样点
sound(x,Fs);

接下来就是将语音信号进行传输的Simulink仿真了，simout就是经过传输的语音信号，但注意这是pcm八进制编码的形式，必须要进行解码，同时和原信号进行对比，验证传输的正确性。

%% 译码
demodata=simout(2:160001);
zz=pcm_decode(demodata,0.8);
figure;
subplot(1,1,1);
plot(zz);
title('译码的语音信号');
sound(zz,Fs);
%x是初始的语音信号矩阵
figure;
plot(x,'b');
hold on
%zz是译码得到的语音信号
plot(zz,'r');
legend('原语音信号','恢复的语音信号');
title('语音信号对比');

五、引用

[1]邢巨伟.模拟信号的数字化传输-PCM编码[J].黑龙江科技信息,2014(33):116.
[2]陈奎.语音PCM编码通信系统的SIMULINK仿真[J].福建电脑,2012,28(10):115-117.
[3]赵守彬,李鸿刚.利用Matlab实现PCM编码的A律13折线法量化[J].科技信息,2009(36):547.