Efficient direction-of-arrival estimation using low-rank approximations

Abstract

A novel and useful Direction of Arrival (DOA) estimation system and method is provided that allows for DOA-dependent calibration with a complexity of order O(NlogN). In one embodiment, the architecture includes several Fast Fourier Transform (FFT) machines operating in parallel with coefficient multiplication before and after FFT operations. These pre- and post-FFT coefficients are computed using an optimal low-rank approximation of the distortion matrix using singular value decomposition. The values ​​of the pre- and post-FFT calibration coefficients before and after each rank are calculated from the singular value decomposition of the distortion matrix C=B / F, where B is the digital beamforming (DBF) matrix and F is the ideal FFT matrix. A method for obtaining the beamforming matrix B is also disclosed. This architecture utilizes the relatively low rank required for full matrix multiplication, N. 2 This achieves K×N×log2N operations, where K is the rank of the approximation and N is the length of the FFT, which is significantly smaller than the approximation of K×N×log2N operations. We disclose a circuit and method for computing the calibration coefficients, as well as a simple proof of approximating this to a general beamforming matrix.

Description

[Technical Field]

[0001] The subject matter disclosed herein relates to the field of imaging radar, sonar, ultrasound, and other sensors that perform distance measurements with FMCW signals and / or angle measurements with digital beamforming and array processing, and more particularly to an efficient mechanism for estimating DOA (Direction of Arrival) using LRA (Low Rank Approximation) techniques. [Background technology]

[0002] In recent years, many industries, such as the automotive and delivery industries, have moved towards autonomous solutions. These autonomous platforms need to operate in their environments, interacting with both stationary and moving objects. To this end, these systems require a sensor suite that allows them to reliably and efficiently sense their surroundings. For example, in order for an autonomous car to plan a route on a road with other vehicles, the trajectory planner must have a 3D map of the environment showing moving objects.

[0003] Visual sensors are not degraded by bad weather or poor visibility due to fog, smoke, sand, storms, etc. They are also limited in estimating radial velocity. LIDAR (Light Detection and Ranging) devices are used to measure the distance to a target by shining a laser beam at the target. However, these devices are expensive, have moving parts, and have very limited range. Radar is an augmentation technology, not a replacement.

[0004] Due to the natural limitations of visual sensors in distance accuracy and reliability issues with optical (e.g. laser) technologies, the best solution to generate this 3D map is through radar technology, which imposes a new set of requirements that modern radars do not comply with.

[0005] Generally, the larger the aperture of a receiving antenna, the more radiation it receives, resulting in a higher sensitivity, or equivalently, a narrower main lobe, and therefore the ability of the receiving antenna to receive weak signals and provide a relatively accurate indication as to their direction.

[0006] On the other hand, vehicle radars, including imaging radars in automobiles, typically have smaller apertures. Automotive radars transmit relatively weak signals due to power consumption and regulatory constraints, but link budgets may require lower sensitivity due to their shorter range and stronger reflected signals from targets. However, while vehicle radars do not need to detect point targets, such as aircraft detected by missiles, they do require high accuracy to provide an image of the environment that serves as input to one or more tracking and post-processing algorithms and / or SLAM (Simultaneous Localization and Mapping) algorithms, which detect the location of nearby obstacles, such as other vehicles or pedestrians, and generate an object list from raw radar detections. A narrower lobe with higher accuracy can provide sharper contours in the target image. Lobe width is determined solely by the equivalent aperture normalized to the wavelength of the transmitted radar signal, rather than by the number of receiving antenna elements within the aperture, which does not adversely affect sensitivity—that is, the ability to detect weak reflected signals—and the ambiguity resolution and sidelobe level.

[0007] Another important performance parameter of imaging radar is the sidelobe level of the antenna array. Sidelobes reflected from strong targets can mask weaker targets or cause false detections. For example, a large object such as a wall located in the direction of the sidelobe will cause reflections from the wall to appear in the mainlobe. This can mask reflections from obstacles such as pedestrians or create phantom obstacles that vehicles can stop at.

[0008] Therefore, it is important to reduce side lobes as much as possible in automotive imaging radar. Furthermore, there is a need for a compact radar switched array antenna that has high azimuth and elevation resolution and accuracy that increases the effective aperture while using a small number of transmit (TX) and receive (RX) elements that meet the requirements of cost, space, power, and reliability.

[0009] Recently, the application of radar in the automotive industry has begun to emerge. High-end cars are already equipped with radar to provide parking assistance and lane departure warning to drivers. Currently, interest in self-driving cars is increasing, and it is expected that self-driving cars will become a major driving force in the automotive industry in the coming years.

[0010] Self-driving cars offer new perspectives for the application of radar technology in automobiles. Automotive radars can not only assist the driver but also play an active role in controlling the vehicle. Therefore, they could become key sensors in the vehicle's autonomous control system.

[0011] Radar is preferred over other alternatives, such as sonar or LIDAR, because it is less affected by weather conditions and the deployed sensors can have a very small negative impact on the vehicle's aerodynamics and appearance. FMCW (Frequency Modulated Continuous Wave) radar is one type of radar that offers several advantages over other radars. For example, it ensures that the distance and speed information of surrounding objects can be detected simultaneously. This information is important for the control system of autonomous vehicles to provide safe and collision-free operation.

[0012] FMCW radar is commonly used for short-range detection, such as automotive radar. Some advantages of FMCW radar in automotive applications include the following: (1) FMCW modulation is relatively easy to generate, offers wide bandwidth, high average power, high accuracy, and low cost due to low-bandwidth processing, enabling very good range resolution and allowing Doppler shift to be used for velocity determination; (2) FMCW radar can operate at short range with good performance; (3) FMCW sensors can be fabricated compactly with a single RF transmit source with an oscillator also used to downconvert the received signal; and (4) because transmission is continuous, modest output power of solid-state components is sufficient.

[0013] A radar system installed in an automobile should be able to provide the information required by the control system in real time. A baseband processing system that can provide enough computing power to meet the real-time system requirements is necessary. The processing system performs digital signal processing on the received signal to extract useful information such as the range and velocity of surrounding objects.

[0014] Vehicles, especially cars, are now increasingly equipped with technology designed to assist drivers in critical situations. In addition to cameras and ultrasonic sensors, automakers are turning to radar as the costs of the related technologies fall. Radar's appeal lies in its ability to quickly and unambiguously measure the speed and distance of multiple objects under all weather conditions. The associated radar signals are frequency modulated and can be analyzed with a spectrum analyzer. In this way, developers of radar components can automatically detect, measure, and display signals in the time and frequency domains, up to frequencies of 500 GHz.

[0015] There is also currently great interest in using radar in the field of autonomous vehicles, which are expected to become even more common in the future. Automotive millimeter-wave radar is suitable for use in collision avoidance and autonomous driving. Millimeter-wave frequencies between 77 and 81 GHz are less susceptible to interference from rain, fog, snow, and other weather factors, dust, and noise than ultrasonic and laser radars. These automotive radar systems typically include a high-frequency radar transmitter that transmits a radar signal in a known direction. The transmitter may transmit the radar signal in a continuous or pulsed mode. These systems further include a receiver connected to an appropriate antenna system that receives echoes or reflections from the transmitted radar signal. Each such reflection or echo represents an object illuminated by the transmitted radar signal.

