ep741.m

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

%《数字信号处理--原理、实现及应用》第七章 例7.4.2 程序：ep742.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 调用remez函数设计带阻滤波器
f=[0.2,0.35,0.65,0.8];
m=[1,0,1];
rp=1;rs=60;
dat1=(10^(rp/20)-1)/(10^(rp/20)+1);dat2=10^(-rs/20);
rip=[dat1,dat2,dat1];
[N,fo,mo,w]=remezord(f,m,rip);
hn=remez(N,fo,mo,w);
%以下为绘图检验部分
hw=fft(hn,512);
n=0:N;
subplot(2,2,1);stem(n,hn,'.');
xlabel('n');ylabel('h(n)');
grid on
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)')