fig741.m

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

%《数字信号处理--原理、实现及应用》第六章 图7.4.1绘图程序:fig741.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月

clear;close all
fc=1/2; fs=11/20;
Rp=2;As=20;Fs=2;
f=[fc,fs];m=[1,0];
dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];
[N,fo,mo,wo]=remezord(f,m,dev,Fs);
N=N+1;
hn=remez(N,fo,mo,wo);
hw=fft(hn,512);
w=0:511;w=2*w/512;
subplot(2,2,1);
plot(w,20*log10(abs(hw)));grid;
axis([0,max(w)/2,-50,10]);
xlabel('ω/π');ylabel('20lg(|Hg(ω)|)')
%line([0,0.4],[hw(1),hw(1)]); %画线检验设计结果
line([fc,fc],[-50,10]);line([fs,fs],[-50,10]);
subplot(2,2,2);
plot(w,abs(hw));grid;
axis([0,max(w)/2,0,1.2]);
xlabel('ω/π');ylabel('|Hg(ω)|')
line([0,0.6],[hw(1),hw(1)]); %画线检验设计结果
line([0,0.6],[1-0.1146,1-0.1146]); %画线检验设计结果
line([0,1],[0.1,0.1]);
line([fc,fc],[0,1.2]);line([fs,fs],[0,1.2]);