[0016] ADAS (Advanced Driver Assistance Systems) are systems developed to automate, adapt, and enhance vehicle systems for safety and better driving. Safety features are designed to avoid collisions and accidents by providing technology that alerts the driver to potential problems, or to implement safeguards and take over control of the vehicle to avoid a collision. Adaptive features may automate lighting, provide adaptive cruise control, automate braking, incorporate GPS / traffic alerts, connect to a smartphone, warn the driver of other vehicles or hazards, keep the driver in the correct lane, or show what is in the blind spot.

[0017] There are many forms of ADAS available, with some features built into the vehicle or available as add-on packages. There are also aftermarket solutions available. ADAS rely on input from multiple data sources, including automotive imaging, LIDAR, radar, image processing, computer vision, and in-vehicle networks. Additional input can come from other sources outside the main vehicle platform, such as other vehicles, known as V2V (vehicle-to-vehicle) or vehicle-to-infrastructure systems (e.g., cellular or Wi-Fi data networks).

[0018] Advanced driver assistance systems are currently one of the fastest growing areas in automotive electronics, with steadily increasing adoption of industry-wide quality standards for vehicle safety systems ISO 26262 and the development of technology-specific standards such as IEEE 2020 for image sensor quality and communication protocols such as vehicle information APIs.

[0019] In recent years, many industries, such as the automotive and delivery industries, have moved towards autonomous solutions. These autonomous platforms operate in their environments, interacting with both stationary and moving objects. For this purpose, these systems require a sensor suite that allows them to reliably and efficiently sense their surroundings. For example, in order for an autonomous vehicle to plan a route on a road with other vehicles, the trajectory planner must have a 3D map of the environment showing moving objects.

[0020] Visual sensors are not degraded by bad weather and poor visibility (e.g., fog, smoke, sand, rainstorms, storms, etc.). They are also limited in estimating radial velocity. LIDAR devices (Light Detection and Ranging) are used to measure the distance to a target by shining a laser beam at it. However, they are expensive, have mostly moving parts, and have a very limited range. Therefore, automotive radar is considered an augmentation technology, not a replacement.

[0021] In the automotive field, radar sensors are key components for realizing comfort and safety, such as adaptive cruise control (ACC) or collision mitigation systems (CMS). As the number of automotive radar sensors operating simultaneously and in close proximity to each other increases, radar sensors may receive signals from other radar sensors. Reception of external signals (interference) can cause problems such as ghost targets or a poor signal-to-noise ratio. Figure 1 illustrates such an automotive interference scenario with direct interference from several surrounding vehicles.

[0022] A well-known way to reduce the number of antenna elements in an array is to use a MIMO technique known as a "virtual array," in which separable (e.g., orthogonal) waveforms are transmitted (usually simultaneously) from different antennas and digitally processed to create a larger effective array. The shape of this "virtual array" is a special convolution of the positions of the transmitting and receiving antennas.

[0023] It is also known that bandpass sampling allows the delamped signal to be sampled at a lower A / D frequency while retaining range information for targets at ranges consistent with the designed bandpass filter.

[0024] Achieving high resolution simultaneously in angle, range and Doppler dimensions is a major challenge, as resolution increases linearly with hardware complexity (among other things).

[0025] Additionally, DOA (direction of arrival) estimation is also a critical component in any radar system. For imaging radar, this is typically performed digitally and is commonly referred to as digital beamforming (DBF). Prior art methods that fall into this category involve linear operations that can be implemented as matrix multiplications. Other nonlinear approaches are typically referred to as super-resolution techniques. However, these techniques require very high computational power and are difficult to implement in low-cost consumer sensors.

[0026] However, a problem with prior art DOA estimation techniques is typically the number of calculations required. Solutions to this problem typically involve either (1) using additional computing power, (2) lowering the frame rate, or (3) reducing the number of range-Doppler bins, or any combination of the above. Therefore, it would be desirable to have a radar system that performs DOA estimation without compromising the above parameters and that presents a relatively low computational load. Furthermore, the radar should achieve comparable SLL (side lobe level) compared to prior art full DBF calculations. Summary of the Invention

[0027] The present invention provides a system and method for direction-of-arrival estimation with order O(NlogN) complexity that enables DOA-dependent calibration. In one embodiment, the architecture includes several FFT (Fast Fourier Transform) machines operating in parallel with coefficient multiplication before and after the FFT operation. These pre- and post-FFT coefficients are computed using an optimal low-rank approximation of the distortion matrix using singular value decomposition. The values ​​of the pre- and post-FFT calibration coefficients before and after each rank are calculated from the singular value decomposition of the distortion matrix C=B / F, where B is the DBF (Digital Beamforming) matrix and F is the ideal FFT matrix. A method for obtaining the beamforming matrix B is also disclosed.

[0028] This architecture has a relatively low rank compared to the N required for full matrix multiplication. 2 This realizes K×N×log2N operations, where K is the rank of the approximation and N is the length of the fast Fourier transform, which is significantly smaller than the approximation of K×N×log2N operations. Circuits and methods for computing the calibration coefficients, as well as a simple proof of approximating this to a general beamforming matrix, are further disclosed.

[0029] Thus, in accordance with the present invention, there is provided a method of estimating DOA (Direction of Arrival) of a signal for use in a radar system, the method comprising: receiving input data; multiplying each element of the input data by a plurality of sets of a priori coefficients to obtain a first plurality of results; performing a plurality of Fast Fourier Transform operations on the first plurality of results to generate a second plurality of results; multiplying each element of the second plurality of results by the plurality of sets of a posteriori coefficients to obtain a third plurality of results; and summing the third plurality of results to obtain an approximate DOA estimate.

[0030] Also according to the present invention, there is provided a method for estimating a direction of arrival (DOA) of a signal for use in a radar system, comprising: receiving input data x; multiplying each element of the input data x by k sets of a priori coefficients Vk to generate k first results diag(Vk)·x; performing k fast Fourier transform operations F on the k first results to generate k second results Fdiag(Vk)·x; and multiplying each element of the k second results by k sets of a posteriori coefficients Uk to generate k third results.

number

[0031] According to the present invention, there is provided an apparatus for estimating the direction of arrival (DOA) of a signal for use in a radar system, operative to receive a receive antenna array response x, multiplying each element of the antenna array response by k sets of a priori coefficients Vk to generate k first results diag(Vk)·x, performing k Fast Fourier Transform operations F on the k first results to generate k second results F·diag(Vk)·x, and multiplying each element of the k second results by k sets of a posteriori coefficients Uk to generate k third results.

number

