v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
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.
- Page 279 and 280: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 281 and 282: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 283 and 284: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 285 and 286: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 287 and 288: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 289 and 290: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 291 and 292: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 293 and 294: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 295 and 296: 4.5. CONSTRAINING CARDINALITY 295 m
- Page 297 and 298: 4.5. CONSTRAINING CARDINALITY 297 m
- Page 299 and 300: 4.5. CONSTRAINING CARDINALITY 299 a
- Page 301 and 302: 4.5. CONSTRAINING CARDINALITY 301 f
- Page 303 and 304: 4.5. CONSTRAINING CARDINALITY 303 n
- Page 305 and 306: 4.5. CONSTRAINING CARDINALITY 305 W
- Page 307 and 308: 4.5. CONSTRAINING CARDINALITY 307 t
- Page 309 and 310: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 311 and 312: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 313 and 314: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 315 and 316: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 317 and 318: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 319 and 320: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 321 and 322: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 323 and 324: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 325 and 326: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 327 and 328: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 329: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 333 and 334: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 335 and 336: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 337 and 338: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 339 and 340: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 341 and 342: 4.7. CONVEX ITERATION RANK-1 341 fi
- Page 343 and 344: 4.7. CONVEX ITERATION RANK-1 343 Gi
- Page 345 and 346: Chapter 5 Euclidean Distance Matrix
- Page 347 and 348: 5.2. FIRST METRIC PROPERTIES 347 co
- Page 349 and 350: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 351 and 352: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 353 and 354: 5.4. EDM DEFINITION 353 The collect
- Page 355 and 356: 5.4. EDM DEFINITION 355 5.4.2 Gram-
- Page 357 and 358: 5.4. EDM DEFINITION 357 D ∈ EDM N
- Page 359 and 360: 5.4. EDM DEFINITION 359 5.4.2.2.1 E
- Page 361 and 362: 5.4. EDM DEFINITION 361 ten affine
- Page 363 and 364: 5.4. EDM DEFINITION 363 spheres: Th
- Page 365 and 366: 5.4. EDM DEFINITION 365 By eliminat
- Page 367 and 368: 5.4. EDM DEFINITION 367 where Φ ij
- Page 369 and 370: 5.4. EDM DEFINITION 369 5.4.2.2.6 D
- Page 371 and 372: 5.4. EDM DEFINITION 371 10 5 ˇx 4
- Page 373 and 374: 5.4. EDM DEFINITION 373 corrected b
- Page 375 and 376: 5.4. EDM DEFINITION 375 by translat
- Page 377 and 378: 5.4. EDM DEFINITION 377 Crippen & H
- Page 379 and 380: 5.4. EDM DEFINITION 379 where ([√
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.