nba.m

来自「利用窗函数法实现FIR滤波器的设计」· M 代码 · 共 34 行

Nwin1=10;   
Nwin2=20;
Nwin3=50;
Nwin4=100;
w1=boxcar(Nwin1);        %矩形窗
[y1,f1]=freqz(w1,1,512);   %频率特性
mag1=abs(y1);             %幅频特性
w2=boxcar(Nwin2);        %矩形窗
[y2,f2]=freqz(w2,1,512);   %频率特性
mag2=abs(y2);             %幅频特性
w3=boxcar(Nwin3);        %矩形窗
[y3,f3]=freqz(w3,1,512);   %频率特性
mag3=abs(y3);             %幅频特性
w4=boxcar(Nwin4);        %矩形窗
[y4,f4]=freqz(w4,1,512);   %频率特性
mag4=abs(y4);             %幅频特性
figure(1);
subplot(221);
plot (f1/pi,20*log10(mag1/max(mag1))); 
xlabel('归一化频率');ylabel('振幅/dB');
subplot(222);
plot (f2/pi,20*log10(mag2/max(mag2))); 
xlabel('归一化频率');ylabel('振幅/dB');
subplot(223);
plot (f3/pi,20*log10(mag3/max(mag3))); 
xlabel('归一化频率');ylabel('振幅/dB');
subplot(224);
plot (f4/pi,20*log10(mag4/max(mag4))); 
xlabel('归一化频率');ylabel('振幅/dB');
grid on;