ms25_1_1.m

《数字信号处理实验(MATLAB版)》这

%频率采样法设计低通
for N=21:40:61                           %由使用者输入N
wc=0.4*pi;                               %输入截止频率
N1=fix(wc/(2*pi/N));                     %样点间隔为2*pi/N,N1为wc的样点数
N2=N-2*N1-1;                             %N2为阻带样点数
A=[ones(1,N1+1),zeros(1,N2),ones(1,N1)]; %建立符幅特性样本序列
theta=-pi*(N-1)/N*[0:N-1];               %建立相位特性样本序列 
Hk=A.*exp(j*theta);                      %建立频率特性样本序列
h=real(ifft(Hk))                        %由反变换求脉冲序列,去掉运算误差造成的虚部
[db,mag,pha,grd,w]=freqz_m(h,1);         %由脉冲序列求得频率特性
m=(N-1)/2;wa=[0:m-1]/m;                       %为作图求0-1π的样点数,建立对应的频率向量
subplot(2,1,1),plot(w/pi,mag),hold
axis([0,1,-0.1,1.2]);title('N=21和N=61时滤波器的幅频响应');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
subplot(2,1,2),plot(w/pi,db),hold
axis([0,1,-50,1.2]);title('N=21和N=61时滤波器的幅频响应(dB)');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
end