数字信号处理实验报告实验四.doc

--------------------------校验:_____________-----------------------日期:_____________数字信号处理实验报告实验四数字信号处理实验报告实验四:FIR滤波器设计实验班级:13080241姓名:徐可学号:1308024121实验四FIR滤波器设计实验一、。、,给定抽样频率为Ωs=3π×104(rad/s),通带截止频率为Ωp=3π×103(rad/s),阻带起始频率为Ωst=6π×103(rad/s),阻带最小衰减大于50dB。,观察其频率响应函数有什么变化三、实验设备计算机,Matlab软件四、;;。五、=2*pi/(30000*pi);rs=50;wp=3000*pi*T;ws=6000*pi*T;wc=(wp+ws)/2/pi;beta=; %ceil(x)为大于等于x的最小整数hdn=kaiser(N,beta); %产生长度为N的凯塞窗函数hn=fir1(N-1,wc,'low',kaiser(N,beta)); %用凯塞窗函数设计低通滤波器subplot(121);stem(0:N-1,hn,'b.'); %绘制滤波器时域波形;xlabel('n');ylabel('h(n)');gridon;omega=linspace(0,pi,512);mag=freqz(hn,[1],omega); %计算单位抽样响应的频率响应magdb=20*log10(abs(mag)); subplot(122);plot(omega/pi,magdb,'b'); %绘制对数幅度特性曲线axis([01-805]);xlabel('\omega/\pi');ylabel('20lg|H(e^j^\omega)|');gridon;=2*pi/(30000*pi);wp=3000*pi*T/pi;ws=6000*pi*T/pi;wc=(wp+ws)/2;delta_w=ws-wp;hn=fir1(N-1,wc,'low',hamming(N));omega=linspace(0,pi,512);mag=freqz(hn,1,omega);magdb=20*log10(abs(mag));subplot(121),stem([0:N-1],hn,'.'xlabel('n');ylabel('h(n)');title('单位抽样响应');subplot(122),plot(omega/pi,magdb);gridon;axis([01-10010]);xlabel('\omega/\pi');ylabel('dB');title('幅值频率响应');=2*pi/(30000*pi);wp=3000*pi*T/pi;ws=6000*pi*T/pi;wc=(wp+ws)/2;delta_w=ws-wp;hn=fir1(N