%利用Hamming窗設(shè)計FIR數(shù)字低通濾波器
ap=3;as=30;
fp=30; fs=50;
Fs=200;
Wp=0.3*pi;Ws=0.5*pi;
%確定濾波器階數(shù)
f=(fs-fp)/Fs;
N=3.3/f;
M=N-1;
fprintf('濾波器階數(shù)=% .0f\n',M);
%理想低通截頻
Wc=(Wp+Ws)/2;
%窗函數(shù)的值
k=0:M;
Wk=0.54-0.46*cos(2*pi*k/M)
%濾波器的單位抽樣響應(yīng)
hd=sin((k-M/2)*Wc)/(pi*(k-M/2))
h=hd.*Wk;
subplot(221)
plot(h)
%設(shè)計FIR低通濾波器的增益響應(yīng)
omega=linspace(0,pi,512);
mag=freqz(h,[1],omega);
magdb=20*log(abs(mag));
subplot(222)
plot(omega/pi,magdb);
xlabel('normalized frequency');
ylabel('gain,dB');grid;
x=load('D:\學習\課程資料\BBB.TXT');
subplot(223);
plot(x(:,1));
title('真實信號');
y=zeros(1,512);
for t=33:512
    for n=1:32
     y(t)=y(t)+h(n).*x(t-n);
    end
end
subplot(224);
plot(y);
title('輸出信號');