MATLAB课程设计论文-FIR数字高通滤波器设计

t(2,1,2);plot(xlabel,X(1:150001))Р原始语音信号(无噪声)采样后X1的频谱图:Р图3原始语音信号采样后X1的频谱图Р加噪声语音信号采样后X的频谱图如图4所示:Р图4加噪声语音信号采样后X的频谱图Р由图分析信号噪声频率为200Hz,滤波器设置参数如下:Рfp = 300; %设置通带截止频率Рfs = 500;%设置阻带最低频率Рdelta_w = 2*pi*(fs-fp)/Fs;%fs-fp是过度带宽;转换为数字频率过渡带delta_w=0.0285Рwc = pi*(fs+fp)/Fs;%数字截止频率wc=0.0570Р3.采用矩形窗设计FIR高通滤波器РN = ceil(1.8*pi/delta_w); %计算滤波器的阶次N=199,矩形:1.8*piРif N/2 ==0Р N = N+1; %确保N为奇数РendРtau = (N-1)/2;%tau=363.5000Рfor n = 1:NРh(n) =(sin(pi*(n-1-tau))-sin(wc*(n-1-tau)))/(pi*(n-1-tau));РendРh((N+1)/2)=1-wc/pi; %根据滤波器类型确定极限值РWn=(boxcar(N));Рh1 = h.*Wn'; %加窗后得到h1,即为设计的FIR滤波器。Рh2 = fir1(N-1,1-wc/pi,boxcar(N));Рfigure;subplot(2,1,1);stem(h1); %滤波器单位冲击响应序列及幅频响应РF=abs(fft(h1)); xlabel1 = 0:Fs/(N-1):Fs/2;Рsubplot(2,1,2);plot(xlabel1,F(1:(N+1)/2));Р滤波器单位冲击响应序列及幅频响应如图5所示:Р图5滤波器单位冲击响应序列及幅频响应图Р高通滤波器频域波形图如图6所示: