真实熵预测图书借阅

最新推荐文章于 2022-12-09 20:24:11 发布
最新推荐文章于 2022-12-09 20:24:11 发布
阅读量253 收藏 1
点赞数
分类专栏： 算法 文章标签： python 真实熵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_48615832/article/details/119617569
版权

计算真实熵

计算真实熵的逻辑可参考我的上一篇文章：真实熵(the actual entropy)
代码参考：真实熵与人类行为可预测性

#真实熵
import pandas as pd
import math
# 横向最多显示多少个字符， 一般80不适合横向的屏幕，平时多用200
pd.set_option('display.width', 500)
# 显示所有列
pd.set_option('display.max_columns', None)

#图书类别
Book_category = ['A','B','C','D','E','F','G','H','I','J','K','N','O','P','Q','R','S','T','U','V','X','Z']
# # 读取数据
# Person_data = pd.read_excel(r'Person_8748847336.xlsx').iloc[0: , 1:]
LOAN_DATE = ['2013/1/1', '2013/1/1', '2013/1/1', '2013/1/1', '2013/1/1', '2013/2/20', '2013/4/19', '2013/4/19', '2013/4/19', '2013/4/24', '2013/4/24', '2013/4/24', '2013/4/24', '2013/4/24', '2013/4/24', '2013/4/24', '2013/4/27', '2013/5/8', '2013/5/15', '2013/5/15', '2013/5/17', '2013/6/6', '2013/6/8', '2013/6/8', '2013/10/9', '2013/10/29', '2013/10/29', '2013/11/7', '2013/11/7', '2013/11/14', '2013/11/14', '2013/11/14', '2013/11/25', '2013/12/2', '2013/12/2', '2013/12/4', '2013/12/6', '2013/12/11', '2013/12/11', '2013/12/18', '2013/12/18', '2014/1/1', '2014/1/2', '2014/2/17', '2014/2/17', '2014/3/13', '2014/3/13', '2014/4/15', '2014/4/29', '2014/4/29', '2014/4/29', '2014/5/5', '2014/9/23', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/10/10', '2014/12/2', '2014/12/2', '2014/12/2', '2014/12/4', '2014/12/4', '2015/1/14', '2015/2/28', '2015/2/28', '2015/3/5', '2015/3/5', '2015/3/5', '2015/3/5', '2015/3/11', '2015/3/11', '2015/3/11', '2015/3/12', '2015/3/12', '2015/3/12', '2015/3/12', '2015/3/12', '2015/3/19', '2015/3/19']
time_number = ['2013/1/1', '2013/2/20', '2013/4/19', '2013/4/24', '2013/4/27', '2013/5/8', '2013/5/15', '2013/5/17', '2013/6/6', '2013/6/8', '2013/10/9', '2013/10/29', '2013/11/7', '2013/11/14', '2013/11/25', '2013/12/2', '2013/12/4', '2013/12/6', '2013/12/11', '2013/12/18', '2014/1/1', '2014/1/2', '2014/2/17', '2014/3/13', '2014/4/15', '2014/4/29', '2014/5/5', '2014/9/23', '2014/10/10', '2014/12/2', '2014/12/4', '2015/1/14', '2015/2/28', '2015/3/5', '2015/3/11', '2015/3/12', '2015/3/19']
ITEM_CALLNO = ['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'J', 'D']

N = len(ITEM_CALLNO)    #样本数量
t_n = len(time_number)    #次数

def contains(small, big):
    for i in range(len(big)-len(small)+1):
        if big[i:i+len(small)] == small:
            return True
    return False

def actual_entropy(l):
    n = len(l)
    sequence = []
    sum_A = 0
    for i in range(0, n):
        for j in range(i+1, n+1):
            s = l[i:j]
            if not contains(list(s), sequence): # s is not contained in previous sequence
                print(f'i={i + 1}， s：{list(s)}， Λ{i+1}={len(list(s))}， 序列：{sequence}')
                sum_A += len(s)    #sum_Λ
                sequence.append(l[i])
                break
            else:
                if j == n:
                   print(f'i={i + 1}， s：不存在这个序列， Λ{i + 1}={n - (i + 1) + 2}， 序列：{sequence}')
                   sum_A += n - (i + 1) + 2
                   sequence.append(l[i])

    ae = 1 / (sum_A / n ) * math.log(n)
    print("真实熵：",ae)
    return ae

actual_entropy(ITEM_CALLNO)

结果：

i=1， s：['H']， Λ1=1， 序列：[]
i=2， s：['H', 'H']， Λ2=2， 序列：['H']
i=3， s：['H', 'H', 'H']， Λ3=3， 序列：['H', 'H']
i=4， s：['H', 'H', 'D']， Λ4=3， 序列：['H', 'H', 'H']
i=5， s：['H', 'D']， Λ5=2， 序列：['H', 'H', 'H', 'H']
i=6， s：['D']， Λ6=1， 序列：['H', 'H', 'H', 'H', 'H']
i=7， s：['F']， Λ7=1， 序列：['H', 'H', 'H', 'H', 'H', 'D']
i=8， s：['F', 'F']， Λ8=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F']
i=9， s：['F', 'D']， Λ9=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F']
i=10， s：['D', 'E']， Λ10=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F']
i=11， s：['E']， Λ11=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D']
i=12， s：['E', 'D']， Λ12=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E']
i=13， s：['D', 'D']， Λ13=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E']
i=14， s：['D', 'D']， Λ14=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D']
i=15， s：['D', 'F', 'F', 'H']， Λ15=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D']
i=16， s：['F', 'F', 'H']， Λ16=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D']
i=17， s：['F', 'H']， Λ17=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F']
i=18， s：['H', 'B']， Λ18=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F']
i=19， s：['B']， Λ19=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H']
i=20， s：['B', 'B']， Λ20=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B']
i=21， s：['B', 'C']， Λ21=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B']
i=22， s：['C']， Λ22=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B']
i=23， s：['D', 'D', 'I']， Λ23=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C']
i=24， s：['D', 'I']， Λ24=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D']
i=25， s：['I']， Λ25=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D']
i=26， s：['F', 'I']， Λ26=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I']
i=27， s：['I', 'F', 'D']， Λ27=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F']
i=28， s：['F', 'D', 'I']， Λ28=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I']
i=29， s：['D', 'I', 'F', 'F']， Λ29=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F']
i=30， s：['I', 'F', 'F']， Λ30=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D']
i=31， s：['F', 'F', 'D', 'F']， Λ31=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I']
i=32， s：['F', 'D', 'F']， Λ32=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F']
i=33， s：['D', 'F', 'F', 'H', 'F']， Λ33=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F']
i=34， s：['F', 'F', 'H', 'F']， Λ34=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D']
i=35， s：['F', 'H', 'F']， Λ35=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F']
i=36， s：['H', 'F']， Λ36=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F']
i=37， s：['F', 'F', 'F', 'F']， Λ37=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H']
i=38， s：['F', 'F', 'F', 'D', 'F']， Λ38=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F']
i=39， s：['F', 'F', 'D', 'F', 'D']， Λ39=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F']
i=40， s：['F', 'D', 'F', 'D']， Λ40=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F']
i=41， s：['D', 'F', 'D']， Λ41=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F']
i=42， s：['F', 'D', 'D']， Λ42=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D']
i=43， s：['D', 'D', 'D', 'K']， Λ43=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F']
i=44， s：['D', 'D', 'K']， Λ44=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D']
i=45， s：['D', 'K']， Λ45=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D']
i=46， s：['K']， Λ46=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D']
i=47， s：['K', 'D']， Λ47=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K']
i=48， s：['D', 'I', 'I']， Λ48=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K']
i=49， s：['I', 'I']， Λ49=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D']
i=50， s：['I', 'I']， Λ50=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I']
i=51， s：['I', 'D']， Λ51=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I']
i=52， s：['D', 'H']， Λ52=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I']
i=53， s：['H', 'D', 'D']， Λ53=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D']
i=54， s：['D', 'D', 'D', 'D']， Λ54=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H']
i=55， s：['D', 'D', 'D', 'D']， Λ55=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D']
i=56， s：['D', 'D', 'D', 'D']， Λ56=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D']
i=57， s：['D', 'D', 'D', 'D']， Λ57=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D']
i=58， s：['D', 'D', 'D', 'D', 'D']， Λ58=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D']
i=59， s：['D', 'D', 'D', 'D', 'D', 'F']， Λ59=6， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D']
i=60， s：['D', 'D', 'D', 'D', 'F']， Λ60=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D']
i=61， s：['D', 'D', 'D', 'F', 'K']， Λ61=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=62， s：['D', 'D', 'F', 'K']， Λ62=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=63， s：['D', 'F', 'K']， Λ63=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=64， s：['F', 'K']， Λ64=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=65， s：['K', 'D', 'D']， Λ65=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F']
i=66， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ66=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K']
i=67， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ67=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D']
i=68， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ68=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D']
i=69， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ69=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D']
i=70， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ70=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D']
i=71， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ71=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D']
i=72， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ72=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D']
i=73， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']， Λ73=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=74， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'J']， Λ74=11， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=75， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'J']， Λ75=10， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=76， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'J']， Λ76=9， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=77， s：['D', 'D', 'D', 'D', 'D', 'D', 'D', 'J']， Λ77=8， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=78， s：['D', 'D', 'D', 'D', 'D', 'D', 'J']， Λ78=7， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=79， s：['D', 'D', 'D', 'D', 'D', 'J']， Λ79=6， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=80， s：['D', 'D', 'D', 'D', 'J']， Λ80=5， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=81， s：['D', 'D', 'D', 'J']， Λ81=4， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=82， s：['D', 'D', 'J']， Λ82=3， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=83， s：['D', 'J']， Λ83=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=84， s：['J']， Λ84=1， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D']
i=85， s：不存在这个序列， Λ85=2， 序列：['H', 'H', 'H', 'H', 'H', 'D', 'F', 'F', 'F', 'D', 'E', 'E', 'D', 'D', 'D', 'F', 'F', 'H', 'B', 'B', 'B', 'C', 'D', 'D', 'I', 'F', 'I', 'F', 'D', 'I', 'F', 'F', 'D', 'F', 'F', 'H', 'F', 'F', 'F', 'F', 'D', 'F', 'D', 'D', 'D', 'K', 'K', 'D', 'I', 'I', 'I', 'D', 'H', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'F', 'K', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'D', 'J']
真实熵： 1.1106628141225792
下一次借书的种类数： 2.1594483556205337

计算可预测性

上面算得真实熵S=1.1106628141225792，用户可能借阅的图书种类有N=22种
计算方法

在这里插入图片描述
计算代码：

import mpmath
N = 22
S = 1.1106628141225792

def getPredictability(N, S):
    f = lambda x: (((1-x)/(N-1)) **(1-x))* x**x - 2**(-S)
    root = mpmath.findroot(f, 1)
    return float(root.real)

result = getPredictability(N, S)
print(result)

计算结果：

0.8724859773218735

因此最大可预测值为87.25%

最大可预测性与真实熵S的关系

import mpmath
import numpy as np
import matplotlib.pyplot as pl

def getPredictability(N, S):
    f = lambda x: (((1-x)/(N-1)) **(1-x))* x**x - 2**(-S)
    root = mpmath.findroot(f, 1)
    return float(root.real)

N  = 22
slist = np.arange(0, 1.6 ,0.1)
plist = []
for S in slist:
    p = getPredictability(N, S)
    plist.append(p)
print(slist)
print(plist)

pl.plot(slist, plist, 'g-o')
pl.xlabel('$S$', fontsize = 20)
pl.ylabel('$\Pi{max}$ ', fontsize = 20)
pl.show()

结果：

[0.  0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.  1.1 1.2 1.3 1.4 1.5]
[1.0, 0.9922076523891302, 0.9828942537046568, 0.9727602729432984, 0.9619949144642034, 0.9506950566016172, 0.9389200034509394, 0.9267095019884533, 0.9140916051169893, 0.9010866935833577, 0.8877097567024888, 0.8739717816424107, 0.8598806436493308, 0.8454416965880155, 0.8306581726777373, 0.8155314542381094]

在这里插入图片描述

最大可预测性与访问目标种类数量N的关系

import mpmath
import numpy as np
import pandas as pd
import matplotlib.pyplot as pl

def getPredictability(N, S):
    f = lambda x: (((1-x)/(N-1)) **(1-x))* x**x - 2**(-S)
    root = mpmath.findroot(f, 1)
    return float(root.real)

nlist  = np.arange(100, 1000, 200)
slist = np.arange(0, 1.1 ,0.1)
nsplist = []
for N in nlist:
    for S in slist:
        p = getPredictability(N, S)
        nsplist.append([N, S, p])

df = pd.DataFrame(nsplist, columns = ['N', 'S', 'P'])
print(df)
groups = df.groupby('N')
for name, group in groups:
    pl.plot(group.S, group.P, label = "N = "+ str(name), marker='o', linestyle='-')
pl.legend()
pl.xlabel('$S$', fontsize = 20)
pl.ylabel('$\Pi{max}$ ', fontsize = 20)
pl.show()

结果：

      N    S         P
0   100  0.0  1.000000
1   100  0.1  0.993476
2   100  0.2  0.985928
3   100  0.3  0.977866
4   100  0.4  0.969422
5   100  0.5  0.960665
6   100  0.6  0.951639
7   100  0.7  0.942371
8   100  0.8  0.932884
9   100  0.9  0.923194
10  100  1.0  0.913312
11  300  0.0  1.000000
12  300  0.1  0.994145
13  300  0.2  0.987482
14  300  0.3  0.980427
15  300  0.4  0.973087
16  300  0.5  0.965516
17  300  0.6  0.957748
18  300  0.7  0.949807
19  300  0.8  0.941710
20  300  0.9  0.933468
21  300  1.0  0.925094
22  500  0.0  1.000000
23  500  0.1  0.994408
24  500  0.2  0.988087
25  500  0.3  0.981415
26  500  0.4  0.974489
27  500  0.5  0.967361
28  500  0.6  0.960060
29  500  0.7  0.952607
30  500  0.8  0.945017
31  500  0.9  0.937303
32  500  1.0  0.929474
33  700  0.0  1.000000
34  700  0.1  0.994569
35  700  0.2  0.988452
36  700  0.3  0.982009
37  700  0.4  0.975332
38  700  0.5  0.968466
39  700  0.6  0.961441
40  700  0.7  0.954276
41  700  0.8  0.946986
42  700  0.9  0.939582
43  700  1.0  0.932072
44  900  0.0  1.000000
45  900  0.1  0.994683
46  900  0.2  0.988710
47  900  0.3  0.982428
48  900  0.4  0.975923
49  900  0.5  0.969241
50  900  0.6  0.962408
51  900  0.7  0.955444
52  900  0.8  0.948362
53  900  0.9  0.941172
54  900  1.0  0.933884

在这里插入图片描述

结论

1、真实熵越大，可预测性越小；
2、在熵不变的条件下，去过的地方（种类）越多，可预测性越强。

Fo*(Bi)
关注
专栏目录
4.2 图书借阅系统数据库设计 --MySQL
天罡gg的专栏
04-11 3万+
大家好，我是天罡gg，一个有十多年丰富经验的高级架构师，参与过很多系统的数据库设计，在数据库设计方面有相当丰富的经验。正赶上这篇实战专栏的数据库设计，所以今天让我们来一起做一下《图书借阅系统的数据库设计》，一篇既有理论知识，又能实战落地的数据库设计！本文从【需求分析】开始，到【概念结构设计】的画ER图, 再到【逻辑结构设计】的ER图转关系模型,最后到【物理结构设计】的MySQL表设计, 再加上【具体业务分析】的SQL实现，你是不是已经怀疑真有这么齐全吗？你没看错，全文HTML统计1万3千多字。
真实熵求解matlab代码
11-27
求解真实熵actual_entropy，即Lempel-Ziv算法所求出的熵。文件里包含两个程序，pcactentory是主程序，containarray是检测序列是否重复的程序。还需要说明，程序里没有计算第一个元素的长度，如果需要再最后k+1即可
真实熵(the actual entropy)
Fo*(Bi)的博客
03-10 1016
真实熵参考文献【Limits of Predictability in Human Mobility】 上文的中文版可以看这篇文章： 熵与空间位置预测 熵的可预测性定义 公式2.6为熵的定义，根据链式法则，可以将熵 S 进一步表示为如下形式： 随机熵、时间无关熵和真实熵 文中的【43】就是【Limits of Predictability in Human Mobility】 对每个用户来说，随机熵包含的信息量最小，仅包含用户使用基站的数量信息，所以Sirand最大，与时间无关的熵Siunc包含
随机熵、时间无关熵、真实熵理论+代码
m0_53961910的博客
12-09 787
最近看了一些和熵值有关的论文，在此总结一下。主要包括随机熵、时间无关熵（信息熵）、真实熵，并附加了他们对应的代码实现。
熵（entropy）、交叉熵（cross-entropy）
Yue_Zengying的博客
05-21 1658
原文：一文搞懂熵(Entropy),交叉熵(Cross-Entropy) 一、熵 1、混乱程度，不确定性，信息量？ 不同的人对熵有不同的解释：混乱程度，不确定性，惊奇程度，不可预测性，信息量等等，面对如此多的解释，第一次接触时难免困惑。本文第一部分，让我们先一起搞明白 熵 究竟是什么？ 信息论中熵的概念首次被香农提出，目的是寻找一种高效/无损地编码信息的方法：以编码后数据的平均长度来衡量高效性，平均长度越小越高效；同时还需满足“无损”的条件，即编码后不能有原始信息的丢失。这样，香农提出了熵的定义：无损编码事
交叉熵、相对熵、KL散度、JS散度
Oliverkehl的专栏
03-25 1020
以前没有从另一个角度理解过熵：用最小的bit位数描述一个分布，这个bit数就是熵 交叉熵就是用分布Q(x)拟合真实分布P(x)，最常见的就是LR的Loss 相对熵就是交叉熵减去真实分布的熵 KL散度就是相对熵，但是KL散度不具有对称性 JS散度是 JSD(P||Q)=1/2*KL(P||(P+Q)/2) +1/2*KL(Q||(P+Q)/2) 可见JS散度是具有对称性的 ...
c语言实现图书借阅系统
YSL
06-12 1万+
#include <stdio.h> #include <string.h> #include <stdlib.h> #include <conio.h> struct book //定义图书结构体 { int num; // 书号 char name[20]; //书名 char publish[20]; // 出版社 s...
html图书借阅源码,图书借阅管理系统代码图书管理系统源代码
weixin_36280317的博客
07-08 4021
《图书借阅管理系统代码图书管理系统源代码》由会员分享，可在线阅读，更多相关《图书借阅管理系统代码图书管理系统源代码(4页珍藏版)》请在人人文库网上搜索。1、include#include #include #define num 7 /全局定义void showmain() /显示菜单printf(nn 欢迎使用图书借阅管理系统nn);printf( 1、查询所有图书nn);printf( 2、归...
基于数据挖掘的高校图书馆图书借阅流量预测.pdf
07-14
在介绍《基于数据挖掘的高校图书馆图书借阅流量预测》这一研究文献中的关键知识点之前，需要明确，这篇文章是由陈桂菊撰写，发表于《微型电脑应用》期刊的2020年第36卷第6期，主要研究了如何通过数据挖掘技术提升...
图书借阅管理系统详细需求分析.doc
11-21
图书借阅管理系统是一种用于图书馆日常运营的信息管理系统，旨在提高图书借阅、归还、续借、预约等操作的效率，并方便管理员对馆内资源进行有效管理和监控。以下是对该系统详细需求的分析： 1 **背景** 图书借阅...
关于图书借阅管理系统的设计思考_图书借阅管理_
10-03
图书借阅管理系统是图书馆信息化建设的核心组成部分，它旨在提高图书借阅、归还、续借等操作的效率，同时提供全面的图书管理和用户服务。在当前数字化时代，图书借阅管理系统已经成为图书馆不可或缺的一部分，它通过...
JAVA图书借阅管理系统源代码+数据库，软件工程课程设计
最新发布
02-08
使用Java开发的图书管理系统，读者可以注册登录，登录时会判断账号类型再分别跳到各自对应的页面，读者可以查找，借阅，还书，查看历史借阅记录，修改个人资料，密码等，管理员可以修改添加图书，修改添加读者，修改...
熵的总结
tzl2093的博客
02-09 5708
知识点：自信息量，熵，交叉熵，相对熵，联合熵， 条件熵，信息增益
牛顿迭代法（Newton’s Method）迭代求根的Python程序
热门推荐
Fo*(Bi)的博客
10-26 4万+
参考链接： 用python算微积分及牛顿迭代求解高阶方程 不用迭代，直接得到高次方程的解： from sympy import * from sympy.abc import x def func(x): return x**3 - 0.165*x**2 + 3.993*10**(-4) result = solve(func(x),x) print(result) 结果： [-0.0437370862084818 + 0.e-24*I, 0.0623775815137495 + 0.e-22*
Python数据分析与挖掘——线性回归预测模型
Fo*(Bi)的博客
02-27 4万+
线性回归模型属于经典的统计学模型，该模型的应用场景是根据已 知的变量（自变量）来预测某个连续的数值变量（因变量）。例如，餐 厅根据每天的营业数据（包括菜谱价格、就餐人数、预定人数、特价菜 折扣等）预测就餐规模或营业额；网站根据访问的历史数据（包括新用 户的注册量、老用户的活跃度、网页内容的更新频率等）预测用户的支 付转化率；医院根据患者的病历数据（如体检指标、药物服用情况、平 时的饮食习惯等）预测某种疾病发生的概率。 站在数据挖掘的角度看待线性回归模型，它属于一种有监督的学习 算法，即在建模过程中必须同时具
Python有限差分法——向前差分，向后差分和中心差分的Python程序
Fo*(Bi)的博客
10-19 2万+
有限差分法 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。 该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 分类 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为
第二届全国大学生算法设计与编程挑战赛（赛题，共10个）
Fo*(Bi)的博客
11-23 2万+
n = 10 #跑道长度 role_style = ['>.wwws..s.','.>.wwws.s.','>>>>>s.m.m','>w.wss..s.'] for role in role_style: v = 1 T = 0 t = 0 tc = 0 r_list = list(role) for item in r_list: if item == '.': #平地
马尔可夫链(Markov chain)的性质
Fo*(Bi)的博客
08-14 2万+
本文内容主要参考：李航老师的《统计学习方法》 以下介绍离散状态马尔可夫链的性质。可以自然推广到连续状态马尔可夫链。 一、不可约
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值