Sign up & Download
Sign in

Convex Optimization

by Stephen Boyd, Lieven Vandenberghe
Optimization Methods and Software (2010)

Abstract

We are developing a dual panel breast-dedicated PET system using LSO scintillators coupled to position sensitive avalanche photodiodes (PSAPD). The charge output is amplified and read using NOVA RENA-3 ASICs. This paper shows that the coincidence timing resolution of the RENA-3 ASIC can be improved using certain list-mode calibrations. We treat the calibration problem as a convex optimization problem and use the RENA-3s analog-based timing system to correct the measured data for time dispersion effects from correlated noise, PSAPD signal delays and varying signal amplitudes. The direct solution to the optimization problem involves a matrix inversion that grows order (n3) with the number of parameters. An iterative method using single-coordinate descent to approximate the inversion grows order (n). The inversion does not need to run to convergence, since any gains at high iteration number will be low compared to noise amplification. The system calibration method is demonstrated with measured pulser data as well as with two LSO-PSAPD detectors in electronic coincidence. After applying the algorithm, the 511keV photopeak paired coincidence time resolution from the LSO-PSAPD detectors under study improved by 57%, from the raw value of 16.30.07 ns FWHM to 6.920.02 ns FWHM (11.520.05 ns to 4.890.02 ns for unpaired photons).

Cite this document (BETA)

Available from www.informaworld.com
Page 189
hidden

Convex Optimization

4.7 Vector optimization 175
objective function. We denote a general vector optimization problem as
minimize (with respect to K) f0(x)
subject to fi(x) ≤ 0, i = 1, . . . ,m
hi(x) = 0, i = 1, . . . , p.
(4.56)
Here x ∈ Rn is the optimization variable, K ⊆ Rq is a proper cone, f0 : Rn → Rq
is the objective function, fi : Rn → R are the inequality constraint functions, and
hi : Rn → R are the equality constraint functions. The only difference between this
problem and the standard optimization problem (4.1) is that here, the objective
function takes values in Rq, and the problem specification includes a proper cone
K, which is used to compare objective values. In the context of vector optimization,
the standard optimization problem (4.1) is sometimes called a scalar optimization
problem.
We say the vector optimization problem (4.56) is a convex vector optimization
problem if the objective function f0 isK-convex, the inequality constraint functions
f1, . . . , fm are convex, and the equality constraint functions h1, . . . , hp are affine.
(As in the scalar case, we usually express the equality constraints as Ax = b, where
A ∈ Rp×n.)
What meaning can we give to the vector optimization problem (4.56)? Suppose
x and y are two feasible points (i.e., they satisfy the constraints). Their associated
objective values, f0(x) and f0(y), are to be compared using the generalized inequal-
ity K . We interpret f0(x) K f0(y) as meaning that x is ‘better than or equal’ in
value to y (as judged by the objective f0, with respect to K). The confusing aspect
of vector optimization is that the two objective values f0(x) and f0(y) need not be
comparable; we can have neither f0(x) K f0(y) nor f0(y) K f0(x), i.e., neither
is better than the other. This cannot happen in a scalar objective optimization
problem.
4.7.2 Optimal points and values
We first consider a special case, in which the meaning of the vector optimization
problem is clear. Consider the set of objective values of feasible points,
O = {f0(x) | ∃x ∈ D, fi(x) ≤ 0, i = 1, . . . ,m, hi(x) = 0, i = 1, . . . , p} ⊆ Rq,
which is called the set of achievable objective values. If this set has a minimum
element (see §2.4.2), i.e., there is a feasible x such that f0(x) K f0(y) for all
feasible y, then we say x is optimal for the problem (4.56), and refer to f0(x) as
the optimal value of the problem. (When a vector optimization problem has an
optimal value, it is unique.) If x? is an optimal point, then f0(x?), the objective
at x?, can be compared to the objective at every other feasible point, and is better
than or equal to it. Roughly speaking, x? is unambiguously a best choice for x,
among feasible points.
A point x? is optimal if and only if it is feasible and
O ⊆ f0(x?) +K (4.57)
Page 303
hidden
Part II
Applications

Sign up today - FREE

Mendeley saves you time finding and organizing research. Learn more

  • All your research in one place
  • Add and import papers easily
  • Access it anywhere, anytime

Start using Mendeley in seconds!

Already have an account? Sign in

Readership Statistics

1116 Readers on Mendeley
by Discipline
 
 
 
by Academic Status
 
47% Ph.D. Student
 
10% Student (Master)
 
8% Post Doc
by Country
 
32% United States
 
8% Germany
 
6% China