1. Introduction
Many articles and reports at conferences in the USA, China, France, Italy, Great Britain, and Russia are devoted to achieving angular super-resolution. The methods of MUSIC (multiple signal classification) [
1,
2,
3], ESPRIT [
4], the deconvolution method [
5,
6], the maximum entropy method (MME) of Berg [
7], Borgiotti–Lagunas [
8], Capon [
9], the maximum likelihood method [
10], and others [
11] have been studied in detail. New methods are being developed for solving problems of forming images of objects in the radio range with inaccurately specified information due to the presence of random components of measurement data [
12]. These methods successfully solve one-dimensional problems. The performances of most algorithms based on them are not enough to work in real-time modes. Most methods (including the well-known MUSIC and ESPRIT) are applicable only for narrow-band signals and are ineffective when using broadband signals. This article proposes a new method for achieving super-resolution. This method allows for increasing the achievable level of the angular resolution.
The developed algebraic methods and their varieties [
13,
14,
15,
16,
17,
18,
19,
20] are free from the listed disadvantages, and allow obtaining solutions with super-resolution when using wide- and ultra-wide-band signals as well as solving two-dimensional problems [
21]. Numerical experiments have shown that algebraic methods have significantly higher noise immunity [
17,
18,
19].
The most important property of these methods involves the use of a priori information of the source signals, which allows increasing the stability of solutions to which the super-resolution problem reduces. Algebraic methods are quite simple, which makes it possible to implement high-speed algorithms even when relatively cheap computing devices are used. This allows, unlike many other methods [
22,
23,
24,
25,
26,
27,
28], to apply them in real-time modes.
After the introduction, the article considers the formulation of the inverse problem of achieving a super-solution. The following is a brief justification of the algebraic method for solving the problem with additional conditions. The next section sets and solves the problem of increasing the degree of super-resolution based on a new mathematical method of digital signal processing for the antenna array.
2. Problem Statement
The achievement of angular super-resolution is based on solving inverse problems described by one- and two-dimensional Fredholm integral equations of the first kind of convolution type
where
is the angular distribution of the reflected signal amplitude,
is the two-dimensional radiation pattern (RP) of the antenna system,
is the envelope of the received signal dependence when scanning the viewing sector
. During an active location,
should be considered the product of the transmitting and receiving radiation patterns, while at passive and semi-active locations,
is the receiving RP. The task is to find an approximate solution of the integral Equation (
1) with respect to
, i.e., it reduces to the solution of the inverse problem.
In direct measuring, signal describes a real distribution with an angular resolution on the angles and not exceeding the Rayleigh criterion, i.e., , where is the wavelength, D is the size of the antenna system. The task is to find an algorithm for digital signal processing that provides with a resolution better than , i.e., with super-resolution.
Considering
as an unknown function and
as the given one, we arrive at (
1) in the form of the Fredholm equation of the first kind of convolution type. The task is inverse and ill-posed according to Hadamard.
3. An Algebraic Method for Finding an Approximate Solution
We give a brief description of an algebraic method for achieving super-resolution in one- and two-dimensional problems of radar, remote sensing, and hydroacoustics [
12,
13,
14,
15,
16,
17] based on an approximate solution (
1). The RP of the antenna system, the received signal
, and the angular distribution of the reflected signal
are described by physically realizable functions from the Hilbert space
of quadratically integrable functions in a bounded two-dimensional domain
. The boundaries of the domain
are selected when constructing a solution for each specific problem, based on the analysis of the values of
, which decreases to almost zero when angles are removed from the values corresponding to the maximum. The values of
at the boundaries of the selected areas are naturally considered zero.
The solution
is represented as a decomposition according to the chosen complete orthonormal system of
functions in the area of the source location
.
An important question is the choice of the system of functions. The
system is determined based on a priori information about the solution, as well as reasonable assumptions about the properties of the solution. For example, in the remote sensing problem, if it is known that reflection occurs from a smoothly inhomogeneous surface with the possible presence of highly reflective small-sized inclusions, it is advisable to use wavelets. If the radar object under study is a collection of individual objects with small angular dimensions that are not resolved by direct observation, step functions should be used, etc. Thus, the problem is mathematically reduced to solving the inverse problem (
1) with an additional condition: an approximate solution is sought on sets of functions defined by sets of introduced basic functions
. To increase the resolution by
N times,
N basic functions should be used. The number of selected functions is limited by the stability of the solutions obtained and may vary with different approaches.
If the area
corresponds to the width of the beam, then the Rayleigh criterion corresponds to one function in (
2), then for
, we obtain a solution with super-resolution. The degree of the desired super-resolution increases with increasing the number of functions used
M in (
2). The described approach allows parameterization of inverse problems and reduces them to the search for expansion coefficients
from (
2). Indeed, using (
2), we obtain
where
The
coefficients can be found by minimizing the standard deviation
of the final sum on the right side (
3) of the signal under study
in the scanning sector, i.e., from the set of equations
or in matrix notation, by entering the
-vector of coefficients
, and
As a result, the inverse problem is reduced to solving a system of linear algebraic equations (SLAE) (
5)–(
7).
The choice of the system of functions used to represent the solution is carried out on the basis of a priori information about the solution or a mathematical model of the object under study is used.
The main feature of Equation (
7) is its poor conditionality because the unstable inverse problem is solved. The conditioning number of matrices
—as experiments on mathematical models show—increase exponentially with the increasing
N. Considering this factor, the solutions of SLAE (
7) are constructed as iterative processes with a sequential increase in the number
M of functions used in (
2). This makes it possible to approach the maximum achievable level of angular super-resolution for a specific task. The iterative process is convergent since it is easy to show that adding a new term in (
3) improves the approximation to
U.
4. Over-Resolution by Digital Antenna Arrays
Numerical experiments have shown that the SLAE of form (
7), due to the poor conditionality, makes it possible to obtain 2–3, under favorable conditions; 5–7
functions in the representation (
2) without significant distortion of the real distribution
, which exceeds the Rayleigh criterion by 2–8 times.
We justify the possibility of further increasing the level of the achieved super-resolution by modern systems based on digital antenna arrays (DAA) using a new modification of the algebraic method.
Consider a linear equidistant DAA of
elements with a distance
d between neighboring emitters. We consider the amplitude RP of individual DAA emitters to be weakly directional and equal to a constant within the scanning area. Then the radiation pattern of the antenna array in the direction
is
where
and
are the complex amplitudes of the current on the
n-th emitter
with accuracy up to a constant equal in the case under consideration to
. The angles
and
lay in the same plane with the linear DAA and are counted from the direction perpendicular to it.
The signal (
1) received when scanning along the angle
for a DAA with a directional diagram in form (
8) is the sum of
with
Here, is the signal received by the individual emitter of the DAA.
It is known that, in this case, solution
up to the normalizing factor can be represented in the form:
where
is an arbitrary function orthogonal to all exponents (
9) on the segment
.
In contrast to the above-known method (
3)–(
7), we will use the fact that in DAA, signals received by each emitter can be represented digitally. Using (
10), the solution representation in form (
2), and taking into account that signal
C is a complex conjugate,
, we obtain
Next, by entering a rectangular matrix
with elements
we come to a new SLAE with respect to the vector of unknown expansion coefficients
With the same dimension, the conditionality matrix numbers
from (
7) are significantly larger than those of matrices
from (
14). Greater stability of the system (
14) solutions allows increasing the number of functions in the solution (
2) representation. This makes it possible to increase the degree of super-resolution, as well as successfully process signals with higher levels of random components.
The solution of the entire problem is constructed in the form of an iterative process with sequentially increasing
M—the number of
functions in (
2) and, hence, in (
12)–(
14).
At the beginning of the process, at
, the system (
14) is reconditioned. In this case, as is usually recommended, a pseudo-solution is sought—a vector that minimizes the norm of the discrepancy of the system of equations. Thus, the problem with solving SLAE (
14) is replaced by the problem of finding the global minimum of the function
, where
is the norm of vector
. The norm is the sum of squares of vector components, and the described pseudo-solution search is an application of the least squares method.
For
, we obtain a normal system. The SLAE solution is carried out by the numerical LU decomposition method (lower–upper decomposition by Tadeusz Banachiewicz, 1938) [
28], which is based on the Gaussian sequential exclusion algorithm. In this case, the results are accurate up to numerical rounding errors.
An important advantage of the method lies in its speed because for the most time-consuming operations—calculating the elements of the matrices —the use of the received signal (vector ) is not required. When solving inverse problems of image reconstruction with the same matrix , but with a different , it is enough to calculate the LU decomposition of . The speed increases even more if there is a database of a specific measuring system in the form of LU decompositions of for the most commonly used systems of functions. The application of this approach makes it possible to ensure the operation of autonomous robotic systems and real-time remote control systems.
The achieved level of super-resolution is limited by the instability of SLAE solutions (
14). With increasing matrix
dimensions, the problems become incorrect, since negligibly small rounding errors in the calculation process lead to unacceptable errors in solutions.
An additional increase in the level of the achieved super-resolution is possible when carrying out the regularization of problems based, for example, on more complete a priori information about the solution, which is often available in real problems [
17,
19].
For
, the system (
14) is underdetermined. To obtain the only adequate solution to the problem, it is necessary to redefine it by adding some a priori data about the solution, i.e., about the unknown vector
.
Usually, the solution of an underdetermined SLAE (
14) is considered to be a normal pseudo-solution, i.e., one solution that has a minimum norm
. This solution is found when minimizing the norm
on a previously obtained family of SLAE solutions.
Generalization of the developed approach (
8)–(
14) to two-dimensional problems does not cause fundamental difficulties.
5. Numerical Results
In a mathematical model, during numerical experiments, the results of image reconstruction with conventional flat antenna arrays by the algebraic method (
3)–(
7) and the DAA of the same dimension were compared using the presented modification of the algebraic method. The classical problem of the angular resolutions of two point objects located within the beam width of DAA
, i.e., not resolved by direct observation, was investigated.
Figure 1 shows the results of the image reconstructions for two identical signal sources shown as a thin solid line (Chart 1). The angular distance between the sources was
. Step functions were chosen to represent solutions as
functions. The legends in the result figures are given in the text. The dashed polylines (Chart 2) here and in the following figures show the solutions found by the usual algebraic methods (
3)–(
7). Two objects were resolved, and the localization of each object was
of the beam width
.
A solid bold polyline (Chart 3) shows the solution obtained on the basis of a modification of the algebraic method (
8)–(
14). The achieved localization of two objects was
with almost accurate determination of their angular coordinates.
When trying to improve localization, i.e., when using more than
M functions in (
2), distortions appear in the solutions in the forms of false sources, the amplitude values and their numbers increase rapidly with the increase in the number of
functions used.
When the objects under study approach, the quality of the approximate solutions obtained begins to deteriorate.
Figure 2 shows the found images
at a distance of
in the same notation as in
Figure 1. For illustration, an upper curve (Chart 4) is added showing the change in the amplitude of a received signal
when scanning the sector. Both methods still allow one to resolve objects. The localization of objects has not changed, however, when using the algebraic method (
3)–(
7), a false source with an amplitude up to
of the amplitude of the true source appeared.
With the further convergence of objects, more functions are required for an adequate representation of the solution, and, consequently, the stability of solutions decreases.
Figure 3 shows the solutions obtained with the highest possible resolution by the two methods compared with a very small angular distance between objects of
(a solid thin polyline, Chart 1). The algebraic method (
3)–(
7) did not allow resolving objects (dotted polyline, Chart 2).
The modification of the method (
8)–(
14), which allows using an underdetermined system, made it possible to resolve objects based on a normal pseudo-solution, providing a ten-fold excess of the Rayleigh criterion (bold polyline, Chart 3). The resulting error in determining the angular position of the objects was
.
Thus, DAA-based systems, due to a new signal processing method, can significantly increase the level of angular super-resolution compared to known methods with satisfactory image quality. The proposed modification makes it possible to successfully restore images of complex objects with super-resolution, including smoothly inhomogeneous regions and areas with large gradients of the amplitude distribution of the reflected signal.
Thus,
Figure 4 shows the solutions with the maximally achieved resolution obtained by the algebraic method and its newly presented modification. The true distribution
is presented in the form of a solid thin curve (Chart 1); in the form of a dotted curve (Chart 2)—the solution found by the algebraic method is found, in the form of a bold curve (Chart 3)—the solution for the DAA based on (
10)–(
14) is found. If the first solution as a whole correctly reflects the general averaged form of the distribution
, then the second allows one to highlight the details of the image, albeit with slight distortions.