電気通信大学 情報理工学研究科 情報学専攻 2026年8月実施 選択問題 確率・オペレーションズリサーチ
Author
shank (co-authored with GPT 6.1 SOL)
Description
この科目(確率・オペレーションズリサーチ)を受験する場合には,問1は必ず解答し,問2と問3はいずれか一方のみを選択して解答すること。
期待値が μ = 100 \mu=100 μ = 100 ,標準偏差が σ = 20 \sigma=20 σ = 20 の分布に従う互いに独立な確率変数 X 1 , X 2 , … , X n X_1,X_2,\ldots,X_n X 1 , X 2 , … , X n を考える。また,これらの標本平均を X ˉ n = 1 n ∑ i = 1 n X i \bar{X}_n=\frac{1}{n}\sum_{i=1}^{n}X_i X ˉ n = n 1 ∑ i = 1 n X i とする。
問 1-1: 標本平均 X ˉ n \bar{X}_n X ˉ n の期待値 E [ X ˉ n ] E[\bar{X}_n] E [ X ˉ n ] を求めよ。
問 1-2: 標本平均 X ˉ n \bar{X}_n X ˉ n の分散 V [ X ˉ n ] V[\bar{X}_n] V [ X ˉ n ] を求めよ。
以下では,標本平均 X ˉ n \bar{X}_n X ˉ n が μ \mu μ から ± 2 \pm 2 ± 2 以上外れる確率 P ( ∣ X ˉ n − μ ∣ ≥ 2 ) P(|\bar{X}_n-\mu|\geq 2) P ( ∣ X ˉ n − μ ∣ ≥ 2 ) を 1 % 1\% 1% 以下に抑える最小の n n n ,すなわち
P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 0.01 (1) P(|\bar{X}_n-\mu|\geq 2)\leq 0.01
\tag{1} P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 0.01 ( 1 )
となる n n n を求めたい。
問 1-3: 中心極限定理により,n n n が十分に大きいとき,X ˉ n \bar{X}_n X ˉ n は近似的に正規分布に従うものとする。この近似の下で,以下の確率変数 Z Z Z を考える。
Z = X ˉ n − μ σ / n Z=\frac{\bar{X}_n-\mu}{\sigma/\sqrt{n}} Z = σ / n X ˉ n − μ
式 (1) を Z Z Z を用いた不等式 P ( ∣ Z ∣ ≥ f ( n ) ) ≤ 0.01 P(|Z|\geq f(n))\leq 0.01 P ( ∣ Z ∣ ≥ f ( n )) ≤ 0.01 に置き換えるとき,f ( n ) f(n) f ( n ) を n n n の式として求めよ。
問 1-4: 標準正規分布表から,P ( ∣ Z ∣ ≥ 2.58 ) P(|Z|\geq 2.58) P ( ∣ Z ∣ ≥ 2.58 ) がおおよそ 0.01 0.01 0.01 であることがわかっている。そこで P ( ∣ Z ∣ ≥ 2.58 ) = 0.01 P(|Z|\geq 2.58)=0.01 P ( ∣ Z ∣ ≥ 2.58 ) = 0.01 と仮定して,問 1-3 の結果から式 (1) を満たす最小の n n n の近似値を求めよ(小数点以下を切り上げて答えよ)。
問 1-5: 以下のチェビシェフの不等式
P ( ∣ X ˉ n − μ ∣ ≥ ε ) ≤ V [ X ˉ n ] ε 2 P(|\bar{X}_n-\mu|\geq\varepsilon)
\leq\frac{V[\bar{X}_n]}{\varepsilon^2} P ( ∣ X ˉ n − μ ∣ ≥ ε ) ≤ ε 2 V [ X ˉ n ]
を用いて,式 (1) を満たす最小の n n n を求めよ。
問 1-6: 問 1-4 と問 1-5 で求めた n n n の差を確認し,なぜこのような差が生じるのかを簡潔に述べよ。
Keywords:期待値:expected value,標準偏差:standard deviation,互いに独立:mutually independent,確率変数:random variable,標本平均:sample mean,分散:variance,確率:probability,最小:smallest,中心極限定理:central limit theorem,近似:approximation,正規分布:normal distribution,標準正規分布表:standard normal distribution table,おおよそ:approximately,切り上げ:round up,チェビシェフの不等式:Chebyshev's inequality
この問題を選択する場合には問3に解答してはいけない。
二項分布 B ( n , p ) B(n,p) B ( n , p ) に従う確率変数 X X X についての問いに答えよ。ただし,確率関数(確率質量関数)は下記で与えられるものとする。
P ( X = k ) = n C k p k ( 1 − p ) n − k , ( k = 0 , 1 , … , n ) P(X=k)={}_n C_k\,p^k(1-p)^{n-k},\qquad(k=0,1,\ldots,n) P ( X = k ) = n C k p k ( 1 − p ) n − k , ( k = 0 , 1 , … , n )
問 2-1: 次の事象 A から E のうち,二項分布に従うモデルとして不適切なものをすべて選べ。
A. 表が出る確率が 0.5 のコインを 10 回投げた時に,表が出る回数。
B. 赤玉 3 個,白玉 7 個が入った袋から,玉を 1 個取り出して色を確認し,それを袋に戻さずに 10 回連続で取り出したときの赤玉の回数。
C. 1 時間に平均約 2 回発生する稀なシステムエラーについて,ある特定の 1 時間におけるエラーの発生回数。
D. ある交差点を通過する車が左折する確率を 10 % 10\% 10% としたとき,次に左折車が現れるまでに通過した直進車の台数。
E. ランダムにシャッフルしたトランプの山札を上から順番にめくって,最初に「キング」が出るまでにめくるカードの枚数。
問 2-2: 二項分布に従う確率変数 X X X のモーメント母関数 E [ e t X ] E[e^{tX}] E [ e tX ] を n , p , t n,p,t n , p , t の式として求めよ。また,1 次モーメント E [ X ] E[X] E [ X ] ,2 次モーメント E [ X 2 ] E[X^2] E [ X 2 ] をそれぞれ n n n と p p p の式として求めよ。ただし n n n は 2 以上とする。
以下では,製造ラインで生産される部品の不良率を 0 < p < 1 0<p<1 0 < p < 1 とする。このラインからランダムに 2 個の部品を取得して,それを 1 組とする。また,1 組に含まれる不良品の数を確率変数 X X X とし,これは二項分布 B ( 2 , p ) B(2,p) B ( 2 , p ) に従うものとする。
管理者は 1 組内の不良品数 X X X が事前に決めておいた基準値 k k k 未満の場合にはその 1 組を「出荷」し,それ以外の場合には「廃棄」する。「出荷」した場合には,1 組内の不良品が 1 個ならば c 1 c_1 c 1 円の損失を,不良品が 2 個ならば c 2 c_2 c 2 円の損失を被るとする。「廃棄」を選択した場合は不良品の個数に関わらず B B B 円の損失を被るとする。ただし c 1 < B < c 2 c_1<B<c_2 c 1 < B < c 2 である。また,基準値を k k k としたときの損失の期待値を R ( k ) R(k) R ( k ) とする。
問 2-3: k = 1 k=1 k = 1 としたときの損失の期待値 R ( 1 ) R(1) R ( 1 ) を p , B p,B p , B の式として求めよ。
問 2-4: k = 2 k=2 k = 2 としたときの損失の期待値 R ( 2 ) R(2) R ( 2 ) を p , c 1 , B p,c_1,B p , c 1 , B の式として求めよ。
問 2-5: 損失の期待値 R ( k ) R(k) R ( k ) が最小になる k k k を求めよ。
Keywords:二項分布:binomial distribution,確率質量関数:probability mass function,事象:event,不適切:inappropriate,モーメント母関数:moment generating function,1 次モーメント:first moment,2 次モーメント:second moment,製造ライン:production line,不良率:defective product rate,組:pair,未満:less than,出荷:ship,廃棄:discard,損失:loss
この問題を選択する場合には,問2に解答してはいけない。
以下の式 (1)–(5) からなる最大化問題をシンプレックス法で解くことを考える:
max 0.4 x 1 + 0.5 x 2 ( 1 ) s.t. 0.3 x 1 + 0.1 x 2 ≤ 2.7 , ( 2 ) 0.5 x 1 + 0.5 x 2 = 6 , ( 3 ) 0.6 x 1 + 0.4 x 2 ≥ 6 , ( 4 ) x 1 , x 2 ≥ 0. ( 5 ) \begin{array}{rll}
\max & 0.4x_1+0.5x_2 & (1) \\
\text{s.t.} & 0.3x_1+0.1x_2\leq 2.7, & (2) \\
& 0.5x_1+0.5x_2=6, & (3) \\
& 0.6x_1+0.4x_2\geq 6, & (4) \\
& x_1,x_2\geq 0. & (5)
\end{array} max s.t. 0.4 x 1 + 0.5 x 2 0.3 x 1 + 0.1 x 2 ≤ 2.7 , 0.5 x 1 + 0.5 x 2 = 6 , 0.6 x 1 + 0.4 x 2 ≥ 6 , x 1 , x 2 ≥ 0. ( 1 ) ( 2 ) ( 3 ) ( 4 ) ( 5 )
このためには,以下の手続きを実行しなければならない。
制約をスラック変数および人為変数を用いて等式制約に変換する。
最初の非基底変数を選んで実行基底解を求める。ここで,基底解とは変換した等式制約からなる連立 1 次方程式の解であり,実行基底解とはその基底変数の値が全て非負となる基底解のことである。
求めた実行基底解の値を改善できるかどうかを考えるために,シンプレックス表を作る。もし改善できなければシンプレックス法を終了する。改善できる場合は,目的関数の値を最も大きく改良する非基底変数を 1 つ選ぶ。
現行の基底変数の中から非基底変数にする変数を選ぶ。ステップ 3 で選んだ非基底変数を増加させると,等式制約を満たすように,現行の基底変数の値を減少させなければならない。したがって,現行の基底変数が負にならない範囲で,選んだ非基底変数を増加させる。
次に基底解に入る基底変数と次に基底解から出る非基底変数の交点係数を 1 とし,その列の他の係数を 0 にするようにシンプレックス表を更新する。
シンプレックス表を更新すると,実行基底解も更新されるので,ステップ 3 に戻ってこれまでの手順を繰り返す。
問 3-1: ステップ 3 で最初に作ったシンプレックス表を解答欄にしたがって作成せよ。ただし,目的関数の係数は解答欄にあるものへ合わせることとし,M M M は十分大きな罰金を表す。また,z j z_j z j の行は c j c_j c j と x j x_j x j の列の各要素の積の和である。
問 3-2: ステップ 2 で選んだ非基底変数を答えよ。
問 3-3: 1 回目のステップ 3 から 5 でシンプレックス表を更新したとき,変数 x 2 x_2 x 2 の列における係数はそれぞれ 1 / 3 , 1 / 3 , 1 / 5 1/3,\ 1/3,\ 1/5 1/3 , 1/3 , 1/5 となる。それに基づいてこのシンプレックス表を完成させなさい。
問 3-4: 式 (1)–(5) の最大化問題の最適解(各変数の値)とその解における目的関数の値をシンプレックス法で求めよ。さらに,解答欄にしたがってシンプレックス表を完成させよ。ただし,1 回目の表は問 3-1,2 回目の表は問 3-3 で解答しているので,ここには 3 回目以降のシンプレックス表を書くこと。また解答欄には 5 回目まで表を用意しているがすべてを埋める必要があるとは限らない。
Keywords:最大化問題:maximization problem,シンプレックス法:simplex method,制約:constraints,スラック変数:slack variable,人為変数:artificial variable,等式制約:equality constraints,非基底変数:non basic variable,実行基底解:basic feasible solution,非負:non-negative,目的関数:objective function,係数:coefficient,表:table,罰金:penalty,最適解:optimal solution,完成:complete
题目描述
报考本科目(概率与运筹学)时,问1为必答题,问2与问3只能选择其中一题作答。
考虑相互独立的随机变量 X 1 , X 2 , … , X n X_1,X_2,\ldots,X_n X 1 , X 2 , … , X n ,它们服从期望为 μ = 100 \mu=100 μ = 100 、标准差为 σ = 20 \sigma=20 σ = 20 的分布。定义它们的样本均值为 X ˉ n = 1 n ∑ i = 1 n X i \bar{X}_n=\frac{1}{n}\sum_{i=1}^{n}X_i X ˉ n = n 1 ∑ i = 1 n X i 。
问 1-1: 求样本均值 X ˉ n \bar{X}_n X ˉ n 的期望 E [ X ˉ n ] E[\bar{X}_n] E [ X ˉ n ] 。
问 1-2: 求样本均值 X ˉ n \bar{X}_n X ˉ n 的方差 V [ X ˉ n ] V[\bar{X}_n] V [ X ˉ n ] 。
以下希望求出最小的 n n n ,使样本均值 X ˉ n \bar{X}_n X ˉ n 与 μ \mu μ 的偏差绝对值至少为 2 的概率 P ( ∣ X ˉ n − μ ∣ ≥ 2 ) P(|\bar{X}_n-\mu|\geq 2) P ( ∣ X ˉ n − μ ∣ ≥ 2 ) 不超过 1 % 1\% 1% ,即满足
P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 0.01. (1) P(|\bar{X}_n-\mu|\geq 2)\leq 0.01.
\tag{1} P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 0.01. ( 1 )
问 1-3: 根据中心极限定理,当 n n n 足够大时,假定 X ˉ n \bar{X}_n X ˉ n 近似服从正态分布。在这一近似下,考虑下列随机变量 Z Z Z :
Z = X ˉ n − μ σ / n . Z=\frac{\bar{X}_n-\mu}{\sigma/\sqrt{n}}. Z = σ / n X ˉ n − μ .
将式 (1) 改写成使用 Z Z Z 的不等式 P ( ∣ Z ∣ ≥ f ( n ) ) ≤ 0.01 P(|Z|\geq f(n))\leq 0.01 P ( ∣ Z ∣ ≥ f ( n )) ≤ 0.01 时,求 f ( n ) f(n) f ( n ) 关于 n n n 的表达式。
问 1-4: 由标准正态分布表可知,P ( ∣ Z ∣ ≥ 2.58 ) P(|Z|\geq 2.58) P ( ∣ Z ∣ ≥ 2.58 ) 约为 0.01 0.01 0.01 。因此,假定 P ( ∣ Z ∣ ≥ 2.58 ) = 0.01 P(|Z|\geq 2.58)=0.01 P ( ∣ Z ∣ ≥ 2.58 ) = 0.01 ,利用问 1-3 的结果,求满足式 (1) 的最小 n n n 的近似值(小数部分向上取整)。
问 1-5: 利用下列切比雪夫不等式,求满足式 (1) 的最小 n n n :
P ( ∣ X ˉ n − μ ∣ ≥ ε ) ≤ V [ X ˉ n ] ε 2 . P(|\bar{X}_n-\mu|\geq\varepsilon)
\leq\frac{V[\bar{X}_n]}{\varepsilon^2}. P ( ∣ X ˉ n − μ ∣ ≥ ε ) ≤ ε 2 V [ X ˉ n ] .
问 1-6: 比较问 1-4 与问 1-5 所求出的 n n n 的差异,并简要说明产生这一差异的原因。
关键词:期望:expected value,标准差:standard deviation,相互独立:mutually independent,随机变量:random variable,样本均值:sample mean,方差:variance,概率:probability,最小:smallest,中心极限定理:central limit theorem,近似:approximation,正态分布:normal distribution,标准正态分布表:standard normal distribution table,大约:approximately,向上取整:round up,切比雪夫不等式:Chebyshev's inequality
选择本题时,不得作答问3。
对于服从二项分布 B ( n , p ) B(n,p) B ( n , p ) 的随机变量 X X X ,回答下列问题。其概率函数(概率质量函数)如下:
P ( X = k ) = n C k p k ( 1 − p ) n − k , ( k = 0 , 1 , … , n ) . P(X=k)={}_n C_k\,p^k(1-p)^{n-k},\qquad(k=0,1,\ldots,n). P ( X = k ) = n C k p k ( 1 − p ) n − k , ( k = 0 , 1 , … , n ) .
问 2-1: 从以下 A 至 E 的情形中,选出所有不适合用二项分布建模的情形。
A. 将正面朝上的概率为 0.5 的硬币投掷 10 次,正面出现的次数。
B. 从装有 3 个红球、7 个白球的袋中,每次取出 1 个球并确认颜色,取出后不放回;连续取球 10 次时,取出的红球个数。
C. 对于平均每小时约发生 2 次的罕见系统错误,在某个特定的 1 小时内发生的错误次数。
D. 设通过某个交叉路口的车辆左转的概率为 10 % 10\% 10% ,在下一辆左转车辆出现之前通过的直行车辆数。
E. 将一副扑克牌随机洗牌后,依次从牌堆顶部翻牌,直到第一次出现“K”为止翻出的牌数。
问 2-2: 求服从二项分布的随机变量 X X X 的矩母函数 E [ e t X ] E[e^{tX}] E [ e tX ] ,用 n , p , t n,p,t n , p , t 表示。此外,分别用 n n n 与 p p p 表示一阶矩 E [ X ] E[X] E [ X ] 和二阶矩 E [ X 2 ] E[X^2] E [ X 2 ] 。其中,n ≥ 2 n\geq 2 n ≥ 2 。
以下设某条生产线所生产零件的不良率为 0 < p < 1 0<p<1 0 < p < 1 。从该生产线随机取得 2 个零件,组成一组。将一组中不良品的个数记为随机变量 X X X ,并假定 X X X 服从二项分布 B ( 2 , p ) B(2,p) B ( 2 , p ) 。
管理者预先确定一个阈值 k k k :当一组中的不良品个数 X X X 小于 k k k 时,将该组“出货”;否则将该组“报废”。出货时,如果一组中有 1 个不良品,则损失 c 1 c_1 c 1 日元;如果有 2 个不良品,则损失 c 2 c_2 c 2 日元。选择报废时,无论不良品个数多少,均损失 B B B 日元。假定 c 1 < B < c 2 c_1<B<c_2 c 1 < B < c 2 。将阈值设为 k k k 时的期望损失记为 R ( k ) R(k) R ( k ) 。
问 2-3: 当 k = 1 k=1 k = 1 时,求期望损失 R ( 1 ) R(1) R ( 1 ) ,用 p , B p,B p , B 表示。
问 2-4: 当 k = 2 k=2 k = 2 时,求期望损失 R ( 2 ) R(2) R ( 2 ) ,用 p , c 1 , B p,c_1,B p , c 1 , B 表示。
问 2-5: 求使期望损失 R ( k ) R(k) R ( k ) 最小的 k k k 。
关键词:二项分布:binomial distribution,概率质量函数:probability mass function,事件:event,不适当:inappropriate,矩母函数:moment generating function,一阶矩:first moment,二阶矩:second moment,生产线:production line,不良率:defective product rate,组:pair,小于:less than,出货:ship,报废:discard,损失:loss
选择本题时,不得作答问2。
考虑用单纯形法求解由下列式 (1)–(5) 构成的最大化问题:
max 0.4 x 1 + 0.5 x 2 ( 1 ) s.t. 0.3 x 1 + 0.1 x 2 ≤ 2.7 , ( 2 ) 0.5 x 1 + 0.5 x 2 = 6 , ( 3 ) 0.6 x 1 + 0.4 x 2 ≥ 6 , ( 4 ) x 1 , x 2 ≥ 0. ( 5 ) \begin{array}{rll}
\max & 0.4x_1+0.5x_2 & (1) \\
\text{s.t.} & 0.3x_1+0.1x_2\leq 2.7, & (2) \\
& 0.5x_1+0.5x_2=6, & (3) \\
& 0.6x_1+0.4x_2\geq 6, & (4) \\
& x_1,x_2\geq 0. & (5)
\end{array} max s.t. 0.4 x 1 + 0.5 x 2 0.3 x 1 + 0.1 x 2 ≤ 2.7 , 0.5 x 1 + 0.5 x 2 = 6 , 0.6 x 1 + 0.4 x 2 ≥ 6 , x 1 , x 2 ≥ 0. ( 1 ) ( 2 ) ( 3 ) ( 4 ) ( 5 )
为此,必须执行以下步骤。
使用松弛变量和人工变量,将约束转换为等式约束。
选定初始的非基变量,求出可行基解。这里,基解是由转换后的等式约束构成的线性方程组的解;可行基解是所有基变量的值均非负的基解。
为判断能否改善所求可行基解的值,建立单纯形表。如果无法改善,则结束单纯形法;如果可以改善,则选择一个使目标函数值改善最大的非基变量。
从当前的基变量中选出一个将变为非基变量的变量。增大步骤 3 所选择的非基变量时,为满足等式约束,必须减小当前基变量的值。因此,在当前基变量均不变为负数的范围内,增大所选择的非基变量。
将下一次进入基解的基变量与下一次离开基解的非基变量所对应的交点系数变为 1,并更新单纯形表,使该列其他系数变为 0。
更新单纯形表后,可行基解也随之更新,因此返回步骤 3,重复以上过程。
问 3-1: 按照答题栏的格式,建立在步骤 3 中首次得到的单纯形表。目标函数的系数应与答题栏中给出的系数一致,M M M 表示足够大的罚值。此外,z j z_j z j 行由 c j c_j c j 与 x j x_j x j 列各对应元素乘积之和构成。
问 3-2: 写出在步骤 2 中所选择的非基变量。
问 3-3: 第一次执行步骤 3 至 5、更新单纯形表后,变量 x 2 x_2 x 2 所在列的系数分别为 1 / 3 , 1 / 3 , 1 / 5 1/3,\ 1/3,\ 1/5 1/3 , 1/3 , 1/5 。据此完成该单纯形表。
问 3-4: 用单纯形法求式 (1)–(5) 所给最大化问题的最优解(各变量的值)及该解对应的目标函数值。此外,按照答题栏的格式完成单纯形表。第 1 张表已在问 3-1 中作答,第 2 张表已在问 3-3 中作答,因此这里应填写第 3 张及之后的单纯形表。答题栏提供了直到第 5 张的表格,但不一定需要全部填写。
关键词:最大化问题:maximization problem,单纯形法:simplex method,约束:constraints,松弛变量:slack variable,人工变量:artificial variable,等式约束:equality constraints,非基变量:non basic variable,可行基解:basic feasible solution,非负:non-negative,目标函数:objective function,系数:coefficient,表:table,罚值:penalty,最优解:optimal solution,完成:complete
Kai
問 1-1
期待値の線形性より,
E [ X ˉ n ] = 1 n ∑ i = 1 n E [ X i ] = 100. E[\bar X_n]=\frac{1}{n}\sum_{i=1}^n E[X_i]=100. E [ X ˉ n ] = n 1 i = 1 ∑ n E [ X i ] = 100.
問 1-2
独立性より共分散は 0 0 0 であるから,
V [ X ˉ n ] = 1 n 2 ∑ i = 1 n V [ X i ] = 20 2 n = 400 n . V[\bar X_n]=\frac{1}{n^2}\sum_{i=1}^n V[X_i]
=\frac{20^2}{n}=\frac{400}{n}. V [ X ˉ n ] = n 2 1 i = 1 ∑ n V [ X i ] = n 2 0 2 = n 400 .
問 1-3
標準偏差が 20 / n 20/\sqrt n 20/ n であることから,
∣ X ˉ n − μ ∣ ≥ 2 ⟺ ∣ Z ∣ ≥ 2 20 / n = n 10 . |\bar X_n-\mu|\geq 2
\iff |Z|\geq\frac{2}{20/\sqrt n}
=\frac{\sqrt n}{10}. ∣ X ˉ n − μ ∣ ≥ 2 ⟺ ∣ Z ∣ ≥ 20/ n 2 = 10 n .
したがって,f ( n ) = n / 10 f(n)=\sqrt n/10 f ( n ) = n /10 である。
問 1-4
問題文で指定された近似の下では,
n 10 ≥ 2.58 ⟺ n ≥ 25.8 2 = 665.64. \frac{\sqrt n}{10}\geq 2.58
\iff n\geq 25.8^2=665.64. 10 n ≥ 2.58 ⟺ n ≥ 25. 8 2 = 665.64.
よって,求める最小の整数は n = 666 \boxed{n=666} n = 666 である。
問 1-5
チェビシェフの不等式より,
P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 400 / n 2 2 = 100 n . P(|\bar X_n-\mu|\geq2)
\leq\frac{400/n}{2^2}=\frac{100}{n}. P ( ∣ X ˉ n − μ ∣ ≥ 2 ) ≤ 2 2 400/ n = n 100 .
この上界を 0.01 0.01 0.01 以下にする条件は n ≥ 10000 n\geq10000 n ≥ 10000 である。したがって,この不等式によって条件を保証できる最小の整数は n = 10000 \boxed{n=10000} n = 10000 である。
問 1-6
差は 10000 − 666 = 9334 10000-666=9334 10000 − 666 = 9334 である。中心極限定理による値は標本平均の分布を正規分布で近似して求めた値であり,有限の n n n に対する厳密な保証ではない。一方,チェビシェフの不等式は平均と分散のみから成り立つ一般的な上界であり,分布の形を利用しないため保守的になる。なお,10000 10000 10000 はこの上界から得られる十分条件であって,元の分布における実際の最小標本数を決定したものではない。
問 2-1
固定回数の独立なベルヌーイ試行における成功回数として不適切なのは,B , C , D , E \boxed{\mathrm{B,C,D,E}} B , C , D , E である。
A:独立なコイン投げを仮定すれば B ( 10 , 0.5 ) B(10,0.5) B ( 10 , 0.5 ) に従う。
B:非復元抽出なので試行が独立ではなく,成功確率も変化する。超幾何分布によるモデルであり,特に本問では全 10 10 10 個を取り出すため赤玉数は必ず 3 3 3 となる。
C:一定時間内の発生回数であり,固定された試行回数がない。一定率で独立に発生するという仮定を加える場合は,ポアソン分布によるモデルが適する。
D:最初の左折車が現れるまでの待ち回数である。各車の左折・非左折を独立な試行として扱う場合は,幾何分布によるモデルになる。
E:非復元抽出で最初のキングが現れるまでの待ち回数であり,固定回数の独立試行の成功回数ではない。
問 2-2
二項定理より,モーメント母関数は
M X ( t ) = E [ e t X ] = ∑ k = 0 n ( n k ) ( p e t ) k ( 1 − p ) n − k = ( 1 − p + p e t ) n . M_X(t)=E[e^{tX}]
=\sum_{k=0}^n {n\choose k}(pe^t)^k(1-p)^{n-k}
=(1-p+pe^t)^n. M X ( t ) = E [ e tX ] = k = 0 ∑ n ( k n ) ( p e t ) k ( 1 − p ) n − k = ( 1 − p + p e t ) n .
これを微分すると,
M X ′ ( t ) = n p e t ( 1 − p + p e t ) n − 1 , M X ′ ′ ( t ) = n p e t ( 1 − p + p e t ) n − 1 + n ( n − 1 ) p 2 e 2 t ( 1 − p + p e t ) n − 2 . \begin{aligned}
M_X'(t)&=npe^t(1-p+pe^t)^{n-1},\\
M_X''(t)&=npe^t(1-p+pe^t)^{n-1}
+n(n-1)p^2e^{2t}(1-p+pe^t)^{n-2}.
\end{aligned} M X ′ ( t ) M X ′′ ( t ) = n p e t ( 1 − p + p e t ) n − 1 , = n p e t ( 1 − p + p e t ) n − 1 + n ( n − 1 ) p 2 e 2 t ( 1 − p + p e t ) n − 2 .
よって,
E [ X ] = n p , E [ X 2 ] = n p + n ( n − 1 ) p 2 = n p ( 1 − p ) + n 2 p 2 . \boxed{E[X]=np},\qquad
\boxed{E[X^2]=np+n(n-1)p^2=np(1-p)+n^2p^2}. E [ X ] = n p , E [ X 2 ] = n p + n ( n − 1 ) p 2 = n p ( 1 − p ) + n 2 p 2 .
問 2-3
不良品が 0 0 0 個の組を出荷した場合の損失を 0 0 0 とする。k = 1 k=1 k = 1 では X = 0 X=0 X = 0 の場合だけ出荷し,X = 1 , 2 X=1,2 X = 1 , 2 の場合に廃棄するので,
R ( 1 ) = B { P ( X = 1 ) + P ( X = 2 ) } = B { 2 p ( 1 − p ) + p 2 } = B ( 2 p − p 2 ) . R(1)=B\{P(X=1)+P(X=2)\}
=B\{2p(1-p)+p^2\}
=\boxed{B(2p-p^2)}. R ( 1 ) = B { P ( X = 1 ) + P ( X = 2 )} = B { 2 p ( 1 − p ) + p 2 } = B ( 2 p − p 2 ) .
問 2-4
k = 2 k=2 k = 2 では X = 0 , 1 X=0,1 X = 0 , 1 の場合に出荷し,X = 2 X=2 X = 2 の場合だけ廃棄する。したがって,
R ( 2 ) = c 1 P ( X = 1 ) + B P ( X = 2 ) = 2 c 1 p ( 1 − p ) + B p 2 . R(2)=c_1P(X=1)+BP(X=2)
=\boxed{2c_1p(1-p)+Bp^2}. R ( 2 ) = c 1 P ( X = 1 ) + BP ( X = 2 ) = 2 c 1 p ( 1 − p ) + B p 2 .
問 2-5
損失額を非負とする通常の解釈では,0 ≤ c 1 < B < c 2 0\leq c_1<B<c_2 0 ≤ c 1 < B < c 2 である。閾値を整数とすると,常に廃棄する k ≤ 0 k\leq0 k ≤ 0 ,k = 1 k=1 k = 1 ,k = 2 k=2 k = 2 ,常に出荷する k ≥ 3 k\geq3 k ≥ 3 の四つの方策を比較すればよい。
R ( k ) = B ( k ≤ 0 ) , R ( k ) = 2 c 1 p ( 1 − p ) + c 2 p 2 ( k ≥ 3 ) . \begin{aligned}
R(k)&=B &&(k\leq0),\\
R(k)&=2c_1p(1-p)+c_2p^2 &&(k\geq3).
\end{aligned} R ( k ) R ( k ) = B = 2 c 1 p ( 1 − p ) + c 2 p 2 ( k ≤ 0 ) , ( k ≥ 3 ) .
0 < p < 1 0<p<1 0 < p < 1 および c 1 < B < c 2 c_1<B<c_2 c 1 < B < c 2 より,
B − R ( 2 ) = B ( 1 − p ) 2 + 2 ( B − c 1 ) p ( 1 − p ) > 0 , R ( 1 ) − R ( 2 ) = 2 ( B − c 1 ) p ( 1 − p ) > 0 , R ( 3 ) − R ( 2 ) = ( c 2 − B ) p 2 > 0. \begin{aligned}
B-R(2)&=B(1-p)^2+2(B-c_1)p(1-p)>0,\\
R(1)-R(2)&=2(B-c_1)p(1-p)>0,\\
R(3)-R(2)&=(c_2-B)p^2>0.
\end{aligned} B − R ( 2 ) R ( 1 ) − R ( 2 ) R ( 3 ) − R ( 2 ) = B ( 1 − p ) 2 + 2 ( B − c 1 ) p ( 1 − p ) > 0 , = 2 ( B − c 1 ) p ( 1 − p ) > 0 , = ( c 2 − B ) p 2 > 0.
したがって,最適な整数の閾値は k = 2 \boxed{k=2} k = 2 である。閾値に実数も許す場合,1 < k ≤ 2 1<k\leq2 1 < k ≤ 2 はすべて同じ出荷・廃棄方策を与える。
以下では,x 3 x_3 x 3 を式 (2) のスラック変数,x 4 x_4 x 4 を式 (3) の人為変数,x 5 x_5 x 5 を式 (4) の余剰変数,x 6 x_6 x 6 を式 (4) の人為変数として表を作る。いずれも非負とし,大 M M M 法の目的関数を
z = 2 5 x 1 + 1 2 x 2 − M x 4 − M x 6 z=\frac25x_1+\frac12x_2-Mx_4-Mx_6 z = 5 2 x 1 + 2 1 x 2 − M x 4 − M x 6
と定める。等式制約は
3 10 x 1 + 1 10 x 2 + x 3 = 27 10 , 1 2 x 1 + 1 2 x 2 + x 4 = 6 , 3 5 x 1 + 2 5 x 2 − x 5 + x 6 = 6. \begin{aligned}
\frac3{10}x_1+\frac1{10}x_2+x_3&=\frac{27}{10},\\
\frac12x_1+\frac12x_2+x_4&=6,\\
\frac35x_1+\frac25x_2-x_5+x_6&=6.
\end{aligned} 10 3 x 1 + 10 1 x 2 + x 3 2 1 x 1 + 2 1 x 2 + x 4 5 3 x 1 + 5 2 x 2 − x 5 + x 6 = 10 27 , = 6 , = 6.
各表の基底変数を x B i x_{B_i} x B i ,その目的関数係数を c B i c_{B_i} c B i とする。z j = ∑ i c B i a i j z_j=\sum_i c_{B_i}a_{ij} z j = ∑ i c B i a ij ,Δ j = c j − z j \Delta_j=c_j-z_j Δ j = c j − z j とし,正の Δ j \Delta_j Δ j が最大の変数を基底に入れる。右辺列の z j z_j z j は現在の目的関数値を表す。以下の表はこの変数名・符号規約に従う。
問 3-1
初期基底を ( x 3 , x 4 , x 6 ) (x_3,x_4,x_6) ( x 3 , x 4 , x 6 ) とすれば,その値は ( 27 / 10 , 6 , 6 ) (27/10,6,6) ( 27/10 , 6 , 6 ) である。c j = ( 2 / 5 , 1 / 2 , 0 , − M , 0 , − M ) c_j=(2/5,1/2,0,-M,0,-M) c j = ( 2/5 , 1/2 , 0 , − M , 0 , − M ) より,初期表は次のようになる。
基底/評価行 x 1 x_1 x 1 x 2 x_2 x 2 x 3 x_3 x 3 x 4 x_4 x 4 x 5 x_5 x 5 x 6 x_6 x 6 右辺 x 3 x_3 x 3 3 / 10 3/10 3/10 1 / 10 1/10 1/10 1 1 1 0 0 0 0 0 0 0 0 0 27 / 10 27/10 27/10 x 4 x_4 x 4 1 / 2 1/2 1/2 1 / 2 1/2 1/2 0 0 0 1 1 1 0 0 0 0 0 0 6 6 6 x 6 x_6 x 6 3 / 5 3/5 3/5 2 / 5 2/5 2/5 0 0 0 0 0 0 − 1 -1 − 1 1 1 1 6 6 6 z j z_j z j − 11 M / 10 -11M/10 − 11 M /10 − 9 M / 10 -9M/10 − 9 M /10 0 0 0 − M -M − M M M M − M -M − M − 12 M -12M − 12 M Δ j \Delta_j Δ j 11 M / 10 + 2 / 5 11M/10+2/5 11 M /10 + 2/5 9 M / 10 + 1 / 2 9M/10+1/2 9 M /10 + 1/2 0 0 0 0 0 0 − M -M − M 0 0 0 —
問 3-2
初期の非基底変数は x 1 , x 2 , x 5 \boxed{x_1,x_2,x_5} x 1 , x 2 , x 5 であり,いずれも 0 0 0 に設定する。
問 3-3
M M M が十分大きいとき,初期表で最大の正の評価値を持つのは x 1 x_1 x 1 である。x 1 x_1 x 1 列の正の係数に対する比は
27 / 10 3 / 10 = 9 , 6 1 / 2 = 12 , 6 3 / 5 = 10. \frac{27/10}{3/10}=9,\qquad
\frac{6}{1/2}=12,\qquad
\frac{6}{3/5}=10. 3/10 27/10 = 9 , 1/2 6 = 12 , 3/5 6 = 10.
最小比は 9 9 9 なので,x 1 x_1 x 1 を入れ,x 3 x_3 x 3 を出す。第1行の 3 / 10 3/10 3/10 をピボットにすると,第2表は次のようになる。
基底/評価行 x 1 x_1 x 1 x 2 x_2 x 2 x 3 x_3 x 3 x 4 x_4 x 4 x 5 x_5 x 5 x 6 x_6 x 6 右辺 x 1 x_1 x 1 1 1 1 1 / 3 1/3 1/3 10 / 3 10/3 10/3 0 0 0 0 0 0 0 0 0 9 9 9 x 4 x_4 x 4 0 0 0 1 / 3 1/3 1/3 − 5 / 3 -5/3 − 5/3 1 1 1 0 0 0 0 0 0 3 / 2 3/2 3/2 x 6 x_6 x 6 0 0 0 1 / 5 1/5 1/5 − 2 -2 − 2 0 0 0 − 1 -1 − 1 1 1 1 3 / 5 3/5 3/5 z j z_j z j 2 / 5 2/5 2/5 2 / 15 − 8 M / 15 2/15-8M/15 2/15 − 8 M /15 4 / 3 + 11 M / 3 4/3+11M/3 4/3 + 11 M /3 − M -M − M M M M − M -M − M 18 / 5 − 21 M / 10 18/5-21M/10 18/5 − 21 M /10 Δ j \Delta_j Δ j 0 0 0 11 / 30 + 8 M / 15 11/30+8M/15 11/30 + 8 M /15 − 4 / 3 − 11 M / 3 -4/3-11M/3 − 4/3 − 11 M /3 0 0 0 − M -M − M 0 0 0 —
x 2 x_2 x 2 列は,題意のとおり 1 / 3 , 1 / 3 , 1 / 5 1/3,1/3,1/5 1/3 , 1/3 , 1/5 となっている。
問 3-4
第2表では x 2 x_2 x 2 を入れる。比は 27 , 9 / 2 , 3 27,9/2,3 27 , 9/2 , 3 であり,最小比 3 3 3 に対応する x 6 x_6 x 6 を出す。第3行の 1 / 5 1/5 1/5 をピボットにすると,第3表を得る。
基底/評価行 x 1 x_1 x 1 x 2 x_2 x 2 x 3 x_3 x 3 x 4 x_4 x 4 x 5 x_5 x 5 x 6 x_6 x 6 右辺 x 1 x_1 x 1 1 1 1 0 0 0 20 / 3 20/3 20/3 0 0 0 5 / 3 5/3 5/3 − 5 / 3 -5/3 − 5/3 8 8 8 x 4 x_4 x 4 0 0 0 0 0 0 5 / 3 5/3 5/3 1 1 1 5 / 3 5/3 5/3 − 5 / 3 -5/3 − 5/3 1 / 2 1/2 1/2 x 2 x_2 x 2 0 0 0 1 1 1 − 10 -10 − 10 0 0 0 − 5 -5 − 5 5 5 5 3 3 3 z j z_j z j 2 / 5 2/5 2/5 1 / 2 1/2 1/2 − 7 / 3 − 5 M / 3 -7/3-5M/3 − 7/3 − 5 M /3 − M -M − M − 11 / 6 − 5 M / 3 -11/6-5M/3 − 11/6 − 5 M /3 11 / 6 + 5 M / 3 11/6+5M/3 11/6 + 5 M /3 47 / 10 − M / 2 47/10-M/2 47/10 − M /2 Δ j \Delta_j Δ j 0 0 0 0 0 0 7 / 3 + 5 M / 3 7/3+5M/3 7/3 + 5 M /3 0 0 0 11 / 6 + 5 M / 3 11/6+5M/3 11/6 + 5 M /3 − 11 / 6 − 8 M / 3 -11/6-8M/3 − 11/6 − 8 M /3 —
ここでは Δ 3 − Δ 5 = 1 / 2 > 0 \Delta_3-\Delta_5=1/2>0 Δ 3 − Δ 5 = 1/2 > 0 なので,x 3 x_3 x 3 を入れる。正の係数を持つ行の比は
8 20 / 3 = 6 5 , 1 / 2 5 / 3 = 3 10 . \frac{8}{20/3}=\frac65,\qquad
\frac{1/2}{5/3}=\frac3{10}. 20/3 8 = 5 6 , 5/3 1/2 = 10 3 .
よって x 4 x_4 x 4 を出し,第2行の 5 / 3 5/3 5/3 をピボットにすると,第4表は次のようになる。
基底/評価行 x 1 x_1 x 1 x 2 x_2 x 2 x 3 x_3 x 3 x 4 x_4 x 4 x 5 x_5 x 5 x 6 x_6 x 6 右辺 x 1 x_1 x 1 1 1 1 0 0 0 0 0 0 − 4 -4 − 4 − 5 -5 − 5 5 5 5 6 6 6 x 3 x_3 x 3 0 0 0 0 0 0 1 1 1 3 / 5 3/5 3/5 1 1 1 − 1 -1 − 1 3 / 10 3/10 3/10 x 2 x_2 x 2 0 0 0 1 1 1 0 0 0 6 6 6 5 5 5 − 5 -5 − 5 6 6 6 z j z_j z j 2 / 5 2/5 2/5 1 / 2 1/2 1/2 0 0 0 7 / 5 7/5 7/5 1 / 2 1/2 1/2 − 1 / 2 -1/2 − 1/2 27 / 5 27/5 27/5 Δ j \Delta_j Δ j 0 0 0 0 0 0 0 0 0 − M − 7 / 5 -M-7/5 − M − 7/5 − 1 / 2 -1/2 − 1/2 1 / 2 − M 1/2-M 1/2 − M —
M > 1 / 2 M>1/2 M > 1/2 ならすべての評価値が非正となるため,ここで終了する。人為変数は x 4 = x 6 = 0 x_4=x_6=0 x 4 = x 6 = 0 なので元の問題に対しても実行可能であり,
( x 1 , x 2 ) = ( 6 , 6 ) , z max = 2 5 ⋅ 6 + 1 2 ⋅ 6 = 27 5 = 5.4 . \boxed{(x_1,x_2)=(6,6)},\qquad
\boxed{z_{\max}=\frac25\cdot6+\frac12\cdot6=\frac{27}{5}=5.4}. ( x 1 , x 2 ) = ( 6 , 6 ) , z m a x = 5 2 ⋅ 6 + 2 1 ⋅ 6 = 5 27 = 5.4 .
補助変数の値は ( x 3 , x 4 , x 5 , x 6 ) = ( 3 / 10 , 0 , 0 , 0 ) (x_3,x_4,x_5,x_6)=(3/10,0,0,0) ( x 3 , x 4 , x 5 , x 6 ) = ( 3/10 , 0 , 0 , 0 ) である。第5表は不要である。
また,元の等式制約から x 2 = 12 − x 1 x_2=12-x_1 x 2 = 12 − x 1 を代入すると,実行可能範囲は 6 ≤ x 1 ≤ 15 / 2 6\leq x_1\leq15/2 6 ≤ x 1 ≤ 15/2 ,目的関数は 6 − x 1 / 10 6-x_1/10 6 − x 1 /10 となる。これは x 1 x_1 x 1 に関して減少するので,x 1 = 6 x_1=6 x 1 = 6 で最大となり,上の結果とも一致する。