v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

330 CHAPTER 4. SEMIDEFINITE PROGRAMMING Φ Figure 87: MRI radial sampling pattern in DC-centric Fourier domain representing 4.1% (10 lines) subsampled data. Only half of these complex samples, in theory, need be acquired for a real image because of conjugate symmetry. Due to MRI machine imperfections, samples are generally taken over full extent of each radial line segment. MRI acquisition time is proportional to number of lines. Radial sampling in the Fourier domain can be simulated by Hadamard product ◦ with a binary mask Φ∈ R n×n whose nonzero entries could, for example, correspond with the radial line segments in Figure 87. To make the mask Nyquist-centric, like DFT matrix F , define a permutation matrix for circular shift 4.46 Θ ∆ = [ 0 I I 0 ] ∈ S n (753) Then given subsampled Fourier domain (MRI k-space) measurements in incomplete K ∈ C n×n , we might constrain F(U) thus: and in vector form, (35) (1664) ΘΦΘ ◦ F UF = K (754) δ(vec ΘΦΘ)(F ⊗F ) vec U = vec K (755) Because measurements K are complex, there are actually twice the number of equality constraints as there are measurements. 4.46 Matlab fftshift().

4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 331 We can cut that number of constraints in half via vertical and horizontal mask Φ symmetry which forces the imaginary inverse transform to 0 : The inverse subsampled transform in matrix form is and in vector form later abbreviated where F H (ΘΦΘ ◦ F UF)F H = F H KF H (756) (F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) vec U = (F H ⊗F H ) vec K (757) P vec U = f (758) P ∆ = (F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) ∈ C n2 ×n 2 (759) Because of idempotence P = P 2 , P is a projection matrix. Because of its Hermitian symmetry [140, p.24] P = (F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) = (F ⊗F ) H δ(vec ΘΦΘ)(F H ⊗F H ) H = P H (760) P is an orthogonal projector. 4.47 P vec U is real when P is real; id est, when for positive even integer n [ ] Φ Φ = 11 Φ(1, 2:n)Ξ ∈ R n×n (761) ΞΦ(2:n, 1) ΞΦ(2:n, 2:n)Ξ where Ξ∈ S n−1 is the order-reversing permutation matrix (1608). In words, this necessary and sufficient condition on Φ (for a real inverse subsampled transform [246, p.53]) demands vertical symmetry about row n +1 and 2 horizontal symmetry 4.48 about column n+1. 2 Define ⎡ ⎤ 1 0 0 −1 1 0 ∆ = ∆ −1 1 . .. ∈ R ... ... ... n×n (762) ⎢ ⎣ ... ⎥ 1 0 ⎦ 0 T −1 1 4.47 (759) is a diagonalization of matrix P whose binary eigenvalues are δ(vec ΘΦΘ) while the corresponding eigenvectors constitute the columns of unitary matrix F H ⊗F H . 4.48 This condition on Φ applies to both DC- and Nyquist-centric DFT matrices.

4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 331<br />

We can cut that number of constraints in half via vertical and horizontal<br />

mask Φ symmetry which forces the imaginary inverse transform to 0 : The<br />

inverse subsampled transform in matrix form is<br />

and in vector form<br />

later abbreviated<br />

where<br />

F H (ΘΦΘ ◦ F UF)F H = F H KF H (756)<br />

(F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) vec U = (F H ⊗F H ) vec K (757)<br />

P vec U = f (758)<br />

P ∆ = (F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) ∈ C n2 ×n 2 (759)<br />

Because of idempotence P = P 2 , P is a projection matrix. Because of its<br />

Hermitian symmetry [140, p.24]<br />

P = (F H ⊗F H )δ(vec ΘΦΘ)(F ⊗F ) = (F ⊗F ) H δ(vec ΘΦΘ)(F H ⊗F H ) H = P H<br />

(760)<br />

P is an orthogonal projector. 4.47 P vec U is real when P is real; id est, when<br />

for positive even integer n<br />

[<br />

]<br />

Φ<br />

Φ = 11 Φ(1, 2:n)Ξ<br />

∈ R n×n (761)<br />

ΞΦ(2:n, 1) ΞΦ(2:n, 2:n)Ξ<br />

where Ξ∈ S n−1 is the order-reversing permutation matrix (1608). In words,<br />

this necessary and sufficient condition on Φ (for a real inverse subsampled<br />

transform [246, p.53]) demands vertical symmetry about row n +1 and<br />

2<br />

horizontal symmetry 4.48 about column n+1.<br />

2<br />

Define ⎡<br />

⎤<br />

1 0 0<br />

−1 1 0<br />

∆ =<br />

∆ −1 1 . ..<br />

∈ R ...<br />

...<br />

...<br />

n×n (762)<br />

⎢<br />

⎣<br />

...<br />

⎥<br />

1 0 ⎦<br />

0 T −1 1<br />

4.47 (759) is a diagonalization of matrix P whose binary eigenvalues are δ(vec ΘΦΘ) while<br />

the corresponding eigenvectors constitute the columns of unitary matrix F H ⊗F H .<br />

4.48 This condition on Φ applies to both DC- and Nyquist-centric DFT matrices.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!