Next Article in Journal
Risk-Sensitive Maximum Principle for Controlled System with Delay
Next Article in Special Issue
Min–Max Dynamic Programming Control for Systems with Uncertain Mathematical Models via Differential Neural Network Bellman’s Function Approximation
Previous Article in Journal
Programming for High-Performance Computing on Edge Accelerators
Previous Article in Special Issue
Dynamics of Evolutionary Differential Equations with Several Spatial Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algebraic Methods for Achieving Super-Resolution by Digital Antenna Arrays

by
Boris A. Lagovsky
1 and
Evgeny Ya. Rubinovich
2,*
1
Department of Applied Mathematics of Russian Technological University, 119454 Moscow, Russia
2
Trapeznikov Institute of Control Sciences of RAS, 117997 Moscow, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(4), 1056; https://doi.org/10.3390/math11041056
Submission received: 3 December 2022 / Revised: 3 February 2023 / Accepted: 10 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue Dynamics and Control Theory with Applications)

Abstract

:
The actual modern problem of developing and improving measurement and observation systems (including robotic ones) is to increase the volume and quality of the information received. Increasing the angle resolution to values significantly exceeding the Rayleigh criterion, i.e. achieving super-resolution is one of important ways to solve the problem. Angular super-resolution which makes it possible to detail images of research objects and their individual fragments, improves the quality of solutions to detection, recognition and identification problems, increases the range of such systems. In many papers methods developed by authors to achieve a super-resolution based on approximate solutions of inverse problems in the form of Fredholm integral equation of the first kind of convolution type called algebraic are presented. The methods used, as well as their varieties, make it possible to reduce solutions of inverse problems posed to solving sets of linear algebraic equations (SLAE). This paper presents results of further improvement of algebraic methods based on intelligent analysis of received signals. It is shown that their use in systems based on digital antenna arrays makes it possible to increase the achieved degree of exceeding the Rayleigh criterion. In the course of numerical experiments with a mathematical model, the stability of the solutions obtained and their adequacy were confirmed. The numerical results obtained open the following possibilities: (1) obtaining images of studied objects with a resolution exceeding the Rayleigh criterion by 4 to 10 times, (2) determining the angular coordinates of individual small-sized objects as part of multi-element complex objects (group targets), (3) clarifying boundaries of extended objects and their individual elements, (4) localizing individual bright objects on a smoothly inhomogeneous reflective background. Applying presented new methods does not require a significant computing power, what allows you to work in a real time mode using relatively simple and inexpensive computing devices. The ways of further improvement of presented algebraic methods for solving applied inverse problems are described.

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
U ( α , φ ) = Ω f ( α α , φ φ ) I ( α , φ ) d α d φ ,
where I ( α , φ ) is the angular distribution of the reflected signal amplitude, f ( α , φ ) is the two-dimensional radiation pattern (RP) of the antenna system, U ( α , φ ) is the envelope of the received signal dependence when scanning the viewing sector Ω . During an active location, f ( α , φ ) should be considered the product of the transmitting and receiving radiation patterns, while at passive and semi-active locations, f ( α , φ ) is the receiving RP. The task is to find an approximate solution of the integral Equation (1) with respect to I ( α , φ ) , i.e., it reduces to the solution of the inverse problem.
In direct measuring, signal U ( α , φ ) describes a real distribution I ( α , φ ) with an angular resolution on the angles α and φ not exceeding the Rayleigh criterion, i.e., θ r = λ / D , where λ is the wavelength, D is the size of the antenna system. The task is to find an algorithm for digital signal processing U ( α , φ ) that provides I ( α , φ ) with a resolution better than θ r , i.e., with super-resolution.
Considering I ( α , φ ) as an unknown function and U ( α , φ ) 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 U ( α , φ ) , and the angular distribution of the reflected signal I ( α , φ ) are described by physically realizable functions from the Hilbert space L 2 Ω 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 U ( α , φ ) , which decreases to almost zero when angles are removed from the values corresponding to the maximum. The values of I ( α , φ ) at the boundaries of the selected areas are naturally considered zero.
The solution I ( α , φ ) is represented as a decomposition according to the chosen complete orthonormal system of g m ( α , φ ) functions in the area of the source location Ω .
I ( α , φ ) = m = 1 b m g m ( α , φ ) m = 1 M b m g m ( α , φ ) .
An important question is the choice of the system of functions. The g m ( α , φ ) 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 g m ( α , φ ) . 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 M > 1 , 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 b m , m = 1 , 2 , , M from (2). Indeed, using (2), we obtain
U ( α , φ ) = m = 1 b m h m ( α , φ ) m = 1 M b m h m ( α , φ ) ,
where
h m ( α , φ ) = Ω f ( α α , φ φ ) g m ( α , φ ) d α d φ .
The b m coefficients can be found by minimizing the standard deviation δ 2 of the final sum on the right side (3) of the signal under study U ( α , φ ) in the scanning sector, i.e., from the set of equations
Ω U ( α , φ ) h j ( α , φ ) h m ( α , φ ) d α d φ m = 1 M b m Ω h j ( α , φ ) h m ( α , φ ) d α d φ , j = 1 , 2 , , M ,
or in matrix notation, by entering the B -vector of coefficients b j , and
V j = Ω U ( α , φ ) h j ( α , φ ) d α d φ , G j m = Ω h j ( α , φ ) h m ( α , φ ) d α d φ , j , m = 1 , 2 , , M ,
V = G B .
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 g m ( α , φ ) 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 G —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 g m ( α , φ ) functions in the representation (2) without significant distortion of the real distribution I ( α , φ ) , 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 2 N + 1 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 α 0 is
f ( α 0 α ) = n = N N J n exp i k d n sin α = n = N N exp i k d n ( sin α 0 sin α ) ,
where k = 2 π / λ and J n are the complex amplitudes of the current on the n-th emitter ( n = N , , N ) with accuracy up to a constant equal in the case under consideration to exp i k d n sin α 0 . The angles α 0 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
U ( α ) = n = N N C n exp i k d n sin α ,
with
C n = Ω exp i k d n sin α I ( α ) d α , n = 0 , ± 1 , , ± N .
Here, C n is the signal received by the individual emitter of the DAA.
It is known that, in this case, solution I ( α ) up to the normalizing factor can be represented in the form:
I ( α ) = n = N N exp i k d n sin α C n + ψ ( α ) ,
where ψ ( α ) is an arbitrary function orthogonal to all exponents (9) on the segment [ λ / d ; λ / d ] .
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, C n = C n * , we obtain
D n = R e ( C n ) m = 1 M b m Ω cos i k d n sin α g n ( α ) d α , n = 0 , , N .
Next, by entering a rectangular matrix H with elements
H n m = Ω cos i k d n sin α g m ( α ) d α , n , m = 0 , , N ,
we come to a new SLAE with respect to the vector of unknown expansion coefficients B
D = H B .
With the same dimension, the conditionality matrix numbers G from (7) are significantly larger than those of matrices H 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 g m ( α ) functions in (2) and, hence, in (12)–(14).
At the beginning of the process, at M < N , 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 f ( B ) = H B D , where A is the norm of vector A . 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 M = N , 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 H —the use of the received signal (vector D ) is not required. When solving inverse problems of image reconstruction with the same matrix H , but with a different D , it is enough to calculate the LU decomposition of H . The speed increases even more if there is a database of a specific measuring system in the form of LU decompositions of H for the most commonly used systems of g m ( α , φ ) 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 H 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 M > N , 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 B .
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 B min . This solution is found when minimizing the norm B 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 0.5 , 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 0.9 θ 0.5 . Step functions were chosen to represent solutions as g m ( α ) 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 1 / 6 of the beam width θ 0.5 .
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 ( 1 / 30 ) θ 0.5 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 g m ( α ) functions used.
When the objects under study approach, the quality of the approximate solutions obtained begins to deteriorate. Figure 2 shows the found images I ( α ) at a distance of 0.2 θ 0.5 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 U ( α ) 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 1 / 3 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 0.1 θ 0.5 (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 0.06 θ 0.5 .
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 I ( α ) 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 I ( α ) , then the second allows one to highlight the details of the image, albeit with slight distortions.

6. Conclusions

This article proposes a new method for the digital processing of measurement results, allowing technical systems based on digital antenna arrays (DAA) to increase the achieved level of angular super-resolution. The method is based on the intelligent analysis of signals received by each element of DAA and solving inverse problems described by the Fredholm integral equation of the first kind of convolution type with additional conditions. The results of numerical studies have shown that, compared with other varieties of algebraic methods, the proposed signal processing has greater noise immunity and allows obtaining solutions with high levels of angular super-resolutions, with good image recovery quality of complex objects. The performances of algorithms based on the presented method allow using them in real-time autonomous and remotely controlled robotic systems.

Author Contributions

Conceptualization, B.A.L. and E.Y.R.; methodology, B.A.L. and E.Y.R.; software, B.A.L. and E.Y.R.; validation, B.A.L. and E.Y.R.; formal analysis, B.A.L. and E.Y.R.; investigation, B.A.L. and E.Y.R.; writing—original draft preparation, B.A.L. and E.Y.R.; writing—review and editing, B.A.L. and E.Y.R.; visualization, B.A.L. and E.Y.R.; supervision, B.A.L. and E.Y.R.; project administration, E.Y.R.; funding acquisition, E.Y.R. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the Russian Science Foundation (RSF), grant no. 23-29-00448.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DAAdigital antenna array(s)
RPradiation pattern
SLAEsets of linear algebraic equations
MUSICmultiple signal classification
ESPRITestimation of signal parameters via rotational invariant techniques
MMEmethod of maximum entropy
LUlower–upper

References

  1. Waweru, N.P.; Konditi, D.B.O.; Langat, P.K. Performance Analysis of MUSIC Root-MUSIC and ESPRIT. DOA Estimation Algorithm. Int. J. Electr. Comput. Energ. Electron. Comm. Eng. 2014, 8, 209–216. [Google Scholar]
  2. Rao, B.D.; Hari, K.V.S. Performance Analysis of Root-MUSIC. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 1939–1949. [Google Scholar] [CrossRef]
  3. Kim, K.T.; Seo, D.K.; Kim, H.T. Ecient radar target recognition using the MUSIC algorithm and invariant features. IEEE Trans. Antennas Propag. 2002, 50, 325–337. [Google Scholar]
  4. Lavate, T.B.; Kokate, V.K.; Sapkal, A.M. Performance Analysis of MUSIC and ESPRIT. DOA Estimation Algorithms for Adaptive Array Smart Antenna in Mobile Communication. In Proceedings of the 2010 Second International Conference on Computer and Network Technology, Bangkok, Thailand, 23–25 April 2010; pp. 308–311. [Google Scholar]
  5. Ahmed, A.; Recht, B.; Romberg, J. Blind Deconvolution Using Convex Programming. IEEE Trans. Inf. Theory 2012, 60, 1711–1732. [Google Scholar] [CrossRef] [Green Version]
  6. Almeida, M.S.; Figueiredo, M.A. Deconvolving images with unknown boundaries using the alternating direction method of multipliers. IEEE Trans. Image Process 2013, 22, 3074–3086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Dudik, M.; Phillips, S.J.; Schapire, R.E. Maximum entropy density estimation with generalized regularization and an application to species distribution modeling. J. Mach. Learn. Res. 2007, 8, 1217–1260. [Google Scholar]
  8. Borgiotti, G.V.; Kaplan, L.J. Superresolution of uncorrected interference sources by using adaptive array technique. IEEE Trans. Antennas Propagat. 1979, AP-27, 842–845. [Google Scholar] [CrossRef]
  9. Tan, W.Q.; Hou, Y.G. Estimation of direction of source arrival based upon MUSIC and Capon. J. Nanchang Inst. Technol. 2008, 27, 20–23. [Google Scholar]
  10. Stoica, P.; Sharman, K.C. Maximum likelihood methods for direction-of-arrival estimation. IEEE Trans. Acoust. Speech Signal Process. 1990, 38, 1132–1143. [Google Scholar] [CrossRef]
  11. Cetin, M.; Karl, W.C. Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization. IEEE Trans. Image Process. 2001, 10, 623–631. [Google Scholar] [CrossRef] [Green Version]
  12. Candes, E.J.; Fernandezgranda, C. Super-Resolution from Noisy Data. J. Fourier Anal. Appl. 2013, 19, 1229–1254. [Google Scholar] [CrossRef] [Green Version]
  13. Morse, P.; Feshbach, H. Methods of Theoretical Physics; McGraw-Hill Book Company: New York, NY, USA; Toronto, ON, Canada; London, UK, 1953; 1978p. [Google Scholar]
  14. Lagovsky, B.A.; Rubinovich, E.Y. Algorithms for digital processing of measurement data providing angular super-resolution. Mechatron. Autom. Control 2021, 22, 349–356. [Google Scholar]
  15. Lagovsky, B.A. Angular super-resolution in two-dimensional radar problems. Radio Eng. Electron. 2021, 66, 853–858. [Google Scholar]
  16. Lagovsky, B.; Rubinovich, E. Achieving Angular Superresolution of Control and Measurement Systems in Signal Processing. Adv. Syst. Sci. Appl. 2021, 21, 104–116. [Google Scholar]
  17. Lagovsky, B.; Samokhin, A.; Shestopalov, Y. Angular Superresolution Based on a Priori Information. Radio Sci. 2021, 5, 1–11. [Google Scholar] [CrossRef]
  18. Lagovsky, B.A.; Samokhin, A.B.; Shestopalov, Y.V. Increasing effective angular resolution measuring systems based on antenna arrays. In Proceedings of the 2016 URSI International Symposium on Electromagnetic Theory (EMTS), Espoo, Finland, 14–18 August 2016; pp. 432–434. [Google Scholar]
  19. Lagovsky, B.A.; Chikina, A.G. Superresolution in signal processing using a priori information. In Proceedings of the 2017 International Conference on Electromagnetics in Advanced Applications (ICEAA), Verona, Italy, 11–15 September 2017; pp. 944–947. [Google Scholar]
  20. Shchukin, A.A.; Pavlov, A.E. Parameterization of user functions in digital signal processing for obtaining angular superresolution. Russ. Technol. J. 2022, 10, 38–43. [Google Scholar] [CrossRef]
  21. Lagovsky, B.A.; Samokhin, A.B.; Shestopalov, Y.V. Creating Two-Dimensional Images of Objects with High Angular Resolution. In Proceedings of the 2014 24th International Crimean Conference Microwave & Telecommunication Technology, Sevastopol, Ukraine, 7–13 September 2014; pp. 114–115. [Google Scholar]
  22. Nickel, U. Angular superresolution with phased array radar: A review of algorithms and operational constraints. IEEE Proc. Commun. Radar Signal Process. 1989, 134, 53–59. [Google Scholar] [CrossRef]
  23. Charbonnier, L.; Blanc-F’eraud, G.; Barlaud, M. Deterministic edge-preserving regularization in computed imaging. IEEE Trans. Image Process. 1997, 6, 298–310. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Zhang, Y.; Huang, Y. Angular super-resolution for scanning radar with improved regularized iterative adaptive approach. IEEE Geosci. Remote Sens. Lett. 2016, 13, 1–5. [Google Scholar] [CrossRef]
  25. Yang, J.; Kang, Y.; Zhang, Y. A Bayesian angular super-resolution method with lognormal constraint for sea-surface target. IEEE Access 2020, 8, 13419–13428. [Google Scholar] [CrossRef]
  26. Ramani, S.; Liu, Z.; Rosen, J.; Nielsen, J.; Fessler, J.A. Regularization parameter for nonlinear iterative image restoration and MRI selection reconstruction using GCV and SURE- based methods. IEEE Trans. Image Process. 2012, 21, 3659–3672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Evdokimov, N.A.; Lukyanenko, D.V.; Yagola, A.G. Application of multiprocessor systems for the solution of Fredholm integral equations. Comput. Methods Program. 2009, 10, 263–267. [Google Scholar]
  28. Schwarzenberg-Czerny, A. On matrix factorization and efficient least squares solution. Astron. Astrophys. Suppl. Ser. 1995, 110, 405–410. [Google Scholar]
Figure 1. Image reconstruction of two point sources: 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method.
Figure 1. Image reconstruction of two point sources: 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method.
Mathematics 11 01056 g001
Figure 2. Image reconstruction of point sources at a distance of 0.2 θ 0.5 : 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method; 4-change in the amplitude of a received signal U ( α ) when scanning the sector.
Figure 2. Image reconstruction of point sources at a distance of 0.2 θ 0.5 : 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method; 4-change in the amplitude of a received signal U ( α ) when scanning the sector.
Mathematics 11 01056 g002
Figure 3. Image reconstruction of point sources at a distance of 0.1 θ 0.5 : 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method.
Figure 3. Image reconstruction of point sources at a distance of 0.1 θ 0.5 : 1-signal sources; 2-solution found by the usual algebraic method; 3-solution obtained on the basis of a modification of the algebraic method.
Mathematics 11 01056 g003
Figure 4. The original and reconstructed image of a complex object: 1-true distribution I ( α ) ; 2-solution found by the algebraic method; 3-solution for the DAA.
Figure 4. The original and reconstructed image of a complex object: 1-true distribution I ( α ) ; 2-solution found by the algebraic method; 3-solution for the DAA.
Mathematics 11 01056 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lagovsky, B.A.; Rubinovich, E.Y. Algebraic Methods for Achieving Super-Resolution by Digital Antenna Arrays. Mathematics 2023, 11, 1056. https://doi.org/10.3390/math11041056

AMA Style

Lagovsky BA, Rubinovich EY. Algebraic Methods for Achieving Super-Resolution by Digital Antenna Arrays. Mathematics. 2023; 11(4):1056. https://doi.org/10.3390/math11041056

Chicago/Turabian Style

Lagovsky, Boris A., and Evgeny Ya. Rubinovich. 2023. "Algebraic Methods for Achieving Super-Resolution by Digital Antenna Arrays" Mathematics 11, no. 4: 1056. https://doi.org/10.3390/math11041056

APA Style

Lagovsky, B. A., & Rubinovich, E. Y. (2023). Algebraic Methods for Achieving Super-Resolution by Digital Antenna Arrays. Mathematics, 11(4), 1056. https://doi.org/10.3390/math11041056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop