以前我在大厅里发表过一个帖子,推导了破产风险的计算公式,但在那个帖子里分3次写的,比较凌乱,使用的代号也不够规范,推导过程也绕了好大弯子,有必要重新整理一下。
使用符号(注意大小写不一样):
p:单次赌博赢的概率 q:单次赌博输的概率 P(k):资金为k时的破产概率 B:资金(美元) b:每次下注量 W:盈利能力(美元) S:盈利能力标准差(美元) N:以前赌博样本数(小时)
使用例子: 设计一个游戏,用100张卡片,红色的55张,白色的45张。 随机抽取一张卡片,如果时红色,甲赢乙输,如果时白色,乙赢甲输。 初始甲有100美元,乙有100亿美元,每次下注1美元。 求甲破产的概率。
P(0)=1,宣布破产 甲的资金有p可能变成B+1,有q可能变成B-1。下面式子不难理解: P(B)=pP(B+1)+qP(B-1)
pP(B)+qP(B)=pP(B+1)+qP(B-1)
也可以写成: P(B+1)-P(B)=(q/p)(P(B)-P(B-1))
那么: P(2)-P(1)=(q/p)(P(1)-P(0)) P(3)-P(2)=(q/p)(P(2)-P(1))=(q/p)^2(P(1)-P(0)) P(4)-P(3)=(q/p)(P(3)-P(2))=(q/p)^3(P(1)-P(0))
…… P(B+1)-P(B)=(q/p)^B(P(1)-P(0))
P(m)-P(B)=(P(m)-P(m-1)) + …+ (P(B+1)-P(B))
求和得 P(m)-P(B)=(P(1)-P(0))((q/p)^B-(q/p)^m)/(1-q/p) 当m=0时: 1-P(B)=(1-P(1))(1-(q/p)^B)(1-(q/p)) 当m趋于无穷大时,也就是钱无穷多是P(m)趋于0,不会破产 P(B)=(1-P(1))(q/p)^B/(1-(q/p)) 得: (1-P(B))/P(B)=(1-(q/p)^B)/(q/p)^B
P(B)=(q/p)^B 如果每次下注量为b美元 显然: P(B)=(q/p)^(B/b)
解释:甲有10000美元每次下注100美元的破产概率,与甲有100美元每次下注1美元的破产概率是一样的。
如果继续玩的EV与现在的盈利能力一样,则有: EV=(p-q)b=W //EV=(55/100-45/100)100=10=W p+q=1 解得: p=1/2+W/2/b //55% q=1/2-W/2/b //45%
考虑盈利标准差: S^2=1/N*Σ(p(b-W)^2+q(-b-W)^2) =1/N*Σ(b^2+W^2-2(p-q)bW) =1/N*Σ(b^2)-W^2 =Σ(b^2)-W^2 得: S^2+W^2=b^2
解释: Σ为从1到N的求和 标准差S公式: S^2=1/NΣ(x-W)^2 其中W为x的平均值 赌博时下注量为b时 赢的概率p,赢时下注量与平均值的差为b-W 输的概率q,输时下注量与平均值的差为-b-W
P(B)=(q/p)^(B/b) =((1/2-W/2/b)/(1/2+W/2/b))^(B/b) =((1-W/b)/(1+W/b))^(B/b) ln(P(B))=(B/b)(ln(1-W/b)-ln(1+W/b)) 在x很小时: ln(1+x)可以近似写成x ln(P(B))=-2W/b*(B/b)=-2WB/b^2 =-2WB/(S^2+W^2)
P(B)=e^(-2WB/(S^2+W^2))
|