close all;
clear all;
[y,Fs,bits] = wavread('G:\113.wav');
s=y(1001:129256);
ss=s';
for m=1:1001;
yy(m,:)=ss((m-1)*128+1:(m-1)*128+256);%分组1000×256重叠128点
end;
yy2=yy.^2;
for m=1:1001;
shortt(m)=sum(yy2(m,:));%时域能量
end;
a3=abs(fft(yy,256,2));%fft变换取模幅度值
a4=a3.^2;%功率谱值，
x=zeros(1001,123);
x1=zeros(1001,123);
x=a3(:,6:128); %滤波
x1=a4(:,6:128);
t=sum(x,2); %求和，幅度能量
t1=sum(x1,2); %功率能量
zs=0;
for m=1:1001 %求过零率
for sub=2:256
zs=zs+abs(sign(yy(m,sub))-sign(yy(m,sub-1)))/2;
end;
guozero(m)=zs;
zs=0;
end;
for m=6:998;%前向3点后向5点
if t1(m)>8000
d(m)=1;
elseif t1(m-1)>8000
d(m)=1;
elseif t1(m-2)>8000;
d(m)=1;
elseif t1(m-3)>8000
d(m)=1;
elseif t1(m-4)>8000
d(m)=1;
elseif t1(m-5)>8000
d(m)=1;
elseif t1(m+1)>8000
d(m)=1;
elseif t1(m+2)>8000
d(m)=1;
elseif t1(m+3)>8000
d(m)=1;
else d(m)=0;
end;
end;
for m=4:998
%if t1(m)>4500
if t1(m-1)>7500
d1(m)=1;
%end;
elseif t1(m-2)>7500
d1(m)=1;
%end;
elseif t1(m-3)>7500
d1(m)=1;
% end;
elseif t1(m+1)>7500
d1(m)=1;
%end;
elseif t1(m+2)>7500
d1(m)=1;
%end;
elseif t1(m+3)>7500
d1(m)=1;
%end;
else d1(m)=0;
end;
end;
for m=1:998;
m
n=s((m-1)*128+1:(m-1)*128+256);%分组重叠128点
if d(m)==1
n=n.*hanning(256);%频域处理加汉明窗
else
n=(1-2*31.25*7/8000)*n.*hamming(256);%加高通滤波器
end;
n_fft=fft(n);%fft变换
n_abs=abs(n_fft);%取模
n_angle=angle(n_fft);%取相位
n_abs2=n_abs.^2;%平方得到功率谱
check(m)=sum(n_abs2(6:128));%功率谱的能量
chechk1(m)=mean(n_abs2);%平均能量
if d(m)==0 %无声检测
noise_abs2=n_abs2;;%更新噪声
n_abs2(5)=0.4*n_abs2(5);
% n_abs2(253)=0.4*n_abs2(253);
n_abs2=0.005*n_abs2;%噪声帧衰减
end;
sum1=mean(noise_abs2);%噪声帧的平均能量
if d(m)==1 %语音帧处理
if t1(m)>4500 %能量较大的帧，做谱减法
for sub=1:256
if(n_abs2(sub)-2.0*sum1)>0
n_abs2(sub)=n_abs2(sub)-2.0*sum1;
else
n_abs2(sub)=0.04*n_abs2(sub);%谱减时遇到负值的处理
end;
end;
elseif t1(m)>3000
n_abs2=0.3*n_abs2; %过度帧整体衰减
else
n_abs2=0.1*n_abs2;
end;
end;
n_abs=n_abs2.^0.5;
z_fft=n_abs.*(cos(n_angle)+i*sin(n_angle));
z=ifft(z_fft); %逆FFT变换
z=real(z);
p((m-1)*128+65:(m-1)*128+192)=z(65:192);%输出文件
end;
wavwrite(p,Fs,'G:\002211.wav');
