Markov Chain Monte Carlo (MCMC) Methods: Example 11 (Matlab)
Markov Chain Monte Carlo (MCMC) Methods: Example 11 (Matlab)
Gibbs Sampler
Example 11 (Matlab)
Consider again deblurring of the constellation image (b) in
of Example 7 in which ` 1 -norm based prior
the case (ii)P
πpr(x) ∝ exp( ni=1 −α|xi |) was used resulting to a posterior of
the form
n
1 X
π(x | y ) ∝ exp − 2 ky − Axk2 − α |xi | .
2σ
i=1
Find using Gibbs sampler a conditional mean (CM) estimate
with σ = 0.005 and α = 20 and compare the result to the
MAP estimate computed in Example 7. In order to speed up
the sampling procedure, restrict the problem first into a
region of interest (ROI), referring here to a set in which the
pixels a priori seem to be nonzero based on the blurred
image.
Markov chain Monte Carlo (MCMC) methods
Gibbs Sampler
1 T 1 ασ 2 − yeiT ai
= −α + 2
yei ai − 2 aiT ai xi = 0, xi = ,
aiT ai
i.e.
σ σ
Example 12 (Matlab)
Consider deblurring of the constellation image in the case of
the conditionally gaussian (hierarchical) prior model
1
σ −1 A
σ −1 y
2
= exp −
x −
2 Dθ−1/2 0
(k−1)
x = (σ −2 AT A + Dθ−1
(k−1) ) (σ A y + σ −1 AT z1 + Dθ−1/2
−1 −1 T
(k−1) z2 )
= (AT A + σ 2 Dθ−1
(k−1) )
−1
(σAT y + σAT z1 + σ 2 Dθ−1/2
(k−1) z2 ),
Gamma Inv. G.
Markov chain Monte Carlo (MCMC) methods
Gibbs Sampler
Gamma Hyperprior
Gamma Hyperprior