%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是為求取一列信號數據的自協方差函數
%      
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=zixiefangcha(data)
n=max(size(data));n=n/50;
for k=1:n-1;
    temp=0;
    for i=k+1:n;
        temp=temp+data(i)*data(i-k);
    end;
    R(k)=temp/(n-k);
end;