srelevent.m

来自「该程序可以方便的求解出钢丝绳漏磁检测信号的自相关性」· M 代码 · 共 38 行

%装载采集的信号wirerope.mat
file5_17=['f:\9_5_19.xls'];
x5_17_1=xlsread(file5_17);
wirerope=x5_17_1';
save wirerope;
load wirerope;
%=============================
%将信号中第4500到第27000个采样点赋给s
indx=1:100;

s=wirerope(indx);
%=============================
%画出原始信号
%figure ;
%plot(s);
%title('原始信号');
%=============================
i=1;
a=s';
sum=0;
for k=1:100,
    sum=sum+a(k);
end
xave=sum/100;
sum2=0;
for k=1:100,
    sum2=sum2+(a(k)-xave)^2;
end
sum1=0;
for k=1:100,
    if (k+i)<=100
      sum1=sum1+(a(k)-xave)*(a(k+i)-xave);
    else 
      sum1=sum1+(a(k)-xave)*(a(k-i)-xave);
    end
end
  srele=sum1/sum2
%en1=log2(100)