function [M,m,df1,f]=T2F(m,ts,df,fs)
%------------------------輸入參數
%m:信號
%ts:系統時域采樣間隔
%df:所需的頻率分辨率
%fs:系統采樣頻率
%---------------------輸出(返回)參數
%M:傅里葉變換后的頻譜序列
%m: 輸入信號參與過傅里葉變換后對應的序列，需要注意的是，該序列與輸入信號m的區別，其長度是不一樣的，輸入的m長度不一定是2的整數次冪，而傅里葉變換要求輸入信號長度為2的整數次冪，故傅里葉變換前需對輸入的m信號進行補零，其長度有所增加，故輸出參數中的m為補零后的輸入信號，其長度與輸入參數m不一樣，但與M,f長度是一樣的，并且，其與時間序列t所對應的序列m(1:length(t))與輸入參數中的m是一致的。
%df1:返回的頻率分辨率
%f:與M相對應的頻率序列
[M,m,df1]=fftseq(m,ts,df);
f = [0:df1:df1*(length(m)-1)] -fs/2;    %頻率向量
M=M/fs;