[0032] According to the present invention, there is provided a PCB (Printed Circuit Board) assembly including a plurality of transmitting antennas fabricated on one side of the PCB assembly, a plurality of receiving antennas fabricated on an opposite side of the PCB assembly, and a transceiver coupled to the plurality of transmitting antennas and the plurality of receiving antennas, the transceiver operative to generate and transmit transmission signals to the one or more transmitting antennas and receive wave signals reflected by the one or more receiving antennas; and a transceiver coupled to the transceiver, the transceiver receiving input data x, multiplying each element of the input data x by a set of k pre-coefficients Vk to generate k first results diag(Vk)x, performing fast Fourier transform operations F on the k first results to generate k second results Fdiag(Vk)x, and multiplying each element of the k second results by a set of k post-coefficients Uk to generate k third results

number

[0033] The present invention will now be described in more detail in the following exemplary embodiments and with reference to the drawings, in which identical or similar elements are in part designated by identical or similar reference numerals and in which features of the various exemplary embodiments are combinable. The present invention is herein described, by way of example only, with reference to the accompanying drawings, in which:

[0034] [Figure 1] FIG. 1 illustrates an example of a street scene incorporating several vehicles equipped with automotive radar sensor units. [Figure 2] FIG. 1 illustrates an exemplary radar system incorporating multiple receivers and transmitters. [Figure 3] FIG. 1 illustrates an exemplary radar transceiver constructed in accordance with the present invention. [Figure 4]FIG. 1 is a high-level block diagram illustrating an exemplary MIMO FMCW radar of the present invention. [Figure 5] FIG. 1 is a block diagram illustrating an exemplary DRP (Digital Radar Processor) IC configured in accordance with the present invention. [Figure 6] FIG. 1 is a high-level block diagram illustrating an exemplary DOA (Direction of Arrival) estimation using an LRA (Low Rank Approximation) technique. [Figure 7] FIG. 1 illustrates an exemplary method for low-rank approximation-based DOA estimation. [Figure 8] FIG. 1 illustrates an exemplary method for calculating a priori and a posteriori coefficients. DETAILED DESCRIPTION OF THE INVENTION

[0035] In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be apparent to those skilled in the art that the present invention may be practiced without these specific details. In other instances, well-known methods, procedures, and components have not been described in detail so as not to obscure the present invention.

[0036] Among these disclosed benefits and improvements, other objects and advantages of the present invention will become apparent from the following description taken in conjunction with the accompanying drawings. Detailed embodiments of the present invention are disclosed herein. However, it will be understood that the disclosed embodiments are merely exemplary of the disclosure, which may be embodied in various forms. Furthermore, each of the examples given in connection with various embodiments of the present invention are illustrative only and are not intended to be limiting.

[0037] The subject matter which is regarded as the invention is particularly pointed out and distinctly claimed in the claims at the end of this specification, but the invention, both as to organization and method of operation, together with objects, features, and advantages thereof, may best be understood by reference to the following detailed description when read in conjunction with the accompanying drawings.

[0038] The drawings constitute a part of this specification and include exemplary embodiments of the present invention, illustrating various objects and features thereof. The drawings are not necessarily to scale, and some features may be enlarged or minimized to show details of particular components. Furthermore, measurements, specifications, and the like shown in the drawings are for illustrative purposes and not limiting. Accordingly, specific structural and functional details disclosed herein are not intended to be limiting, but should merely be interpreted as a representative basis for teaching those skilled in the art to variously employ the present invention. Furthermore, where considered appropriate, reference numerals may be repeated among the drawings to indicate corresponding or similar elements.

[0039] The illustrated embodiments of the present invention may be implemented using electronic components and circuits that are for the most part known to those skilled in the art, and therefore will not be described in more detail than is necessary for the understanding and appreciation of the concepts underlying the present invention and in order not to be obvious from the knowledge of the present invention.

[0040] Any reference herein to a method should be applied mutatis mutandis to a system capable of carrying out the method. Any reference herein to a system should be applied mutatis mutandis to a method that may be carried out by the system.

[0041] Throughout this specification and claims, the following terms take the meanings expressly associated therewith, unless the context clearly dictates otherwise. As used herein, the phrases "in one embodiment," "in an exemplary embodiment," and "in some embodiments" do not necessarily refer to the same embodiment, although they may. Furthermore, as used herein, the phrases "in another embodiment," "in an alternative embodiment," and "in some other embodiments" do not necessarily refer to different embodiments, although they may. Thus, as described below, various embodiments of the invention can be readily combined without departing from the scope or spirit of the invention.

[0042] Also, as used herein, unless the context clearly dictates otherwise, the term "or" is an inclusive "or" operator and is equivalent to the term "and / or." The term "based on" is non-exclusive and allows for based on additional unlisted factors, unless the context clearly dictates otherwise. Also, throughout this specification, the meanings of "a," "an," and "the" include plural referents. The meaning of "in" includes "in" and "on."

[0043] Frequency Modulated Continuous Wave (FMCW) radar is a radar that uses frequency modulation. The operating principle of FMCW radar is to transmit a continuous wave with increasing (or decreasing) frequency. Such a wave is called a chirp. Figure 4 shows an example of a chirp waveform 10. After being reflected by an object, the transmitted wave is received by a receiver. An example of a transmitted chirp waveform 12 and a received (i.e., reflected) chirp waveform 14 at the receiver is shown in Figure 5.

[0044] Considering the use of radar in automotive applications, four frequency bands of varying bandwidths are currently available to automakers between 24 GHz and 77 GHz. The 24 GHz ISM band has a maximum bandwidth of 250 MHz, while the ultra-wideband (UWB) band between 76 and 81 GHz offers up to 5 GHz. Bands with bandwidths up to 4 GHz are located between the 77 and 81 GHz frequencies. Currently, many applications use this band. Note that other frequencies allocated for this use include 122 GHz and 244 GHz, which are only 1 GHz wide. Having sufficient bandwidth is important for radar applications because the signal bandwidth determines the range resolution.

[0045] Conventional digital beamforming FMCW radars are characterized by very high resolution across the radial, angular, and Doppler dimensions. Imaging radars are based on the well-known technique of phased arrays using ULAs (Uniformly Linearly Distributed Arrays). It is well known that the far-field beam pattern of a linear array architecture is obtained using a Fourier transform. Range measurements are obtained by performing a Fourier transform on the delamp signal generated by multiplying the conjugate of the transmitted and received signals. Radar range resolution is determined by the radar's RF bandwidth and is equal to the speed of light, c, divided by twice the RF bandwidth. Doppler processing is performed by performing a Fourier transform across the slow time dimension, and its resolution is limited by the coherent processing interval (CPI), i.e., the total transmission time used for Doppler processing.

[0046] When using radar signals in automotive applications, it is desirable to simultaneously determine the velocity and distance of multiple objects in one measurement cycle. Conventional pulse radars cannot easily handle such a task, as they can only determine distance based on the timing offset between the transmitted and received signals within one cycle. When velocity is also to be determined, frequency modulated signals, such as linear FMCW (frequency modulated continuous wave) signals, are used. Pulse Doppler radars can also measure the Doppler offset directly. The frequency offset between the transmitted and received signals is also known as the beat frequency. The beat frequency is determined by the Doppler frequency component f D and the delay component f T The Doppler component contains information about velocity, and the delay component contains information about range. Because there are two unknowns, range and velocity, determining the desired parameter requires two beat frequency measurements. A second signal, whose frequency is linearly modified immediately after the first, is incorporated into the measurement.

