ep651.m

来自「高西全2009年最新版<数字信号处理>源码程序」· M 代码 · 共 13 行

%《数字信号处理》第六章 例6.5.1 设计程序：ep651.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用双线性变换法设计DF
wpz=0.8;wsz=0.44;rp=3;rs=15;
[N,wc]=buttord(wpz,wsz,rp,rs);      %调用buttord和butter直接设计数字滤波器   
[Bz,Az]=butter(N,wc,'high');
wk=0:pi/512:pi;
Hz=freqz(Bz,Az,wk);
subplot(2,2,1);
plot(wk/pi,20*log10(abs(Hz)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-100,10]);