概率论与数理统计浙大第五版 第一章 部分习题+R代码

习题一

5.

(1)

设事件A为从中任意抽取5片，至少是2片安慰剂

$P(A)=\frac{C_5^2{C}_5^3+C_5^3{C}_5^2+C_5^4{C}_5^1+C_5^5}{C_{10}^5}=1-P(\overline{A})=1-\frac{C_5^1{C}_5^4+C_5^5}{C_{10}^5}=\frac{226}{252}=\frac{113}{126}$

t=choose(5,2)*choose(5,3)+choose(5,3)*choose(5,2)+choose(5,4)*choose(5,1)+choose(5,5)
s=choose(10,5)
p=t/s
p

(2)

设事件A为第一次取到安慰剂，事件B为第二次取到安慰剂，事件C为第三次取到安慰剂

$P(ABC)=P(C|AB)P(B|A)P(A)=\frac{1}{2}\ast \frac{4}{9}\ast \frac{3}{8}=\frac{1}{12}$

6.

(1)

设事件A为最小号码为5号

$P(A)=\frac{C_5^2}{C_{10}^3}=\frac{1}{12}$

(2)

设事件A为最大号码为5号

$P(A)=\frac{C_4^2}{C_{10}^3}=\frac{1}{20}$

7.

设事件A为顾客恰能按所订颜色如数得到订货

$P(A)=\frac{C_{10}^4{C}_4^3{C}_3^2}{C_{17}^9}=\frac{252}{2431}$

8.

(1)

设事件A为恰有90件次品

$P(A)=\frac{C_{400}^{90}{C}_{1100}^{110}}{C_{1500}^{200}}$

(太大了没算，见R代码及结果)

(2)

设事件A为至少有两件次品

$P(A)=1-P(\overline{A})=1-\frac{C_{400}^1{C}_{1100}^{199}+C_{1100}^{200}}{C_{1500}^{200}}$

9.

设事件A为四只鞋子中至少有两只配成一双

$P(A)=1-P(\overline{A})=1-\frac{C_5^4{C_2^1}^5}{C_{10}^4}=\frac{13}{21}$

10.

设事件A为其排列结果为ability

$P(A)=\frac{C_1^1{C}_2^1{C}_2^1{C}_1^1{C}_1^1{C}_1^1{C}_1^1}{A_{11}^7}=\frac{1}{415800}$

11.

设事件A为杯子中球的最大个数为1

$P(A)=\frac{C_4^3}{3^4}$

设事件B为杯子中球的最大个数为2

$P(B)=\frac{C_4^1{C}_3^2}{3^4}$

设事件C为杯子中球的最大个数为3

$P(C)=\frac{C_4^1{C}_3^1}{3^4}$

以上做法是错误的，再此感谢热心网友指出错误，下面说明下正解，以及错误原因

---

原题意思应该是每个杯子的都是本质不同的，例如：1 1 1 0 表示一二三号杯子各有一个球，但不等同于 0 1 1 1，同理得证及求球的个数的排列。

设事件A为杯子中球的最大个数为1

$P(A)=\frac{A_4^3}{4^3}$

设事件B为杯子中球的最大个数为2

（三个球中选两个球捆绑一起，在对应投入杯中，下同理

$P(B)=\frac{C_3^2{A}_4^2}{4^3}$

设事件C为杯子中球的最大个数为3

$P(C)=\frac{C_3^3{A}_4^1}{4^3}$

R代码让我偷个懒 在此略过 (～￣▽￣)～

12.

设事件A为发生一个部件强度太弱

$P(A)=\frac{C_{10}^1{C}_{47}^3...C_{26}^3{C}_{23}^3}{C_{50}^3{C}_{47}^3...C_{26}^3{C}_{23}^3}=\frac{C_{10}^1}{C_{50}^3}=\frac{1}{1960}$

13.

(1)

设事件A为任选4名学生，一、二、三、四年级的学生各占一名

$P(A)=\frac{C_5^1{C}_2^1{C}_3^1{C}_2^1}{C_{12}^5}=\frac{4}{33}$

(2)

设事件A为任选5名学生，一、二、三、四年级的学生均包含在内

$P(A)=\frac{C_5^2{C}_2^1{C}_3^1{C}_2^1+C_5^1{C}_2^2{C}_3^1{C}_2^1+C_5^1{C}_2^1{C}_3^2{C}_2^1+C_5^1{C}_2^1{C}_3^1{C}_2^2}{C_{12}^5}=\frac{10}{33}$

R

代码

详细解释在代码注释中

# 定义一个函数，该函数为所求的积分函数
# 因此在计算上述定积分和计算sqrt(2)同理
# 在此用例为sqrt(2)
f = function(x)
  sqrt(2)
# 生成x
x = seq(1.5, 3.5, length=100)
y = rep(0, length(x))
j = 1
# 计算每个x对应的y值
for (i in x) {
  y[j] = f(i)
  j = j + 1
}
# 根据函数划出积分曲线
plot(x, y, type='o')
# 确定积分边界
abline(v = 2)
abline(v = 3)
# 简单设定随机点 (x,y)|2<=x<=3,0<=y<=100
# 即随机点的分布面积为100
s = 100
a = 0
cnt <- 1
# 有兴趣的朋友可以适当的多循环几次，这样结果更为精确，但是运行时间太长，并没有跑太久
while (cnt <= 100000) {
    # 随机点 xx，yy
    xx = 2 + runif(1)
    yy = 100 * runif(1)
    # 点在积分面积内则点为红色
    if( f(xx) > yy ) {
        a = a + 1
        points(x = xx, y = yy, pch = 20, cex = 1, col = "red")
    }
    # 否则为绿色
    else
        points(x = xx, y = yy, pch = 20, cex = 1, col = "green")
  cnt = cnt + 1
}
# 计算积分
print((a / 10000) * s)