数字信号处理

R数字高通滤波器滤波后序列2频域'); Р%窗函数法(哈明窗)设计线性相位带阻FIR数字滤波?Рwlp=0.1*pi;Рwls=0.15*pi;Рwus=0.4*pi;Рwup=0.45*pi;Рbt=wls-wlp;РN0=ceil(6.2*pi/bt);РN=N0+mod(N0+1,2);%产生滤波函数序列长度Рwp=[(wlp+wls)/2/pi,(wus+wup)/2/pi];Рhn=fir1(N-1,wp,'stop',hamming(N));Р%傅里叶变换Рfh=fft(hn);Рfigure;Рsubplot(2,1,1),plot(abs(fh)),grid;Рaxis([0,N,0,1.5]);Рtitle('哈明窗,带阻FIR数字高通滤波器幅频特性');Рsubplot(2,1,2),plot(phase(fh)),grid;Рaxis([0,N,0,4]);Рtitle('哈明窗,带阻FIR数字高通滤波器相频特性'); Р%滤波效果Рfigure;Рsubplot(2,2,1),plot(filter(hn,1,y1)),grid;Рtitle('哈明窗,带阻FIR数字高通滤波器滤波后序列1时域');Рsubplot(2,2,2),plot(filter(hn,1,y2)),grid;Рtitle('哈明窗,带阻FIR数字高通滤波器滤波后序列2时域');Рbtft1=fft(filter(hn,1,y1));Рbtft2=fft(filter(hn,1,y2));Рsubplot(2,2,3),plot(abs(btft1)),grid;Рtitle('哈明窗,带阻FIR数字高通滤波器滤波后序列1频域'); Рsubplot(2,2,4),plot(abs(btft2)),grid;Рtitle('哈明窗,带阻FIR数字高通滤波器滤波后序列2频域');