filter1.m

来自「自己编写GUI界面实现语音增强」· M 代码 · 共 16 行

h=fir1(20,0.15,'high',kaiser(21,3.6));
[Ha,f]=freqz(h,1,256,8000);
Ha_db=20*log10(abs(Ha));
figure(1);
plot(f,Ha_db,1.27,-6,'r*');
Ha=abs(Ha);
[y,fs,bit]=wavread('g:\ex3.wav');
for i=1:256:256*100,
Y=abs(fft(y(i:i+255),256));
Y=Y.*Ha;
Ha(10)
y(i:i+255)=ifft(Y,256);
end;
wavwrite(y,fs,bit,'g:\filter1.wav');
%[Y,f]=psd(y,256,8000);
%plot(f,Y,'r-');