%《MATLAB在電子信息課程中的應(yīng)用》(第三版)第七章例7.27程序q727
% 用remez函數(shù)設(shè)計FIR低通濾波器
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear;close all
fc=1/4; fs=5/16; 			% 輸入給定指標
Rp=3;As=60;Fs=2;
f=[fc,fs];m=[1,0]; 		% 計算remezord函數(shù)所需參數(shù)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函數(shù)所需參數(shù)
hn=remez(N,fo,mo,W);		% 調(diào)用remez函數(shù)進行設(shè)計
hw=fft(hn,512);			% 求設(shè)計出的濾波器頻率特性
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]); %畫線檢驗設(shè)計結(jié)果
line([1/4,1/4],[-90,5]);line([5/16,5/16],[-90,5]);