[0047] Using an FM chirp sequence, both parameters can be determined in one measurement cycle. Because a single chirp is very short compared to the full measurement cycle, each beat frequency is dominated by a delay component f T In this way, the range can be ascertained immediately after each chirp. Determining the phase shift between several consecutive chirps in the sequence allows the Doppler frequency to be determined using a Fourier transform, allowing the vehicle's speed to be calculated. Note that the velocity resolution improves as the length of the measurement cycle increases.

[0048] MIMO (multiple-input, multiple-output) radar is a type of radar that transmits and receives signals using multiple TX and RX antennas. Each transmit antenna in the array independently radiates a waveform signal that is different from the signals radiated from the other antennas. Alternatively, the signals may be identical but transmitted at non-overlapping times. The reflected signals belonging to each transmitter antenna can be easily separated at the receiver antenna because (1) orthogonal waveforms are used for transmission or (2) they are received at non-overlapping times. A virtual array containing the information from each transmit antenna to each receive antenna is fabricated. Thus, if there are M transmit antennas and N receive antennas, the virtual array has M·N independent transmit and receive antenna pairs using only M+N physical antennas. This characteristic of MIMO radar systems provides several advantages, including improved spatial resolution, a larger antenna aperture, and improved sensitivity for detecting slow-moving objects.

[0049] As mentioned above, the signals transmitted from different TX antennas are orthogonal. The orthogonality of the transmitted waveforms can be obtained using TDM (Time Division Multiplexing), Frequency Division Multiplexing or Spatial Coding. In the examples and explanations presented herein, TDM is used, which allows only a single transmitter to transmit at a time.

[0050] The radar of the present invention operates to reduce complexity, cost, and power consumption by implementing a time-multiplexed MIMO FMCW radar, as opposed to full MIMO FMCW. The time-multiplexed approach to automotive MIMO imaging radar offers significant cost and power advantages compared to full MIMO radar. Full MIMO radar simultaneously transmits multiple separable signals from multiple transmit array elements. These signals must typically be separated in each receive channel using a bank of matched filters. In this case, the entire virtual array is loaded at once.

[0051] In time-multiplexed MIMO, only one TX (transmit) array element transmits at a time. The transmit side is greatly simplified and a bank of matched filters per RX (receive) channel is not required. The virtual array is gradually loaded over time, taking the time it takes to transmit from all TX elements in the array.

[0052] 2 shows a high-level block diagram illustrating an exemplary radar system incorporating multiple receivers and transmitters. The radar system, generally referenced 280, includes a DRP (Digital Radar Processor) / signal processor 282 that performs signal processing functions, including, among other things, DOA (Direction of Arrival) estimation utilizing the inventive LRA (Low Rank Approximation) mechanism; a plurality of N transmitters TX1 through TXN 284, each coupled to a transmit antenna 288; and a plurality of M receivers RX1 through RXM 286, each coupled to a receive antenna 290. TX data lines 292 connect the DRP to the transmitters, RX lines 294 connect the receivers to the DRP, and control signals 296 are provided by the DRP to each of the transmitters 284, 286. Note that N and M can be any positive integers greater than 1.

[0053] 3 illustrates an exemplary radar transceiver constructed in accordance with the present invention. The radar transceiver, generally referenced 80, includes a transmitter 82, a receiver 84, and a controller 83. The transmitter 82 includes a nonlinear frequency hopping sequencer 88, an FMCW chirp generator 90, an LO (local oscillator) 94, a mixer 92, a power amplifier (PA) 96, and an antenna 98.

[0054] The receiver 84 includes an antenna 100, an RF front end 101, a mixer 102, an IF block 103, an ADC 104, fast time range processing 106, slow time processing (Doppler and precision range) 108, and azimuth and elevation processing 110.

[0055] In operation, the nonlinear frequency hopping sequencer 88 generates a nonlinear starting frequency hopping sequence. The starting frequency of each chirp is input to the FMCW chirp generator 90, which functions to generate a chirp waveform at the specific starting frequency. The chirps are upconverted via mixer 92 to the appropriate band (e.g., the 80 GHz band) according to an LO 94. The upconverted RF signal is amplified via a PA 96 and output to an antenna 98, which may include an antenna array in the case of a MIMO radar.

[0056] On the receive side, echo signals received by the antenna 100 are input to the RF front-end block 101. In a MIMO radar, the receive antenna 100 includes an antenna array. The signal from the RF front-end circuit is mixed with the transmit signal via a mixer 102 to generate a beat frequency that is input to an IF filter block 103. The output of the IF block is converted to digital form via an ADC 104 and input to a fast time processing block 106 to generate coarse range data. A slow processing block 108 functions to generate both fine range data and Doppler velocity data. Then, azimuth and elevation data are calculated via an azimuth / elevation processing block 110. 4D image data 112 is input to downstream image processing and detection. Note that in one embodiment, the azimuth / elevation processing block 110 performs DOA (direction of arrival) estimation using the inventive LRA (low rank approximation) mechanism.

[0057] FIG. 4 shows a high-level block diagram illustrating an exemplary MIMO FMCW radar of the present invention. The radar transceiver sensor, generally referenced 40, includes multiple transmit circuits 66, multiple receive circuits 58, a ramp or chirp generator 60 including an LO (local oscillator) 61, a nonlinear frequency hopping sequencer 62, an optional TX element sequencer 75 (dashed line), and a DRP (digital radar processor) / signal processing block 44 including block 45 that, in one embodiment, provides direction of arrival (DOA) estimation utilizing the LRA (low rank approximation) mechanism of the present invention. In operation, the radar transceiver sensor typically communicates with, and may be controlled by, a host 42. Each transmit block includes a power amplifier 70 and an antenna 72. The transmitter receives the transmit signal output of the chirp generator 60, which is fed to a PA within each transmit block. The optional TX element sequencer (dashed line) generates multiple enable signals 64 that control the transmit element sequence. It will be appreciated that DOA estimation may be performed in radar systems with or without TX element sequencing, and with or without MIMO operation. Furthermore, DOA estimation is not limited to implementation in MIMO FMCW radars, but may also be performed using other types of radar systems.

[0058] Each receive block includes an antenna 58, an LNA (Low Noise Amplifier) ​​50, a mixer 52, an IF (Intermediate Frequency) block 54, and an ADC (Analog-to-Digital Converter) 56. The signal processing block 44 may include any suitable electronic device capable of processing, receiving, or transmitting data or instructions. For example, the processing unit may include one or more of a microprocessor, a CPU (Central Processing Unit), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), a DSP (Digital Signal Processor), a GPU (Graphical Processing Unit), or a combination of these devices. As described herein, the term "processor" refers to a single processor. of Processing Unit, Multiple Processors , complex It is intended to encompass a number of processing units, or other suitably configured computing elements.

