muin.asv

计算两个离散时间序列变量间的信息熵和互信息

clear all
close all
clc
f=2;m=50;x=0;
for fs=1:5:50
    x=x+1;
    for n=1:10000
        y(n)=sin(2*pi*f*n/fs);
        u(n)=0.02*sin(2*pi*f*n/fs);
    end
s=y;d=y+m*u;
h1(x)=entropy(s);h2(x)=entropy(d);
MI(x)=mutualinfo(s,d);
ECC(x)=2*MI(x)/(h1(x)+h2(x));
end
disp(h1);disp(h2);disp(MI);disp(ECC);
subplot(2,2,1);plot(h1);
xlabel('Frequency');
ylabel('entropy');
subplot(2,2,2);plot(h2);
xlabel('Frequency');
ylabel('entropy');
subplot(2,2,3);plot(MI);xlabel('Frequency');
ylabel('entropy');
subplot(2,2,4);plot(ECC);