Decision Making Under Uncertainty: Solutions for ECE 313 Problem Set 7, Assignments of Statistics

Solutions to problem set 7 of the university of illinois at urbana-champaign's ece 313: decision-making under uncertainty course. The solutions cover various topics such as maximum-likelihood decision rules, bayesian decision rules, and likelihood ratios for different detection problems.

Typology: Assignments

Pre 2010

Uploaded on 02/24/2010

koofers-user-iva
koofers-user-iva 🇺🇸

10 documents

1 / 3

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
University of Illinois Fall 2009
ECE 313: Problem Set 7: Solutions
Decision-Making Under Uncertainty
1. [A warm-up exercise]
(a) The likelihood matrix Lis as shown below and the maximum-likelihood decision rule is indicated
shading.
Hypothesis X=0 X=1 X=2 X=3
H10.4 0.3 0.2 0.1
H00.1 0.2 0.3 0.4
It is easy to get PFA = sum of unshaded entries on H0row = 0.1+0.2=0.3
and PMD = sum of unshaded entries on H1row = 0.1+0.2=0.3 also.
(b) By the law of total probability, P(E)=π0PFA +π1PMD =0.7×0.3+0.3×0.3=0.3.
(c) The joint probability matrix J=LD where D= diag[π1,π0] is as shown below, and the Bayesian
decision rule is as indicated by the shading.
Hypothesis X=0 X=1 X=2 X=3
H10.12 0.09 0.06 0.03
H00.07 0.14 0.21 0.28
P(E) = sum of all the unshaded entries in the Jmatrix = 0.25 is smaller than the average error
probability for the ML rule, as it should be.
(d) In general, the Jmatrix can be written as
Hypothesis X=0 X=1 X=2 X=3
H10.4π10.3π10.2π10.1π1
H00.1π00.2π00.3π00.4π0
Now, if 0.1π0>0.4π1=0.4(1 π0), that is, π0>0.8, then the Bayesian decision is in favor of H0
whenever X= 0. But, if π0>0.8, then the same decision (favoring H0) holds whenever Xequals
1, 2, or 3 also. (Work it out!) Hence, if π0>0.8, the Bayesian decision always is in favor of H0.
From the symmetry of the problem, it should be obvious that if π1>0.8, the Bayesian decision
always is in favor of H1. The incredulous are invited to work out the details for themselves.
Alternatively, the likelihood ratio takes on values, 4, 3
2,2
3and 1
4, and hence always exceeds
π0
π1=π0
1π0if π0<0.2, that is, if π1>0.8, and can never exceed π0
π1=π0
1π0if π0>0.8. The
result to be proved follows from this.
2. [Which route did she fly?]
(a) If H0is the true hypothesis, then Xtakes on values 0, 1, 2 with probabilities (1 p)2,2p(1 p),
and p2respectively. If H1is the true hypothesis, then Xtakes on values 0, 2, 4 with probabilities
(1 q)2,2q(1 q), and q2respectively. The likelihood matrix is thus
Hypothesis X=0 X=1 X=2 X=4
H1(1 q)20 2q(1 q)q2
H0(1 p)22p(1 p)p20
(b) When X= 0, the ML decision is in favor of H1or H0according as (1 q)2(1 p)2, that is,
according as 1 q1p, or pq.
When X= 1, the ML decision is always in favor of H0regardless of the values of pand q.
When X= 2, the ML decision is in favor of H1or H0according as 2q(1 q)p2. Note that the
maximum value of 2q(1 q) is 1
2(prove this!) and so if p> 1
20.707 . . ., then the ML decision
will always be in favor of H0regardless of the value of q.
When X= 4, the ML decision is always in favor of H1regardless of the values of pand q.
pf3

Partial preview of the text

Download Decision Making Under Uncertainty: Solutions for ECE 313 Problem Set 7 and more Assignments Statistics in PDF only on Docsity!

University of Illinois Fall 2009

ECE 313: Problem Set 7: Solutions

Decision-Making Under Uncertainty

  1. [A warm-up exercise]

