



Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Solutions to problem 1 and 2 of the ece544 homework 2. It includes the estimation of mean and correlation matrices for given images, determination of appropriate models for each image based on correlation matrices, and the design of an energy function for favoring specific configurations in problem 2. The document also mentions the use of the gibbs sampler for simulating ising and 8-level potts models in problems 3 and 4.
Typology: Assignments
1 / 6
This page cannot be seen from the preview
Don't miss anything!




Problem 1
(a) The mean and 4 × 4 correlation matrices can be estimated as follows:
μˆ =
n
x(n)
Rˆ(k) =
n
x(n)x(n + k), k ∈ { 0 , 1 , 2 , 3 }
2 .
We obtain the following values for the five images considered. In each case, a least-
squares fit could be used to determine whether ρ
‖k‖ (isotropic model) or ρ
|k 1 |+|k 2 | (sepa-
rable model) is more appropriate.
Elaine: ˆμ = 136.3565, Rˆ =
Man: ˆμ = 89.0075, ˆR =
Aerial: ˆμ = 180.5718, Rˆ =
Clock: ˆμ = 185.9803, Rˆ =
Couple: ˆμ = 123.1772, Rˆ =
(b) This is a separable AR(1) process:
x(n 1 , n 2 ) = 100 + y(n 1 , n 2 )
where
y(n 1 , n 2 ) = 1 − ρy(n 1 − 1 , n 2 ) − ρy(n 1 − 1 , n 2 ) + ρ
2 y(n 1 − 1 , n 2 − 1) + e(n 1 , n 2 )
and e(n 1 , n 2 ) is white noise with mean 0 and variance 30(1 − ρ
2 )
2
. Realizations of this
process may be obtained by initializing the boundary of the image x (e.g., with zeroes)
and recursively applying the formula above. A more efficient implementation would
exploit the separability nature of the process by implementing horizontal filtering in a
first step, and vertical filtering in a second step:
z(n 1 , n 2 ) = ρz(n 1 − 1 , n 2 ) + e(n 1 , n 2 )
y(n 1 , n 2 ) = ρy(n 1 , n 2 − 1) + z(n 1 , n 2 )
(c) For ρ = 0. 1 , 0. 5 , 0 .95 we obtain the following realizations:
ρ = 0. 1 ρ = 0. 5 ρ = 0. 95
Problem 3: Simulation of Ising model.
After 100 iterations of the Gibbs sampler, we obtain the following configurations:
− 1 = 0. 1 T
− 1 = 0. 5 T
− 1 = 0. 88
− 1 = 1. 5 T
− 1 = 10
Observe that clusters become larger as the temperature parameter T decreases.
Problem 4: Simulation of 8-level Potts model.
After 100 iterations of the Gibbs sampler, we obtain the following configurations:
− 1 = 0. 1 T
− 1 = 0. 5 T
− 1 = 0. 88
− 1 = 1. 5 T
− 1 = 10
Observe that clusters become larger as the temperature parameter T decreases.