my_filter.m

用matlab编写的关于设计滤波因子的源代码

%%%%%%%%%%%%%%%%%%%  fir滤波因子的设计  %%%%%%%%%%%%%%%
T=0.002;   % 采样间隔
f0=60;     % 截止频率
N=129;     % 采样点为（2*N+1）
for k=1:N-1
    h(N-k)=sin(-k*2*pi*f0*T)/(-pi*k);
    h(N)=0;
    h(k+N)=sin(-k*2*pi*f0*T)/(pi*k);
end
for i = 1:N 
    h(i) = h(i) * (0.54 - 0.46 * cos(2 * pi * i / (N - 1)));   % 用汉明窗加权
end

% %%%%%%%%%%%%%%%% 对滤波因子加窗函数 %%%%%%%%%%%%%%%
% for k=1:2*N-1    
%     if k<=N-64||k>=N+64
%         h(k)=0;
%     end
% end
%%%%%%%%%%%%%%%%
H=abs(fft(h)).^2;
figure(1);
plot(h);