Computer Science Theory for the Information Age-2: 高维空间中的正方体和Chernoff Bounds

时间:2022-08-31 11:12:37

高维空间中的正方体和Chernoff Bounds

本文将介绍高维空间中正方体的一些性质,以及一个非常常见也是非常有用的概率不等式——Chernoff Bounds。

考虑$d$维单位正方体$C=\{x|0\leq x_i\leq 1,i=1,\cdots,d\}$,其中心点为$(\frac{1}{2},\cdots,\frac{1}{2})$,体积为1。现在我们将其半径收缩到$1-\frac{c}{d}$,其体积为$(1-\frac{c}{d})^d\leq e^{-c}$,所以当$d$很大时,高维正方体的体积总是分布在其边缘地带。

定义超平面$H=\{x|\sum_{i=1}^dx_i=\frac{d}{2}\}$,即过中心点但不过原点的对角面。现在我们从正方体$C$中均匀随机的产生观察点$x$(相当与从$[0,1]$独立均匀的选取$x_1,\cdots,x_d$),$x=(x_1,x_2,\cdots,x_d)$到$H$的距离为:

\begin{equation} L=\frac{1}{\sqrt{d}}|(\sum_{i=1}^dx_i-\frac{d}{2}|\end{equation}

这个距离平方的期望为:

\begin{equation}\mathbb{E}(L^2)=\frac{1}{d}\mathbb{E}[(\sum_{i=1}^dx_i-\frac{d}{2})^2]=\frac{1}{d}\mathop{Var}[\sum_{i=1}^dx_i]=\frac{1}{d}\frac{d}{12}=\frac{1}{12}\end{equation}

其中$\mathbb{E}(\sum_{i=1}^dx_i)=\frac{d}{2}, \mathop{Var}(\sum_{i=1}^dx_i)=\frac{d}{4}$。所以根据Markov不等式$\mathbb{P}(|x|\geq a)\leq\frac{\mathbb{E}(|x|)}{a}$可得:

$$\mathbb{P}(L\geq t)=\mathbb{P}(L^2\geq t^2)\leq\frac{\mathbb{E}(L^2)}{t^2}=\frac{1}{12t^2}$$

因此我们可以得到如下引理:

引理一 在$C$内随机均匀的选一点,则该点到超平面的距离在$t$以内的概率至少为$1-\frac{1}{12t^2}$,即$\mathbb{P}(L\leq t)\geq1-\frac{1}{12t^2}$。

接下去,我们将证明一个比引理一更一般的引理,这个引理在证明Chernoff Bounds时会用到。

引理二 令$x_1,x_2,\cdots,x_d$为独立的随机变量,且$0\leq x_i \leq 1$,$\mathbb{E}(x_i)=p_i$。令$y_i=x_i-p_i$,且记$\mu=\sum_{i=1}^dp_i$。那么对任意的正整数$n$有:

\begin{equation}\mathbb{E}[(\sum_{i=1}^dy_i)^n]\leq \mathop{Max}\{(2n\mu)^\frac{n}{2},n^n\}\end{equation}

证明:首先,我们将$(y_1+y_2+\cdots+y_d)^n$写成单项式的求和形式,即$(y_1+y_2+\cdots+y_d)^n=\sum_{I\in S}\prod_{i\in I}y_i^{r_i}$,其中$r_i$表示在每一个单项式中$y_i$出现的次数,$I$表示非零$r_i$对应的下标集合,$S=\{I|\sum_{i\in I}r_i=n\}$。所以$\mathbb{E}[(y_1+y_2+\cdots+y_d)^n]=\mathbb{E}[\sum_{I\in S}\prod_{i\in I}y_i^{r_i}]$。

现在我们先计算其中单个单项式的期望。由于随机变量之间的相互独立性,所以$\mathbb{E}(\prod_{i\in I}y_i^{r_i})=\prod_{i\in I}\mathbb{E}(y_i^{r_i})$,另外又因为$\mathbb{E}(y_i)=0$,所以这里我们可以只考虑$r_i\geq 2$,所以每个集合$I$的大小将小于等于$\frac{n}{2}$,即$|I|\leq\frac{d}{2}$。由于$y_i\in [-p_i,1-p_i]$,所以:

\begin{align*}\mathbb{E}[|y_i^{r_i}|]&\leq\mathbb{E}(y_i^2)=\mathbb{E}[(x_i-p_i)^2]\\&=\mathbb{E}(x_i^2)-p_i^2\leq \mathbb{E}(x_i^2)\leq\mathbb{E}(x_i)=p_i\end{align*}

因此,

$$\prod_{i\in I}\mathbb{E}(y_i^{r_i})\leq \prod_{i\in I}\mathbb{E}(|y_i^{r_i}|)\leq\prod_{i\in I}p_i\triangleq p(I)$$

也就是每个单项式的期望不会超过$p(I)$,所以$\mathbb{E}[(\sum_{i=1}^dy_i)^r]\leq\sum_{I,|I|\leq\frac{n}{2}}p(I)N(I)$,其中$N(I)$表示此单项式出现的次数。且$I$对应的单项式数量不会超过按如下方式产生的单项式数量,即每次从$|I|$(因为该单项式只有$|I|$个因子可供选择)中选择一个因子,然后选择$n$次,故$n(I)\leq |I|^n$。

同时,

\begin{align}\sum_{I:|I|=t}p(I)&=\sum_{I:|I|=t}(\prod_{i\in I}p_i)\leq (\sum_{i=1}^dp_i)^t\frac{1}{t!}\label{equ:exp1}\\&=\frac{\mu^t}{t!}\approx\frac{\mu^t}{\sqrt{2\pi t}(\frac{t}{e})^t}\label{equ:exp2}\end{align}

其中等式\ref{equ:exp1}成立的原因:所有$t$个不同的$p_i$相乘的和必定小于从全部的$d$个$p_i$中选$t$次,并且把重复的$t!$个排列当成相同的单项式。等式\ref{equ:exp2}成立是因为$t!\approx \sqrt{2\pi t}(\frac{t}{e})^t$。所以:

\begin{equation}\mathbb{E}[(\sum_{i=1}^d)^r]\leq\sum_{t=1}^\frac{r}{2}\frac{\mu^tt^n}{\sqrt{2\pi}t^te^{-t}}\leq\frac{\mathop{Max}_{t=1}^\frac{n}{2}f(t)}{\sqrt{2\pi}}\sum_{t=1}^\frac{r}{2}t^n\end{equation}

这里$f(t)=\frac{(e\mu)^t}{t^t}$。对$f(t)$求导可知,在$t<\mu$时,$f(t)$为增函数;在$t>\mu$时,$f(t)$为减函数。故我们可以分两种情况讨论:1)当$\mu<\frac{n}{2}$时,$\mathop{Max}_{t=1}^\frac{n}{2}f(t)=f(\mu)=e^\mu\leq e^\frac{n}{2}$;2)当$\mu>\frac{n}{2}$时,$\mathop{Max}_{t=1}^\frac{n}{2}f(t)=f(\frac{n}{2})\leq\frac{(2e\mu)^\frac{n}{2}}{n^\frac{n}{2}}$。所以:

