清华大学随机过程答案1
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
{(
0
)( 2
)}
Cξ (t1, t2) = E
η
cos
t1
−
1 2
cos
t1
η
cos
t2
−
1 2
cos
t2
= E {η2}
cos t1 cos t2 −
E
{η} cos t1
cos t2
+
1 4
cos t1
cos t2
=
1 3
cos t1 cos t2
−
1 2
cos t1 cos t2
+
1 4
cos t1
其中 n = 1, 2 . . . , k = 1, 2 . . .
3. 质点在直线上做随机运动,即在 t = 1, 2, 3, · · · 时质点可以在 x 轴上往右或往左做一个 单位距离的随机游动。若往右移动一个单位距离的概率为 p,往左移动一个单位距离的 概率为 q,即 P {ξ (i) = +1} = p,P {ξ (i) = −1} = q,p + q = 1,且各次游动是相互统 ∑n 计独立的。经过 n 次游走,质点所处的位置为 ηn = η (n) = ξi。
即 {ξ (t)} 为严平稳随机过程。
7. 定义随机过程 x(t) = a − bt,其中,a~N (0, σ) 和 b~N (0, σ) 为独立的高斯随机变量,证
明其样本轨道与 t 轴的交点位于区间 (0, T ) 的概率为 (arctan T ) /π。
参考答案:
样本轨道与
t
轴的交点为
a ,则所求概率为
由全概率公式,有
() ( P (Br) = P (A) P (Br | A) + P A P Br
) A=
m ( 1 )r +
n ×1r
m+n 2 m+n
P (A | Br) =
P (A) P (Br P (Br)
| A)
=
m ( 1 )r
m m+n
(
12m)+rn+
2
n m+n
×1r
=
m m + n · 2r
质点反向移动,即有 n − m 次 ξi = −1。 ∑n
则 ηn = ξi = m + (n − m) (−1) = 2m − n = k
i=1
又各次游走是相互统计独立的,则
()
(
)
P (ηn = k) =
n pmqn−m = m
n
n+k 2
n+k n−k
p2q2
()
∑n 则 E [η (n) = k] = (2m − n)
参考答案:
(1) V = a 时,一条样本轨道为典型的正弦曲线。 2
(2) ξ (0) = 0,fξ(0)(x) = δ(x);ξ (π/2ω) = V ,其概率密度同 V 一样。
(π) ξ
4ω
=
V
√ 2
,
fξ(
π 4ω
)
(x)
=
√ 2 a
,
0
<
0, 其他
x
<
√a 2
() 5π
ξ 4ω
=
V
−
√ 2
清华大学电子工程系版权所有
1
概率论与随机过程 (2) ,homework1_intro_solution © 清华大学电子工程系
1. 袋中装有 m 只正品硬币,n 只次品硬币(次品硬币两面都有国徽)。在袋中任取一只,将 它掷 r 次。已知每次都得到国徽,问取得的硬币是正品的概率。
参考答案:
设 Br =“出现 r 次国徽面”,A =“任取一只是正品”。
(1) 随机过程 {X (n)} 的一条典型轨道如图 1.1 所示:
图 1.1
清华大学电子工程系版权所有
2
(2)
n
为奇数时,X
(n)
=
−1, 1,
硬币正面 ,即
P
(X (n)
=
1)
=
p (X(n)
=
−1)
=
1,
硬币反面
2
同理,n
为偶数时,X
(n)
=
1, −1,
硬币正面 ,即
P
(X (n)
=
1)
=
p (X(n)
(1) 试画出随机过程 {X (n)} 的典型样本轨道。 (2) 求随机过程 {X (n)} 的一维概率分布列。 (3) 对两时刻 n, n + k,求 X (n) 和 X (n + k) 的两维联合分布列,n = 1, 2, · · · , k =
1, 2, · · · 。
概念复习: 样本轨道,随机过程有限维分布列 参考答案:
T
常见问题: 忘了写范围。
6. 设有随机过程 ξ (t) = A cos(ωt + Θ),其中相位 Θ 是一个均匀分布于 (−π, π) 间的随机 变量,判断 ξ (t) 是否为严平稳过程。 概念复习:宽平稳与严平稳概念 参考答案:
P (ξ(t1) ⩽ x1, ξ(t2) ⩽ x2, · · · , ξ(tn) ⩽ xn)
b
P
{ 0
⩽
a b
⩽
T }。
利用二元随机变量函数的已有结论 ∫∞
fz(z) = |y| fxy(yz, y)dy
−∞
其中 z = x/y
以及 {fab(a, b) =
1
−
a2 +b2 2σ
} e 2πσ
经过运算可得 FC(c) =
1 2
+
1 π
arctan
(c)
故P
0
⩽
a b
⩽
T
= FC(T ) − FC(0) = (arctan T ) /π
7
,
fξ(
5π 4ω
)
(x)
=
√ 2 a
,
−
√a 2
0, 其他
<
x
<
0
清华大学电子工程系版权所有
6
常见问题:没有弄清楚随机过程中样本空间和参数 t 的关系,一旦 t 固定,得到了是一 个随机变量 X,其概率密度的写法为 p(X = x) = fX (x)。
清华大学电子工程系版权所有
参考文献
[1] 陆大纟金. 随机过程及其应用. 清华大学出版社, 1986. [2] 陆大纟金,张灏. 随机过程及其应用(第二版). 清华大学出版社, 2012.
参考答案:
考虑 t 时刻的随机变量 ξ(t),设其所处的锯齿波的起始点时刻为 t0,则随机变量 t0 ∼
U (t
−
T, t),注意到
ξ(t)
=
t−t0 T
A,ξ
(t)
是
t0
的单调函数,取值范围为
[0, A]:
1
11 1
pξ(t)(x) = p(t0) dξ(t)
=
T
A
=
, x ∈ [0, A] A
dt
=
−1)
=
1,
硬币反面
2
综上,P
(X (n)
=
1)
=
p (X(n)
=源自文库
−1)
=
1。 2
(3) 由于多次抛硬币的实验相互之间独立,则随机变量 x(n), x(n + k) 也相互独立,则
P (X(n) = 1, X(n + k) = 1) = P (X(n) = 1) P (X(n + k) = 1) = 1/4; P (X(n) = 1, X(n + k) = −1) = P (X(n) = 1) P (X(n + k) = −1) = 1/4; P (X(n) = −1, X(n + k) = 1) = P (X(n) = −1) P (X(n + k) = 1) = 1/4; P (X(n) = −1, X(n + k) = −1) = P (X(n) = −1) P (X(n + k) = −1) = 1/4;
cos t2
=
1 12
cos t1 cos t2
Comments:此题考察的是自相关函数和自协方差函数的定义。
5. ([1] 第一章习题 3) 设有一随机过程 ξ(t),它的样本函数为周期性的锯齿波。图 1.2 画出 了两个样本函数图。各样本函数具有同一形式的波形,其区别仅在于锯齿波的起点位置
清华大学电子工程系版权所有
=
1
∫π f (ωh + θ)dθ
2π −π
注意,周期函数 f (θ) 在任一个周期的积分都相等。所以
1
∫π f (θ)dθ =
1
∫π f (ωh + θ)dθ
2π −π
2π −π
所以,对 ∀n, ∀t1, t2, · · · , tn, ∀h 都有
P (ξ (t1) ⩽ x1, ξ (t2) ⩽ x2, · · · , ξ (tn) ⩽ xn) = P (ξ (t1 + h) ⩽ x1, ξ (t2 + h) ⩽ x2, · · · , ξ (tn + h) ⩽ xn)
= P (A cos (ωt1 + Θ) ⩽ x1, A cos (ωt2 + Θ) ⩽ x2, · · · , A cos (ωtn + Θ) ⩽ xn)
= =
1 2π 1
∫π ∫−ππ
p (A cos (ωt1 f (θ) dθ
+
θ)
⩽
x1, A cos (ωt2
+
θ)
⩽
x2,
···
, A cos (ωtn
Rξ (t1, t2) = E {η cos t1η cos t2} = E {η2} cos t1 cos t2
∫1
= η2dη cos t1 cos t2
0
E
{η}
= =
1 ∫3
cos t1
1
ηdη
cos t2 1
=
0
2
E {ξ} = E {η cos t} =
∫
1
η cos tdη
=
1
cos t
当 n1 > n2 时可得到类似的结论。
(3) 二维概率分布列
ξ(n1)
1
-1
ξ(n2)
1
p2
pq
-1
pq
q2
Rξξ (n1, n2) = p2 + q2 − 2pq
常见错误: 把随机过程 {η (n)} 的自相关函数和 {ξ(n)} 的搞混淆。认为 η (n1) 和 η (n2) 独立。
4. ([1] 第一章习题 7) 设有随机过程 {ξ (t) , −∞ < t < ∞},ξ (t) = η cos (t),其中 η 为均匀分 布于 (0, 1) 间的随机变量,求 {ξ (t)} 的自相关函数 Rξ (t1, t2),自协方差函数 Cξ (t1, t2)。 参考答案:
+
θ)
⩽
xn) dθ
2π −π
清华大学电子工程系版权所有
5
其中 f (θ) ≜ p (A cos (ωt1 + θ) ⩽ x1, A cos (ωt2 + θ) ⩽ x2, · · · , A cos (ωtn + θ) ⩽ xn),是 一个以 2π 为周期的函数。
∀h, P (ξ (t1 + h) ⩽ x1, ξ (t2 + h) ⩽ x2, · · · , ξ (tn + h) ⩽ xn)
i=1
(1) 求 {η (n)} 的均值函数。 (2) 求 {η (n)} 的自相关函数 Rηη(n1, n2)。 (3) 给定时刻 n1, n2,求随机过程 {ξ(n)} 的二维概率密度函数及相关函数。
概念复习:
离散状态离散时间随机过程数字特征
参考答案:
(1) 解法一:设在 n 次游走中有 m 次质点正向移动,即有 m 次 ξi = +1,有 n − m 次
8.([2] 第一章习题 1) 设随机过程 ξ(t) = V sin ωt,其中 ω 为常数,V 为服从 (0, a) 内均匀 分布的随机变量。
(1) 画出 ξ(t) 的某一条样本轨道。 (2) 求 ξ(0),ξ (π/4ω),ξ (π/2ω),ξ (5π/4ω) 的概率密度。
概念复习:
连续状态连续时间随机变量的一维分布
n
pmqn−m = pn − qn。
m=0
m
∑n
解法二:因各次游走是相互统计独立的,则 E [η (n)] = E[ξi] = (p − q)n。
i=1
清华大学电子工程系版权所有
3
(2) 假设 n1 < n2,
Rηη (n1, n2) = E [η (n1) η (n2)] = E {η (n1) [η (n1) + η (n2) − η (n1)]} = E[η (n1)]2 + E [η (n1)] E [η (n2) − η (n1)] = {E [η (n1)]}2 + V ar [η (n1)] + (p − q)2n1 (n2 − n1) = (p − q)2n21 + n1V ar [ξi] + (p − q)2n1 (n2 − n1) = (p − q)2n1n2 + n1[1 − (p − q)2]
Commets:此题考察的是条件概率的定义以及乘法公式。
2. 考虑一个如下定义的离散时间随机过程 X (n) , n = 1, 2, · · · 。无限次抛掷一枚硬币,对 n = 1, 2, · · · ,如果第 n 次抛掷结果为正面,则 X (n) = (−1)n;如果第 n 次抛掷结果为 反面,则 X (n) = (−1)n+1。
4
x W
$
W
D
t
7
x W
$
W
E
t
7
图 1.2
不同。设在 t = 0 后的第一个零值点位于 τ0,τ0 是一个随机变量,它在 (0, T ) 内均匀分
布,即
1/T fτ0 (t) = 0
(0 ⩽ t ⩽ T ) (其他值)
若锯齿波的幅度为 A,求随机过程 ξ(t) 的一维概率密度。