[0059] For example, a processor may include one or more general-purpose CPU cores and, optionally, one or more special-purpose cores (e.g., DSP cores, floating point, gate arrays, etc.) The one or more general-purpose cores execute general-purpose opcodes, while the special-purpose cores perform functions specific to their purpose.

[0060] The connected or embedded memory may include DRAM (Dynamic Random Access Memory) or EDO (Extended Data Output) memory, or other types of memory such as ROM, static RAM, flash, NVSRAM (Non-Volatile Static Random Access Memory), removable memory, bubble memory, or any combination of the above. The memory stores electronic data that can be used by the device. For example, the memory may store electronic data or content such as radar-related data, audio and video files, documents and applications, device settings and user settings, timing and control signals, or data, data structures, or databases for various modules. The memory may be configured as any type of memory.

[0061] The transmitted and received signals are mixed (i.e., multiplied) to produce a signal that is processed by the signal processing unit 44. The multiplication process produces two signals, one with a phase equal to the difference of the multiplied signals and the other with a phase equal to the sum of the phases. The sum signal is filtered and the difference signal is processed by the signal processing unit. The signal processing unit performs all necessary processing of the received digital signals and also controls the transmitted signals. Some functions performed by the signal processing block include determining coarse range, velocity (i.e., Doppler), fine range, elevation and azimuth angles, performing interference detection, mitigation and avoidance, and performing SLAM (Simultaneous Localization and Mapping).

[0062] 5 is a block diagram illustrating an exemplary digital radar processor IC of the present invention. The radar processor IC, generally referenced 390, includes several chip service functions 392, including a temperature sensor circuit 396, a watchdog timer 398, a POR (power-on reset) circuit 400, etc.; a PLL system 394, including a power domain circuit 402; an RPU (radar processing unit) 404, including a parallel FFT engine 406, a data analyzer circuit 408, a DMA (direct memory access) circuit 410, and a DOA estimation / LRA mechanism 411; a CPU block 412, including a TX / RX control block 414, a safety core block 418, an L1 and L2 cache memory circuit 424, and a DOA estimation / LRA mechanism 425; a memory system 426; and an I / F (interface) circuit 428. In one embodiment, the RPU is configured to perform DOA (direction of arrival) estimation utilizing the LRA (low-rank approximation) mechanism of the present invention in either the RPU or the CPU, or partially in both.

[0063] The TX / RX control circuit 414 may incorporate settling time control to eliminate frequency source settling time, mutual interference detection, mitigation and avoidance block 416. The safety core block 418 includes a system watchdog timer circuit 420 and an RFBIST circuit adapted to perform continuous testing of RF elements in the radar system. The I / F circuit includes interfaces for radar output data 430, TX control 432, RX control 434, external memory 436, and RF clock 438.

[0064] It is noted that the digital radar processor circuitry 390 can be implemented on monolithic silicon or across several integrated circuits, depending on the particular implementation. Similarly, the transmit and receive circuits can be implemented on a single IC or across multiple ICs, depending on the particular implementation.

[0065] In one embodiment, the DRP390 is used in an automotive radar FMCW MIMO-based system. Such systems require multiple transmit and receive channels to achieve the desired range, azimuth, elevation, and velocity. The more channels, the better the resolution performance. Depending on the implementation, multiple transmit channels may be integrated into a single chip, and multiple receive channels may be integrated into a single chip. The system may include multiple TX and RX chips. Each TX and RX chip can operate as part of a larger system adapted to achieve maximum system performance. In one embodiment, the system also includes at least one control channel. The control channel operates to configure both the TX and RX devices.

[0066] The present invention provides a compact radar switched array antenna with high azimuth and elevation angular resolution and accuracy and increased effective aperture while using a small number of TX and RX elements. The present invention also provides a compact radar antenna array with high azimuth and elevation angular resolution and accuracy and increased effective aperture while reducing undesirable side lobes.

[0067] One embodiment of the present invention relates to a method for increasing the effective aperture of a radar switch / MIMO antenna array using a small number of transmit and receive array elements, where the array of physical radar receive / transmit elements is arranged into at least two opposing RX rows and at least two opposing TX columns, whereby each row contains a plurality of receive elements equally spaced from one another and each column contains a plurality of transmit elements equally spaced from one another, and the array forms a rectangular physical aperture.

[0068] Used as a switching array, the first TX element from one column is activated to transmit a radar pulse during a given time slot, and reflections of the first transmission are received by all RX elements, creating virtual duplicates of two opposing RX rows centered on an origin determined by the position of the first TX element within the rectangular physical aperture.

[0069] This process is repeated for all remaining TX elements during different time slots, creating two virtual, opposing RX rows centered on the origin determined by the position of each activated TX element within the rectangular physical aperture. During each period, transmission reflections from each TX element are received by all RX elements. In this way, a rectangular virtual aperture with dimensions twice the dimensions of the rectangular physical aperture is realized by the two opposing RX rows. This virtual aperture determines the radar beamwidth and sidelobes.

[0070] Note that the above replication method works equally well in MIMO or hybrid switch / MIMO designs where some signals are transmitted simultaneously by multiple TX array elements using orthogonal waveforms that are subsequently separated at the receiver.

[0071] Direction of arrival (DOA) estimation It should be noted that the DOA estimation / LRA mechanism of the present invention is applicable to many types of radars and is not intended to be limited to the exemplary radar system disclosed herein. For example, LRA beamforming can be applied to radars incorporating a uniform linear array (ULA), in which all antenna sensors are arranged in a line and the distance between adjacent sensors is the same as for any two adjacent sensors. MIMO FMCW radar is presented herein for illustrative purposes only to help explain the principles of the DOA estimation mechanism of the present invention.

[0072] Digital Beamforming (DBF) is a well-known technique for determining the direction in which a target is located, also known as the Direction of Arrival (DOA). An array antenna with multiple antenna elements is used to receive reflected waves from the target. The direction of the target is determined by applying a DOA estimation method, such as the well-known beamforming method.

[0073] In a direction-of-arrival (DOA) estimation method using an array antenna, the beamforming method scans the array antenna's main lobe in multiple directions and determines the direction where the output power is maximized as the DOA. Note that the angular resolution is determined by the width of the main lobe. Therefore, if it is desired to increase the resolution so that the directions of multiple targets can be determined, it is desirable to increase the number of antenna elements and thereby the aperture length of the array. Similar methods are minimum-norm methods that determine DOA based on extended algorithms such as eigenvalues ​​and eigenvectors of the correlation matrix of the array's received signals via signal parameter estimation using rotation-invariant techniques (ESPRIT), MUSIC (Multiple Signal Classification), and estimation of signal parameters. Considering these techniques, the order of the correlation matrix, i.e., the number of antenna elements, determines the number of targets that can be detected. Therefore, to be able to determine the directions of many targets, it is desirable to increase the number of antenna elements.

