我是靠谱客的博主 腼腆白云,最近开发中收集的这篇文章主要介绍【OR】Robust Optimization(8): Probability and tail boundsProb and Tail BoundsExamples,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

Navigator

  • Prob and Tail Bounds
  • Examples
    • example 1
    • example 2
    • example 3

Prob and Tail Bounds

机会约束的概率形式为
P r o b [ ∑ i = 1 n z i U i ≥ t ] ≤ ε Prob[sum_{i=1}^n z_iU_igeq t]leq varepsilon Prob[i=1nziUit]ε
其中, U i U_i Ui是随机变量, z i z_i zi为固定参数. 可以使用很多方法来控制概率,其中一些tail bounds的方法证明是有有效的.
最简单的tail boundMarkov inequality,即对于非负的随机变量 U U U a ≥ 0 ageq 0 a0
P r o b ( U ≥ a ) ≤ E [ U ] a Prob(Ugeq a)leq frac{mathbb{E}[U]}{a} Prob(Ua)aE[U]
根据此可以推出Chebyshev bound
P r o b ( U > E U + t ) ≤ P r o b ( ( U − E U ) 2 ≥ t 2 ) ≤ E ( U − E U ) 2 t 2 Prob(U>mathbb{E}U+t)leq Prob((U-mathbb{E}U)^2geq t^2)leq frac{mathbb{E}(U-mathbb{E}U)^2}{t^2} Prob(U>EU+t)Prob((UEU)2t2)t2E(UEU)2
为了描述极端低概率事件的bounds,一种常用的技术是chernoff bound.
P r o b ( U ≥ t ) = P r o b ( e λ U ≥ e λ U ) ≤ E e λ U e λ t Prob(Ugeq t)=Prob(e^{lambda U}geq e^{lambda U})leq frac{mathbb{E}e^{lambda U}}{e^{lambda t}} Prob(Ut)=Prob(eλUeλU)eλtEeλU
或者等价的
P r o b ( U ≥ t ) ≤ inf ⁡ λ ≥ 0 E e λ U e − λ U Prob(Ugeq t)leq inf_{lambdageq 0}mathbb{E}e^{lambda U}e^{-lambda U} Prob(Ut)λ0infEeλUeλU

Examples

example 1

Z Z Z为一个零均值高斯随机变量,方差为 σ 2 sigma^2 σ2
E exp ⁡ ( λ Z ) = exp ⁡ ( λ 2 σ 2 2 ) mathbb{E}exp(lambda Z)=exp(frac{lambda^2sigma^2}{2}) Eexp(λZ)=exp(2λ2σ2)
并且对于 t ≥ 0 tgeq 0 t0
P r o b ( Z ≥ t ) ≤ inf ⁡ λ ≥ 0 exp ⁡ ( λ 2 σ 2 2 − λ t ) = exp ⁡ ( − t 2 2 σ 2 ) Prob(Zgeq t)leq inf_{lambdageq 0}exp(frac{lambda^2sigma^2}{2}-lambda t)=exp(-frac{t^2}{2sigma^2}) Prob(Zt)λ0infexp(2λ2σ2λt)=exp(2σ2t2)

Whenever there is a parameter σ ≥ 0 sigmageq 0 σ0 such that a mean-zero random variable satisfies
E e λ U ≤ exp ⁡ ( λ 2 σ 2 2 ) mathbb{E}e^{lambda U}leq exp(frac{lambda^2sigma^2}{2}) EeλUexp(2λ2σ2)
we call the random variable sub-Gaussian with parameter σ 2 sigma^2 σ2, this is an extremely useful concept.

example 2

Indeed, we have Hoeffding’ s inequality which gives sharp bounds on sums of sub-Gaussian random variables.

U 1 , U 2 , … , U n U_1, U_2, dots, U_n U1,U2,,Un为独立的mean-zero sub Gaussian random variables参数为 σ = ( σ 1 , … , σ n ) sigma=(sigma_1, dots, sigma_n) σ=(σ1,,σn)
P r o b ( ∑ i = 1 n U i ≥ t ) ≤ exp ⁡ ( − t 2 2 ∥ σ ∥ 2 2 ) Prob(sum_{i=1}^n U_igeq t)leq exp(-frac{t^2}{2lVertsigmarVert_2^2}) Prob(i=1nUit)exp(2σ22t2)
使用Chernoff bound,令 Z = 1 T U Z=mathbf{1}^TU Z=1TU,可以得到
P r o b ( Z ≥ t ) ≤ exp ⁡ ( λ 2 ∑ i = 1 n σ i 2 2 − λ t ) Prob(Zgeq t)leq exp(frac{lambda^2sum_{i=1}^nsigma_i^2}{2}-lambda t) Prob(Zt)exp(2λ2i=1nσi2λt)
λ ≥ 0 lambdageq 0 λ0区间最小化可以得到tail bound.

