zoukankan      html  css  js  c++  java
  • Maximum Entropy Model(最大熵模型)初理解

    0,熵的描述

    熵(entropy)指的是体系的混沌的程度(可也理解为一个随机变量的不确定性),它在控制论概率论数论天体物理生命科学等领域都有重要应用,在不同的学科中也有引申出的更为具体的定义,是各领域十分重要的参量。熵由鲁道夫·克劳修斯(Rudolf Clausius)提出,并应用在热力学中。后来在,克劳德·艾尔伍德·香农(Claude Elwood Shannon)第一次将熵的概念引入到信息论中来。----baidu

     下面我们将从随机变量开始一步一步慢慢理解熵。

    1,随机变量(random variable)

    1.1 随机变量(random variable)

    什么事随机变量?

    表示随机现象(在一定条件下,并不总是出现相同结果的现象称为随机现象)各种结果的实值函数(一切可能的样本点)。如掷一颗骰子,它的所有可能结果是出现1点、2点、3点、4点、5点和6点 ,若定义X为掷一颗骰子时出现的点数,则X为一随机变量。

    随机变量   X{1,2,3,4,5,6}

          图(1)

     

    1.2 随机变量概率(The probability of a random variable)

    什么是随机变量的概率?

    要全面了解一个随机变量,不但要知道它取哪些值,而且要知道它取这些值的规律,即要掌握它的概率分布。概率分布可以由分布函数刻画。若知道一个随机变量的分布函数,则它取任何值和它落入某个数值区间内的概率都可以求出。所以我们可以P(X=x)其中一种情况出现的概率。而P(X)我们叫它为概率分布函数。

    如上述掷一颗骰子,X是均匀分布 X~U[1,6]。而P(X)的分布函数如下,也可以看出P(X=1)=1/6

                  图(2)

    又如某一地区的大学生身高为正态分布,若定义X为男性身高可能出现的值,则X也是一个随机变量,服从X~N(172.70, 8.01)。用P(X)表示随机变量的概率分布。

    分布函数

                                      

    概率分布图如下,而每个学生的身高都对应了一个概率,如P(X=1.7)就能得到相应的概率

                  图(3)

     1.3 随机变量的期望(Expected value)

    期望又如何表示,表示什么?

    假设随机变量X有值x1概率为p1,X有值x2概率为p2,..X有值xk概率为pk。则离散随机变量的期望可以定义为:

    如骰子点数的期望E[X]=1/6(1+2+3+4+5+6)

    其实期望就是我们生活中常常遇到的平均值,相当于用一个值简单的描述一个随机变量的分布情况与取值情况。

    1.4 随机变量的随机性(randomness of a random variable)

    如何刻画一个随机变量的随机性(不确定性)?

    a:韦小宝用灌铅作弊骰子投掷的随机变量X更随机 还是我们用正常骰子投掷更随机呢?

    明显正常骰子投掷更为随机。正常骰子投掷时有6种可能性{1,2,3,4,5,6},而韦小宝投掷时只有一种可能性6.

    b:那么随机性又是受什么影响呢?

    X的值?明显不是!X取{1,2,3,4,5,6}与X取{2,3,4,5,6,7}可能性都是6种。

    那么什么影响了X的随机性呢?

                    图(4)

    从上面对比图可以看出P(X)不同,它们的随机性不同。正常骰子的分布更为随机。

    但是图(5)哪个随机变量分布哪个更随机呢?

    第一个图服从均匀分布X~U(1,6)

    第二个图服从高斯分布X~N(3.5,1)

    两组都是X都是{1,2,3,4,5,6}只是它们取得概率不同。 

                    图(5)              

     因此我们急需引入一个概念来表述随机变量的随机性,这也就是我们将要说的熵(Entropy)。

    2, 最大熵原理(Principle of maximum entropy)

    2.1 熵的定义

    熵的计算公式:(为什么要将熵计算公式定义成这样? 香农这样定义肯定有他的道理哈。在后面推导以及应用的时候就能感受到香农这么定义的强大。

     

    其中,x表示随机变量,与之相对应的是所有可能输出的集合,定义为符号集,随机变量的输出用x表示。P(x)表示输出概率函数。变量的不确定性越大,熵也就越大,把它搞清楚所需要的信息量也就越大.  

    现在我们可以用熵的公式来比较图4与5中到底哪个更随机了。

    P(X)均匀分布,正常骰子投掷时的熵:      =2.584963                           

    P(X)特殊分布,韦小宝骰子投掷时的熵:=0

    P(X)正态分布,正态骰子投掷的熵:H(X)=2.028845

    代入正态分布函数,可以得到 6个点的概率分别为{ 0.0175283 0.1295176 0.3520653 0.3520653 0.1295176 0.0175283}

    代入公式H(X)可以推出H(X)=

    -(0.0175283*log2(0.0175283)+
       0.1295176*log2(0.1295176)+
       0.3520653*log2(0.3520653)+
       0.3520653*log2(0.3520653)+
       0.1295176*log2(0.1295176)+
       0.0175283*log2(0.0175283)
      )=2.028845

    故得出结论:对于投掷骰子这个事件,随机变量X,当X概率分布P(X)是均匀分布时,熵H(X)值最大,是最随机的。

    猜测:对于一个随机变量X,当它的分布是均匀分布时,它的熵H(X)是最大的(X是最随机的)。       这也就是我们说的最大熵。

    2.2 单约束 最大熵推导

    单约束最大熵的基本想法就是在一定条件下(概率和为1),找到一个分布p*(X),使熵H(X)的值达到最大。可以写成:

            

    在约束下求最大值,使用拉格朗日乘子法。设

                    

    推出拉格朗日方程

    要求最大值,则考虑对求偏到导得到:

          可以改写成:   

    =>  

    从而得到方程组

               得到   

    得到

    证明猜想是正确的,一个随机变量X,当p(X)为均匀分布时,熵H(X)最大。

    2.3 多约束 最大熵推导

    多约束最大熵的基本想法就是在多个条件下(共m+1个约束),找到一个分布p*(X),使熵H(X)的值达到最大。可以写成:

    =>拉格朗日方程

    对所有求偏导

    =>   

    =>

    又因为:

    => 

    =>

    3,最大熵模型(Maximum Entropy Model)

     上面讲的都是关于一个随机变量的熵H(X),对于条件概率模型P(Y|X)的条件熵H(Y|X)也常常被用到,定义为

     

     其中集合Z包含变量X,Y的所有范围。(x,y都是向量)

    具体推导如下:

     

    求最大熵就相当于求一个条件分布p(y|x)使得条件熵H(y|x)最大,其中x,y表示向量

     

    已知存在m+1个约束:

      

    p(x,y)表示边缘概率,f(x,y)表示x与y的函数,~p(x,y)表示经验分布

    使用拉格朗日乘子法可以得到拉格朗日方程:

    于各个部分对p(x1),p(x2)...求偏导数,

    最后解方程可得

    具体推导请看 http://pan.baidu.com/s/1i31hnEX 

      

    5,引用

    http://pan.baidu.com/s/1i31hnEX

    http://pan.baidu.com/s/1ntBO2pj

    http://pan.baidu.com/s/1nt9M7ln

    http://pan.baidu.com/s/1o6v7vfW

    http://pan.baidu.com/s/1hqvJ9lE

    http://pan.baidu.com/s/1qWHhYSO

  • 相关阅读:
    M3U8文件
    线程间的通信方式以及线程与进程的区别
    string类的实现
    进程间的通信方式
    GDB调试多线程
    面向对象与面向过程的优缺点
    同步与异步,阻塞与非阻塞的区别
    计算机cpu、寄存器、内存区别
    Linux下Makefile中动态链接库和静态链接库的生成与调用
    使用session来存储用户的登录信息
  • 原文地址:https://www.cnblogs.com/qianwen/p/4001574.html
Copyright © 2011-2022 走看看