自适应干扰对消程序 matlab

%导入冲击信号和参考信号
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