%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是為求取一列信號數據的Fourier級數展開系數
%   輸出數據中R(1,1)相當于a0，R(1,2)無意義，R(i,1)相當于a(i-1)，R(i,2)相當于b(i-1)
%      
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=fuliyejishuxishu(t,data);
T=max(t);
N=13;                                          %N為測點數目
for i=1:N;
    dataan=data.*cos(2*(i-1)*pi*data/T);
    databn=data.*sin(2*(i-1)*pi*data/T);
    R(i,1)=2/T*xinhaojifen(t,dataan);
    R(i,2)=2/T*xinhaojifen(t,databn);
end;