example 3

Bound random variables(Hoeffding’s lemma) U ∈ [ a , b ] Uin[a, b] U[a,b]为一个mean-zero random variable,对于 a ≤ 0 ≤ b aleq 0leq b a0b U U U是一个sub-Gaussian,并且
E [ exp ⁡ ( λ U ) ] ≤ exp ⁡ ( λ 2 ( b − a ) 2 8 ) mathbb{E}[exp(lambda U)]leq exp(frac{lambda^2(b-a)^2}{8}) E[exp(λU)]exp(8λ2(ba)2)
Proof:
由之前的方程可以得到
P r o b ( ∑ t = 1 n U i ≥ t ) ≤ exp ⁡ ( − 2 t 2 ∑ i = 1 n ( b i − a i ) 2 ) Prob(sum_{t=1}^nU_igeq t)leq exp(-frac{2t^2}{sum_{i=1}^n(b_i-a_i)^2}) Prob(t=1nUit)exp(i=1n(biai)22t2)
由指数函数的凸性可以得到
e λ x ≤ e λ a b − x b − a + e λ b x − a b − a e^{lambda x}leq e^{lambda a}frac{b-x}{b-a}+e^{lambda b}frac{x-a}{b-a} eλxeλababx+eλbbaxa
关于随机变量 U U U取期望,可以得到
E [ e λ U ] ≤ b b − a e λ a + − a b − a e λ b mathbb{E}[e^{lambda U}]leq frac{b}{b-a}e^{lambda a}+frac{-a}{b-a}e^{lambda b} E[eλU]babeλa+baaeλb
z = λ ( b − a ) , p = b b − a , p z = λ b z=lambda(b-a), p=frac{b}{b-a},pz=lambda b z=λ(ba),p=babpz=λb ( 1 − p ) z = − λ a (1-p)z=-lambda a (1p)z=λa所以
E [ e λ U ] ≤ p λ a + ( 1 − p ) e λ b = e ( p − 1 ) z [ p + ( 1 − p ) e z ] mathbb{E}[e^{lambda U}]leq p^{lambda a}+(1-p)e^{lambda b}=e^{(p-1)z}[p+(1-p)e^z] E[eλU]pλa+(1p)eλb=e(p1)z[p+(1p)ez]
f ( z ) = ( p − 1 ) z + log ⁡ ( p + ( 1 − p ) e z ) f(z)=(p-1)z+log(p+(1-p)e^z) f(z)=(p1)z+log(p+(1p)ez),求出一阶导数和二阶导数得到
{ f ′ ( z ) = ( p − 1 ) + 1 − p p e − z + 1 − p f ′ ′ ( z ) = p ( 1 − p ) e − z ( p e − z + 1 − p ) 2 ≤ 1 4 begin{cases} f'(z)=(p-1)+frac{1-p}{pe^{-z}+1-p}\ f''(z)=frac{p(1-p)e^{-z}}{(pe^{-z}+1-p)^2}leq frac{1}{4} end{cases} {f(z)=(p1)+pez+1p1pf(z)=(pez+1p)2p(1p)ez41
且可以计算出 f ( 0 ) = 0 , f ′ ( 0 ) = 0 f(0)=0, f'(0)=0 f(0)=0,f(0)=0,根据Taylor公式可以得到 f ( z ) ≤ 1 8 z 2 f(z)leq frac{1}{8}z^2 f(z)81z2,当 b = − a b=-a b=a时,可以得到以下简单的证明
E [ e λ U ] ≤ b 2 b e − λ b + b 2 b e λ b = ∑ k ∈ Z + λ k b k k ! = 1 + ∑ k ≥ 1 λ 2 k b 2 k ( 2 k ) ! ≤ 1 + ∑ k ≥ 1 ( λ b 2 ) 2 1 k ! = exp ⁡ ( λ 2 b 2 2 ) mathbb{E}[e^{lambda U}]leq frac{b}{2b}e^{-lambda b}+frac{b}{2b}e^{lambda b}=sum_{kinmathbb{Z}_+}frac{lambda^kb^k}{k!}=1+sum_{kgeq 1}frac{lambda^{2k}b^{2k}}{(2k)!}leq 1+sum_{kgeq 1}(frac{lambda b}{2})^2frac{1}{k!}=exp(frac{lambda^2b^2}{2}) E[eλU]2bbeλb+2bbeλb=kZ+k!λkbk=1+k1(2k)!λ2kb2k1+k1(2λb)2k!1=exp(2λ2b2)

最后

以上就是腼腆白云为你收集整理的【OR】Robust Optimization(8): Probability and tail boundsProb and Tail BoundsExamples的全部内容,希望文章能够帮你解决【OR】Robust Optimization(8): Probability and tail boundsProb and Tail BoundsExamples所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(40)

评论列表共有 0 条评论

立即
投稿
返回
顶部