ms25_3.m

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

%频率采样法设计低通(2点过渡带)
N=61;T1=0.5886;T2=0.1065;                
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),T1,T2,zeros(1,N2-2),T2,T1,ones(1,N1-1)]; 
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份
wp=2*pi/N*fix(wc/(2*pi/N)-1);            %确定wp的位置
Rp=-(min(db(1:fix(wp/dw)+1)))            %检验通带波动
ws=wp+3*2*pi/N;         %确定ws的位置
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.1]);title('理想幅频、样点序列及实际滤波器幅频响应');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
set(gca,'XTickMode','manual','XTick',[0,wp/pi,wc/pi,ws/pi,1]);
set(gca,'YTickMode','manual','YTick',[0,T2,T1,1]);grid
subplot(2,1,2),plot(w/pi,db);
axis([0,1,-100,1.2]);title('实际滤波器幅频响应(dB)');
xlabel('频率(单位:\pi)');ylabel('G(dB)');
set(gca,'XTickMode','manual','XTick',[0,wp/pi,wc/pi,ws/pi,1]);
set(gca,'YTickMode','manual','YTick',[-60,-40,-20,0]);grid