[0074] Typically, during the radar calibration process, a beamforming matrix is ​​computed, and during DOA estimation, this matrix is ​​multiplied by the array response vector. The DBF matrix is ​​then calculated using the DOA-dependent calibration vectors {b1,...,b N}, where N is the number of DOA angles to be estimated (i.e., the number of angles to scan in azimuth and / or elevation), b i ∈C P , and P is the number of antenna elements. Therefore, the DBF matrix B can be written as:

number

number

[0075] In the particular case of an ideal ULA (Uniform Linear Array), the calibration vector is a matched filter of the corresponding steering vector a(θ). Thus,

number

number

number

[0076] One advantage of using FFT operations over multiplication using a DBF matrix B is that it reduces the computational complexity. A typical matrix multiplication is 2While the FFT requires NlogN operations, the FFT requires only NlogN operations. For high-resolution radars with a relatively large number of antenna elements (e.g., P ≈ 100) and a large number of range-Doppler bins over which DOA estimation must be performed, this difference in computational complexity is very significant and poses a significant challenge for real-time implementation. However, a problem arises in that a low-complexity FFT is only effective for an ideal uniform array without any impairments.

[0077] Therefore, for an ideal uniform linear antenna array, a relatively simple FFT operation can be used with NlogN operations for DOA estimation. However, in the real world, when the antenna patterns of different elements are not identical due to manufacturing tolerances and other influences, different complex correction vectors are used for each DOA, and the different complex correction vectors function to correct the antenna impairments as much as possible using, for example, DBF (digital beamforming) matrices. This operation is N 2 This requires operations (i.e., matrix multiplication) that, for high-resolution radars that can exceed N100 elements, must be performed for each range-Doppler (twice for azimuth and elevation), making standard calculations infeasible for such radar systems.

[0078] LRA (Low Rank Approximation) for DBF DOA Estimation Solutions to the above problems include applying additional computational power, lowering the frame rate, and / or reducing the number of range-Doppler bins. Alternatively, a low-rank approximation-based DOA mechanism can be used, which attempts to solve the above problems while requiring a relatively low computational load and does not compromise other parameters. Furthermore, the LRA-based DOA mechanism achieves SLL comparable to full DBF computation.

[0079] In practice, when ULAs are used (i.e., not necessarily dense), and even when considering impaired arrays, the DBF matrix resembles an FFT matrix with an additional constant phase-gain calibration. Thus, it can be expressed as follows:

number

number

[0080] Figure 6 shows a high-level block diagram illustrating an exemplary low-rank approximation for DBF DOA estimation. The LRA circuit, generally referenced 120, includes a number of pre-multipliers 122, an FFT computation block 124, a post-multiplier 126, and an adder 128. Figure 7 shows an exemplary method for low-rank approximation-based DOA estimation.

[0081] The LRA approximation method of the present invention takes advantage of the fact that B is similar to the FFT matrix F. In one embodiment, K FFT machines 124 operating in parallel are used, where K is the rank of the approximation. Before each FFT operation, the antenna array response x 121 is multiplied by a set of coefficients called pre-coefficients Vk via multiplier 122 (S150). An FFT is then performed on each multiplication result X·V (S152). Similarly, following each FFT operation, each FFT bin is multiplied by a different set of coefficients called post-coefficients Uk via multiplier 126 (S154). The K results of the FFT operation are

number

number

[0082] Therefore, Figure 6 shows a K-rank approximation to matrix multiplication of DBF matrices. We assume that this approximation is close to the FFT matrix with some deviation even when the antenna array is not ideal. In one embodiment, the values ​​of the calibration coefficients before and after the FFT operation at each rank are calculated from the singular value decomposition of the distortion matrix C = B / F (i.e., element-wise or Hadamard division), where B is the DBF matrix and F is the ideal FFT matrix. This architecture is advantageous because the low rank required for full matrix multiplication is N 2 This achieves K×N×log2N operations, which is significantly smaller than the previous K×N×log2N operations. In practice, a rank of 4 turns out to be sufficient to compensate for 3D phase center misalignment and weak leakage between antenna channels. Therefore, the largest four values ​​are used and the remaining values ​​are zero. This is the optimal low-rank approximation of the distortion matrix C.

[0083] Note that when K=P, the approximation is perfect and any desired DBF matrix can be implemented. However, this is more computationally intensive than the matrix multiplication described above. In one embodiment, it is possible to select K<<P while still achieving a very accurate approximation. For an ideal ULA antenna, K=1 is sufficient, which is effectively equivalent to a single FFT operation where the input and output are multiplied by a constant vector. In an exemplary embodiment, a value of 4 is selected for K, which yields satisfactory results (i.e., four largest singular values). As shown in FIG. 6, the exemplary circuit utilizes four pre-multipliers, four FFT computation blocks, and four post-multipliers. It will be appreciated that the LRA mechanism of the present invention can be implemented using any desired rank depending on the particular application.

[0084] It is noted that in high-resolution radar, the amount of data that needs to be processed is typically enormous. In one embodiment, a matched filter, i.e., ML (Maximum Likelihood), is used to estimate the DOA. Other techniques are typically more computationally intensive. In the case of a ULA (Uniform Linear Array), the matched filter reduces to an FFT, especially in the ULA case. It is noted that due to the computational efficiency of the FFT, using other alternative methods may require more computations.

[0085] For example, consider a radar with 256 range bins, 1024 Doppler bins, for a total of 256 × 1024 = 262,144 range-Doppler bins, N = 128 azimuth bins, and M = 32 elevation bins. For an ideal ULA, it is possible to perform ML spatial processing (i.e., azimuth and elevation) via FFT with a complexity of N × M × log2(N × M) ≈ 50e3. However, for a non-ULA array, full matrix multiplications are required, resulting in a complexity of N × M × (N + M) ≈ 655e3, a more than 13-fold increase. Considering spatial processing for all range-Doppler bins, full matrix multiplications require approximately 158 billion or more operations per CPI. However, such a large number of calculations is impractical for low-cost consumer radar sensors.

[0086] Existing radars on the market typically have relatively small array sizes, such as 3x4, 6x8, or 12x16. For radars with small arrays, the processing load of estimating DOA is feasible and can even be performed in software. However, even for arrays of only 12x16 size, the number of calculations (i.e., multiplications) becomes so large that a hardware solution is required. Therefore, N 2 Conventional techniques using arithmetic are impractical for radars with larger array sizes, such as 48x48. The LRA method described herein, with a complexity of KxNxlog2N, comes very close in efficiency to the "pure" or "full" FFT DBF mechanism, but obtains performance reasonably close to the full matrix multiplication method.

[0087] To overcome these problems, in one embodiment, the DOA estimation mechanism of the present invention uses a non-ideal MLE (Maximum Likelihood Estimation) matched filter. The MLE matched filter provides a metric for the amount of energy coming from a particular direction. A vector "b" is generated specifically as in equation (1) for a particular azimuth or elevation angle, e.g., 25 degrees. These can be thought of as FFT coefficients, which are complex numbers corresponding to particular Fourier frequencies. These coefficients actually measure spatial frequency. In terms of radar antennas, spatial frequency infers the direction of the phased array. This is how the amount of energy received from the 25 degree direction is calculated. This is done for all the desired DOA angles X1...X being scanned. P , i.e., a linear combination of targets. This standard model is for phased arrays and is called the steering vector. The matched filter is just the complex conjugate.

