改进后的互信息法

%确定仿真数据spcv2（空簧失气），速度200，1通道（车体前部加速度）数据的互信息
clear
Data1=load （‘t.mat‘）;
Data=Data1（1:20487）;
c=2048*0;
h=c+1:c+2000;
x=Data（h）;%取2048个数据
for t=0:20
a=t+1:length（x）+t;
y=Data（c+a）;%x延迟个数据后的y数据
drelx=drawelem2（x）;%对序列x进行不同元素的提取并按从小到大进行排序
drely=drawelem2（y）;%对序列y进行不同元素的提取并按从小到大进行排序
stx=statx2（x）;
sty=statx2（y）;
MI（t+1）=mutualinformation22（xdrelxstxydrelysty）;
end
MI
for i=1:19
    if MI（i+1）        t=i;
        break
    end
end
t
figure（1）
plot（0:20MI‘--b.‘）
grid
xlabel（‘延迟时间‘）;
ylabel（‘互信息‘）;

%调用函数1
function drelx=drawelem2（x）
%提取数组中不同元素，按从小到大排序 
%es=[1 2 3 3];%需提取的矩阵
es=sort（x）;
%es=fix（esx*100）/100;
j=1;%设置初值
for i=1:length（es）-1%矩阵各元素