m5high.m

来自「基于matlab7.0仿真编程 - 巴特沃茲低通及高通滤波器- Butterwo」· M 代码 · 共 12 行

clf;
fp=6500;fs=1500
wp=2*pi*fp;ws=2*pi*fs
[N,wn]=buttord(wp,ws,0.5,40,'s')
[b,a]=butter(N,wn,'high','s')
wa=0:(3*ws)/511:3*ws
h=freqs(b,a,wa)
H=zplane(h,1)
plot(wa/(2*pi),20*log10(abs(h)))
xlabel('Frequency, Hz');ylabel('Gain, dB')
title('Gain Response of Butterworth''s High-pass Filter')
disp('Transfer function=');disp(H);disp('Factorial Order N=');disp(N);