[0088] Each vector "b" corresponds to a θ1,...θN direction. As shown previously, using a large number of virtual antennas (e.g., 128) makes the number of calculations very large. Radar resolution is related to the size of the aperture. A radar with wide element spacing has a large aperture, which means it has a narrow beamwidth, which improves resolution. However, regardless of the aperture size, the number of DOA calculations required is related to the quantity of elements.

[0089] It is noted that although the radar does not transmit from all elements simultaneously as in prior art phased array radar, the mathematics are the same. For example, signals transmitted simultaneously from multiple antennas are combined in the air. In DBF, signals received by multiple receive antenna elements are digitally added. Mathematically, this represents the same thing. Nevertheless, it is noted that the techniques described in this invention can be applied to any DBF setup, be it TD-MIMO, simultaneous transmit MIMO (such as OFDM), or even full ULA with a receiver having a single transmitter.

[0090] In one embodiment, the receive antenna elements are organized as ULAs. There is a different set of coefficients for each direction. In a standard beamforming scheme, using ideal ULAs, this matrix is k}, resulting in a DFT (Discrete Fourier Transform) matrix. A set of angles is selected and computed very efficiently for a given set of frequencies. In the FFT, it is desirable to take advantage of symmetries such as even and odd, positive and negative. A matched filter is computed for a specific set of frequencies. DOA estimation typically covers a specific set of frequencies. Therefore, for ULA and selected DOA groups, standard FFT computation can be used.

[0091] In one embodiment, the data is multiplied by a calibration window vector and then subjected to a DFT as a matrix multiplication, which forms the output of the system before SLAM.

[0092] Thus, in general, it is desired to compute the quantity B·X, which usually requires many computations, of the order O(N 2 ) Instead, the quantity B·X is approximated as before using an FFT operation. Note that the equation that approximates the calculation of B is not the FFT matrix, but is close enough. The distortion matrix is ​​computed using B, which is computed using any desired well-known technique. We assume that B is close to the FFT, and the distortion matrix is ​​computed using B. Next, we compute the SVD (singular value decomposition) to determine the pre- and post-coefficients.

[0093] In one embodiment, the DOA estimation mechanism can be elegantly and efficiently implemented in hardware, where X represents the virtual element array. After range / Doppler processing, the DOA estimate is calculated. For each range, Doppler, and row of the virtual array, azimuth processing is performed, representing 256 x 1024 x 128 = 33.5 million calculations. If elevation is also estimated, an additional 33.5 million calculations are performed, assuming the same resolution for azimuth and elevation. Note that the ranks of azimuth and elevation do not have to be the same and can be different, so different coefficients and a different rank may be used for elevation. The higher the rank, the better the approximation. The rank and other related parameters can be dynamically programmed and selected.

[0094] Beamforming matrix determination B In the following, we disclose a method for determining the beamforming matrix B by inverting the array response matrix A.

[0095] The beamforming matrix has the following relationship: By Remember to be defined I want to .

number

[0096] In general, each direction N can represent any chosen angle, and in Fourier beamforming (or, as in the present case, beamforming closely related to Fourier beamforming), the N possible values ​​of the angle depend on the input wavelength λ, the array spacing d, and the FFT length N. FFT is determined as follows:

number

[0097] The array response matrix is ​​defined by the inverse relationship of equation (7) as follows:

number

number

number

number

number

number

[0098] In index notation, the pth element is given by the following formula:

number

[0099] Note that the dimensions of A are P x ​​N.

[0100] We propose that if we construct a controlled measurement setup (e.g., an anechoic chamber, or an outdoor setup with minimal clutter) such that the environment has only a single point target at angle q with RCS R (set to one without loss of generality), then

number

[0101] Next, the inventor has the following:

number

number

[0102] The required setup, Equation (12), is an idealization that requires infinite SNR (signal-to-noise ratio), but for practical purposes is achievable in a typical anechoic chamber within a ~3 dB FOV (field of view) around the boresight. As we move further away from this FOV, the SNR decreases until Equation (12) no longer holds. The maximum angle at which Equation (12) holds is tθ FOV We define this set of angles within the FOV as

number

number

[0103] We define the total number of angles measured within this FOV as Q. Therefore, the measured array response A has post-measurement dimensions P × Q.

number

[0104] These Q measurement angles are not restricted to belong to any angular grid. We interpolate the rows of the measured array response to a grid of angles defined by the length P of the array, i.e., we define the Q angles

number

number

number

number

number

number

number

number

[0105] We assume that A, like the beamforming matrix B, is slightly different from the inverse Fourier transform matrix F.

number

[0106] therefore,

number

[0107] This completes the inversion because:

number

[0108] This results in

number

number

number

[0109] Calibration coefficient calculation 8 shows a diagram illustrating an exemplary method for calculating the pre- and post-coefficients. In one embodiment, the coefficients of the LRA architecture are computed in an optimal manner using the assumption that the DBF matrix is ​​similar to the FFT matrix. To do this, the following residual or distortion matrix C is first calculated (S130) as follows:

number

number

number

number

number

number

number

[0110] The result of this operation is an approximate residual matrix (S138) given as follows:

number

number

[0111] To show that the coefficients in equations (27) and (28) are the coefficients in the architecture and circuit of Figure 6, we compare the prior art DBF and LRA-based DOA estimation operations. On the one hand, the DBF approach calculates the next DOA estimate y for the scan direction θi from the antenna array response x as follows: i Calculate the following.

number

[0112] Next, the inventors K Approximating C by, we get:

number

[0113] Inserting equation (31) into equation (30) gives:

number

[0114] therefore,

number

number

number

number

[0115] Any arrangement of components that achieve that same functionality is effectively "associated" such that the desired functionality is achieved. Thus, any two components combined herein to achieve a particular functionality can be considered to be "associated" with one another such that the desired functionality is achieved, regardless of architecture or intermediate components. Similarly, any two such associated components can also be considered to be "operably coupled" or "operably connected" with one another such that the desired functionality is achieved.

[0116] Furthermore, those skilled in the art will recognize that the above boundaries between operations are merely exemplary. Multiple operations may be combined into a single operation, a single operation may be distributed among additional operations, and operations may be performed in at least partial overlapping fashion. Furthermore, alternative embodiments may include multiple instances of a particular operation, and the order of operations may be changed in various other embodiments.

