找回密码
 欢迎注册
楼主: mathe

[转载] 抛硬币出现连续正面的概率

[复制链接]
发表于 2012-9-2 08:25:44 | 显示全部楼层
问一下,有一人,抛硬币直到出现连续7个正面为止,问他抛了n>=7次出现以上情况的概率
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-9-23 10:39:05 | 显示全部楼层
好几年前就关注这个问题,根据楼主推导,用matlab编出来了,谢谢楼主

评分

参与人数 1金币 +20 收起 理由
gxqcn + 20 首贴奖励,欢迎常来。

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-9-23 10:40:45 | 显示全部楼层
问题可做进一步推广,比如单次正面概率p不等于0.5的情况又如何?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-9-24 13:14:00 | 显示全部楼层
楼主在第4楼的发言提供了很好的基本思路!
推广:
以p代替概率1/2,
直接编程运算Plianxu(n,t,p),概率转移矩阵M可以稍微改动,
1<=j<=t,     M(1,j)=1-p;
t+1<=j<=n+1,    M(1,j)=0;
2<=i<=t,    M(i,i-1)=p;
t+1<=i<=n+1,    M(i,i)=1;   运用MTALAB运算,可以比较容易得出n<=1000下的概率
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-10-26 10:56:19 | 显示全部楼层
楼主求的是至少5次连续出现正面,如果要求连续5次正面,应该是p(至少5次)-p(至少6次),这样减出的结果有可能分别求出连续2、3次的概率,有可能p(2)<p(3),谁能解释下原因

评分

参与人数 1金币 +20 收起 理由
gxqcn + 20 首帖奖励,欢迎常来。

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-20 15:17:01 | 显示全部楼层
问一下,m次实验中连续出现n次正面的概率。有公式吗

评分

参与人数 1金币 +20 收起 理由
gxqcn + 20 首帖奖励,欢迎常来。

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-21 09:37:43 | 显示全部楼层
问一下,m次实验中连续出现n次正面的概率。有公式吗
lianghr 发表于 2013-1-20 15:17


不知你要的是怎样的公式?
简单公式肯定是没有的,现在世上最好的公式是楼主在31#楼给出的公式。其中,n对应你的m;t对应你的n。如果m和n都不大,求解高次方程解有难处,也可按楼主在4#楼给出的矩阵公式计算:
p(m,n)=uM^mb
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-12-12 20:36:11 | 显示全部楼层
一般情况的公式Run
只不过上面给出的是生成函数的方法。我有一种思路大家看看:
问题  试验中事件A发生的概率设为p, 那么n次独立重复试验中,至少连续发生r次的概率P(r,n)为多少?
由于至少连续发生r次,也就是说,至少存在一个不少于r次的连续事件发生,且剩余中可以存在一系列分隔开的连续发生事件(次数可以小于r,也可以大于或等于r)。因此,原始问题就相当于至多连续发生r次,r+1次,r+2次,...,n次的概率之和。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2015-8-7 16:42:57 | 显示全部楼层
这个r应该怎么算呢?

评分

参与人数 1金币 +20 收起 理由
gxqcn + 20 首帖奖励,欢迎常来。

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2018-12-20 17:45:03 | 显示全部楼层
好老的帖子。发生在我进入论坛之前,盘古开天地的时候。
这个完全是被sheng_jianguo在帖子 https://bbs.emath.ac.cn/thread-15625-1-1.html 带逛过来的,

我来计算一下,答案是精确值 $6993823047305143749226306585/158456325028528675187087900672 = 0.0441372286404216877406785935788813088029671991207867937263842...$,与6#,7#的计算完全吻合。
附注Mathematica代码:
  1. t=10;m=DiagonalMatrix[ConstantArray[1/2,t],-1,t+1];
  2. m[[1]]=ConstantArray[1/2,t+1];m[[1]][[t+1]]=0;m[[t+1]][[t+1]]=1;
  3. MatrixPower[m,100][[-1,1]]
复制代码


毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-4-23 23:24 , Processed in 0.056148 second(s), 22 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表