untitled.m

理想带通滤波器的逼近,在某种程度上可以代替理想带通滤波器.

input_filename='e:\B012_Ori.wav';
>> N=240;
>> Y=wavread(input_filename,[1 18000]);
>> L=length(Y);
>> LL=length(Y)/N;
>> Em=zeros(1,(LL-1)*240);
>> for ii=1:(LL-1)*240
temp=Y(ii:ii+240-1);
temp=temp.*hamming(240);
Em(ii)=sum(temp.*temp);
end
plot(jj,Em);
jj=[1:(LL-1)*240];
subplot(4,1,1);
YY=fft(Em,256);
>> f=500*(0:128)/256;
>> plot(f,YY(1:129));