基于matlab语音信号分析的低通高通滤波源程序

function hfdata=shighpassf(f,fdata,fc) N=length(f); df=f(2)-f(1); fs=N*df; lowpa=zeros(1,N); lfdata=zeros(1,N); wc=fc*2*pi for h=0:N-1 w=h*df*2*pi; if w>2*pi*fs/2 w=w-2*pi*fs; end lowpa(h+1)=(-1*w^2)/(-1*w^2+sqrt(2)*i*w*wc+wc^2); end hfdata=lowpa.*fdata; w1=[0:df:fs-df]; subplot(2,1,1); plot(w1,abs(lowpa)); xlabel('频率'); ylabel('高通滤波器响应'); subplot(2,1,2); plot(w1,abs(hfdata)); xlabel('频率'); ylabel('高通滤波后频谱'); pause;