数字语音处理之中心削波

%读入数据 采样fs=8000Hz 采样位数16bit
fid=fopen('voice.txt','rt');
[a,count]=fscanf(fid,'%f',[1,inf]);
L=length(a);
m=max(a);
for i=1:L;
  a(i)=a(i)/m;
end

%找到归一化后数据的最大值和最小值
m=max(a);
n=min(a);
%为保证幅度之和横坐标轴对称，采用计算公式是n+(m-n)/2，合并为(m+n)/2
ht=(m+n)/2;
for i=1:L;
  a(i)=a(i)-ht;
end
figure(1);
subplot(2,1,1);
plot(a,'k');
axis([0,1711,-1,1]);
title('中心削波前语音波形');
xlabel('样点数');
ylabel('幅度');

coeff=0.7;
th0=max(a)*coeff;
for k=1:L;
   if a(k)>=th0
      a(k)=a(k)-th0;