\begin{align}\mathbb{E}[(\sum_{i=1}^dy_i)^r]&\leq\frac{2}{\sqrt{2\pi}}\mathop{Max}[(\frac{2e\mu}{n})^\frac{n}{2},e^\frac{n}{2}](\frac{n}{2})^n\label{equ:exp3}\\&\leq\mathop{Max}[(\frac{en\mu}{2})^\frac{n}{2},(\frac{en^2}{4})^\frac{n}{2}]\nonumber\\&\leq\mathop{Max}[(2n\mu)^\frac{n}{2},n^n]\nonumber\end{align}

其中利用了不等式$\sum_{t=1}^\frac{n}{2}t^n\leq\int_{0}^{\frac{n}{2}}x^ndx\leq\frac{n}{2(n+1)}(\frac{n}{2})^n\leq\frac{1}{2}(\frac{n}{2})^n$。

好了,有了上面的这个引理后,我们就可以证明这个有用的Chernoff Bounds。

定理一 Chernoff Bounds

假设$x_i,y_i,\mu$与引理二中的一样,那么:

\begin{equation}\mathbb{P}(|\sum_{i=1}^dy_i|\geq t)\leq 3e^{-\frac{t^2}{12\mu}},\quad\text{for } 0<t\leq 3\mu\label{equ:cher1}\end{equation}

