采用MATLAB对A律13折线PCM信源编码进行仿真，基本要求： 1. 输入一个单音信号：x(n)=sin(2*pi*f0*t) ，该单音信号的频率f0=800hz ，对其进行采样获取x(n)=sin(2*pi*f0*n/fs) ，其中，采样频率为fs=32000hz ，n的取值范围为0～255，绘制单音的波形，横轴为时间。 2. 把单音的浮点数据转化为量化单位(即整数)。 3. 对输入量化数据进行13折线PCM编码，将编码信号直接进行译码输出，绘制译码信号的时域波形。

以下是MATLAB代码实现：

% 采样频率
fs = 32000;
% 单音频率
f0 = 800;
% 采样点数
n = 0:255;
% 采样信号
x = sin(2*pi*f0*n/fs);
% 绘制单音波形
plot(n, x);
xlabel('采样点');
ylabel('幅值');
title('单音波形');

% 量化位数
q = 13;
% 量化电平数
L = 2^q;
% 量化电平区间
delta = 2/2^q;
% 量化数据
xq = round(x/delta);
% 绘制量化信号波形
figure;
plot(n, xq);
xlabel('采样点');
ylabel('量化电平');
title('量化信号波形');

% PCM编码
code = zeros(size(xq));
for i = 1:length(xq)
    % 找到量化电平在量化区间内的位置
    [~, index] = min(abs((xq(i) + 0.5)*delta - (-1:2/3:L/3)));
    % 编码为4位二进制数
    code(i) = index - 1;
end
% PCM译码
xq_hat = (code + 0.5)/L*2 - 1;
xq_hat = xq_hat(2:end); % 去掉第一个无效码字

% 绘制译码信号波形
figure;
plot(n(2:end), xq_hat);
xlabel('采样点');
ylabel('幅值');
title('译码信号波形');

运行后会依次绘制单音波形、量化信号波形和译码信号波形。