|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
3 X9 M1 @) a/ x% g7 s* r通用MATLAB代码
, I$ o R# a' H: N% g% ?5 @/ }主函数main.m3 e5 C, g/ K) `3 Z4 z$ Y! t# P
' a$ O- ]0 I5 N+ ~2 ~
- clc
- u1 = rand(4,1);
- u2 = [2;32;6666;5];
- wind_size = size(u1,1);
- mi = calmi(u1, u2, wind_size);
! a' I/ h4 o9 Q- w5 Q8 m$ e. B+ p
) i' ]8 ]9 J, d7 O; t$ y" O
: K) B" v* ]- ?, I2 r5 [calmi.m
2 H1 u- D; Q3 \! N: u+ j5 z$ z' U
- %计算两列向量之间的互信息
- %u1：输入计算的向量1
- %u2：输入计算的向量2
- %wind_size：向量的长度
- function mi = calmi(u1, u2, wind_size)
- x = [u1, u2];
- n = wind_size;
- [xrow, xcol] = size(x);
- bin = zeros(xrow,xcol);
- pmf = zeros(n, 2);
- for i = 1:2
- minx = min(x(:,i));
- maxx = max(x(:,i));
- binwidth = (maxx - minx) / n;
- edges = minx + binwidth*(0:n);
- histcEdges = [-Inf edges(2:end-1) Inf];
- [occur,bin(:,i)] = histc(x(:,i),histcEdges,1); %通过直方图方式计算单个向量的直方图分布
- pmf(:,i) = occur(1:n)./xrow;
- end
- %计算u1和u2的联合概率密度
- jointOccur = accumarray(bin,1,[n,n]); %（xi，yi）两个数据同时落入n*n等分方格中的数量即为联合概率密度
- jointPmf = jointOccur./xrow;
- Hx = -(pmf(:,1))'*log2(pmf(:,1)+eps);
- Hy = -(pmf(:,2))'*log2(pmf(:,2)+eps);
- Hxy = -(jointPmf(：）)'*log2(jointPmf(：）+eps);
- MI = Hx+Hy-Hxy;
- mi = MI/sqrt(Hx*Hy);8 m$ h0 B0 |0 A+ \9 `8 ]$ B
. d, F+ O: Q& d/ t; J5 v& g7 ?
|
|
/1 
发表于 2020-11-11 14:52
收藏
淘帖
支持！
反对！