fir.m

FIR滤波器

%利用Hamming窗设计FIR数字低通滤波器
ap=3;as=30;
fp=30; fs=50;
Fs=200;
Wp=0.3*pi;Ws=0.5*pi;
%确定滤波器阶数
f=(fs-fp)/Fs;
N=3.3/f;
M=N-1;
fprintf('滤波器阶数=% .0f\n',M);
%理想低通截频
Wc=(Wp+Ws)/2;
%窗函数的值
k=0:M;
Wk=0.54-0.46*cos(2*pi*k/M)
%滤波器的单位抽样响应
hd=sin((k-M/2)*Wc)/(pi*(k-M/2))
h=hd.*Wk;
subplot(221)
plot(h)
%设计FIR低通滤波器的增益响应
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('输出信号');