[程序分享] matlab求K-L信息量值的代码 [推广有奖]

1关注
0粉丝

等待验证会员

已卖：3份资源

本科生

27%

还不是VIP/贵宾

威望: 0 级
论坛币: 27 个
通用积分: 2.0000
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 470 点
帖子: 52
精华: 0
在线时间: 89 小时
注册时间: 2019-2-9
最后登录: 2019-10-6

楼主

谢谢谢222 发表于 2019-3-4 23:06:54 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

clear all
close all
load('p1.txt');
load('p2.txt');
syms a1;
syms a2;
n=length(p1);%数据的个数
m=length(p2);%数据的个数
k1=0;k2=0;kL1=[];kL2=[];
a1=sum(p1);%存储p1数据的和
a2=sum(p2);%存储p2数据的和
p=p1./a1;
q=p2./a2;
%%%%滞后循环%%%%%%
for L=1:10
for t=L:n-L
k1=k1-p(t)*(log(p(t)/q(t+L)));
end
k1=abs(k1);
kL1=[kL1,k1];
end
%%%%%%先行循环%%%%%%%%%

double kL2;
for L=-10:-1
for t=abs(L)+1:n
k2=k2-p(t)*(log(p(t)/q(t+L)));
end
k2=abs(k2);
kL2=[kL2,k2];
end
kL=[kL2,kL1];
c=find(kL==min(min(kL)));%找出最小值位置L行c列对应的L值
minkL=kL(c);%K1中的最小值
请问已经求出最小的K-L信息量了，怎么用代码实现求出最小K-L信息量对应的滞后阶次呢，谢谢