基于MATLAB的A律非均匀量化

（1）A律压缩特性：

其中，A 为压缩系数；y 为归一化的压缩器输出电压；x 为归一化的压缩器输入电压
（2）MATLAB实现

function y=a_pcm(x,a)   %x为一个序列，值在0到1，A律量化将得到的结果存在序列y中
t=1/a;                  %a为一个正实数，大于1
for i=1:length(x)
    if x(i)>=0          %判断输入序列值是否大于0
        if (x(i)<=t)
            y(i)=(a*x(i))/(1+log(a)); %若值小于1/a，用此算法
        else
            y(i)=(1+log(a*x(i)))/(1+log(a)); %若值大于1/a，用此算法
        end
    else
        if (x(i)>=-t)
            y(i)=-(a*-x(i))/(1+log(a));     %若值小于0，用此算法
        else
            y(i)=-(1+log(a*-x(i)))/(1+log(a));
        end
    end
end

上述为A律量化函数

x=[0:0.001:1];%定义幅度序列
y1=a_pcm(x,1);%参数为1的A律曲线
y2=a_pcm(x,10);
y3=a_pcm(x,87.65);
pl