MUSIC as described so far is what I would call an ``imaging method''
in the sense that it is attempting to locate objects, but not to
quantify their properties (such as the scattering strength *q*).
We will show now that, with only a little more computational effort, we
can also determine the scattering coefficients. The resulting
algorithm will then be distinguished as an ``inversion method,''
rather than imaging.

Our data are again the matrix elements of the response matrix . We know that can be represented equally in two ways by its scattering expansion and its SVD so that

**K** = _n=1^N q_n**H**_n **H**_n^T
= _m=1^M _m**V**_m **V**_m^T,
recalling that there are *M* sensors in our array, *N* scatterers,
and that the singular vectors are assumed to be normalized to unity.

One formulation of the inversion problem that is very similar in spirit to MUSIC is therefore to find a matrix

#<299#>**K** = _r=1^L^2 q_r**H**_r **H**_r^T,
whose elements agree with those of the measured according to
some criterion such as least-squares. The sum in (Ktilde) is
taken over grid points in an (pixel units) model space in
2D. Generalization to model spaces having other
shapes and/or to 3D is straightforward, as there is no use made here
of any symmetries in the choice of sensor array arrangement, or of
the locations of the test pixels chosen as possible targets and
indexed by *r*. We will want to take advantage
of the orthogonality of the singular vectors . So we
suppose that it is possible to solve an equation of the form

_r=1^L^2 q_r**H**_r **H**_r^T
_m=1^M _m**V**_m **V**_m^T,
in the least-squares sense for the coefficients associated with
some set of points (pixel or voxel centers) in the model space.

Applying the singular vectors to both the right and left sides of (KeqK), we find

_r=1^L^2 q_r(**H**_r^T**
V**^*_m)^2
= _m,
which can then be rewritten in matrix notation as

(**H**_1^T**V**^*_1)^2 & ...& (**
H**_L^2^T**V**^*_1)^2
(**H**_1^T**V**^*_2)^2 & ...& (**
H**_L^2^T**V**^*_2)^2
& &
(**H**_1^T**V**^*_M)^2 & ...& (**
H**_L^2^T**V**^*_M)^2
q_1 q_2 q_L^2=
_1 _2 _M.
We see that estimates of the coefficients can now be found
by solving this rather large () and
underdetermined least-squares system. This form
is similar to various kinds of tomography, with the singular values
taking the role of the data, and the singular vectors taking the
role of the discrete scanning modality. (For comparison,
in traveltime tomography the data are the traveltimes, the matrix
elements are the ray-path lengths through the cells, and
the model parameters are the slownesses in the cells.) Note that the elements of
the present matrix are in general complex, as are the scattering coefficients
, while the singular values are all
nonnegative real numbers.

The computational situation can be improved significantly by noting the similarity of the matrix elements to the direction cosines already treated in MUSIC. In fact, if we take the magnitude of each matrix element, we have the square of a direction cosine. Furthermore, if we sum these magnitudes along the columns of the matrix, then we have obtained exactly the quantity

**H**_r^2_m=1^M^2(**V**_m,**H**_r),
which -- except for the normalization factor --
is the sum over all the direction cosines associated with the
vector . This column sum is therefore a measure of
what we might call the ``coverage'' of the pixel centered at .If there is no coverage, this column sum is zero and the pixel
does not contain a measurable target.
If the normalized coverage is close to unity, then this pixel is
one that MUSIC will clearly identify as a target. For intermediate
values, the situation is somewhat ambiguous, as it is in the normal
MUSIC algorithm, but clearly the closer the normalized sum is to
unity, the more likely it is that the pixel contains a target.

Now we have a clear pixel selection criterion available, based on these column
sums. Thus, we can reduce the size of the inversion problem posed in
(bigmatrixeqn) very quickly by throwing out all pixels that have
small column sums, and/or, equivalently, those whose values are
. We can classify the remaining pixels
by the magnitudes of the normalized sums, choosing to keep only (say)
those pixels having the *M* largest normalized column sums. Or if there
are clearly only *N* << *M* singular values that are nonzero, then we could
reduce the problem still further, and keep only those pixels having the
*N* largest normalized column sums. If we have an matrix when we
have finished this selection process, then we can simply invert the
matrix to find the values of the remaining 's.
If on the other hand, we have an matrix that is not
square remaining after this process of elimination, then we will again have
to solve the inversion problem, by using overdetermined least-squares.

Another possibility when the singular values have a gradual fall off in
magnitude as , but no precipitous drop to zero, is to
multiply (bigmatrixeqn) on the left by a diagonal matrix whose
nonzero elements are the singular values to some power *p*.
Then, the resulting equation is

_r=1^L^2 q_r_m^p(**H**_r^T**
V**^*_m)^2
= _m^p+1,
or, in vector/matrix notation,

_1^p(**H**_1^T**V**^*_1)^2 & ...& _1^p(**
H**_L^2^T**V**^*_1)^2
_2^p(**H**_1^T**V**^*_2)^2 & ...& _2^p(**
H**_L^2^T**V**^*_2)^2
& &
_M^p(**H**_1^T**V**^*_M)^2 & ...& _M^p(**
H**_L^2^T**V**^*_M)^2
q_1 q_2 q_L^2=
_1^p+1 _2^p+1 _M^p+1.
We can apply a MUSIC-like processing scheme to this matrix (that
could then be called ``weighted MUSIC,'' which is similar to some methods
suggested by other authors).
This approach permits the singular values to determine in a natural
way the cutoff (if any) in the contributions from those singular vectors
of negligible importance to the inversion. For example, when *p*=2,
these column sums of the magnitudes of these elements
are just , which are the matrix elements
of the time-reversal operator with each vector .

9/18/2001