二、设计一个高通Chebyshow型数字滤波器
wp=100;ws=80;Fs=300;rp=1;rs=45; %数字滤波器的各项指标;
WP=100*2*pi; %把数字滤波器的频率特征转换成模拟滤波器的频率特征;WS=300*2*pi;
[N,Wn]=cheb2ord(WP,WS,rp,rs,'s'); %Chebyshev II型滤波器参数计算(模拟域);[Z,P,K]=cheb2ap(N,rs); %创建Chebyshev滤波器原型;
[A,B,C,D]=zp2ss(Z,P,K); %表达式从零极点增益形式转换成状态方程形式;[AA,BB,CC,DD]=lp2hp(A,B,C,D,Wn); %实现低通到高通滤波器类型的转换;[a,b,c,d]=bilinear(AA,BB,CC,DD,Fs); %采用双线性变换法,从模拟高通到数字高通;[P,Q]=ss2tf(a,b,c,d); %表达式从状态方程形形式转换成传输函数形式;figure(1);
freqz(P,Q); %绘出频率响应;
[H,W]=freqz(P,Q);
figure(2);
plot(W*Fs/(2*pi),abs(H));grid;
xlabel('频率/Hz');
ylabel('幅值');