exe060301_2.m

来自「Matlab程序为清华大学胡广书版信号处理的上机习题答案」· M 代码 · 共 30 行

%----------------------------------------------------
% exe060301_2.m
% to design Chebyshev high-pass DF with s=2/Ts*(z-1)/(z+1)
%----------------------------------------------------
clear all;

fp=400;fs=300;Fs=1000;
rp=3;rs=35;
% 
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
% Firstly to finish frequency prewarping ;
wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2);
[n,wn]=cheb1ord(wap,was,rp,rs,'s');
[z,p,k]=cheb1ap(n,rp);    %
[bp,ap]=zp2tf(z,p,k); %
[bs,as]=lp2hp(bp,ap,wap);   %

[bz,az]=bilinear(bs,as,Fs) %
[h,w]=freqz(bz,az,256,Fs);  %

figure;
plot(w,20*log10(abs(h)));grid on;
ylabel('High-pass DF: dB')

hp=unwrap(angle(h));
figure;
plot(w,hp);grid on;
ylabel('Phase of TranFunction')