(a) The likelihood matrix L is as shown below and the maximum-likelihood decision rule is indicated shading. Hypothesis X = 0 X = 1 X = 2 X = 3 H 1 0. 4 0. 3 0. 2 0. 1 H 0 0. 1 0. 2 0. 3 0. 4 It is easy to get P (^) FA = sum of unshaded entries on H 0 row = 0.1 + 0.2 = 0. 3 and P (^) MD = sum of unshaded entries on H 1 row = 0.1 + 0.2 = 0.3 also. (b) By the law of total probability, P (E) = π 0 P (^) FA + π 1 P (^) MD = 0. 7 × 0 .3 + 0. 3 × 0 .3 = 0.3. (c) The joint probability matrix J = LD where D = diag[π 1 , π 0 ] is as shown below, and the Bayesian decision rule is as indicated by the shading. Hypothesis X = 0 X = 1 X = 2 X = 3 H 1 0. 12 0. 09 0. 06 0. 03 H 0 0. 07 0. 14 0. 21 0. 28 P (E) = sum of all the unshaded entries in the J matrix = 0.25 is smaller than the average error probability for the ML rule, as it should be. (d) In general, the J matrix can be written as Hypothesis X = 0 X = 1 X = 2 X = 3 H 1 0. 4 π 1 0. 3 π 1 0. 2 π 1 0. 1 π (^1) H 0 0. 1 π 0 0. 2 π 0 0. 3 π 0 0. 4 π (^0) Now, if 0. 1 π 0 > 0. 4 π 1 = 0.4(1 − π 0 ), that is, π 0 > 0 .8, then the Bayesian decision is in favor of H (^0) whenever X = 0. But, if π 0 > 0 .8, then the same decision (favoring H 0 ) holds whenever X equals 1, 2, or 3 also. (Work it out!) Hence, if π 0 > 0 .8, the Bayesian decision always is in favor of H 0. From the symmetry of the problem, it should be obvious that if π 1 > 0 .8, the Bayesian decision always is in favor of H 1. The incredulous are invited to work out the details for themselves. Alternatively, the likelihood ratio takes on values, 4, 32 , 23 and 14 , and hence always exceeds π (^0) π 1 =^

π (^0) 1 −π 0 if^ π^0 <^0 .2, that is, if^ π^1 >^0 .8, and can never exceed^

π 0 π 1 =^

π 0 1 −π 0 if^ π^0 >^0 .8.^ The result to be proved follows from this.

  1. [Which route did she fly?] (a) If H 0 is the true hypothesis, then X takes on values 0, 1, 2 with probabilities (1 − p)^2 , 2p(1 − p), and p 2 respectively. If H 1 is the true hypothesis, then X takes on values 0, 2, 4 with probabilities (1 − q)^2 , 2q(1 − q), and q 2 respectively. The likelihood matrix is thus Hypothesis X = 0 X = 1 X = 2 X = 4 H 1 (1 − q)^2 0 2 q(1 − q) q 2 H 0 (1 − p)^2 2 p(1 − p) p 2 0 (b) When X = 0, the ML decision is in favor of H 1 or H 0 according as (1 − q)^2 ≷ (1 − p)^2 , that is, according as 1 − q ≷ 1 − p, or p ≷ q. When X = 1, the ML decision is always in favor of H 0 regardless of the values of p and q. When X = 2, the ML decision is in favor of H 1 or H 0 according as 2q(1 − q) ≷ p 2. Note that the maximum value of 2q(1 − q) is 12 (prove this!) and so if p > √^12 ≈ 0. 707.. ., then the ML decision will always be in favor of H 0 regardless of the value of q. When X = 4, the ML decision is always in favor of H 1 regardless of the values of p and q.

