程序中涉及的自定義函數
定義傅里葉變換函數：
fuction[f,sf]=T2F(t,st)
dt=t(2)-t(1);
T=t(end);
df=1/T;
N=length(st);
f=-N/2*df:df:N/2*df-df;
sf=fft(st);
sf=T/N*fftshift(sf);

定義傅里葉反變換函數：
fuction[t,st]=F2T(f,sf)
df=f(2)-f(1);
Fmx=(f(end)-f(1)+df);
dt=1/Fmx;
T=dt*N;
%t=-T/2:dt:T/2-dt;
t=0:dt:T-dt;
sff=fftshift(sf);
st=Fmx*ifft(sff);


定義低通濾波函數：
fuction[t st]=lpf(f,sf,B);
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);


定義一個新序列將輸入的序列擴展成間隔為N-1個0的序列：

fuction[out]=sigexpand(d,M)
N=length(d);
out=zeros(M,N);
out(1, : )=d;
out=reshape(out,1,M*N);