%導入沖擊信號和參考信號
clear;
A=load（‘E:\碩士畢業設計\實驗數據\2012.6.21?未進梁\10g\1000HZ\大油門轉\模擬采樣數據.lvm‘）;
len1=20;?%沖擊梁相關運算子序列長度len1
len2=10;?%參考梁相關運算子序列超前和滯后沖擊梁各len2
n1=50;?%截取總數據所需子序列len1的段數，num1>len2
num1=201;?%沖擊梁總數據截取起始點
num2=num1+n1*len1-1;?%沖擊梁總數據截取終止點
num3=num1-len2;?%參考梁總數據截取起始點?
num4=num2+len2;?%參考梁總數據截取終止點
B=A（num1:num22）;?%取沖擊信號列總數據
C=A（num3:num43）;?%取參考信號列總數據

%從沖擊梁總數據中需截取n1段長度為len1的子序列
X=0;?%X為對消后總數據
for?p=0:n1-1
????M=B（len1*p+1:len1*p+1+len1-1）;?%沖擊梁相關運算子序列len1個數據
????N=C（len1*p+1:len1*p+1+2*len2+len1-1）;?%參考梁相關運算子序列（2*len2+len1）個數據
????
%以len1為單位進行自適應對消，未平移和每次平移后的相關序列
????for?i=1:2*len2+1?%正好進行（2*len2+1）次相關運算
??????s=0;?%s為每次相關運算得到的相關系數
????????fo