IIR高通、带通和带阻数字滤波器设计 源代码

Fs=10000; Wp=[1500,4000]*2/Fs; Ws=[2000,3000]*2/Fs; Ap=-2;As=-13; Nn=256; [N,Wn]=buttord(wp,ws,Ap,As); %计算巴特沃斯滤波器阶次和截止频率 [B,A]=butter(N,Wn,'stop'); %频率变换法设计巴特沃斯带通滤波器 [h,f]=freqz(B,A,256,Fs); %数字滤波器响应 subplot(211); plot(f,abs(h)); axis([0,5000,0,1]); xlabel('频率/Hz'); ylabel('振幅/dB'); Title('带阻巴特沃斯滤波器') grid on; subplot(212); plot(f,180/pi*unwrap(angle(h))); xlabel('频率/KHz'); ylabel('相位/dB'); grid on;