ms25_1_2.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 16 行

%频率采样法设计低通(0-2*pi)
N=21;                           %由使用者输入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];               %建立相位特性样本序列 
wa=[0:N-1]/N*2;                          %为作图建立对应的频率向量
subplot(2,1,1),plot(wa,A,'.-');
axis([0,2,-0.1,1.2]);title('理想幅频响应及样点序列(N=21)');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
subplot(2,1,2),plot(wa,theta,'.-');
axis([0,2,-90,1]);title('理想相频响应及样点序列(N=21)');
xlabel('频率(单位:\pi)');ylabel('\phi(\omega)');