shiyan4.m

来自「数字信号处理课程的有关实验」· M 代码 · 共 18 行

wp1=3000*pi;                  %通带截止角频率   
ws1=6000*pi;                  %阻带截止角频率
wsam1=30000*pi;               %采样角频率
fsam1=wsam1/(2*pi);           %采样截止频率
passrad=(wp1+ws1)/2/fsam1;    %截止频率
A=3.3;
wdelta=(ws1-wp1)/fsam1; 
N=ceil(2*pi/wdelta*A);         %滤波器的阶数
w1=hamming(N+1);               %用海明窗实现
L=N/2+1;
n=1:1:N+1;
hd=sin(passrad*(n-L))./(pi*(n-L));            %理想低通滤波器
if(N==ceil(N/2)*2)
   hd(L)=passrad/pi;
end
h1=hd.*w1';                            %加窗
[mag1,rad]=freqz(h1);
plot(rad,20*log10(abs(mag1)));         %频谱显示