function [t,st]=lpf(f,sf,B)
% 本函數為通帶為B的低通濾波器
% 入口參數：f:頻率樣本向量
%           sf:頻譜向量
%           B:低通濾波器帶寬
% 出口參數：t:時間向量
%           st:輸出信號向量
df=f(2)-f(1);
T=1/df;
hf=zeros(1,length(f));
bf=[-floor(B/df):floor(B/df)]+floor(length(f)/2);
hf(bf)=1;
yf=hf.*sf;
[t,st]=F2T(f,yf);
st=real(st);