hc861.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第八章例8.6.1程序hc861
% 由低通数字变换为高通数字频带的变换示例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
ws = [0.4*pi, 0.7*pi];                  % 数字阻带边缘频率
wp = [0.25*pi, 0.8*pi];                 % 数字通带边缘频率
Rp = 1;	                                % 通带波动(dB) 
As = 40;                                % 阻带波动(dB) 
Ripple = 10 ^ (-Rp/20);                 % 通带波动
Attn = 10 ^ (-As/20);                   % 阻带波动
% 计算切比雪夫-2型滤波器参数
[N,wn] = cheb2ord(wp/pi,ws/pi,Rp,As);
% 数字切比雪夫-2型带阻滤波器设计:
[b,a] = cheby2(N,As,ws/pi,'stop');		% 调用综合设计函数
% 级联型式实现:
[sos,G] = tf2sos(b,a)                   % 变换为级联结构
[H,w]=freqz(b,a);                       % 求频率特性
plot(w,abs(H))                          % 绘图
grid on
set(gcf,'color','w')                 % 置图形背景色为白