ep742.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理--原理、实现及应用》第七章 例7.4.3 程序：ep743.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 调用remez函数设计低通滤波器
Fs=10000;   %
f=[1500,2500];  %边界频率为模拟频率(Hz)
m=[1,0];
rp=1;rs=40;
dat1=(10^(rp/20)-1)/(10^(rp/20)+1);dat2=10^(-rs/20);
rip=[dat1,dat2];
[M,fo,mo,w]=remezord(f,m,rip,Fs);M=M+1; %边界频率为模拟频率(Hz)时必须加入采样频率Fs
hn=remez(M,fo,mo,w);
%以下为绘图检验部分
hw=fft(hn,512);
n=0:M;
subplot(2,2,1);stem(n,hn,'.');
xlabel('n');ylabel('h(n)');
w=0:511;w=2*w/512;
subplot(2,2,2);
plot(w,20*log10(abs(hw)));grid on;
axis([0,max(w)/2,-80,10]);
xlabel('ω/π');
ylabel('幅度(dB)')