comp2.m

来自「多变量的互信息的matlab的函数的源码」· M 代码 · 共 30 行

clear all
close all
clc
f=2;m=50;x=0;
for fs=1: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,1,1);
plot(h1,'b');
hold on
plot(h2,'y');
xlabel('Frequency');
ylabel('entropy(s/b,d/y)');
axis([1,50,0,5]);
subplot(2,1,2);
plot(MI,'r');
hold on
plot(ECC,'g');
xlabel('Frequency');
ylabel('MI/r,ECC/g');
axis([1,50,0,5]);