(c) The entry 2p(1 − p) is always shaded on the H 0 row and is not included in the sum P (^) FA = sum of unshaded entries on H 0 row. Hence, P (^) FA ≤ 1 − 2 p(1 − p). Similarly, the entry q 2 is always shaded on the H 1 row and is not included in the sum P (^) MD = sum of unshaded entries on H 1 row. Hence, P (^) MD ≤ 1 − q 2. (d) The joint probability matrix is Hypothesis X = 0 X = 1 X = 2 X = 4 H 1 π 1 (1 − q)^2 0 2 q(1 − q)π 1 π 1 q 2 H 0 π 0 (1 − p)^2 2 p(1 − p)π 0 π 0 p 2 0 When X = 0, the Bayesian decision is in favor of H 1 or H 0 according as π 1 (1 − q)^2 ≷ π 0 (1 − p)^2 , that is, according as (^11) −−qp ≷

π 0 /π 1 or equivalently, p ≷ 1 − √^1 π− q 0 /π^1

When X = 1, the Bayesian decision is always in favor of H 0 regardless of the values of p and q. When X = 2, the Bayesian decision is n favor of H 1 or H 0 according as π 1 2 q(1 − q) ≷ π 0 p 2. When X = 4, the Bayesian decision is always in favor of H 1 regardless of the values of p and q.

  1. [Detection problem with Poisson distributed observations]

(a) The ML rule is to decide a ONE is sent if Λ(X) ≥ 1, where Λ is the likelihood ratio function, defined by

Λ(k) =

P (X = k | one is sent) P (X = k | zero is sent)

e−λ^1 λ k 1 /k! e−λ^0 λ k 0 /k!

λ (^1) λ (^0)

) (^) k e−(λ^1 −λ^0 )^ = 3k^ e−^4 ≈

3 k

  1. 6 Therefore, the ML decision rule is to decide a ONE is sent if X ≥ 4. (b) The MAP rule is to decide a ONE is sent if Λ(X) ≥ π π^0 1 , where Λ is the likelihood ratio found in part (a). So the MAP rule decides a ONE is sent if 3

X

  1. 6 ≥^5 ,^ or equivalently, if^ X^ ≥^5.^ Note that when X = 4, the ML rule decides that a ONE was transmitted, not a ZERO, but because ZEROes are so much more likely to be transmitted than ONEs, the MAP rule decides in favor of a ZERO in this case.

  2. [Detection problem with a geometric model]

(a) If X = n, the likelihood ratio has value

Λ(n) =

p 1 (1 − p 1 )n−^1 p 0 (1 − p 0 )n−^1

p (^1) p (^0)

1 − p (^1) 1 − p (^0)

) (^) n− 1

1 if (n − 1) ln

1 − p (^1) 1 − p (^0)

ln

p (^0) p (^1)

Since p 1 < p 0 , we have that 1−p 1 > 1 −p 0 and ln((1−p 1 )/(1−p 0 )) > 0. Therefore, the maximum likelihood decision rule is

“Decide that H 1 is the true hypothesis if X > 1 +

ln

p 0 p 1

ln

1 −p 1 1 −p 0

Note that there is less chance of a successful transmission on Route 1 than on Route 0, and hence large values of X should lead to the decision that Route 1 was used. (b) The minimum-error-probability (MEP) decision rule decides that H 1 is the true hypothesis if the likelihood ratio exceeds the threshold π 0 /π 1. Now Λ(1) = p 1 /p 0 < 1. Since 1 − p 1 > 1 − p 0 , we see that

Λ(n) =

p 1 (1 − p 1 )n−^1 p 0 (1 − p 0 )n−^1

p 1 (1 − p 1 )n−^2 p 0 (1 − p 0 )n−^2

1 − p (^1) 1 − p (^0)

= Λ(n − 1)

1 − p (^1) 1 − p (^0)

Λ(n − 1),

and so Λ(1) = p 1 /p 0 is the smallest value of the likelihood ratio. It follows that if π 0 /π 1 = π 0 /(1 − π 0 ) < p 1 /p 0 , that is, if π 0 < p 1 /(p 0 + p 1 ), the MEP decision rule will always decide that H 1 is the true hypothesis regardless of the observed value of X. On the other hand, since Λ(n) increases monotonically without bound as n increases, there is no value of π 0 < 1 for which π 0 /π 1 can be guaranteed to be larger than the likelihood ratio no matter what value X takes on. 2