zixiefangcha.m

来自「自己编写的比较全面的故障诊断matlab函数程序」· M 代码 · 共 20 行

%**************************************************************************************************************************************************************
%   本程序为函数功能，目的是为求取一列信号数据的自协方差函数
%      
%   中国北方发动机研究所试验测试中心  苏志霄   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;