matlab打开pcm文件格式,用MATLAB实现PCM编码

PCM编码程序

%修改抽样间隔和程序最后两行

%可以观察量化输入输出关系

n=input('请输入量化级数,k=[8]');

if isempty(n),

n=8;

end

%产生题目要求的函数

endtime = 2; %结束时间

dt = 0.1; %抽样间隔

t = 0:dt:endtime;

a = zeros(size(t)); %产生一个长度为N的全0序列

for tt=1:1/dt

a(tt+[0:endtime/2-1]*(2/dt))=t(tt); %t when 0 <= t < 1

end

for tt=1/dt+1:2/dt

a(tt+[0:endtime/2-1]*(2/dt))=2-t(tt); % -t+2 when 1 <= t < 2 end;

%A率非线性

A = 87.56;

amax=max(abs(a));

c = zeros(size(a));

for i = 1:length(a);

if((a(i)/amax)<=1/A)

c(i) = A*a(i)/(1+log(A));

end

if((a(i)/amax)>1/A)

c(i) = (1+log(A*a(i)))/(1+log(A));

end

end

%均匀量化

c_quan=c;

b_quan=c_quan;

d=1/n;%量化间隔

q=d.*[1:n];

q=q-d/2;%量化电平

for i=1:n

%定位第i个量化间隔码子

c_quan(find((q(i)-d/2<=c_quan) & (c_quan<=q(i)+d/2)))=...

q(i).*ones(1,length(find((q(i)-d/2<=c_quan)&(c_quan<=q(i)+d/2)))); %赋值为相应的量化电平