q727.m

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

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

clear;close all
fc=1/4; fs=5/16; 			% 输入给定指标
Rp=3;As=60;Fs=2;
f=[fc,fs];m=[1,0]; 		% 计算remezord函数所需参数f,m,dev
dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];
[N,fo,mo,W]=remezord(f,m,dev,Fs);	% 确定remez函数所需参数
hn=remez(N,fo,mo,W);		% 调用remez函数进行设计
hw=fft(hn,512);			% 求设计出的滤波器频率特性
w=[0:511]*2/512;
plot(w,20*log10(abs(hw)));grid;
axis([0,max(w)/2,-90,5]);
xlabel('w/pi');ylabel('Magnitude(dB)')
line([0,0.4],[-3,-3]); %画线检验设计结果
line([1/4,1/4],[-90,5]);line([5/16,5/16],[-90,5]);