q728.m

这是matlab在电子信息课程中的应用

%《MATLAB在电子信息课程中的应用》(第三版)第七章例7.28程序q728
% 用remez函数设计FIR数字高通滤波器
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear; close all
f=[0,23/32,3/4,1];m=[0,0,1,1];
N1=30;W1=[1,1]; hn1=remez(N1,f,m,W1);		% 情况①
k=[0:1023]*2/1024;
Hw1=fft(hn1,1024); 
subplot(3,3,1)
plot(k,20*log10(abs(Hw1))); %求出其幅频特性
title('N1=30,W1=[1,1]')		% 标注
axis([0,1,-40,5]),grid on,pause 			% 只画出正半轴频谱
N2=30;W2=[1,5]; hn2=remez(N2,f,m,W2);		% 情况②
Hw2=fft(hn2,1024); subplot(3,3,2)
plot(k,20*log10(abs(Hw2)));%求出其幅频特性
axis([0,1,-40,5]),grid on,pause 			% 只画出正半轴频谱
title('N2=30,W2=[1,5]')		% 标注
N3=60;W3=[1,1]; hn3=remez(N3,f,m,W3);		% 情况③
Hw3=fft(hn3,1024); subplot(3,3,3)
plot(k,20*log10(abs(Hw3)));%求出其幅频特性
axis([0,1,-40,5]),grid on,pause 			% 只画出正半轴频谱
title('N3=60,W3=[1,1]')		% 标注