\begin{equation}\mathbb{P}(|\sum_{i=1}^dy_i|\geq t)\leq 2\times 2^{-\frac{t}{3}},\quad\text{for } t>3\mu\label{equ:cher2}\end{equation}

证明:令$r$为正偶数,$y=\sum_{i=1}^dy_i$,所以$y^r$是非负的。根据Markov不等式有:$\mathbb{P}(|y|\geq t)=\mathbb{P}(y^r\geq t^r)\leq\frac{\mathbb{E}(y^r)}{t^r}$。根据引理二,有$\mathbb{P}(|y|\geq t)\leq\mathop{Max}[\frac{(2r\mu)^\frac{r}{2}}{t^r},\frac{r^r}{t^r}]$,对所有$r$为偶数均成立。

经过简单的计算(求导),我们可以知道$\frac{(2r\mu)^\frac{r}{2}}{t^r}$的最小值在点$r_{min}=\frac{r^2}{2e\mu}$处取得。由于$r_{min}$不一定会是偶数,所以我们取不超过$r_{min}$的最大偶数$r$,且:

1)对所有的$t$:

\begin{align} (\frac{2r\mu}{t^2})^{-\frac{r}{2}}&\leq e^{-\frac{r}{2}}\label{equ:exp4}\\&\leq e^{1-\frac{t^2}{4e\mu}}\label{equ:exp5}\\&\leq 3e^{-\frac{t^2}{12\mu}}\label{equ:1}\end{align}

其中不等式\ref{equ:exp4}是由于$r\leq\frac{t^2}{2e\mu}$,不等式\ref{equ:exp5}是由于$\frac{r_{min}-r}{2}\leq 1\Longrightarrow -\frac{r}{2}\leq 1-\frac{r_{min}}{2}$.

2)当$0<t\leq 3\mu$时,

\begin{align}\frac{r^r}{t^r}&\leq(\frac{t}{3e\mu})^r\leq(\frac{3\mu}{2e\mu})^r=(\frac{2e}{3})^{-r}\nonumber\\&\leq(\sqrt{e})^{-r}=e^{-\frac{r}{2}}<3e^{-\frac{t^2}{12\mu}}\label{equ:2}\end{align}

综合不等式\ref{equ:1}和\ref{equ:2}可知,不等式\ref{equ:cher1}成立。

对于第二个不等式,选择$r$为不超过$\frac{2}{3}t$的最大偶数,即$r\leq\frac{2}{3}t$。又$t>3\mu$,故有:

\begin{equation}\frac{r^r}{t^r}\leq(\frac{4}{9})^\frac{r}{2}\leq(\frac{1}{2})^\frac{r}{2}=2^{-\frac{r}{2}}\end{equation}

\begin{equation}(\frac{2\mu r}{t^2})^\frac{r}{2}\leq(\frac{2tr}{3t^2})^\frac{r}{2}=(\frac{2}{3}\frac{r}{t})^\frac{r}{2}\leq\frac{1}{2}^\frac{r}{2}=2^{-\frac{r}{2}}\end{equation}

所以$\mathop{Max}[\frac{(2r\mu)^\frac{r}{2}}{t^r},\frac{r^r}{t^r}]\leq 2^{-\frac{r}{2}}$。由于$\frac{\frac{2}{3}t-r}{2}\leq 1\Longrightarrow -\frac{r}{2}\leq 1-\frac{t}{3}$,所以$2^{-\frac{r}{2}}\leq 2^{1-\frac{t}{3}}=2\times 2^{-\frac{t}{3}}$,所以不等式\ref{equ:cher2}成立