[0117] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting. As used herein, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprise" and / or "comprising," when used herein, specify the presence of stated features, integers, steps, operations, elements, and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0118] In the claims, reference signs placed between parentheses should not be construed as limiting the claims. The use of introductory phrases such as "at least one" and "one or more" in a claim should not be construed as implying that the introduction of another claim with the indefinite article "a" or "an" limits any particular claim containing an element from such an introduced claim to an invention containing only one such element, even when that same claim contains the introductory phrase "one or more" or "at least one" and an indefinite article such as "a" or "an." The same applies to the use of definite articles. Unless otherwise specified, terms such as "first," "second," etc. are used to arbitrarily distinguish between the elements described by such terms. Thus, these terms are not necessarily intended to indicate a temporal or other priority of such elements. The mere fact that certain measures are recited in mutually different claims does not indicate that a combination of these measures cannot be advantageously used.

[0119] Corresponding structure, materials, acts, and equivalents of all means-plus-function or step-plus-function elements within the scope of the following claims are intended to include any structure, material, or act for performing a function in combination with other claimed elements as explicitly claimed. The description of the present invention has been presented for purposes of illustration and description, but is not intended to be exhaustive or to limit the invention to the form disclosed. Since numerous variations and modifications will readily occur to those skilled in the art, it is intended that the present invention not be limited to the limited number of embodiments described herein. It is therefore to be understood that all suitable variations, modifications, and equivalents may be employed within the spirit and scope of the present invention. The embodiments were chosen and described in order to best explain the principles of the invention and its practical application, and to enable those skilled in the art to understand the invention in various embodiments with various modifications as suited to the particular use intended.

Claims

1. 1. A method for estimating a direction of arrival (DOA) of a signal for use in a radar system, the method comprising: receiving input data including an antenna array response vector incorporating the energy of reflected signals received from a plurality of directions; multiplying each element of the input data by a set of k a priori coefficients to obtain a first plurality of results diag(Vk) x; performing a plurality of Fast Fourier Transform (FFT) operations on the first plurality of results to generate a second plurality of results F·diag(Vk)·x, the FFT operations being performed in a spatial domain over a channel space of the antenna array response vector; multiplying each element of the second plurality of results by k a posteriori coefficients to generate a third plurality of results; [Equation 1] and summing the third plurality of results to obtain an approximate DOA estimate representing the energy of the reflected signal received from each direction; Including, The value k is the rank of the approximation.

2. The method of claim 1 , further comprising computing a different approximate DOA estimate for each desired azimuth and / or elevation angle scanned.

3. The pre- and post-coefficients are expressed as a distortion matrix [Equation 2] is calculated from the singular value decomposition of [Equation 3] 2. The method of claim 1, wherein: H represents the transpose; N represents the length of the fast Fourier transform; and P represents the number of virtual array elements in the processing direction.

4. The method of claim 3 , wherein the distortion matrix is ​​determined by inverting the array response.

5. The method of claim 4 , wherein the array response is determined by making a series of measurements over multiple angles with a single target in a controlled environment.

6. 1. A method of estimating a direction of arrival (DOA) of a signal for use in a radar system, the method comprising: receiving input data x including an antenna array response vector incorporating the energy of reflected signals received from a plurality of directions; multiplying each element of the input data x by k sets of pre-coefficients Vk to obtain k first results diag(Vk) x; performing k Fast Fourier Transform operations F on the k first results to generate k second results F·diag(Vk)·x; Multiplying each element of the k second results by k sets of a posteriori coefficients Uk to obtain k third results. [Equation 4] and summing the k third results to obtain an approximate DOA estimate y representing the energy of the reflected signal received from each direction; Including, The value k is the rank of the approximation.

7. The method of claim 6 , further comprising computing a different approximate DOA estimate y for each desired azimuth and / or elevation angle scanned.

8. The method of claim 6 , wherein the approximation rank k is different from the azimuth and elevation DOA estimates.

9. The pre-coefficients Vk and the post-coefficients Uk are expressed as the distortion matrix [Equation 5] is calculated from the singular value decomposition of [Equation 6] 7. The method of claim 6, wherein: H denotes the transpose; N denotes the length of the fast Fourier transform; and P denotes the number of virtual array elements in the processing direction.

10. The method of claim 9 , wherein the distortion matrix is ​​determined by inverting the array response.

11. The method of claim 10 , wherein the array response is determined by making a series of measurements over multiple angles with a single target in a controlled environment.

12. 1. An apparatus for estimating a direction of arrival (DOA) of a signal for use in a radar system, the apparatus comprising: Radar signal processing circuit for receiving a receive antenna array response vector x incorporating the energy of reflected signals received from multiple directions Including, The radar signal processing circuit includes: multiplying each element of the antenna array response by k sets of pre-coefficients V to obtain k first results diag(V)·x; performing k Fast Fourier Transform operations F on the k first results to generate k second results F·diag(Vk)·x; Multiplying each element of the k second results by k sets of a posteriori coefficients Uk to obtain k third results. [Equation 7] and summing the k third results to obtain an approximate DOA estimate y representing the energy of the reflected signal received from each direction; It works like this: The value k is the rank of the approximation.

13. 13. The apparatus of claim 12, wherein the radar signal processing circuitry is operative to compute a different approximate DOA estimate y for each desired azimuth and / or elevation angle scanned.

14. The apparatus of claim 12 , wherein the approximation rank k is distinct from the azimuth and elevation DOA estimates.

15. The pre-coefficients Vk and the post-coefficients Uk are expressed as the distortion matrix [Equation 8] is calculated from the singular value decomposition of [Equation 9] 13. The apparatus of claim 12, wherein: Ĥ is a diagonal matrix of singular values ​​of C, H represents the transpose, N represents the length of the fast Fourier transform, and P represents the number of virtual array elements in the processing direction.

16. The apparatus of claim 15 , wherein the distortion matrix is ​​determined by inverting an array response.

17. 17. The apparatus of claim 16, wherein the array response is determined by making a series of measurements over multiple angles with a single target in a controlled environment.

18. A radar sensor for an automobile, comprising:

1. A printed circuit board (PCB) assembly comprising: a plurality of transmitting antennas fabricated on one side of the PCB assembly; a plurality of receiving antennas fabricated on opposite sides of the PCB assembly; a transceiver coupled to the plurality of transmitting antennas and the plurality of receiving antennas, the transceiver operative to generate and transmit transmit signals to the plurality of transmitting antennas and to receive reflected wave signals to the plurality of receiving antennas; a PCB assembly including: a radar signal processing circuit coupled to the transceiver, comprising: receiving input data x including an antenna array response vector incorporating the energy of reflected signals received from a plurality of directions; multiplying each element of the input data x by k sets of pre-coefficients Vk to obtain k first results diag(Vk) x; performing k Fast Fourier Transform operations F on the k first results to generate k second results F·diag(Vk)·x; Multiplying each element of the k second results by k sets of a posteriori coefficients Uk to obtain k third results. [Equation 10] and summing the k third results to obtain an approximate DOA estimate y representing the energy of the reflected signal received from each direction; a radar signal processing circuit operative to perform Including, The value k is the rank of the approximation.

19. The pre-coefficients Vk and the post-coefficients Uk are expressed as the distortion matrix [0011] is calculated from the singular value decomposition of [0012] 20. The sensor of claim 18, wherein: ∑ i = ...

Images