function [t,omg,FT,IFT]=prefourier(Trg,N,OMGrg,K)
%用矩陣左乘來(lái)實(shí)現(xiàn)傅里葉變換前的準(zhǔn)備工作
%輸入：TRG，時(shí)域起止范圍，N，時(shí)域抽樣點(diǎn)數(shù)，OMGrg，頻域起止范圍和頻域抽樣點(diǎn)數(shù)
%輸出： 時(shí)域和頻域抽樣點(diǎn)數(shù)，傅里葉變換和逆變換
F=Trg(2)-Trg(1); %時(shí)域范圍
t=linspace(Trg(1),Trg(2)-T/N,N); %生成抽樣時(shí)間點(diǎn)
OMG=OMGrg(2)-OMGrg(1);          %頻域范圍
omg=linsapce(OMGrg(1),OMGrg(2)-OMG/K,K); %生成抽樣頻率點(diǎn)
FT=F/N*exp(-j*kron(omg,t.'));            %構(gòu)造帶有系數(shù)的傅里葉變換矩陣
IFT=OMG/2/pi/K*exp(j*kron(t,omg.'));     %構(gòu)造帶有系數(shù)的的傅里葉逆變換矩陣