Computer Science Theory for the Information Age-2: 高维空间中的正方体和Chernoff Bounds的更多相关文章

  1. Computer Science Theory for the Information Age-4&colon; 一些机器学习算法的简介

    一些机器学习算法的简介 本节开始,介绍<Computer Science Theory for the Information Age>一书中第六章(这里先暂时跳过第三章),主要涉及学习以 ...

  2. Computer Science Theory for the Information Age-1&colon; 高维空间中的球体

    高维空间中的球体 注:此系列随笔是我在阅读图灵奖获得者John Hopcroft的最新书籍<Computer Science Theory for the Information Age> ...

  3. Computer Science Theory for the Information Age-6&colon; 学习理论——VC定理的证明

    VC定理的证明 本文讨论VC理论的证明,其主要内容就是证明VC理论的两个定理,所以内容非常的枯燥,但对于充实一下自己的理论知识也是有帮助的.另外,VC理论属于比较难也比较抽象的知识,所以我总结的这些证 ...

  4. Computer Science Theory for the Information Age-3&colon; 高维空间中的高斯分布和随机投影

    高维空间中的高斯分布和随机投影 (一)在高维球体表面产生均匀分布点的方法 我们来考虑一个采样问题,就是怎样在高维单位球体的表面上均匀的采样.首先,考虑二维的情况,就是在球形的周长上采样.我们考虑如下方 ...

  5. Computer Science Theory for the Information Age-5&colon; 学习理论——VC维的定义以及一些例子

    学习理论——VC维的定义以及一些例子 本文主要介绍一些学习理论上的东西.首先,我们得明确,从训练集上学习出来的分类器的最终目标是用于预测未知的样本,那么我们在训练的时候该用多少的样本才能使产生的分类器 ...

  6. Discovering the Computer Science Behind Postgres Indexes

    This is the last in a series of Postgres posts that Pat Shaughnessy wrote based on his presentation ...

  7. &lbrack;转载&rsqb; A set of top Computer Science blogs

    This started out as a list of top Computer Science blogs, but it more closely resembles a set: the o ...

  8. CSCI 1100 — Computer Science 1 Homework

    CSCI 1100 — Computer Science 1 Homework 8CS1 Multiverse: ClassesOverviewThis homework is worth 100 p ...

  9. 2012年Elsevier旗下Computer Science期刊最新SCI影响因子排名

    Latest Impact Factor figures from Elsevier's Computer Science Journals Medical Image Analysis Impact ...

随机推荐

  1. Tensorflow serving的编译

    Tensorflow serving提供了部署tensorflow生成的模型给线上服务的方法,包括模型的export,load等等. 安装参考这个 https://github.com/tensorf ...

  2. Hadoop 裡的 fsck 指令

    Hadoop 裡的 fsck 指令,可檢查 HDFS 裡的檔案 (file),是否有 corrupt (毀損) 或資料遺失,並產生 HDFS 檔案系統的整體健康報告.報告內容,包括:Total blo ...

  3. Debian安装Apache2&plus;MySQL5&plus;PHP5(zz)

    转载:http://hi.baidu.com/lostdays/item/1d5e7e4833b4d20fc116134b 终于在Debian用apt-get安装好LAMP了,之前在CentOS使用编 ...

  4. 《Java并发编程实战》第七章 取消与关闭 读书笔记

        Java没有提供不论什么机制来安全地(抢占式方法)终止线程,尽管Thread.stop和suspend等方法提供了这种机制,可是因为存在着一些严重的缺陷,因此应该避免使用. 但它提供了中断In ...

  5. 两款【linux字符界面下】显示【菜单】,【选项】的powershell脚本模块介绍

    两款[linux字符界面下]显示[菜单],[选项]的powershell脚本模块介绍 powershell linux  ps1 menu choice Multiselect 传教士 菜单 powe ...

  6. day18 python之re模块与正则表达式

    正则表达式 正则表达式,就是匹配字符串内容的一种规则. 官方定义:正则表达式是对字符串操作的一种逻辑公式,就是用事先定义好的一些特定字符.及这些特定字符的组合,组成一个“规则字符串”,这个“规则字符串 ...

  7. Android 安全机制

    转:http://www.cnblogs.com/GnagWang/archive/2011/03/21/1990507.html 1 Android 安全机制概述 Android 是一个权限分离的系 ...

  8. MyEclipse连接CVS,如果中间经过一层代理,就没法直接联接到CVS了,哪位知道怎么办?

  9. POJ2456 Aggressive cows 2017-05-11 17&colon;54 38人阅读 评论&lpar;0&rpar; 收藏

    Aggressive cows Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 13993   Accepted: 6775 ...

  10. 在window 8 或windows2012 上用命令行安装framework3&period;5 方法

    找到对应操作系统安装目录的sources文件夹下的sxs文件夹,拷贝到本地电脑,如F:盘 根目录下 CMD(管理员身份)命令: dism.exe /online /enable-feature /fe ...