ms25_1.m

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

%频率采样法设计低通
N=input('N=');                           %由使用者输入N
n=0:N-1;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);         %由脉冲序列求得频率特性
dw=2*pi/1000;                            %dw为频率分辨率，将0～2π分为1000份
Rp=-(min(db(1:fix(wc/dw)+1)))            %检验通带波动
ws=wc+2*pi/N;                            %确定ws的位置,比wc右移2*pi/N
As=-round(max(db(fix(ws/dw)+1:501)))     %检验最小阻带衰减
m=(N-1)/2;wa=[0:m-1]/m;                  %为作图求0-1π的样点数,建立对应的频率向量

subplot(2,1,1),plot(wa,A(1:m),'.-',w/pi,mag);
axis([0,1,-0.1,1.2]);title('理想幅频、样点序列及实际滤波器幅频响应');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
subplot(2,1,2),stem(n,h);
title('滤波器脉冲响应');
xlabel('n');ylabel('h(n)');