PCM编解码的MATLAB实现.pdf

clc

% clear all

close all

fs=100;

t=1/fs:1/fs:1;

n=8;

s1=cos(2*pi*t);

s=s1/max(abs(s1)); %信号归一化

y=Alaw(s); %A 律量化

plot(t,x,'r')

figure(2)

function y=Alaw(s0)

A=87.6;

A0=1/A;

s=abs(s0);

y=zeros(1,length(s));

y(find(s<=A0))=(A*s(find(s<=A0))/(1+log(A))).*sign(s0(s<=A0));

y(find(s>A0))=((1+log(A*s(find(s>A0))))/(1+log(A))).*sign(s0(s>A0));

function x=invAlaw(y)