The present disclosure relates to a radar apparatus.
Radar apparatuses that use radar transmission signals with short wavelengths, including microwaves or millimeter-waves that realize high resolution, have been under study in recent years. To improve safety in outdoor spaces, a demand has arisen for the development of radar apparatuses (referred to as, for example, wide-angle radar apparatuses) for sensing vehicles, as well as small objects such as pedestrians, within wide-angle ranges.
Examples of the configuration of a radar apparatus having a wide-angle sensing range include a configuration that uses a method called direction of arrival (DOA) estimation. In DOA estimation, reflected waves from a target are received by an array antenna constituted by a plurality of antennas (or also referred to as antenna elements), and the direction of arrival (or referred to as the angle of arrival) of the reflected waves is estimated using a signal processing algorithm based on reception phase difference as to element spacing (antenna spacing). Techniques for DOA estimation include, for example, the Fourier method, or examples of the technique for realizing high resolution include the Capon method, multiple signal classification (MUSIC), and estimation of signal parameters via rotational invariance techniques (ESPRIT).
Further, for example, a configuration of a radar apparatus (also referred to sometimes as a multiple input multiple output (MIMO) radar) has been proposed that includes a plurality of antennas (array antenna) on the transmitting side in addition to the receiving side such that the transmit and receive array antennas are used to perform signal processing to perform beam scanning (see, for example, NPL 1).
PTL 1
NPL 1
However, methods for sensing a target object (or target) using a radar apparatus (for example, a MIMO radar) have not been fully studied.
A non-limiting exemplary embodiment of the present disclosure provides a radar apparatus with improved target-object sensing accuracy.
A radar apparatus according to one example of the present disclosure includes: a plurality of transmit antennas that transmit a plurality of transmission signals using a multiplexing transmission; and a transmission circuit that applies phase rotation amounts corresponding to combinations of Doppler shift amounts and code sequences to the plurality of transmission signals, wherein each of the combinations of the Doppler shift amounts and the code sequences has at least one different from other combination, and wherein the number of multiplexes of the code sequence corresponding to at least one of the Doppler shift amounts in the combinations is different from the number of multiplexing of code sequences corresponding to the remaining Doppler shift amount.
It should be noted that general or specific embodiments may be implemented as a system, an apparatus, a method, an integrated circuit, a computer program or a storage medium, or any selective combination of the system, the apparatus, the method, the integrated circuit, the computer program, and the storage medium.
According to an exemplary embodiment of the present disclosure, the target-object sensing accuracy of a radar apparatus can be improved.
Additional benefits and advantages of the disclosed exemplary embodiments will become apparent from the specification and drawings. The benefits and/or advantages may be individually obtained by the various embodiments and features of the specification and drawings, which need not all be provided in order to obtain one or more of such benefits and/or advantages.
A MIMO radar transmits signals (radar transmission waves) multiplexed by, for example, time division, frequency division, or code division from a plurality of transmit antennas (or referred to as a transmit array antenna), receives signals (radar reflected waves) reflected from a nearby object using a plurality of receive antennas (or referred to as a receive array antenna), and separates and receives the multiplexed transmission signals from the respective reception signals. Through the process described above, the MIMO radar can acquire the propagation channel response, which is given by the product of the number of transmit antennas and the number of receive antennas, and the reception signals are used as a virtual receive array to perform array signal processing.
In the MIMO radar, furthermore, an arrangement with appropriate element spacing in the transmit and receive array antennas can extend the antenna aperture in a virtual fashion and improve the angular resolution.
For example, PTL 1 discloses a MIMO radar (hereinafter referred to as a “time-division multiplexing MIMO radar”) that uses, as a multiplexing transmission method for the MIMO radar, time-division multiplexing transmission to transmit a signal from each transmit antenna at a shifted transmission time. Time-division multiplexing transmission can be implemented with a simpler configuration than frequency multiplexing transmission or code multiplexing transmission. In time-division multiplexing transmission, furthermore, a sufficient width of the intervals between transmission times can maintain proper orthogonality between transmission signals. The time-division multiplexing MIMO radar outputs transmission pulses, which are an example of transmission signals, while sequentially switching the transmit antennas in a defined period. The time-division multiplexing MIMO radar receives, at a plurality of receive antennas, signals that are the transmission pulses reflected by an object, performs a correlation process between the reception signals and the transmission pulses, and then performs, for example, spatial fast Fourier transform (FFT) processing (processing for direction of arrival estimation of the reflected waves).
The time-division multiplexing MIMO radar sequentially switches the transmit antennas, from which transmission signals (for example, transmission pulses or radar transmission waves) are to be transmitted, in a defined period. In time-division multiplexing transmission, therefore, it may take a longer time to complete transmission of transmission signals from all the transmit antennas than in frequency division transmission or code division transmission. Accordingly, for example, as in PTL 2, when transmission signals are transmitted from the respective transmit antennas and the Doppler frequency (i.e., the relative velocity of the target) is detected from their reception phase changes, the time interval for observing the reception phase changes (for example, sampling interval) is increased to apply Fourier frequency analysis to detect the Doppler frequency. This reduces the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing (folding) (i.e., the range of detectable relative velocities of the target).
If a reflected wave signal outside the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing (in other words, the range of relative velocities) is assumed to come from the target, the radar apparatus is unable to identify whether the reflected wave signal is the aliasing component, causing the ambiguity (uncertainty) of the Doppler frequency (in other words, the relative velocity of the target).
For example, when the radar apparatus transmits transmission signals (transmission pulses) while sequentially switching Nt transmit antennas in period Tr, it takes a transmission time given by Tr×Nt to complete transmission of the transmission signals from all the transmit antennas. As a result of repeating this time-division multiplexing transmission operation Nc times and applying Fourier frequency analysis to detect the Doppler frequency, the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing is given by ±1/(2Tr×Nt) from the sampling theorem. Accordingly, the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing decreases as the number Nt of transmit antennas increases, and the ambiguity of the Doppler frequency is likely to occur even for lower relative velocities.
In the time-division multiplexing MIMO radar, the ambiguity of the Doppler frequency, described above, is likely to occur for lower relative velocities. In the following, focus is on a method for simultaneously transmitting transmission signals from a plurality of transmit antennas in a multiplexed manner, as an example.
Examples of the method for simultaneously transmitting transmission signals from a plurality of transmit antennas in a multiplexed manner include a method (hereinafter referred to as Doppler multiplexing transmission) for transmitting signals so that a plurality of transmission signals can be separated in the Doppler frequency domain on the receiving side (see, for example, NPL 3).
In Doppler multiplexing transmission, on the transmitting side, transmission signals are simultaneously transmitted from a plurality of transmit antennas in such a manner that, for example, with respect to a transmission signal to be transmitted from a reference transmit antenna, transmission signals to be transmitted from transmit antennas different from the reference transmit antenna are given Doppler shift amounts greater than the Doppler frequency bandwidth of reception signals. In Doppler multiplexing transmission, on the receiving side, filtering is performed in the Doppler frequency domain to separate and receive the transmission signals transmitted from the respective transmit antennas.
In Doppler multiplexing transmission, simultaneous transmission of transmission signals from a plurality of transmit antennas can reduce the time interval for observing reception phase changes to apply Fourier frequency analysis to detect the Doppler frequency (or relative velocity), compared with time-division multiplexing transmission. In Doppler multiplexing transmission, however, since filtering is performed in the Doppler frequency domain to separate the transmission signals of the respective transmit antennas, the effective Doppler frequency bandwidth per transmission signal is restricted.
For example, Doppler multiplexing transmission in which a radar apparatus transmits transmission signals from Nt transmit antennas in period Tr will be described. As a result of repeating this Doppler multiplexing transmission operation Nc times and applying Fourier frequency analysis to detect the Doppler frequency (or relative velocity), the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing is given by ±1/(2×Tr) from the sampling theorem. That is, in Doppler multiplexing transmission, the Doppler frequency range over which the Doppler frequency can be detected without causing aliasing is increased by Nt times compared with time-division multiplexing transmission (for example, ±1/(2Tr×Nt)).
In Doppler multiplexing transmission, as described above, filtering is performed in the Doppler frequency domain to separate transmission signals. Accordingly, the effective Doppler frequency bandwidth per transmission signal is restricted to 1/(Tr×Nt), and thus a Doppler frequency range similar to that in time-division multiplexing transmission is obtained. In Doppler multiplexing transmission, furthermore, in a Doppler frequency band exceeding the effective Doppler frequency range per transmission signal, interference with signals in the Doppler frequency band of any other transmission signal different from the transmission signal may lead to failure to correctly separate the transmission signals.
Accordingly, an exemplary embodiment of the present disclosure describes a method for extending the range of Doppler frequencies at which ambiguity does not occur in Doppler multiplexing transmission. With this method, a radar apparatus according to an exemplary embodiment of the present disclosure can improve target-object sensing accuracy over a wider Doppler frequency range.
Embodiments of the present disclosure will be described in detail with reference to the drawings. In the embodiments, the same constituent elements are identified with the same numerals, and a description thereof is omitted because of redundancy.
The following describes a configuration of a radar apparatus (in other words, MIMO radar configuration) having a transmitting branch in which multiplexed different transmission signals are simultaneously sent from a plurality of transmit antennas, and a receiving branch in which the transmission signals are separated and subjected to reception processing.
The following also describes, as an example, a configuration of a radar scheme (also referred to as, for example, chirp pulse transmission (fast chirp modulation)) that uses frequency-modulated pulse waves such as chirp pulses. Note that the modulation scheme is not limited to that for frequency modulation. For example, an exemplary embodiment of the present disclosure is also applicable to a radar scheme that uses a pulse compression radar configured to transmit a pulse train after performing phase modulation or amplitude modulation.
Further, the radar apparatus performs Doppler multiplexing transmission. In addition, in Doppler multiplexing transmission, the radar apparatus performs coding (for example, code division multiplexing (CDM)) on signals (hereinafter referred to as “Doppler-multiplexed transmission signals”) with different phase rotations (in other words, phase shifts) applied, the number of which corresponds to the number of Doppler multiplexes, and transmits the coded signals in a multiplexed manner (hereinafter referred to as “coded Doppler multiplexing”).
[Configuration of Radar Apparatus]
Radar apparatus 10 includes radar transmitter (transmitting branch) 100 and radar receiver (receiving branch) 200.
Radar transmitter 100 generates radar signals (radar transmission signals) and transmits the radar transmission signals in a defined transmission period using a transmit array antenna made up of plural (for example, Nt) transmit antennas 108.
Radar receiver 200 receives reflected wave signals, which are radar transmission signals reflected by a target object (target) (not illustrated), using a receive array antenna made up of plural receive antennas 202-1 to 202-Na. Radar receiver 200 performs signal processing on the reflected wave signals received by respective receive antennas 202 to, for example, detect the presence or absence of the target object or estimate the distance of arrival, the Doppler frequency (in other words, relative velocity), and the direction of arrival of the reflected wave signals, and outputs information about the estimation results (in other words, position measurement information).
The target object is an object to be detected by radar apparatus 10. Examples of the target object include vehicles (including four-wheel and two-wheel vehicles), a person, and a block or a curb.
[Configuration of Radar Transmitter 100]
Radar transmitter 100 includes radar transmission signal generator 101, phase rotation amount setter 104, phase rotators 107, and transmit antennas 108.
Radar transmission signal generator 101 generates a radar transmission signal. Radar transmission signal generator 101 includes, for example, modulated signal emitter 102 and voltage controlled oscillator (VCO) 103. The following describes each of the components of radar transmission signal generator 101.
For example, as illustrated in
VCO 103 outputs, based on the radar transmission signals (modulated signals) output from modulated signal emitter 102, frequency-modulated signals (hereinafter referred to as, for example, frequency chirp signals or chirp signals) to phase rotators 107 and radar receiver 200 (mixer section 204 described below).
Phase rotation amount setter 104 sets phase rotation amounts for phase rotators 107 (in other words, phase rotation amounts corresponding to coded Doppler multiplexing transmission). Phase rotation amount setter 104 includes, for example, Doppler shift setter 105 and coder 106.
Doppler shift setter 105 sets, for example, a phase rotation amount corresponding to a Doppler shift amount to be applied to each radar transmission signal (for example, chirp signal).
Coder 106 sets a phase rotation amount corresponding to coding. Coder 106 calculates phase rotation amounts for phase rotators 107, based on, for example, the phase rotation amounts output from Doppler shift setter 105 and the phase rotation amount corresponding to coding, and outputs the phase rotation amounts to phase rotators 107. Further, coder 106 outputs, for example, information about code sequences used for coding (for example, elements of orthogonal code sequences) to radar receiver 200 (for example, output switching section 209).
Phase rotators 107 applies the phase rotation amounts input from coder 106 to the chirp signals input from VCO 103 and outputs signals subjected to phase rotation to transmit antennas 108. For example, each phase rotator 107 includes a phase shifter, a phase modulator, and so on (not illustrated). The output signals of phase rotators 107 are amplified to defined transmission power and are then radiated into a space from the respective transmit antennas 108. In other words, phase rotation amounts corresponding to Doppler shift amounts and orthogonal code sequences are applied to radar transmission signals, which are then transmitted from plural transmit antennas 108 in a multiplexed manner.
Next, an example method for setting phase rotation amounts using phase rotation amount setter 104 will be described.
Doppler shift setter 105 sets phase rotation amount ϕndm for applying Doppler shift amount DOPndm and outputs phase rotation amount ϕndm to coder 106. Here, ndm=1, . . . , NDM. NDM denotes a set number of different Doppler shift amounts and is hereinafter referred to as the “number of Doppler multiplexes”.
In Radar apparatus 10, since coding performed by coder 106 is also used, the number of Doppler multiplexes NDM may be set to be smaller than the number Nt of transmit antennas 108 used for multiplexing transmission. The number of Doppler multiplexes NDM is greater than or equal to 2.
Doppler shift amounts DOP1, DOP2, . . . , and DOPDM are assigned different phase rotation amounts by, for example, dividing a phase rotation range greater than or equal to 0 and less than 2π. For example, phase rotation amount ϕndm for applying Doppler shift amount DOPndm is assigned, as given by following equation 1. In the following, the angle is expressed in radian.
In equation 1, for example, in a case where the number of Doppler multiplexes NDM is equal to 2, phase rotation amount ϕ1 for applying Doppler shift amount DOP1 is equal to 0, and phase rotation amount ϕ2 for applying Doppler shift amount DOP2 is equal to π. Likewise, in equation 1, for example, in a case where the number of Doppler multiplexes NDM is equal to 4, phase rotation amount ϕ1 for applying Doppler shift amount DOP1 is equal to 0, phase rotation amount ϕ2 for applying Doppler shift amount DOP2 is equal to π/2, phase rotation amount ϕ3 for applying Doppler shift amount DOP3 is equal to ϕ, and phase rotation amount ϕ4 for applying Doppler shift amount DOP4 is equal to 3π/2. In other words, intervals of phase rotation amounts ϕndm for applying Doppler shift amounts DOPndm are equal.
The assignment of phase rotation amounts for applying Doppler shift amounts DOP1, DOP2, . . . , and DOPDM is not limited to that in this assignment method. For example, the assignment of phase rotation amounts given by equation 1 may be shifted. For example, phase rotation amounts may be assigned such that ϕndm=2π(ndm)/NDM. Alternatively, phase rotation amounts ϕ1, ϕ2, . . . , and ϕDM may be randomly assigned to Doppler shift amounts DOP1, DOP2, . . . , and DOPDM using a phase-rotation-amount assignment table.
Coder 106 sets a phase rotation amount based on one or a plurality of orthogonal code sequences less than or equal to NCM for each of phase rotation amounts ϕ1, . . . , and ϕNDM for applying NDM Doppler shift amounts output from Doppler shift setter 105. Further, coder 106 sets phase rotation amounts based on both the Doppler shift amounts and the orthogonal code sequences, that is, “coded Doppler phase rotation amounts” for generating coded Doppler multiplexed signals, and outputs the phase rotation amounts to phase rotators 107.
The following describes an example of the operation of coder 106.
For example, coder 106 uses NCM orthogonal code sequences, the number of which is equal to the number of codes (in other words, the number of code multiplexes) with code length Loc.
In the following, NCM orthogonal code sequences with code length Loc are represented by Codencm={OCncm(1), OCncm(2), . . . , OCncm(Loc)}. OCncm(noc) denotes the noc-th code element in the ncm-th orthogonal code sequence Codencm. Here, noc denotes the index of a code element, and noc=1, . . . , Loc.
The orthogonal code sequences used in coder 106 are, for example, codes that are orthogonal (uncorrelated) to each other. For example, orthogonal code sequences may be Walsh-Hadamard codes. In this case, code length Loc used to generate NCM orthogonal code sequences, the number of which is equal to the number of codes, is given by following equation 2.
[2]
Loc=2ceil[log
Here, ceil[x] denotes the operator (ceiling function) that outputs a minimum integer greater than or equal to real number x.
For example, in a case where NCM=2, code length Loc of Walsh-Hadamard codes is equal to 2, and orthogonal code sequences are represented by Code1={1, 1} and Code2={1, −1}. When a code element in an orthogonal code sequence is 1, 1=exp(j0), with the phase thereof being 0. When a code element in an orthogonal code sequence is −1, −1=exp(jπ), with the phase thereof being π.
Further, for example, in a case where NCM=4, code length Loc is equal to 4, and orthogonal code sequences are represented by Code1={1, 1, 1, 1}, Code2={1, −1, 1, −1}, Code3={1, 1, −1, −1}, and Code4={1, −1, −1, 1}.
Code elements in an orthogonal code sequence are not limited to real numbers and may include complex number values. For example, an orthogonal code sequence Codencm given by following equation 3 may be used. Here, ncm=1, . . . , NCM. In this case, a code length used to generate NCM orthogonal code sequences, the number of which is equal to the number of codes, is represented by Loc=NCM.
For example, in a case where NCM=3, code length Loc is equal to 3 (=NCM), and coder 106 generates orthogonal code sequences represented by Code1={1, 1, 1}, Code2={1, exp(j2π/3), exp(j4π/3)}, and Code3={1, exp(−j2π/3), exp(−j4π/3)}.
For example, in a case where NCM=4, code length Loc is equal to 4 (=NCM), and coder 106 generates orthogonal code sequences represented by Code1={1, 1, 1, 1}, Code2={1, j, −1, −j}, Code3={1,−1,1,−1}, and Code4={1, −j, −1, j}. Here, j is the imaginary unit.
In Coder 106, the number of code multiplexes (hereinafter referred to as the number of coded Doppler multiplexes) for coding a Doppler multiplexed signal that uses the ndm-th Doppler shift amount DOPndm output from Doppler shift setter 105 is represented by “NDOP_CODE(ndm)”. Here, ndm=1, . . . , NDM.
Coder 106 sets the number of coded Doppler multiplexes NDOP_CODE(ndm) so that, for example, the sum of the numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) for coding Doppler multiplexed signals is equal to the number Nt of transmit antennas 108 used for multiplexing transmission. In other words, coder 106 sets the number of coded Doppler multiplexes NDOP_CODE(ndm) so as to satisfy following equation 4. This enables radar apparatus 10 to perform multiplexing transmission in the Doppler domain and the code domain (hereinafter referred to as coded Doppler multiplexing transmission) using Nt transmit antennas 108.
Here, coder 106 sets, for example, the numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) so as to include different numbers of coded Doppler multiplexes in the range greater than or equal to 1 and less than or equal to NCM. For example, coder 106 sets the numbers of coded Doppler multiplexes such that not all of the numbers of coded Doppler multiplexes are set to NCM, which is equal to the number of codes, but at least one of the numbers of coded Doppler multiplexes is set to be smaller than NCM. In other words, coder 106 sets the numbers of coded Doppler multiplexes for Doppler multiplexed signals to be non-uniform. With this setting, for example, radar apparatus 10 can individually separate and receive signals transmitted from plural transmit antennas 108 in a coded Doppler multiplexed manner through reception processing described below.
Coder 106 sets, in the m-th transmission period Tr, coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) given by following equation 5 for phase rotation amount ϕndm for applying the ndm-th Doppler shift amount DOPndm and outputs coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) to phase rotator 107.
Here, the subscript “ndop_code(ndm)” represents an index less than or equal to the number of coded Doppler multiplexes NDOP_CODE(ndm) for phase rotation amount ϕndm for applying Doppler shift amount DOPndm. For example, ndop_code(ndm)=1, . . . , NDOP_CODE(ndm). Further, angle[x] denotes the operator that outputs a radian phase of real number x. For example, angle[1]=0, angle[−1]=π, angle[j]=π/2, and angle[−j]=−π/2. Further, floor[x] denotes the operator that outputs a maximum integer not greater than real number x, where j is the imaginary unit.
For example, as given by equation 5, coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) provides a constant phase rotation amount for applying Doppler shift amount DOPndm (for example, the first term in equation 5) for the duration of Loc transmission periods, the number of which is equal to the code length used for coding, and applies a corresponding phase rotation amount to each of Loc code elements OCndop_code(ndm)(1), . . . , and OCndop_code(ndm)(Loc) of code Codendop_code(ndm) used for coding (the second term in equation 5).
Further, coder 106 outputs, for each transmission period (Tr), orthogonal code element index OC_INDEX to radar receiver 200 (output switching section 209 described below). OC_INDEX represents an orthogonal code element index that indicates the elements of orthogonal code sequence Codendop_code(ndm). OC_INDEX cyclically varies in the range of 1 to Loc for each transmission period (Tr), as given by following equation 6.
[6]
OC_INDEX=mod(m−1,Loc)+1 (Equation 6)
Here, mod(x, y) denotes a modulo operator and is a function that outputs the remainder of x divided by y. Further, m=1, . . . , Nc. Nc denotes the number of transmission periods used for radar position determination (hereinafter referred to as the “number of transmissions of radar transmission signals”). The number Nc of transmissions of radar transmission signals is set to an integer multiple (Ncode times) of Loc. For example, Nc=Loc×Ncode.
Next, an example method for setting the numbers of coded Doppler multiplexes NDOP_CODE(ndm) for Doppler multiplexed signals to be non-uniform using coder 106 will be described.
For example, coder 106 sets the number of orthogonal code sequences (in other words, the number of code multiplexes or the number of codes) NCM satisfying the condition below. For example, the number of orthogonal code sequences NCM and the number of Doppler multiplexes NDM satisfy the following relationship for the number Nt of transmit antennas 108 used for multiplexing transmission.
(Number of orthogonal code sequences NCM)×(number of Doppler multiplexes NDM)>number Nt of transmit antennas used for multiplexing transmission
For example, among the numbers of orthogonal code sequences NCM and the numbers of Doppler multiplexes NDM satisfying the above-described condition, the use of a combination having a smaller value of the product (NCM×NDM) is desirable in terms of both characteristics and complexity of circuit configuration. Note that among the numbers of orthogonal code sequences NCM and the numbers of Doppler multiplexes NDM satisfying the above-described condition, a combination having a smaller value of the product (NCM×NDM) is not restrictive, and any other combination may be applied.
For example, in a case where Nt=3, the combination of NDM=2 and NCM=2 is desirable.
In this case, the assignment of Doppler shift amounts DOP1 and DOP2 and orthogonal codes Code1 and Code2 is determined in accordance with, for example, as illustrated in
For example,
In
For example, in a case where Nt=4 or 5, the combination of NDM=3 and NCM=2 or the combination of NDM=2 and NCM=3 is desirable.
For example,
In
For example,
In
For example, in a case where Nt=4, NDM=2, and NCM=3, if, as illustrated in
In this embodiment, in contrast, for example, in a case where Nt=4, NDM=2, and NCM=3, the numbers of coded Doppler multiplexes NDOP_CODE are set to be non-uniform for Doppler shift amounts DOP1 and DOP2, such as NDOP_CODE(1)=3 and NDOP_CODE(2)=1, or NDOP_CODE(1)=1 and NDOP_CODE(2)=3, as illustrated in
For example, in a case where Nt=6 or 7, the combination of NDM=4 and NCM=2 or the combination of NDM=2 and NCM=4 is desirable.
For example,
In
For example, as illustrated in
Likewise, for example, in a case where Nt=8, the combination of NDM=3 and NCM=3 or the combination of nDM=5 and NCM=2 is desirable. For example, in a case where Nt=9, the combination of NDM=5 and NCM=2 is desirable. For example, in a case where Nt=10, the combination of NDM=6 and NCM=2 or the combination of NDM=4 and NCM=3 is desirable. The number Nt of transmit antennas 108 is not limited to that in the example described above, and an exemplary embodiment of the present disclosure is also applicable to Nt=11 or more.
Next, an example of how coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) is set will be described.
For example, a description will be given of a case where in coder 106, the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NCM=2, and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used. In this case, for example, if the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1 and NDOP_CODE(2)=2, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 7 to 9 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) to phase rotators 107.
[7]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,ϕ1,ϕ1,2ϕ1,2ϕ1,3ϕ1,3ϕ1, . . . } (Equation 7)
[8]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,ϕ2,ϕ2,2ϕ2,2ϕ2,3ϕ2,3ϕ2, . . . } (Equation 8)
[9]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }={0,π,ϕ2,ϕ2+π,2ϕ2,2ϕ2+π,3ϕ2,3ϕ2+π, . . . } (Equation 9)
Here, as an example, the phase rotation amount for applying Doppler shift amount DOPndm is given by ϕndm=2π(ndm−1)/NDM in equation 1, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, and phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π, are used. In this case, coder 106 sets coded Doppler phase rotation amounts Φ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 10 to 12 and outputs coded Doppler phase rotation amounts ψ1, 1(m), Φ1, 2(m), and ψ2, 2(m) to phase rotators 107. Here, m=1, . . . , Nc. Here, a modulo operation for 2π is performed, and results are expressed in radians ranging from 0 or more to less than 2π (the same applies to the following description).
[10]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,0,0,0,0,0,0, . . . } (Equation 10)
[11]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,π,π,0,0,π,π, . . . } (Equation 11)
[12]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }={0,π,π,0,0,π,π,0, . . . } (Equation 12)
As given by equations 10 to 12, when a phase rotation amount is set to ϕndm=2π(ndm−1)/NDM, into which 2π is equally divided, coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) are changed in transmission periods given by NDM×NCM=2×2=4.
As another example, the phase rotation amount for applying Doppler shift amount DOPndm may be set to ϕndm=2π(ndm)/NDM, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to π, and phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to 0, may be used. In this case, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 13 to 15 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) to phase rotators 107. Here, m=1, . . . , Nc.
[13]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,π,π,0,0,π,π, . . . } (Equation 13)
[14]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,0,0,0,0,0,0, . . . } (Equation 14)
[15]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }={0,π,0,π,0,π,0,π, . . . } (Equation 15)
As given by equations 10 to 12 or equations 13 to 15, the number of phases (for example, two, namely, 0 and π) used for a phase rotation amount (for example, a phase rotation amount for applying a Doppler shift amount) is smaller than the number of transmit antennas 108 used for multiplexing transmission, namely, Nt=3. In other words, as given by equations 10 to 12 or equations 13 to 15, the number of phases (for example, two, namely, 0 and π) used for a phase rotation amount for applying a Doppler shift amount is equal to the number of Doppler shift amounts used for multiplexing transmission (in other words, the number of Doppler multiplexes) NDM=2.
Further, for example, a description will be given of a case where in coder 106, the number of transmit antennas used for multiplexing transmission Nt=6, the number of Doppler multiplexes NDM=4, NCM=2, and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used. In this case, for example, if the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1, NDOP_CODE(2)=1, NDOP_CODE(3)=2, and NDOP_CODE(4)=2, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 16 to 21 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 107. Here, m=1, . . . , Nc.
[16]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,ϕ1,ϕ1,2ϕ1,2ϕ1,3ϕ1,3ϕ1, . . . } (Equation 16)
[17]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,ϕ2,ϕ2,2ϕ2,2ϕ2,3ϕ2,3ϕ2, . . . } (Equation 17)=
[18]
{ψ1,3(1),ψ1,3(2),ψ1,3(3),ψ1,3(4),ψ1,3(5),ψ1,3(6),ψ1,3(7),ψ1,3(8), . . . }={0,0,ϕ3,ϕ3,2ϕ3,2ϕ3,3ϕ3,3ϕ3, . . . } (Equation 18)
[19]
{ψ2,3(1),ψ2,3(2),ψ2,3(3),ψ2,3(4),ψ2,3(5),ψ2,3(6),ψ2,3(7),ψ2,3(8), . . . }={0,π,ϕ3,ϕ3+π,2ϕ3,2ϕ3+π,3ϕ3,3ϕ3+π, . . . } (Equation 19)
[20]
{ψ1,4(1),ψ1,4(2),ψ1,4(3),ψ1,4(4),ψ1,4(5),ψ1,4(6),ψ1,4(7),ψ1,4(8), . . . }={0,0,ϕ4,ϕ4,2ϕ4,2ϕ4,3ϕ4,3ϕ4, . . . } (Equation 20)
[21]
{ψ2,4(1),ψ2,4(2),ψ2,4(3),ψ2,4(4),ψ2,4(5),ψ2,4(6),ψ2,4(7),ψ2,4(8), . . . }={0,π,ϕ4,ϕ4+π,2ϕ4,2ϕ4+π,3ϕ4,3ϕ4+π, . . . } (Equation 21)
Here, as an example, the phase rotation amount for applying Doppler shift amount DOPndm is given by ϕndm=2π(ndm−1)/NDM, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π/2, phase rotation amount ϕ3 for applying Doppler shift amount DOP3, which is equal to π, and phase rotation amount ϕ4 for applying Doppler shift amount DOP4, which is equal to 3π/2, are used. In this case, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 22 to 27 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 107. Here, m=1, . . . , Nc.
As given by equations 22 to 27, when a phase rotation amount is set to ϕndm=2π(ndm−1)/NDM, into which 2π is equally divided, coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) are changed in transmission periods given by NDM×NCM=4×2=8.
As given by equations 22 to 27, furthermore, the number of phases (for example, four, namely, 0, π/2, π, and 3π/2) used for a phase rotation amount (for example, a phase rotation amount for applying a Doppler shift amount) is smaller than the number of transmit antennas 108 used for multiplexing transmission, namely, Nt=6. In other words, as given by equations 22 to 27, the number of phases (for example, four, namely, 0, π/2, π, and 3π/2) used for a phase rotation amount for applying a Doppler shift amount is equal to the number of Doppler shift amounts used for multiplexing transmission (in other words, the number of Doppler multiplexes) NDM=4.
While the description has been given of, as an example, the setting of phase rotation amounts in a case where the number Nt of transmit antennas 108 is equal to 3 and the number of Doppler multiplexes NDM is equal to 2 and in a case where the number Nt of transmit antennas 108 is equal to 6 and the number of Doppler multiplexes NDM is equal to 4, the number Nt of transmit antennas 108 and the number of Doppler multiplexes NDM are not limited to the values described above. For example, the number of phases used for a phase rotation amount may be set to be smaller than the number Nt of transmit antennas 108 used for multiplexing transmission, regardless of the number Nt of transmit antennas 108. Further, the number of phases used for a phase rotation amount for applying a Doppler shift amount may be equal to the number NDM of Doppler shift amounts used for multiplexing transmission.
The foregoing description has been given of a method for setting phase rotation amounts using phase rotation amount setter 104.
In
The sum of the numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) is set to be equal to the number Nt of transmit antennas 108, and Nt coded Doppler phase rotation amounts are respectively input to Nt phase rotators 107.
Each of Nt phase rotators 107 applies coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) input thereto to a chirp signal output from radar transmission signal generator 101 for each transmission period Tr. The outputs of Nt phase rotators 107 (referred to as, for example, coded Doppler multiplexed signals) are amplified to defined transmission power and are then radiated into a space from Nt transmit antennas 108 in a transmit array antenna section.
In the following, phase rotator 107 that applies coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) is as represented by “phase rotator PROT #[ndop_code(ndm), ndm]”. Likewise, transmit antenna 108 that radiates the output of phase rotator PROT #[ndop_code(ndm), ndm] into a space is represented by “transmit antenna Tx #[ndop_code(ndm), ndm]”. Here, ndm=1, . . . , NDM, and ndop_code(ndm)=1, . . . , NDOP_CODE(ndm).
For example, a description will be given of a case where when the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NCM4=2, orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, and the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1 and NDOP_CODE(2)=2. In this case, coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) are output from coder 106 to phase rotators 107 in respective transmission periods.
For example, phase rotator PROT #[1, 1] applies, for each transmission period, phase rotation amount ψ1, 1(m) as given by following equation 28 to a chirp signal generated by radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[1, 1] is output from transmit antenna Tx #[1, 1]. Here, cp(t) denotes a chirp signal for each transmission period.
[28]
exp[jψ1,1(1)]cp(t),exp[jψ1,1(2)]cp(t),exp[jψ1,1(3)]cp(t), . . . ,exp[jψ1,1(Nc)]cp(t) (Equation 28)
Likewise, phase rotator PROT #[1, 2] applies, for each transmission period, phase rotation amount ψ1, 2(m) as given by following equation 29 to a chirp signal generated by radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[1, 2] is output from transmit antenna Tx #[1, 2].
[29]
exp[jψ1,2(1)]cp(t),exp[jψ1,2(2)]cp(t),exp[jψ1,2(3)]cp(t), . . . ,exp[jψ1,2(Nc)]cp(t) (Equation 29)
Likewise, phase rotator PROT #[2, 2] applies, for each transmission period, phase rotation amount ψ2, 2(m) as given by following equation 30 to a chirp signal generated by radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[2, 2] is output from transmit antenna Tx #[2, 2].
[30]
exp[jψ2,2(1)]cp(t),exp[jψ2,2(2)]cp(t),exp[jψ2,2(3)]cp(t), . . . ,exp[jψ2,2(Nc)]cp(t) (Equation 30)
The foregoing description has been given of an example of how coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) is set.
In this embodiment, accordingly, plural transmit antennas 108 are associated with combinations (in other words, assignment) of Doppler shift amounts DOPndm and orthogonal code sequences Codencm such that in each of the combinations, at least one of Doppler shift amount DOPndm or orthogonal code sequence Codencm is different. In this embodiment, furthermore, the number of multiplexes of orthogonal code sequence Codencm (in other words, the number of coded Doppler multiplexes NDOP_CODE(ndm)) corresponding to each Doppler shift amount DOPndm in combinations of Doppler shift amounts DOPndm and orthogonal code sequences Codencm is different.
For example, in this embodiment, as illustrated in
[Configuration of Radar Receiver 200]
In
Each receive antenna 202 receives a reflected wave signal that is a radar transmission signal reflected from a target object (target), and outputs the received reflected wave signal to the corresponding one of antenna channel processors 201 as a reception signal.
Each antenna channel processor 201 includes receiving radio section 203 and signal processor 206.
Receiving radio section 203 includes mixer section 204 and low pass filter (LPF) 205. Receiving radio section 203 mixes, using mixer section 204, a chirp signal input from radar transmission signal generator 101, which is a transmission signal, with the received reflected wave signal, and passes the resulting mixed signal through LPF 205. As a result, a beat signal having a frequency corresponding to the delay time of the reflected wave signal is acquired. For example, as illustrated in
In each antenna channel processor 201-z (where z is equal to any of 1 to Na), signal processor 206 includes analog-to-digital (AD) converter 207, beat frequency analyzer 208, output switching section 209, and Doppler analyzers 210.
The signal (for example, beat signal) output from LPF 205 is passed to signal processor 206 and is converted into a discretely sampled data by AD converter 207.
Beat frequency analyzer 208 performs FFT processing of Ndata pieces of discretely sampled data, which are obtained in a defined time range (range gate), for each transmission period Tr. As a result, signal processor 206 outputs a frequency spectrum in which a peak appears at a beat frequency corresponding to the delay time of the reflected wave signal (radar reflected wave). In the FFT processing, for example, beat frequency analyzer 208 may perform multiplication by a window function coefficient such as of the Han window or the Hamming window. The use of a window function coefficient can suppress sidelobes around the beat frequency peak.
Here, a beat frequency response obtained by the m-th chirp pulse transmission, which is output from beat frequency analyzer 208 in the z-th signal processor 206, is represented by RETz(fb, m). Here, fb denotes the beat frequency index and corresponds to FFT index (bin number). For example, fb=0, . . . , Ndata/2, z=0, Na, and m=1, . . . , NC. A beat frequency having smaller beat frequency index fb indicates that the delay time of a reflected wave signal is shorter (in other words, the distance to the target object is shorter).
Beat frequency index fb can be converted into distance information R(fb) using following equation 31. For this reason, beat frequency index fb is hereinafter referred to as the “distance index fb”.
Here, Bw denotes a frequency-modulation bandwidth in a range gate for a chirp signal, and C0 denotes the speed of light.
Output switching section 209 selectively switches and outputs the output of beat frequency analyzer 208 for each transmission period to the OC_INDEX-th Doppler analyzer 210 among Loc Doppler analyzers 210, based on orthogonal code element index OC_INDEX output from coder 106 in phase rotation amount setter 104. In other words, output switching section 209 selects the OC_INDEX-th Doppler analyzer 210 in the m-th transmission period Tr.
Signal processor 206 includes Loc Doppler analyzers 210-1 to 210-Loc. For example, data is input to the noc-th Doppler analyzer 210 for every Loc transmission periods (Loc×Tr) by using output switching section 209. Accordingly, the noc-th Doppler analyzer 210 performs Doppler analysis for each distance index fb using data of Ncode transmission periods among Nc transmission periods (for example, beat frequency response RFTz(fb, m) output from beat frequency analyzer 208). Here, noc denotes the index of a code element, and noc=1, . . . , Loc.
For example, when Ncode is a power of 2, FFT processing is applicable in Doppler analysis. In this case, the FFT size is Ncode, and a maximum Doppler frequency without causing aliasing derived from the sampling theorem is ±1/(2Loc×Tr). Further, the Doppler frequency interval for Doppler frequency index fs is 1/(Ncode×Loc×Tr), and the range of Doppler frequency index fs is given by fs=−Ncode/2, . . . , 0, . . . , Ncode/2−1.
The following description will be given of a case where Ncode is a power of 2, as an example. When Ncode is not a power of 2, for example, data with zero padding can be used to perform FFT processing with a number of data sizes (FFT sizes) equal to a power of 2. In the FFT processing, Doppler analyzer 210 may perform multiplication by a window function coefficient such as of the Han window or the Hamming window. The application of a window function can suppress sidelobes around the beat frequency peak.
For example, output VFTznoc(fb, fs) of Doppler analyzers 210 in the z-th signal processor 206 is given by following equation 32, where j is the imaginary unit and z=1 to Na.
The foregoing description has been given of the processing performed by the components of signal processor 206.
In
For example, CFAR section 211 performs two-dimensional CFAR processing with the distance axis and the Doppler frequency axis (corresponding to the relative velocity) or CFAR processing that is a combination of one-dimensional CFAR processing operations by power addition of outputs VFTznoc(fb, fs) of Doppler analyzers 210 in the first to Na-th signal processors 206, for example, as given by following equation 33. As the two-dimensional CFAR processing or the CFAR processing that is a combination of one-dimensional CFAR processing operations, for example, processing disclosed in NPL 2 may be applied.
CFAR section 211 adaptively sets a threshold and outputs distance index fb_cfar that provides received power greater than the threshold, Doppler frequency index fs_cfar, and received-power information PowerFT(fb_cfar, fs_cfar) to coded Doppler demultiplexer 212.
For example, when phase rotation amount ϕndm for applying Doppler shift amount DOPndm is determined using equation 1, intervals ΔFD of Doppler shift amounts in the Doppler frequency domain, which are output from Doppler analyzers 210, are equal, where ΔFD=Ncode/NDM. Accordingly, in the outputs of Doppler analyzers 210, peaks are detected at intervals of ΔFD for Doppler-shift multiplexed signals in the Doppler frequency domain. When phase rotation amount ϕndm is determined using equation 1, ΔFD may not sometimes be an integer depending on Ncode and NDM. In this case, equation 50 described below may be used to achieve ΔFD having an integer value. The following describes a reception processing operation using ΔFD having an integer value.
Accordingly, CFAR section 211 may divide the respective outputs of Doppler analyzers 210 into ranges at Doppler shift amount intervals ΔFD and CFAR processing (referred to as, for example, “Doppler domain compression CFAR processing”) after, as given by following equation 34, providing power addition (referred to as, for example, “Doppler domain compression”) for the respective ranges while matching peak positions of Doppler-shift multiplexed signals. Here, fs_comp=−Ncode/2, . . . , −Ncode/2+ΔFD−1.
This can compress the Doppler frequency range for the CFAR processing to 1/NDM to reduce the amount of CFAR processing and can simplify the circuit configuration. In addition, CFAR section 211 enables power addition for NDM Doppler-shift multiplexed signals, resulting in signal to noise ratio (SNR) being improved by about (NDM)1/2. As a result, the radar sensing performance of radar apparatus 10 can be improved.
As a result of Doppler domain compression, the range of Doppler frequency index fs_comp in the Doppler frequency domain is reduced to greater than or equal to −Ncode/2, . . . , and less than or equal to −Ncode/2+ΔFD−1, and the range for the CFAR processing is compressed, leading to a reduction in the operation amount of CFAR processing. In
For example, CFAR section 211, which uses Doppler domain compression CFAR processing, adaptively sets a threshold and outputs distance index fb_cfar that provides received power greater than the threshold, Doppler frequency index fs_comp_cfar, and received-power information PowerFT(fb_cfar, fs_comp_cfar+(nfd−1)×ΔFD) for Doppler frequency indices (fs_comp_cfar+(nfd−1)×ΔFD) of NDM Doppler multiplexed signals, where nfd=1, . . . , NDM, to coded Doppler demultiplexer 212.
Note that phase rotation amount ϕndm for applying Doppler shift amount DOPndm is not limited to that given in equation 1. CFAR section 211 can apply Doppler domain compression CFAR processing, for example, if Doppler-shift multiplexed signals have phase rotation amounts ϕndm such that peaks are detected at constant intervals in the Doppler frequency domain output from Doppler analyzers 210.
Next, an example of the operation of coded Doppler demultiplexer 212 illustrated in
Coded Doppler demultiplexer 212 separates signals transmitted in a coded Doppler multiplexed manner, using the output of Doppler analyzers 210, based on the outputs of CFAR section 211, namely, distance index fb_cfar, Doppler frequency index fs_comp_cfar, and received power information PowerFT(fb_cfar, fs_comp_cfar+(nfd−1)×ΔFD) for Doppler frequency indices (fs_comp_cfar+(nfd−1)×ΔFD) of NDM Doppler multiplexed signals, where nfd=1, . . . , NDM, and performs discrimination (in other words, also referred to as determination or identification) of transmit antennas 108 and discrimination of Doppler frequencies (in other words, Doppler velocities or relative velocities).
As described above, for example, coder 106 in phase rotation amount setter 104 does not set all of NDM numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) to NCM, but sets at least one of the numbers of coded Doppler multiplexes to a value smaller than NCM. For example, coded Doppler demultiplexer 212 performs (1) code separation process to detect a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM and performs (2) discrimination of transmit antennas 108 and discrimination of Doppler frequencies of the target based on the detected coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM.
The following describes the processes (1) and (2) described above, which are performed by coded Doppler demultiplexer 212.
<(1) Code Separation Process (Process for Detecting Coded Doppler Multiplexed Signal for which Number of Coded Doppler Multiplexes is Set to be Smaller than NCM)>
There are NDM candidate correspondences between NDM coded Doppler multiplexed signals and the outputs of the respective Doppler analyzers 210 for Doppler frequency indices (fs_comp_cfar), (fs_comp_cfar+ΔFD), (fs_comp_cfar+2ΔFD), . . . , and (fs_comp_cfar+(NDM−1)ΔFD) of NDM coded Doppler multiplexed signals for distance index fb_cfar output from CFAR section 211.
For example, if Doppler shift amount DOPndm set in Doppler shift setter 105 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, there are NDM candidate correspondences with cyclically shifted elements, as follows, with consideration given to Doppler aliasing. Here, patterns of the candidate correspondences are numbered DopCase=1 to NDM.
DopCase=1: {DOP1, DOP2, . . . , DOPDM-1, DOPDM}
DopCase=2: {DOPDM, DOP1, DOP2, . . . , DOPDM-1}
. . . ,
DopCase=NDM: {DOP2, . . . , DOPDM-1, DOPDM, DOP1}
For example, DopCase=1 indicates a correspondence among Doppler shift amounts in the initial state (when the relative velocity to the target is zero). For example, more aliasing components are included as the relative velocity of the target increases in a direction in which the distance to the target decreases, and the resulting correspondences are associated with DopCase=2, 3, . . . , NDM. In other words, DopCase=NDM, NDM-1, . . . , 2 is associated as the relative velocity of the target increases in a direction in which the distance to the target increases.
Here, a table indicating the position of DOPndm counting from the beginning of each DopCase (the position (or order) of DOPndm in DopCase) can be created in advance, based on the Doppler shift amounts set in Doppler shift setter 105. In the following, DOPposi(DOPndm, DopCase) denotes the operator that outputs the position of DOPndm counting from the beginning of each DopCase. For example, in the above-described example of DopCase, DOPposi(DOP1, 1)=1, DOPposi(DOP1, 2)=2, DOPposi(DOP1, NDM)=NDM, DOPposi(DOP2, 1)=2, DOPposi(DOP2, 2)=3, and DOPposi(DOP2, NDM)=1.
In Doppler analyzers 210, the range of Doppler frequencies without causing aliasing is greater than or equal to −1/(2×Loc×Tr) and less than 1/(2×Loc×Tr), and Doppler aliasing occurs outside this range. For example, in a case where Loc=2, in the outputs of Doppler analyzers 210, the range of Doppler frequencies without causing aliasing is greater than or equal to 1/(4×Tr) and less than 1/(4×Tr).
In
DopCase=1: when −1/(2×Tr)≤fTARGET<−3/(8×Tr), DOP1<DOP2<DOP3<DOP4
DopCase=2: when −3/(8×Tr)≤fTARGET<−1/(4×Tr), DOP4<DOP1<DOP2<DOP3
DopCase=3: when −1/(4×Tr)≤fTARGET<−1/(8×Tr), DOP3<DOP4<DOP1<DOP2
DopCase=4: when −1/(8×Tr)≤fTARGET<0, DOP2<DOP3<DOP4<DOP1
DopCase=1: when 0≤fTARGET<1/(8×Tr), DOP1<DOP2<DOP3<DOP4
DopCase=2: when 1/(8×Tr)≤fTARGET<1/(4×Tr), DOP4<DOP1<DOP2<DOP3
DopCase=3: when 1/(4×Tr)≤fTARGET<3/(8×Tr), DOP3<DOP4<DOP1<DOP2
DopCase=4: when 3/(8×Tr)≤fTARGET<1/(2×Tr), DOP2<DOP3<DOP4<DOP1
Here, for Doppler frequency fTARGET of the target in the range of −1/(2×Tr)≤fTARGET<1/(2×Tr), the candidate patterns of the order of coded Doppler multiplexed signals are numbered DopCase=1 to 4 (=NDM). There are 4 (=NDM) candidate patterns to be associated with the order of coded Doppler multiplexed signals.
In the outputs of Doppler analyzers 210, furthermore, the range of Doppler frequencies without causing aliasing is greater than or equal to −1/(2×Loc×Tr) and less than 1/(2×Loc×Tr). Thus, for Doppler frequency fTARGET of the target in the range of −1/(2×Tr)≤fTARGET<1/(2×Tr), the outputs of Doppler analyzers 210 include (Loc−1) occurrences of aliasing. Therefore, for example, in the example (Loc=2) illustrated in
Coded Doppler demultiplexer 212 performs a code separation process on, for example, the outputs of Doppler analyzers 210 in the z-th signal processor 206 indicated by Doppler frequency indices (fs_comp_cfar±(nfd−1)×ΔFD) of NDM coded Doppler multiplexed signals for distance index fb_cfar output from CFAR section 211.
The code separation process may be performed on, for example, all of nfd=1, . . . , NDM for all the candidates of DopCase=1, . . . , NDM. Note that coded Doppler demultiplexer 212 detects a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM, and performs discrimination of transmit antennas 108 and determination of a target Doppler frequency. Accordingly, for example, coded Doppler demultiplexer 212 performs a code separation process, as given by following equation 35, to reduce the operation amount of the separation process.
The superscript asterisk (*) indicates the complex conjugate operator. Further, nfd=1, . . . , NDM, ncm=1, . . . , NCM, and DopCase=1, . . . , NDM.
For example, when a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM is Doppler shift amount “DOPndm1”, coded Doppler demultiplexer 212 performs a code separation process using candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1 in each DopCase.
Here, code separation signal DeMULzncm(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) obtained in equation 35 represents a code separation signal that uses the ncm-th orthogonal code sequence Codencm for the output of Doppler analyzer 210 in the z-th signal processor 206 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD).
In the exp term of equation 35, since the sampling times for the outputs of Doppler analyzers 210 for each code element OCncm(noc) are shifted, phase correction is performed in accordance with Doppler frequency index (fs_comp_cfar±(DopCase−1))×ΔFD-DOP1).
Further, Doppler shift amount DOPndm set in Doppler shift setter 105 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, and DOP1 falls within the range of fs_comp=−Ncode/2, −Ncode/2+ΔFD−1 in the initial state (when the relative velocity to the target is zero). Accordingly, for example, coded Doppler demultiplexer 212 calculates an amount of phase correction using DOP1 as a reference.
When a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM is Doppler shift amount DOPndm1, signals separated based on NCM orthogonal codes are obtained by equation 35 for candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1, where DopCase=1, . . . , NDM. Accordingly, coded Doppler demultiplexer 212 obtains the outputs of a total of NDM×NCM code separation signals. For example, coded Doppler demultiplexer 212 calculates code separation signals for all receive antennas z=1, Na in accordance with equation 35 and calculates code separation signal power sum Pow_DeMULncm(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) using following equation 36.
Here, a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM does not use all the NCM orthogonal codes. In other words, a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM uses some of the orthogonal codes (for example, one orthogonal code). Thus, code separation signal power sum Pow_DeMULncm(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) corresponding to a coded Doppler multiplexed signal (Doppler shift amount: DOPndm1) for which the number of coded Doppler multiplexes is set to be smaller than NCM, where ncm=1, . . . , NCM and DopCase=1, . . . , NDM, contains a component whose received power has as low power value as about the noise level.
As an example, a description will be given of a case where the number of coded Doppler multiplexes of only a coded Doppler multiplexed signal that uses DOPndm1 is set to (NCM−1), which is smaller than NCM (see, for example, the examples of the assignment of Doppler shift amounts and orthogonal codes illustrated in
In this case, coded Doppler demultiplexer 212 detects, among Pow_DeMULncm(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD), where ncm=1, . . . , NCM and DopCase=1, . . . , NDM, Pow_DeMULncm_min(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase_min)−1)×ΔFD) with the minimum received power. Here, “ncm_min” and “DopCase_min” represent the index numbers of ncm and DopCase for which Pow_DeMULncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) has the minimum received power.
As another example, a description will be given of a case where the number of coded Doppler multiplexes of only a coded Doppler multiplexed signal that uses DOPndm1 is set to ((NCM−2), which is smaller than NCM (see, for example, the examples of the assignment of Doppler shift amounts and orthogonal codes illustrated in
In this case, among Pow_DeMULncm(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD), where ncm=1, . . . , NCM and DopCase=1, . . . , NDM, a certain DopCase does not use two codes (for example, the orthogonal codes indicated by crosses (“x”) in
Accordingly, coded Doppler demultiplexer 212 detects Pow_DeMULncm_min1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase_min)−1)×ΔFD)+Pow_DeMULncm_min2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) for which, for a combination of two different orthogonal codes Codencm1 and Codencm2 (where ncm1≠ncm2) among the NCM orthogonal codes, power sum Pow_DeMULncm1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD)+Pow_DeMULncm2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) has the minimum received power. Here, “ncm_min1”, “ncm_min2”, and “DopCase_min” represent the index numbers of ncm1, ncm2, and DopCase for which Pow_DeMULncm1 (fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD)+Pow_DeMULncm2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) has the minimum received power.
Also when the number of coded Doppler multiplexes is set to be smaller than (NCM−2), coded Doppler demultiplexer 212 detects a combination of candidate codes with the minimum code separation signal power sum among combinations of candidate codes that are not used among the NCM orthogonal codes. Accordingly, a combination of DopCase_min and codes (ncm_min1, ncm_min2, . . . ) can be detected.
Coded Doppler demultiplexer 212 determines, based on the detection result described above, that a correspondence pattern of NDM coded Doppler multiplexed signals and the outputs of the respective Doppler analyzers 210 for Doppler frequency indices (fs_comp_cfar), (fs_comp_cfar+ΔFD), (fs_comp_cfar+2ΔFD), . . . , and (fs_comp_cfar±(NDM−1)ΔFD) of the NDM coded Doppler multiplexed signals for distance index fb_cfar output from CFAR section 211 is DopCase_min among DopCase=1 to NDM.
Further, coded Doppler demultiplexer 212 determines that a coded Doppler multiplexed signal transmitted using DOPndm1 for which the number of coded Doppler multiplexes is set to be smaller than NCM in the pattern corresponding to DopCase_min is the output of Doppler analyzer 210 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar±(DOPposi(DOPndm1, DopCase_min)−1)×ΔFD).
A reduction in the SNR of a reception signal that is a reflected wave may make it difficult to discriminate between a signal power level and a noise power level. To address this difficulty, coded Doppler demultiplexer 212 may introduce a determination condition and perform a process of, for example, adopting a determination result (in other words, detection result) when the determination condition is satisfied and removing (in other words, not adopting) a determination result when the determination condition is not satisfied. This can reduce the probability of erroneous detection of the noise component and the like. For example, coded Doppler demultiplexer 212 may adopt, as a determination result, detection value PMIN with the minimum received power when satisfying PowerFT(fb_cfar, fs_comp_cfa)>LEVDETECT×PMIN. Here, LEVDETECT is a determination threshold. LEVDETECT is a real number satisfying 0<LEVDETECT<1.
<(2) Process for Discrimination of Transmit Antennas 108 and Discrimination of Doppler Frequencies of Target>
For example, coded Doppler demultiplexer 212 detects Doppler frequency fTARGET of the target in the range of −1/(2×Tr)≤fTARGET<1/(2×Tr).
For example, coded Doppler demultiplexer 212 determines the presence or absence of aliasing as in Case (a) and Case (b) below, based on the determination result of DopCase (for example, DopCase_min) and the determination result of the code whose received power is minimum (for example, ncm_min). Then, coded Doppler demultiplexer 212 discriminates transmit antennas 108 and Doppler frequencies of the target.
A plurality of codes whose received power is minimum (for example, ncm_min1, ncm_min2, . . . ) are also applicable by replacing ncm_min in the following description with (ncm_min1, ncm_min2, . . . ).
[Case (a): Case without Aliasing]
For example, when ncm_min matches the index of an orthogonal code sequence that is not used (not assigned) for coded Doppler multiplexing in DOPndm1 for which the number of coded Doppler multiplexes is set to be smaller than NCM, coded Doppler demultiplexer 212 determines that no Doppler aliasing occurs. In other words, when a coded Doppler multiplexed signal corresponding to DOPndm1 is not coded using orthogonal code sequence Codencm_min, coded Doppler demultiplexer 212 determines that no Doppler aliasing occurs.
If it is determined that no Doppler aliasing occurs, Coded Doppler demultiplexer 212 determines the Doppler frequency of the target and a transmit antenna in the following way.
Doppler Frequency Determination for Target:
Coded Doppler demultiplexer 212 determines that the Doppler frequency index of the target is (fs_comp_cfar(DopCase_min−1)×ΔFD−DOP1). For example, the Doppler frequency interval for Doppler frequency index fs_comp_cfar is 1/(Ncode×Loc×Tr). Therefore, coded Doppler demultiplexer 212 determines that Doppler frequency fTARGET of the target is (fs_comp_cfar (DopCase_min−1)×ΔFD−DOP1)/(Ncode×Loc×Tr).
Transmit Antenna Determination:
Coded Doppler demultiplexer 212 determines that for the output of Doppler analyzer 210 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) in the z-th signal processor 206, code separation signal DeMULzncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Codencm is the reception signal corresponding to transmit antenna Tx #[ncm, ndm].
[Case (b): Case with Aliasing]
For example, when ncm_min does not match the index of an orthogonal code sequence that is not used for coded Doppler multiplexing in DOPndm1 for which the number of coded Doppler multiplexes is set to be smaller than NCM, coded Doppler demultiplexer 212 determines that Doppler aliasing occurs. In other words, when a coded Doppler multiplexed signal corresponding to DOPndm1 includes a multiplexed signal coded using orthogonal code sequence Codencm_min, coded Doppler demultiplexer 212 determines that the output of Doppler analyzer 210 for distance index fb_cfar, Doppler frequency index (fs_comp_cfar±(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) in the z-th signal processor 206 exceeds the range of maximum Doppler frequencies without causing aliasing derived from the sampling theorem.
If it is determined that Doppler aliasing occurs, coded Doppler demultiplexer 212 decides (in other words, fixes) a Doppler frequency index of the target and a transmit antenna.
Doppler Frequency Determination for Target:
Coded Doppler demultiplexer 212 determines that the Doppler frequency index of the target is (fs_comp_cfar(DopCase_min−1)×ΔFD)−DOP1−Ncode×Sign(fs_comp_cfar+(DopCase_min−1)×ΔFD−DOP1)). For example, the Doppler frequency interval for Doppler frequency index fs_comp_cfar is 1/(Ncode×Loc×Tr). Therefore, coded Doppler demultiplexer 212 determines that Doppler frequency fTARGET of the target is (fs_comp_cfar+(DopCase_min−1)×ΔFD−DOP1−Ncode×Sign(fs_comp_cfar+(DopCase_min−1)×ΔFD−DOP1))/(Ncode×Loc×Tr). Sign(x) is a sign function and is a function for real numbers x, providing as output 1 when x>0, 0 when x=0, and −1 when x<0.
In this manner, when the outputs of Doppler analyzers 210 contain Doppler aliasing, coded Doppler demultiplexer 212 determines a Doppler frequency of the target, with consideration given to the aliasing component (for example, Ncode×Sign(fs_comp_cfar+(DopCase_min−1)×ΔFD−DOP1)).
Transmit Antenna Determination:
When the outputs of Doppler analyzers 210 contain Doppler aliasing, incorrect phase correction is performed by phase correction used in the code separation process (for example, the exp term in equation 35). This equivalently indicates that a code separation process is performed in a case where the elements of an orthogonal code sequence used for the code separation process are code elements given by following equations 37 and 38.
For example, in a case where NCM=2 and orthogonal code sequence Code1={1, 1} with code length Loc=2 is used, a process equivalent to a separation process using Code2={1, −1}, as given by following equation 39, is performed, regardless of Sign(fest).
In a case where Code2={1, −1} is used, on the other hand, a process equivalent to a separation process using Code1={1, 1}, as given by following equation 40, is performed, regardless of Sign(fest).
Accordingly, in a case where NCM=2 and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, coded Doppler demultiplexer 212 determines that for the output of Doppler analyzer 210 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) in the z-th signal processor 206, code separation signal DeMULz1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code1 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[2, ndm] and code separation signal DeMULz2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code2 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[1, ndm].
For example, in a case where NCM=3 and orthogonal code sequences Code1={1, 1, 1}, Code2={1, exp(j2π/3), exp(j4π/3)}, and Code3={1, exp(−j2π/3), exp(−j4π/3)} with code length Loc=3 (=NCM) are used, determination is made as following.
For example, in a case where Code1={1, 1, 1} is used, as given by following equation 41, a process equivalent to a separation process using Code3 or Code2 in accordance with positive or negative Sign(fest) is performed.
For example, in a case where Code2={1, exp(j2π/3), exp(j4π/3)} is used, as given by following equation 42, a process equivalent to a separation process using Code1 or Code3 in accordance with positive or negative Sign(fest) is performed.
For example, in a case where Code3={1, exp(−j2π/3), exp(−j4π/3)} is used, as given by following equation 43, a process equivalent to a separation process using Code2 or Code1 in accordance with positive or negative Sign(fest) is performed.
Accordingly, in a case where NCM=3 and orthogonal code sequences Code1={1, 1, 1}, Code2={1, exp(j2π/3), exp(j4π/3)}, and Code3={1, exp(−j2π/3), exp(−j4π/3)} with code length Loc=3 (=NCM) are used, coded Doppler demultiplexer 212 determines transmit antenna 108 in the following way for the output of Doppler analyzer 210 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) in the z-th signal processor 206.
Case where Sign(fest)>0:
Coded Doppler demultiplexer 212 determines that code separation signal DeMULz1(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code1 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[3, ndm], that code separation signal DeMULz2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code2 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[1, ndm], and that code separation signal DeMULz3(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code3 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[2, ndm].
Case where Sign(fest)<0:
Coded Doppler demultiplexer 212 determines that code separation signal DeMULz1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code1 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[2, ndm], that code separation signal DeMULz2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code2 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[3, ndm], and that code separation signal DeMULz3(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Code3 is the reception signal of the reflected wave transmitted from transmit antenna Tx #[1, ndm].
As described above, the correspondence between a code used for coding during transmission and a code separated by the code separation process during aliasing determination is decided in advance based on orthogonal code sequence Codencm (ncm=1, . . . , NCM) and Sign(fest). This enables coded Doppler demultiplexer 212 to use, for example, code conversion function AliasConv[ncm, Sign(fest)]=Tx_ncm (that outputs orthogonal code sequence CodeTx_ncm index Tx_ncm separated in the code separation process using Codencm during aliasing determination).
Further, inverse code conversion function AliasConv−1[Tx_ncm, Sign(fest)]=ncm (that outputs index ncm of code sequence Codencm used in the code separation process, where the index of an orthogonal code sequence separated in the code separation process during aliasing determination is Tx_ncm) can also be defined in a similar manner.
For example, in a case where NCM=2 and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, code conversion function AliasConv[ncm, Sign(fest)] is decided in advance such that AliasConv[1, Sign(fest)]=2 and AliasConv[2, Sign(fest]=1. Further, the inverse code conversion function is decided in advance such that AliasConv−1[1, Sign(fest)]=2 and AliasConv−1[2, Sign(fest)]=1.
For example, in a case where NCM=3 and orthogonal code sequences Code1={1, 1, 1}, Code2={1, exp(j2π/3), exp(j4π/3)}, and Code3={1, exp(−j2π/3), exp (−j4π/3)} with code length Loc=3 (=NCM) are used, code conversion function AliasConv[ncm, Sign(fest)] is decided in advance in the following way.
Case where Sign(fest)>0:
Code conversion function AliasConv[ncm, Sign(fest)] is decided in advance such that AliasConv[1, Sign(fest)]=3, AliasConv[2, Sign(fest)]=1, and AliasConv[3, Sign(fest)]=2. Further, the inverse code conversion function is decided in advance such that AliasConv−1[1, Sign(fest)]=2, AliasConv−1[2, Sign(fest)]=3, and AliasConv−1[3, Sign(fest)]=1.
Case where Sign(fest)<0:
Code conversion function AliasConv[ncm, Sign(fest)] is decided in advance such that AliasConv[1, Sign(fest)]=2, AliasConv[2, Sign(fest)]=3, and AliasConv[3, Sign(fest)]=1. Further, the inverse code conversion function is decided in advance such that AliasConv−1[1, Sign(fest)]=3, AliasConv−1[2, Sign(fest)]=1, and AliasConv−1[3, Sign(fest)]=2.
Accordingly, for example, with the use of code conversion function AliasConv[ncm, Sign(fest)], coded Doppler demultiplexer 212 performs transmit antenna determination during aliasing determination in the following way. For example, coded Doppler demultiplexer 212 determines that for the output of Doppler analyzer 210 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) in the z-th signal processor 206, code separation signal DeMULzncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm, DopCase_min)−1)×ΔFD) using Codencm is the reception signal corresponding to transmit antenna Tx #[AliasConv[ncm, Sign(fest)], ndm].
Alternatively, coded Doppler demultiplexer 212 may determine that the reception signal corresponding to transmit antenna Tx #[Tx_ncm, ndm] is a code separation signal given by following equation 44 that uses Codencm (where ncm=AliasConv−1[Tx_ncm, Sign(fest)]).
[44]
DeMULzAliasConv
Also in a case where NCM>3, the code conversion function and the inverse code conversion function can be determined in advance, and the determination of a transmit antenna can be performed in a similar manner.
In this embodiment, accordingly, radar receiver 200 can perform aliasing determination (for example, determination of the presence or absence of aliasing and determination of correspondence DopCase) based on the result of code separation of a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM. This enables radar apparatus 10 to discriminate transmit antennas 108 corresponding to coded Doppler multiplexed signals and to discriminate Doppler frequencies of the target even if aliasing occurs. According to this embodiment, therefore, for example, the range over which a Doppler frequency is detectable without ambiguity can be extended to the range greater than or equal to −1/(2Tr) and less than 1/(2Tr).
For example, the range over which a Doppler frequency is detectable without ambiguity when one transmit antenna 108 is used for transmission is the range greater than or equal to −1/(2Tr) and less than 1/(2Tr). In this embodiment, therefore, even when plural transmit antennas 108 are used, the range over which a Doppler frequency is detectable without ambiguity can be achieved in a manner similar to that when a single antenna is used for transmission.
Next, processing for a plurality of coded Doppler multiplexed signals for which the numbers of coded Doppler multiplexes are set to be smaller than NCM (see, for example,
For example, a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM is not only Doppler shift amount DOPnd but there is a plurality of coded Doppler multiplexed signals for which the numbers of coded Doppler multiplexes are set to be smaller than NCM. In this case, coded Doppler demultiplexer 212 detects ncm_min1, ncm_min2, and DopCase_min in the following way.
For example, when coded Doppler multiplexed signals for which the numbers of coded Doppler multiplexes are set to (NCM−1), which is smaller than NCM, are Doppler shift amounts DOPndm1 and DOPndm2, coded DOP Doppler demultiplexer 212 calculates, based on equation 36, code separation signal power sums Pow_DeMULncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) and Pow_DeMULncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm2, DopCase)−1)×ΔFD). Here, ncm=1, . . . , NCM, and DopCase=1, . . . , NDM.
Here, in a coded Doppler multiplexed signal that uses Doppler shift amount DOPndm1, an orthogonal code that is not used for coded Doppler multiplexing is represented by Codencm1, and in a coded Doppler multiplexed signal that uses Doppler shift amount DOPndm2, an orthogonal code that is not used for coded Doppler multiplexing is represented by Codencm2.
In this case, in code separation signal power sum Pow_DeMULncm1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD)+Pow_DeMULncm2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm2, DopCase)−1)×ΔFD) or code separation signal power sum Pow_DeMULAliasConv[ncm1, Sign(fest)](fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD)+Pow_DeMULAliasConv[ncm2, Sign(fest)](fb_cfar, fs_comp_cfar±(DOPposi(DOPndm2, DopCase)−1)×ΔFD) during Doppler aliasing determination, coded Doppler demultiplexer 212 detects Pow_DeMULncm_min1(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase_min)−1)×ΔFD)+Pow_DeMULncm_min2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm2, DopCase_min)−1)×ΔFD) with the minimum received power when DopCase is equal to any of 1, . . . , and NDM.
Here, “ncm_min1”, “ncm_min2”, and “DopCase_min” represent the index numbers of ncm1, ncm2, and DopCase for which code separation signal power sum Pow_DeMULncm1(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD)+Pow_DeMULncm2(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm2, DopCase)−1)×ΔFD) has the minimum received power.
The number of coded Doppler multiplexed signals for which the numbers of coded Doppler multiplexes are set to be smaller than NCM is not limited to two and may be three or more.
The foregoing description has been given of an example of the operation of coded Doppler demultiplexer 212.
In
For example, direction estimator 213 generates, based on the output of coded Doppler demultiplexer 212, virtual receive array correlation vector h(fb_cfar, fs_comp_cfar) given by following equation 45 or 46 and performs a direction estimation process.
Virtual receive array correlation vector h(fb_cfar, fs_comp_cfar) includes Nt×Na elements, the number of which is equal to the product of the number Nt of transmit antennas and the number Na of receive antennas. Virtual receive array correlation vector h(fb_cfar, fs_comp_cfar) is used in a process for performing direction estimation on reflected wave signals from the target, the direction estimation being based on a phase difference among receive antennas 202. Here, z=1, . . . , Na.
Case without Aliasing:
Case with Aliasing:
For example, direction estimator 213 calculates a spatial profile, with azimuth direction θ in direction estimation evaluation function value PH(θ, fb_cfar, fs_comp_cfar) being variable within a defined angular range. Direction estimator 213 extracts a predetermined number of local maximum peaks in the calculated spatial profile in order from the largest and outputs the azimuth directions of the local maximum peaks as direction-of-arrival estimation values (for example, position measurement output).
There are various methods with direction estimation evaluation function value PH(θ, fb_cfar, fs_comp_cfar) depending on the direction-of-arrival estimation algorithm. For example, an estimation method using an array antenna, as disclosed in NPL 3, may be used.
For example, when a virtual receive array includes Nt×Na elements arranged in a straight line at equal intervals dH, a beamformer method can be given by following equations 47 and 48. Any other technique, such as Capon or MUSIC, may be applicable.
In equation 47, the superscript H denotes the Hermitian transpose operator. Further, a(θu) denotes the direction vector of the virtual receive array relative to an incoming wave in azimuth direction θu.
Azimuth direction θu is a vector that is changed at azimuth interval β1 in an azimuth range over which direction-of-arrival estimation is performed. For example, θu is set as follows.
θu=θmin+uβ1,u=0, . . . ,NU
NU=floor[(θmax−θmin)/β1]+1
Here, floor(x) is a function that returns the largest integer value not greater than real number x.
In equation 47, furthermore, Dcal is an (Nt×Na) order matrix including an array correction coefficient for correcting phase deviations and amplitude deviations across the transmit array antenna and across the receive array antenna, and a coefficient for reducing the influence of coupling of elements across the antennas. If the coupling between antennas in the virtual receive array is negligible, Dcal is a diagonal matrix with diagonal components including an array correction coefficient for correcting phase deviations and amplitude deviations across the transmit array antenna and across the receive array antenna.
For example, direction estimator 213 may output, as position measurement results, distance information based on distance index fb_cfar and Doppler velocity information of the target based on the Doppler frequency determination result for the target, together with the direction estimation result, to, for example, a vehicle control apparatus (not illustrated) in the case of an in-vehicle radar or to an infrastructure control apparatus (not illustrated) in the case of an infrastructure radar.
Doppler frequency information may be converted into the relative velocity component and then output. Doppler frequency index four in the Doppler frequency determination result for the target can be converted into relative velocity component vd(fout) using following equation 49. Here, λ denotes the wavelength of the carrier frequency of an RF signal output from a transmission radio section (not illustrated). Further, Δf denotes the Doppler frequency interval in FFT processing performed in Doppler analyzer 210. For example, in this embodiment, Δf=1/{Ncode×Loc×Tr}.
As described above, in this embodiment, radar apparatus 10 applies phase rotation amounts corresponding to Doppler shift amounts and orthogonal code sequences to radar transmission signals to transmit radar transmission signals (in other words, coded Doppler multiplexed signals) from plural transmit antennas 108 in a multiplexed manner. In this embodiment, plural transmit antennas 108 are associated with combinations of Doppler shift amounts (DOPndm) and orthogonal code sequences (DOPncm) such that in each of the combinations, at least one of Doppler shift amount (DOPndm) or orthogonal code sequence (DOPncm) is different. In this embodiment, furthermore, the number of multiplexes of the orthogonal code sequence (in other words, the number of codes) corresponding to each Doppler shift amount in combinations of Doppler shift amounts and orthogonal code sequences is different. In other words, the numbers of coded Doppler multiplexes for the respective Doppler-multiplexed transmission signals are set to be non-uniform.
Radar apparatus 10 can determine, based on, for example, the received power of a code-separated signal for each coded Doppler multiplexed signal, transmit antenna 108 associated with the coded Doppler multiplexed signal (in other words, the combination of Doppler shift amount and orthogonal code sequence) and the presence or absence of Doppler aliasing (for example, DopCase and the like). This enables radar apparatus 10 to appropriately determine a Doppler frequency of the target even in the presence of Doppler aliasing.
According to this embodiment, therefore, radar apparatus 10 can extend the effective Doppler frequency bandwidth to 1/(Tr) and can extend the Doppler frequency (relative velocity) detection range with no ambiguity. Accordingly, radar apparatus 10 can improve target-object sensing accuracy over a wider Doppler frequency range.
In this embodiment, furthermore, coded Doppler multiplexing, which is performed using both Doppler multiplexing and coding, can reduce the number of Doppler multiplexes compared with the use of only Doppler multiplexing in multiplexing transmission. This can increase the intervals of phase rotation amounts for applying Doppler shifts, thereby, for example, relieving the accuracy requirements (phase modulation accuracy) for the phase shifters and achieving the cost reduction effect of an RF section, including reduction of the man-hours required for the adjustment of the phase shifters.
In this embodiment, furthermore, since coded Doppler multiplexing is performed using both Doppler multiplexing and coding, radar apparatus 10 performs, for each code element, Fourier frequency analysis (FFT processing) to detect the Doppler frequency (detect the relative velocity). Accordingly, for example, compared with Fourier frequency analysis (FFT processing) to detect the Doppler frequency (detect the relative velocity) using only Doppler multiplexing in multiplexing transmission, the FFT size corresponds to (1/code length) and the number of times FFT processing is performed is increased by (code length) times. For example, when the FFT operation amount with FFT size Nc is roughly estimated to be Nc×log2(Nc), coded Doppler multiplexing according to this embodiment has an operation amount ratio of about {Loc×Nc/Loc×log2(Nc/Loc)}/{Nc×log2(Nc)}=1−log2(Loc)/log2(Nc) relative to FFT operation with only Doppler multiplexing. For example, in a case where Loc=2 and Nc=1024, the operation amount ratio is 0.9. The operation reduction effect of FFT processing can be achieved, and the effect of simplification of the circuit configuration and cost reduction can also be achieved.
(Variation 1 of Embodiment 1)
Phase rotation amount ϕndm for applying Doppler shift amount DOPndm is not limited to, for example, the value given in equation 1 and so on. For example, phase rotation amount ϕndm may be a value given by following equation 50. Here, round(x) represents the round function that outputs a rounded integer value for real number x. The round(Ncode/NDM) term is introduced so that the phase rotation amount is an integer multiple of the Doppler frequency interval in Doppler analyzer 210. In equation 50, the angle is expressed in radian.
(Variation 2 of Embodiment 1)
Embodiment 1 has described a case where in coder 106 in phase rotation amount setter 104, the numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) are set to be different (in other words, non-uniform) in a range greater than or equal to 1 and less than or equal to NCM; however, this is not required. For example, at least one of the numbers of coded Doppler multiplexes NDOP_CODE(1), NDOP_CODE(2), . . . , and NDOP_CODE(NDM) may be greater than or equal to 1 and less than NCM, and the number of coded Doppler multiplexes, whose value is 0, may be included.
For example, coded Doppler demultiplexer 212 detects, for combinations of candidate codes that are not used among NCM orthogonal codes, a combination with the minimum code separation signal power sum. Accordingly, combinations of DopCase_min and codes (ncm_min1, ncm_min2, . . . ) can be detected. Further, coded Doppler demultiplexer 212 performs Doppler aliasing determination by utilizing the inclusion of at least one DOPndm1 for which the number of coded Doppler multiplexes is set to be greater than or equal to 1 and less than NCM.
This enables radar apparatus 10 to detect Doppler frequency fTARGET of the target in the range of −1/(2×Tr)≤fTARGET<1/(2×Tr) and to perform discrimination of transmit antennas 108 and discrimination of Doppler frequencies of the target.
For example,
(Variation 3 of Embodiment 1)
Variation 3 describes the sub-array configuration of transmit antennas of a radar apparatus.
A combination of some of the transmit antennas may be used as a sub-array to narrow the beam width of a transmission directional beam pattern to improve transmission directional gain. This narrows the sensible angular range, but increases the sensible distance range. In addition, a beam weight factor for generating a directional beam can be made variable to control the beam direction to be variable.
In
In
Here, Index_TxSubArray represents the element index of a sub-array, and Index_TxSubArray=1, . . . , NSA. Further, λ denotes the wavelength of a radar transmission signal, and dSA denotes sub-array antenna spacing.
Each beam weight multiplier 110 multiplies the output of the corresponding one of phase rotators 107 by beam weight factor WTx(Index_TxSubArray, θTxBF) input from beam weight generator 109. Transmission signals multiplied by beam weight WTx(Index_TxSubArray, θTxBF) are transmitted from NSA sub-array antennas. Here, Index_TxSubArray=1, . . . , NSA.
With the operation described above, radar transmitter 100a can perform transmission such that transmission directional beams can be directed in a predetermined direction using sub-arrays for the outputs of phase rotators 107. This can improve transmission directional gain in a predetermined direction and extend the sensible distance range. Additionally, the SNR can be improved.
Further, radar transmitter 100a can variably set a beam weight factor for generating a transmission directional beam to control the beam direction to be variable.
The configuration for sub-array transmission described in Variation 3 is also applicable to other variations or embodiments.
In the configuration illustrated in
(Variation 4 of Embodiment 1)
Variation 4 describes a case where the Doppler shift amount (or transmission period Tr) varies for each transmission frame. In other words, the intervals of Doppler shift amounts used for coded Doppler multiplexing transmission are each set to vary for each frame (for example, every Nc transmission periods (Nc×Tr)) in which a radar transmission signal is transmitted.
For example, there is a plurality of targets having substantially equal reception levels for the same distance index fb_cfar. In this case, if the intervals of Doppler peaks of the plurality of targets match the intervals of Doppler shift amounts in Doppler multiplexing, coded Doppler demultiplexer 212 may be unable to discriminate transmit antennas 108 and the Doppler frequencies of the targets.
Since the Doppler frequencies of a plurality of targets may differ, the relative movement speeds between the targets and radar apparatus 10 may differ. Radar apparatus 10 continuously performs radar observation such that even if transmit antennas 108 and the Doppler frequencies of the targets fail to be discriminated in radar position measurement output at a certain point in time, the distance between the plurality of targets is different in radar position measurement output at the subsequent point in time, resulting in it being likely that measurement is successful. It is therefore assumed that separate output of each of the plurality of targets is likely to be obtained.
In addition, to more reliably separate signals corresponding to the plurality of targets, for example, radar apparatus 10 may continuously perform radar position measurement by setting at least one of transmission period Tr and the Doppler shift amount to vary for each radar position measurement (for example, every Nc transmission periods (Nc×Tr)).
Accordingly, for example, even if, for the same distance index fb_cfar, the reception levels of Doppler peaks of a plurality of targets are substantially equal, the intervals of the Doppler peaks match the intervals of Doppler shift amounts, and coded Doppler demultiplexer 212 is unable to discriminate transmit antennas 108 and to discriminate the Doppler frequencies of the targets, the intervals of the Doppler shift amounts may be more likely to be different in the subsequent radar position measurement. This enables radar apparatus 10 to more reliably separate the signals corresponding to the plurality of targets.
For example, phase rotation amount ϕndm for applying Doppler shift amount DOPndm may be set as given by following equation 52 (the angle is expressed in radian).
For example, in equation 52, radar apparatus 10 varies δ for each radar position measurement, thereby variably setting the intervals of the Doppler shift amounts. For example, radar apparatus 10 may periodically vary δ for each radar position measurement to 0, 1, 0, 1, . . . .
Further, for example, radar apparatus 10 sets transmission period Tr to vary for each radar position measurement, thereby changing the intervals of the Doppler shift amounts. As a result, an effect equivalent to that when the Doppler shift amounts are variably set is achieved.
(Variation 5 of Embodiment 1)
Variation 5 describes, for example, a method for reducing the influence of interference from a plurality of radar apparatuses for which the same or some frequency bands overlap.
In
Pseudo random noise (PN) codes, M-sequence codes, or Gold codes may be used as pseudo-random codes, for example. Pseudo-random code sequence RCode is constituted by, for example, LRC code elements, as given by following equation 53.
RCode={RC(1),RC(2), . . . ,RC(NLRC)} (Equation 53)
Code elements of a pseudo-random code sequence include, for example, values of 1 to −1, where a phase rotation of 0 is applied for 1 and a phase rotation of π is applied for −1. A pseudo-random code sequence has code length NLRC less than or equal to Nc. Further, for example, random code applier 111 cyclically varies code element index RC_INDEX(m) for the pseudo-random code sequence every m-th transmission period, as given by following equation 54.
[53]
RC_INDEX(m)=mod(m−1,NLRC)+1 (Equation 54)
Further, random code applier 111 outputs random code element RC(RC_INDEX(m)) of pseudo-random code sequence RCode to random code multiplier 214.
In radar transmitter 100b, each of Nt phase rotators 107 applies, for each transmission period Tr, coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) to which code element RC(RC_INDEX(m)) of pseudo-random code sequence RCode is applied, namely, phase rotation amount ψndop_code(ndm), ndm(m) angle[RC(RC_INDEX(m))], to a chirp signal output from radar transmission signal generator 101.
For example, a description will be given of a case where when the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NCM=2, orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, and the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1 and NDOP_CODE(2)=2.
In this case, random code applier 111 applies random code RC(RC_INDEX(m)) to coded Doppler phase rotation amounts 111, ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) and outputs the resulting phase rotation amounts to phase rotators 107 for each transmission period.
Phase rotator PROT #[1, 1] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 55 for each transmission period. The output of phase rotator PROT #[1, 1] is output from transmit antenna Tx #[1, 1]. Here, cp(t) denotes a chirp signal for each transmission period.
[54]
exp[j{ψ1,1(1)+angle[RC(1)]}]cp(t),exp[j{ψ1,1(2)+angle[RC(2)]}]cp(t), . . . , (Expression 55)
Likewise, phase rotator PROT #[1, 2] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 56 for each transmission period. The output of phase rotator PROT #[1, 2] is output from transmit antenna Tx #[1, 2].
[55]
exp[j{ψ1,2(1)+angle[RC(1)]}]cp(t),exp[j{ψ1,2(2)+angle[RC(2)]}]cp(t), . . . , (Expression 56)
Likewise, phase rotator PROT #[2, 2] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 57 for each transmission period. The output of phase rotator PROT #[2, 2] is output from transmit antenna Tx #[2, 2].
[56]
exp[j{ψ2,2(1)+angle[RC(1)]}]cp(t),exp[j{ψ2,2(2)+angle[RC(2)]}]cp(t), . . . , (Expression 57)
In radar receiver 200b, random code multiplier 214 multiplies output signal RFTz(fb, m) of beat frequency analyzer 208 in transmission period m by random code element RC(RC_INDEX(m)) input from random code applier 111. Random code multiplier 214 outputs a signal represented by RC(RC_INDEX(m))×RFTz(fb, m) to output switching section 209. Here, z=1, . . . , Na.
With the operation described above, even in the presence of interference from a plurality of radar apparatuses for which the same or some frequency bands overlap, in radar apparatus 10b, the interference signals can be converted into pseudo-random signals by random code multiplier 214 before being input to Doppler analyzers 210. Accordingly, in the outputs of Doppler analyzers 210, the effect of diffusing the signal power of interference waves over the Doppler frequency domain is achieved. For example, multiplication of a pseudo-random code sequence can reduce peak power of an interference wave to about 1/Ncode. This can greatly reduce the probability of CFAR section 211 in the subsequent stage erroneously detecting peaks of interference waves.
If phase rotators 107 do not have sufficient phase rotation accuracy, a phase rotation error caused by the application of a random code during a transmission period of an orthogonal code may cause interference between orthogonal codes. To address the interference, for example, code length NLRC of the pseudo-random code sequence is set to be less than or equal to Ncode, and random code applier 111 may apply the same random code element during the transmission period of orthogonal code length Loc. For example, random code applier 111 may set pseudo-random code element index RC_INDEX(m) given by following equation 58 for each transmission period m.
This allows the same random code element to be applied during the transmission period of an orthogonal code. Even when phase rotators 107 do not have sufficient phase rotation accuracy, phase rotation errors caused by random codes are constant, and interference between orthogonal codes can be reduced.
For example, a description will be given of a case where when RC_INDEX(m) given by equation 58 is used, the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NCM=2, orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, and the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1 and NDOP_CODE(2)=2.
In this case, random code applier 111 applies random code RC(RC_INDEX(m)) to coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) and outputs the resulting phase rotation amounts to phase rotators 107 for each transmission period.
Phase rotator PROT #[1, 1] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 59 for each transmission period. The output of phase rotator PROT #[1, 1] is output from transmit antenna Tx #[1, 1]. Here, cp(t) denotes a chirp signal for each transmission period.
[58]
exp[j{ψ1,1(1)+angle[RC(1)]}]cp(t),exp[j{ψ1,1(2)+angle[RC(1)]}]cp(t),exp[j{ψ1,1(3)+angle[RC(2)]}]cp(t),exp[j{ψ1,1(4)+angle[RC(2)]}]cp(t), . . . , (Expression 59)
Likewise, phase rotator PROT #[1, 2] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 60 for each transmission period. The output of phase rotator PROT #[1, 2] is output from transmit antenna Tx #[1, 2].
[59]
exp[j{ψ1,2(1)+angle[RC(1)]}]cp(t),exp[j{ψ1,2(2)+angle[RC(1)]}]cp(t),exp[j{ψ1,2(3)+angle[RC(2)]}]cp(t),exp[j{ψ1,2(4)+angle[RC(2)]}]cp(t), . . . , (Expression 60)
Likewise, phase rotator PROT #[2, 2] applies, for each transmission period, phase rotation to a chirp signal generated in radar transmission signal generator 101 in a manner given by following expression 61 for each transmission period. The output of phase rotator PROT #[2, 2] is output from transmit antenna Tx #[2, 2].
[60]
exp[j{ψ2,2(1)+angle[RC(1)]}]cp(t),exp[j{ψ2,2(2)+angle[RC(1)]}]cp(t),exp[j{ψ2,2(3)+angle[RC(2)]}]cp(t),exp[j{ψ2,2(4)+angle[RC(2)]}]cp(t), . . . , (Expression 61)
(Variation 6 of Embodiment 1)
Embodiment 1 has described the case where the output of phase rotator PROT #[ndop_code(ndm), ndm] is transmitted from transmit antenna Tx #[ndop_code(ndm), ndm]; however, this is not required.
For example, the association between each of plural transmit antennas 108 and the assignment of a coded Doppler multiplexed signal (in other words, a combination of Doppler shift amount and code sequence) may be set to vary for each frame in which a radar transmission signal is transmitted.
For example, when continuously performing radar position measurement, radar apparatus 10 may vary transmit antenna 108, which transmits the output of phase rotator PROT #[ndop_code(ndm), ndm], for each radar position measurement (for example, every Nc transmission periods (Nc×Tr)). Here, ndm=1, . . . , NDM, and ndop_code(ndm)=1, . . . , NDOP_CODE(ndm).
For example, radar apparatus 10 may hold a plurality of assignment tables for assignment as to from which of Nt transmit antennas #1, . . . , and #Nt to transmit the output of each of Nt phase rotators PROT #[ndop_code(mdm), ndm]. For example, radar apparatus 10 can modify an assignment table for each radar position measurement (for example, every Nc transmission periods (Nc×Tr)) to variably set transmit antenna 108 to be used for transmission for each radar position measurement.
Accordingly, when continuously performing radar position measurement, radar apparatus 10 variably sets, for the output of each of Nt phase rotators 107 in each radar position measurement, the corresponding one of transmit antennas 108. Accordingly, in cases such as when signals affected by interference (for example, inter-code interference) different depending on transmit antenna 108 are received, the effect of randomizing the influence of interference can be achieved by varying transmit antennas 108.
(Variation 7 of Embodiment 1)
In Embodiment 1, coder 106 applies phase rotation based on one or a plurality of orthogonal code sequences less than or equal to NCM to each of phase rotation amounts ϕ1, . . . , and ϕNDM for applying NDM Doppler shift amounts output from Doppler shift setter 105 to set coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) and outputs coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) to phase rotators 107. However, the processing performed by coder 106 is not limited to this.
Variation 7 describes a case where associations between transmit antennas 108 and Doppler multiplexed signals are set to vary for each transmission period.
For example, the assignment of Doppler multiplexing (in other words, Doppler shift amounts) to transmit antennas 108 may vary for each transmission period without changing the number of Doppler multiplexes used for coded Doppler multiplexing transmission.
For example, coder 106 sets, for each code element of an orthogonal code sequence (for example, each transmission period Tr), a coded Doppler phase rotation amount using a phase rotation amount for applying a different Doppler shift amount. In other words, coder 106 may make the value of Doppler shift amount DOPndm different for each code element of an orthogonal code sequence (for example, each transmission period Tr) in coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) applied to a radar transmission signal to be transmitted from each transmit antenna 108.
The setting of coded Doppler phase rotation amounts in the manner described above achieves effects similar to those of Embodiment 1. Further, for example, when signals affected by interference (for example, inter-code interference) different depending on transmit antenna 108 are received, the effect of randomizing the influence of interference can be achieved by varying transmit antennas 108.
For example, coder 106 may set coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) using following equation 62 instead of equation 5.
In coded Doppler phase rotation amount ψndop_code(ndm), ndm(m) given by equation 62, a phase rotation amount for applying a Doppler shift amount is set to vary for transmission periods in such manner as to be ϕmod(ndm+OC_INDEX-2,NDM)+1 for the duration of Loc transmission periods, the number of which is equal to the code length used for coding (the first term in equation 62), and the phase rotation amounts of Loc code elements OCndop_code(ndm)(1), . . . , and OCndop_code(ndm)(Loc) of code Codendop_code(ndm) used for coding are applied (the second term in equation 62). In equation 62, phase rotation amount ϕmod(ndm+OC_INDEX-2,NDM)+1 for applying a Doppler shift amount is set to vary for each code element (for example, OC_INDEX).
As an example, a description will be given of a case where when Nt=3, NDM=2, and NCM=2, coded Doppler phase rotation amount ψndop_code(ndm), ndm1(m) given by equation 62 is set.
In this case, the assignment of Doppler shift amounts DOP1 and DOP2 and orthogonal codes Code1 and Code2 is determined in accordance with, for example, the setting of NDOP_CODE(1) and NDOP_CODE(2), as illustrated in
As illustrated in
For example, a description will be given of a case where in coder 106, using equation 62, when the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NCM=2, and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1 and NDOP_CODE(2)=2 (see, for example,
[62]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }{0,0,ϕ1,ϕ2,2ϕ1,2ϕ2,3ϕ1,3ϕ2, . . . } (Equation 63)
[63]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }{0,0,ϕ2,ϕ1,2ϕ2,2ϕ1,3ϕ2,3ϕ1, . . . } (Equation 64)
[64]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }{0,π,ϕ2,ϕ1+π,2ϕ2,2ϕ1+π,3ϕ2,3ϕ1+π, . . . } (Equation 65)
In each of coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by equations 63 to 65, phase rotation amounts ϕ1 and ϕ2 for applying Doppler shift amounts DOP1 and DOP2 are alternately used in a period with code length Loc=2.
For ψ1, 2(m) and ψ2, 2(m) given by equations 64 and 65, phase rotations using orthogonal code sequences Code1 and Code2 are applied, with a relationship maintained in which phase rotation amounts ϕ1 and ϕ2 for applying Doppler shift amounts DOP1 and DOP2 are the same phase rotation amount in a period with code length Loc=2. In other words, for ψ1, 2(m) and ψ2, 2(m), phase rotation amounts for applying Doppler shift amounts are changed in a similar manner in a period with code length Loc, and a plurality of orthogonal code sequences are used to perform code multiplexing.
For ψ1, 1(m) and ψ1, 2(m) given by equations 63 and 64, phase rotation amounts ϕ1 and ϕ2 for applying Doppler shift amounts DOP1 and DOP2 are different phase rotation amounts in a period with code length Loc=2. Also for ψ1, 1(m) and ψ2, 2(m) given by equations 63 and 65, phase rotation amounts ϕ1 and ϕ2 for applying Doppler shift amounts DOP1 and DOP2 are different phase rotation amounts in a period with code length Loc=2.
Here, for example, a description will be given of a case where a phase rotation amount for applying Doppler shift amount DOPndm is ϕndm=2π(ndm−1)/NDM given by equation 1, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, and phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π, are used. In this case, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 66 to 68 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) to phase rotators 107. Here, m=1, . . . , Nc.
[65]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }{0,0,0,π,0,0,0,π, . . . } (Equation 66)
[66]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }{0,0,π,0,0,0,π,0, . . . } (Equation 67)
[67]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }{0,π,π,π,0,π,π,π, . . . } (Equation 68)
As another example, a description will be given of a case where when Nt=6, NDM=4, and NCM=2, coded Doppler phase rotation amount ψndop_code(ndm), ndm1(m) given by equation 62 is set.
In this case, the assignment of Doppler shift amounts DOP1, DOP2, DOP3, and DOP4 and orthogonal codes Code1 and Code2 is determined in accordance with, for example, the setting of NDOP_CODE(1), NDOP_CODE(2), NDOP_CODE(3), and NDOP_CODE(4), as illustrated in
As illustrated in
For example, a description will be given of a case where in coder 106, using equation 62, when the number of transmit antennas used for multiplexing transmission Nt=6, the number of Doppler multiplexes NDM=4, NCM=2, and orthogonal code sequences Code1={1, 1} and Code2={1, −1} with code length Loc=2 are used, the numbers of coded Doppler multiplexes are set such that NDOP_CODE(1)=1, NDOP_CODE(2)=1, NDOP_CODE(3)=2, and NDOP_CODE(4)=2. In this case, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 69 to 74 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 107. Here, m=1, . . . , Nc.
[68]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }{0,0,ϕ1,ϕ2,2ϕ1,2ϕ2,3ϕ1,3ϕ2, . . . } (Equation 69)
[69]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }{0,0,ϕ2,ϕ3,2ϕ2,2ϕ3,3ϕ2,3ϕ3, . . . } (Equation 70)
[70]
{ψ1,3(1),ψ1,3(2),ψ1,3(3),ψ1,3(4),ψ1,3(5),ψ1,3(6),ψ1,3(7),ψ1,3(8), . . . }{0,0,ϕ3,ϕ4,2ϕ3,2ϕ4,3ϕ3,3ϕ4, . . . } (Equation 71)
[71]
{ψ2,3(1),ψ2,3(2),ψ2,3(3),ψ2,3(4),ψ2,3(5),ψ2,3(6),ψ2,3(7),ψ2,3(8), . . . }{0,π,ϕ3,ϕ4+π,2ϕ3,2ϕ4+π,3ϕ3,3ϕ4+π, . . . } (Equation 72)
[72]
{ψ1,4(1),ψ1,4(2),ψ1,4(3),ψ1,4(4),ψ1,4(5),ψ1,4(6),ψ1,4(7),ψ1,4(8), . . . }{0,0,ϕ4,ϕ1,2ϕ4,2ϕ1,3ϕ4,3ϕ1, . . . } (Equation 73)
[73]
{ψ2,4(1),ψ2,4(2),ψ2,4(3),ψ2,4(4),ψ2,4(5),ψ2,4(6),ψ2,4(7),ψ2,4(8), . . . }{0,π,ϕ4,ϕ1+π,2ϕ4,2ϕ1+π,3ϕ4,3ϕ1+π, . . . } (Equation 74)
In each of coded Doppler phase rotation amounts ψndop_code(ndm), ndm(m) given by equations 69 to 74, phase rotation amounts ϕndm and ϕmod(ndm, NDM)+1 for applying Doppler shift amounts DOPndm and DOPmod(ndm, NDM)+1 are used in a period with code length Loc=2. In coded Doppler phase rotation amounts ψndop_code(ndm), ndm(m) with ndm being different, phase rotation amounts for applying Doppler shift amounts are different.
Further, for coded Doppler phase rotation amounts ψ1, ndm(m), . . . , and ψn2dop_code(ndm), ndm(m), phase rotations using orthogonal code sequences Code1, . . . , and Codendop_code(ndm) are applied, with a relationship maintained in which phase rotation amounts ϕndm and ϕmod(ndm, NDM)+1 for applying Doppler shift amounts DOPndm and DOPmod(ndm,NDM)+1 are the same phase rotation amount in a period with code length Loc=2.
Here, for example, a description will be given of a case where a phase rotation amount for applying Doppler shift amount DOPndm is ϕndm=2π(ndm−1)/NDM given by equation 1, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π/2, phase rotation amount ϕ3 for applying Doppler shift amount DOP3, which is equal to π, and phase rotation amount ϕ4 for applying Doppler shift amount DOP4, which is equal to 3π/2, are used. In this case, coder 106 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 75 to 80 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 107. Here, m=1, . . . , Nc.
Next, an example of the operation of radar receiver 200 when coded Doppler phase rotation amounts are set by coder 106 described above will be described. In radar receiver 200, a code separation process performed by coded Doppler demultiplexer 212 is different from that in Embodiment 1.
Coded Doppler demultiplexer 212 performs a code separation process on the outputs of Doppler analyzers 210 in the z-th signal processor 206 indicated by Doppler frequency indices (fs_comp_cfar±(nfd−1)×ΔFD) of NDM coded Doppler multiplexed signals for distance indices fb_cfar output from CFAR section 211.
The code separation process may be performed on, for example, all of nfd=1, . . . , NDM for all the candidates of DopCase=1, . . . , NDM. Note that coded Doppler demultiplexer 212 detects a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM, and performs discrimination of transmit antennas 108 and determination of a target Doppler frequency. Accordingly, coded Doppler demultiplexer 212 performs a code separation process, as given by following equation 81, to reduce the operation amount of the separation process.
The superscript asterisk (*) indicates the complex conjugate operator. Further, nfd=1, . . . , NDM, ncm=1, . . . , NCM, and DopCase=1, . . . , NDM.
For example, when a coded Doppler multiplexed signal for which the number of coded Doppler multiplexes is set to be smaller than NCM uses coded Doppler phase rotation amount ψndop_code(ndm), ndm1(m) given by equation 62, coded Doppler demultiplexer 212 performs the code separation process given by equation 81, with consideration given to the use of DOPmod(ndm1+noc-2, NDM)+1 as the Doppler shift amount for the noc-th code element.
In equation 81, coded Doppler demultiplexer 212 performs a code separation process using candidate DOPposi(DOPmod(ndm1+noc-2, NDM)+1, DopCase) including Doppler shift amount DOPmod(ndm1+noc-2, NDM)+1 for the output of Doppler analyzer 210 for every noc-th code element within each DopCase. Here, code separation signal DeMULzncm(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) given by equation 81 represents a code separation signal that uses Codencm for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) in the z-th signal processor 206. The Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) represents a Doppler frequency index in the output of Doppler analyzer 210 for the first code element when the coded Doppler phase rotation amount given by equation 62 is used.
In equation 81, in the exp term, since the sampling times for the outputs of Doppler analyzers 210 for each code element are shifted, phase correction is performed in accordance with Doppler frequency index (fs_comp_cfar+(DopCase−1)×ΔFD-DOP1).
Further, Doppler shift amount DOPndm set in Doppler shift setter 105 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, and DOP1 falls within the range of fs_comp=−Ncode/2, −Ncode/2+ΔFD−1 in the initial state (when the relative velocity to the target is zero). Accordingly, coded Doppler demultiplexer 212 calculates an amount of phase correction using, for example, DOP1 as a reference.
The subsequent processing performed by radar receiver 200 is similar to that in Embodiment 1, and a description thereof is omitted here.
(Variation 8 of Embodiment 1)
Embodiment 1 has described a case where when radar apparatus 10 repeatedly transmits chirp pulses Nc times as radar transmission signals, the center frequencies of the chirp signals are made constant (see, for example,
Variation 8 describes a case where the center frequencies of chirp signals are variably set.
[Configuration of Radar Apparatus]
In the following, for example, radar apparatus 10c transmits radar transmission signals such that center frequencies fc of chirp signals are changed by Δf (for example, increased in a case where Δf>0 and decreased in a case where Δf<0) for every transmission period Tr.
In radar transmitter 100c, radar transmission signal generator 101c includes modulated signal emitter 102, VCO 103c, and transmission frequency controller 112.
For example, modulated signal emitter 102 periodically emits sawtooth-shaped modulated signals for VCO control. Here, the transmission period is represented by Tr.
Transmission frequency controller 112 controls, for each transmission period Tr, center frequency fc of a frequency-modulated signal (chirp signal) to be output from VCO 103c. For example, transmission frequency controller 112 may change center frequency fc of a frequency-modulated signal by Δf for each transmission period Tr.
VCO 103c outputs frequency-modulated signals to phase rotators 107 and radar receiver 200c (for example, mixer section 204), based on the output of transmission frequency controller 112 and the output of modulated signal emitter 102.
In
That is, in
The respective chirp signals may be, for example, chirp signals having the same frequency-modulation bandwidth Bw in time width TA of a range gate. In the example illustrated in
Other operation of radar transmitter 100c illustrated in
Next, an example of the operation of radar receiver 200c of radar apparatus 10c will be described.
In radar receiver 200c, the processing performed by antenna channel processors 201 on signals received by receive antennas 202, and the operation of subsequent CFAR section 211 and coded Doppler demultiplexer 212 are similar to the operation in Embodiment 1. In radar receiver 200c, furthermore, a direction estimation process performed by direction estimator 213c using the output of coded Doppler demultiplexer 212 is also similar to the operation in Embodiment 1.
In radar receiver 200c, for example, a conversion process for Doppler velocity information of the target based on the Doppler frequency determination result for the target (for example, Doppler frequency fTARGET of the target), which is performed by direction estimator 213c, is different from that in Embodiment 1.
The conversion of distance information R(fb) based on the beat frequency index is similar to that in Embodiment 1, and direction estimator 213c may output, for example, based on equation 31, distance information R(fb) using beat frequency index (or distance index) fb.
Direction estimator 213c may output Doppler velocity information vd of the target detected in the following way using, for example, Doppler frequency fTARGET of the target and distance index fb_cfar.
For example, when radar transmission signals for which center frequencies fc of chirp signals are changed by Δf for every transmission period Tr are used, even if the relative velocity of the target is zero, center frequencies fc of chirp signals are changed for every transmission period Tr. Accordingly, a reception signal of radar apparatus 10c includes phase rotation caused by a change in the center frequency of a chirp signal for each transmission period Tr.
Center frequency fc in the m-th transmission period Tr for target distance Rtarget is changed by (m−1)Δf relative to the center frequency in the first transmission period Tr as a reference. Phase rotation amount Δη(m, Rtarget) caused by the change in center frequency fc is given by equation 81-1, with consideration given to reflected-wave arrival time (2Rtarget/Co) from target distance Rtarget. Equation 81-1 indicates the relative phase rotation amount obtained using the phase of the first transmission period Tr as a reference. C0 denotes the speed of light
Therefore, as given by following equation 81-2, direction estimator 213c calculates Doppler velocity information vd(fTARGET, fb_cfar), based on a transformation equation taking into account Δf that is the amount of change in center frequency fc of a chirp signal for each transmission period Tr.
In equation 81-2, the first term corresponds to equation 49 and is a relative Doppler velocity component represented by Doppler frequency fTARGET. The second term in equation 81-2 is a Doppler velocity component that is generated by changing center frequency fc of a chirp signal by Δf for each transmission period Tr. For example, as given by equation 81-2, direction estimator 213c can calculate the true relative Doppler velocity vd(fTARGET, fb_cfar) of the target by removing the Doppler component in the second term from the first term. Here, R(fb_cfar) denotes distance information (distance estimation value) calculated from beat frequency index fb_cfar in accordance with equation 31.
The Doppler range of the target is assumed to be up to ±1/(2×Tr). If vd(fTARGET, fb_cfar) satisfies vd(fTARGET, fb_cfar)<−C0/(4f0Tr), direction estimator 213c may output detected Doppler velocity information vd(fTARGET, fb_cfar) of the target in accordance with following equation 81-3.
Also, since the Doppler range of the target is assumed to be up to ±1/(2×Tr), if vd(fTARGET, fb_cfar) satisfies vd(fTARGET, fb_cfar)>C0/(4f0 Tr), direction estimator 213c may output detected Doppler velocity information vd(fTARGET, fb_cfar) of the target in accordance with following equation 81-4.
As described above, in Variation 8, in radar apparatus 10c, center frequencies fc of chirp signals are changed based on transmission periods Tr of radar transmission signals. For example, radar apparatus 10c transmits radar transmission signals such that center frequencies fc of chirp signals are changed by Δf (for example, increased in a case where Δf>0 and decreased in a case where Δf<0) for every transmission period Tr. Even in this case, as in Embodiment 1, radar apparatus 10c (for example, MIMO radar) can extend the effective Doppler frequency bandwidth to 1/(Tr) and can extend the Doppler frequency (relative velocity) detection range with no ambiguity. In addition, radar apparatus 10c can reduce mutual interference between multiplexed signals to about the noise level. According to Variation 8, therefore, radar apparatus 10c can improve target-object sensing accuracy over a wider Doppler frequency range.
Variation 8 is applicable not only to Embodiment 1 but also to Variations 1 to 7 of Embodiment 1, and respective similar effects can be achieved. For example, in any of Variations 1 to 7 of Embodiment 1, radar transmission signals including chirp signals whose center frequencies fc are changed in the manner illustrated in
In Variation 8, furthermore, for example, radar apparatus 10c transmits radar transmission signals in such a manner that center frequencies fc of chirp signals are changed by Δf for every transmission period Tr, and can thus improve the distance resolution by a change width of the center frequencies of chirp signals (see, for example, NPL 4). According to Variation 8, since the distance resolution can be improved by a change width of the center frequencies of chirp signals, the chirp sweep bandwidth (for example, Bw) can be reduced compared with transmission with the center frequencies of chirp signals kept constant. The reduction in chirp sweep bandwidth can reduce, for example, the transmission period Tr while improving the distance resolution. As a result, in code multiplexing transmission, the detectable Doppler range without ambiguity can further be extended.
(Variation 9 of Embodiment 1)
The periods in which the center frequencies of chirp signals are changed are not limited to transmission periods Tr as in Variation 8. Variation 9 describes a case where the center frequencies of chirp signals are set to vary in increments of transmission periods of a plurality of chirp signals.
For example, Variation 9 describes a case where the center frequencies of chirp signals are set to vary for every Loc transmission periods (Loc×Tr) (in other words, a transmission period of an orthogonal code sequence, hereinafter referred to as “code transmission period”), the number of which is equal to the code length of a single orthogonal code used for Doppler multiplexing transmission.
[Configuration of Radar Apparatus]
In Variation 9, for example, radar apparatus 10d transmits radar transmission signals such that center frequencies fc of chirp signals are changed by Δf (for example, increased in a case where Δf>0 and decreased in a case where Δf<0) for every code transmission period (Loc×Tr).
In radar transmitter 100d, radar transmission signal generator 101d includes modulated signal emitter 102, VCO 103d, and transmission frequency controller 112d. In radar transmitter 100d, phase rotation amount setter 104d includes Doppler shift setter 105 and coder 106d.
For example, in radar transmission signal generator 101d, modulated signal emitter 102 periodically emits sawtooth-shaped modulated signals for VCO control. Here, the transmission period is represented by Tr.
Transmission frequency controller 112d controls, based on orthogonal code element index OC_INDEX output from coder 106d, center frequency fc of a frequency-modulated signal (chirp signal) to be output from VCO 103d for each code transmission period (Loc×Tr).
For example, transmission frequency controller 112d may change center frequency fc of a frequency-modulated signal to be output from VCO 103d by Δf in transmission period Tr corresponding to OC_INDEX=1. In other words, transmission frequency controller 112d controls, in transmission period Tr corresponding to OC_INDEX≠1, center frequency fc of a frequency-modulated signal to be output from VCO 103d to be the same as center frequency fc in previous transmission period Tr. With this control, transmission frequency controller 112d can perform control such that center frequencies fc are changed by Δf for every code transmission period (Loc×Tr).
VCO 103d outputs frequency-modulated signals to phase rotators 107 and radar receiver 200d (for example, mixer section 204), based on the output of transmission frequency controller 112d and the output of modulated signal emitter 102.
In
VCO 103d outputs a chirp signal having center frequency fc(Loc+1) equal to f0+Δf during the (Loc+1)-th transmission period Tr #(Loc+1). Further, VCO 103d outputs chirp signals having center frequencies fc(Loc+2) to fc(2Loc) equal to f0+Δf during the (Loc+2)-th transmission period Tr #(Loc+2) to the (2Loc)-th transmission period Tr #(2Loc), respectively.
Likewise, VCO 103d outputs a chirp signal having center frequency fc(m) equal to f0+floor[(m−1)/Loc]Δf during the m-th transmission period Tr #m.
That is, in
The respective chirp signals may be, for example, chirp signals having the same frequency-modulation bandwidth Bw in time width TA of a range gate. In the example illustrated in
Other operation of radar transmitter 100d illustrated in
Next, an example of the operation of radar receiver 200d of radar apparatus 10d will be described.
In radar receiver 200d, the processing performed by antenna channel processors 201 on signals received by receive antennas 202, and the operation of subsequent CFAR section 211 are similar to the operation in Embodiment 1. In radar receiver 200d, furthermore, a direction estimation process performed by direction estimator 213d using the output of coded Doppler demultiplexer 212d is also similar to the operation in Embodiment 1.
In radar receiver 200d, for example, the operation of coded Doppler demultiplexer 212d and a conversion process for Doppler velocity information of the target, which is performed by direction estimator 213d, is different from those in Embodiment 1.
The conversion of distance information R(fb) based on the beat frequency index is similar to that in Embodiment 1, and direction estimator 213d may output, for example, based on equation 31, distance information R(fb) using beat frequency index (or distance index) fb.
The following describes an example of the operation of coded Doppler demultiplexer 212d, which is different from that in Embodiment 1.
For example, when radar transmission signals for which center frequencies fc of chirp signals are changed by Δf for every code transmission period (Loc×Tr) are used, even if the relative velocity of the target is zero, center frequencies fc of chirp signals are changed for every code transmission period (Loc×Tr). Accordingly, the output of each of Loc Doppler analyzers 210 of radar apparatus 10d includes phase rotation caused by a change in the center frequency of a chirp signal for each code transmission period (Loc×Tr).
For example, center frequency fc in the m-th transmission period Tr for target distance Rtarget is changed by floor[(m−1)/Loc]Δf relative to center frequency fc in the first transmission period Tr as a reference. Accordingly, phase rotation amount Δη(m, Rtarget) caused by the change in center frequency is given by equation 81-5, with consideration given to reflected-wave arrival time (2Rtarget/C0) from target distance Rtarget Equation 81-5 indicates the relative phase rotation amount obtained using the phase of the first transmission period Tr as a reference. C0 denotes the speed of light.
Each of Loc Doppler analyzers 210 performs Doppler analysis that takes into account the phase rotation given by equation 81-5. Here, code transmission period (Loc×Tr) in which center frequency fc of a chirp signal is changed by Δf matches the period for switching Doppler analyzers 210 for every Loc code elements. Accordingly, to correct Doppler phase rotation caused by the time difference of Doppler analysis among Loc Doppler analyzers 210 to separate code-multiplexed signals, coded Doppler demultiplexer 212d uses following equation 81-6 instead of equation 35 described in <(1) code separation process> in Embodiment 1.
Equation 81-6 is different from equation 35 in the additional term for correcting the phase rotation caused by the change in center frequency:
Here, R(fb_cfar) denotes a distance estimation value calculated from beat frequency index fb_cfar in accordance with equation 31.
In equation 81-6, due to the change of Δf in reflected-wave arrival time (2R(fb_cfar)/Co) from R(fb_cfar), the phase rotation amount is given by 2πΔf×(2R(fb_cfar)/Co) in code transmission period (Loc×Tr). Therefore, phase rotation caused by the time difference of Doppler analysis among Loc Doppler analyzers 210 is, for example, for the noc-th Doppler analyzer 210, (noc−1)/Loc times relative to that for the first Doppler analyzer 210. Equation 81-6 is derived from this situation. Note that noc=1, . . . , Loc.
In this manner, the period in which center frequency fc of a chirp signal is changed by Δf is made equal to code transmission period (Loc×Tr) (or transmission periods, the number of which is equal to a divisor of Loc). Accordingly, the period in which the center frequency of a chirp signal is changed matches the period for switching Doppler analyzers 210 for each code element. For example, the period in which center frequency fc of a chirp signal is changed by Δf matches the period for switching Doppler analyzers 210 for each code element, thereby allowing coded Doppler demultiplexer 212d to uniquely decide a phase correction value for each of Loc Doppler analyzers 210. This enables coded Doppler demultiplexer 212d to easily perform phase correction in (1) code separation process. Therefore, coded Doppler demultiplexer 212d can easily perform the code separation process.
Next, an example of the operation of direction estimator 213d that is different from that in Embodiment 1 (for example, a conversion process for Doppler velocity information of the target) will be described.
Direction estimator 213d may output Doppler velocity information vd of the target detected in the following way using, for example, Doppler frequency fTARGET of the target and distance index fb_cfar.
When radar transmission signals for which center frequencies fc of chirp signals are changed by Δf for every code transmission period (Loc×Tr) are used, even if the relative velocity of the target is zero, center frequencies fc of chirp signals are changed for every code transmission period (Loc×Tr). Accordingly, a reception signal of radar apparatus 10d includes phase rotation caused by a change in the center frequency of a chirp signal for each code transmission period (Loc×Tr).
Center frequency fc in the m-th transmission period Tr for target distance Rtarget is changed by floor[(m−1)/Loc]Δf relative to the center frequency in the first transmission period Tr as a reference. Phase rotation amount Δη(m, Rtarget) caused by the change in center frequency fc is given by equation 81-7, with consideration given to reflected-wave arrival time (2Rtarget/Co) from target distance Rtarget. Equation 81-7 indicates the relative phase rotation amount obtained using the phase of the first transmission period Tr as a reference. C0 denotes the speed of light
Therefore, as given by following equation 81-8, direction estimator 213d calculates Doppler velocity information vd(fTARGET, fb_cfar), based on a transformation equation taking into account Δf that is the amount of change in center frequency fc of a chirp signal for each code transmission period (Loc×Tr).
In equation 81-8, the first term corresponds to equation 49 and is a relative Doppler velocity component represented by Doppler frequency fTARGET. The second term in equation 81-8 is a Doppler velocity component that is generated by changing center frequency fc of a chirp signal by Δf for each code transmission period (Loc×Tr). For example, as given by equation 81-8, direction estimator 213d can calculate the true relative Doppler velocity vd(fTARGET, fb_cfar) of the target by removing the Doppler component in the second term from the first term. Here, R(fb_cfar) denotes distance information (distance estimation value) calculated from beat frequency index fb_cfar in accordance with equation 31.
The Doppler range of the target is assumed to be up to ±1/(2×Tr). If vd(fTARGET, fb_cfar) satisfies vd(fTARGET, fb_cfar)<−C0/(4f0Tr), direction estimator 213d may output detected Doppler velocity information vd(fTARGET, fb_cfar) of the target in accordance with following equation 81-9.
Also, since the Doppler range of the target is assumed to be up to ±1/(2×Tr), if vd(fTARGET, fb_cfar) satisfies vd(fTARGET, fb_cfar)>C0/(4f0 Tr), direction estimator 213d may output detected Doppler velocity information vd(fTARGET, fb_cfar) of the target in accordance with following equation 81-10.
As described above, in Variation 9, in radar apparatus 10d, center frequencies fc of chirp signals are changed based on code transmission periods (Loc×Tr). For example, radar apparatus 10d transmits radar transmission signals such that center frequencies fc of chirp signals are changed by Δf (for example, increased in a case where Δf>0 and decreased in a case where Δf<0) for every code transmission period (Loc×Tr). Even in this case, as in Embodiment 1, radar apparatus 10d (for example, MIMO radar) can extend the effective Doppler frequency bandwidth to 1/(Tr) and can extend the Doppler frequency (relative velocity) detection range with no ambiguity. In addition, radar apparatus 10d can reduce mutual interference between code-multiplexed signals to about the noise level. According to Variation 9, therefore, radar apparatus 10d can improve target-object sensing accuracy over a wider Doppler frequency range.
Variation 9 is applicable not only to Embodiment 1 but also to Variations 1 to 7 of Embodiment 1, and respective similar effects can be achieved. For example, in any of Variations 1 to 7 of Embodiment 1, radar transmission signals including chirp signals whose center frequencies fc are changed in the manner illustrated in
According to Variation 9, furthermore, when the period in which center frequency fc of a chirp signal is changed by Δf is a plurality of transmission periods Tr, for example, the period in which center frequency fc of a chirp signal is changed is made to match code transmission period (Loc×Tr). Accordingly, the period in which center frequency fc of a chirp signal is changed is also made to match the period for switching Doppler analyzers 210 for each code element (in other words, the transmission period of a code used for code multiplexing). Therefore, coded Doppler demultiplexer 212d can easily perform phase correction in a code demultiplexing process.
In Variation 9, furthermore, radar apparatus 10d transmits, for example, radar transmission signals in such a manner that center frequencies fc of chirp signals are changed by Δf for every code transmission period (Loc×Tr), with the change width of the center frequencies of chirp signals being given by Δf×Ncode and the distance resolution being given by 0.5C0/(Δf×Ncode). For example, as the value given by Δf×Ncode increases, the distance resolution can be more improved by the change width of the center frequencies of chirp signals. According to Variation 9, since the distance resolution can be improved by a change width of the center frequencies of chirp signals, the chirp sweep bandwidth (for example, Bw) can be reduced compared with transmission with the center frequencies of chirp signals kept constant. The reduction in chirp sweep bandwidth can reduce, for example, the transmission period Tr while improving the distance resolution. As a result, in code multiplexing transmission, the detectable Doppler range without ambiguity can further be extended.
While Variation 9 has described the case where radar transmission signals are used for which center frequencies fc of chirp signals are changed by Δf for every code transmission period (Loc×Tr), the period in which center frequency fc of a chirp signal is changed is not limited to this. For example, radar apparatus 10d may use radar transmission signals for which center frequencies fc of chirp signals are changed by Δf for every transmission periods, the number of which is equal to a divisor of code length Loc (divisor of (Loc×Tr)). When 1 is used among divisors of code length Loc, as in Variation 8, center frequencies fc are changed by Δf for every transmission period Tr.
When radar transmission signals for which center frequencies fc of chirp signals are changed by Δf for every transmission periods, the number of which is equal to divisor E of code length Loc, that is, for every E transmission periods (E×Tr), are used, center frequencies fc of chirp signals are changed by an amount corresponding to Δf×Loc/ε for every code transmission period (Loc×Tr). Thus, in this case, Δf in equations 81-8, 81-9, and 81-10 may be replaced by Δf×Loc/ε. Further, coded Doppler demultiplexer 212d uses following equation 81-11 instead of equation 81-6 in the code separation process. Here, E denotes a divisor of Loc.
(Variation 10 of Embodiment 1)
The method for controlling center frequency fc of a frequency-modulated signal (chirp signal) is not limited to the methods according to Variation 8 (
The configuration of a radar apparatus according to Variation 10 may be similar to, for example, the configuration of radar apparatus 10d according to Variation 9.
In Variation 10, for example, radar apparatus 10d periodically changes center frequency fc of a chirp signal over a plurality of transmission periods within code transmission period (Loc×Tr). In this case, the timing at which center frequencies fc of chirp signals are changed in a round is made to match code transmission period (Loc×Tr), thereby allowing radar apparatus 10d to reduce the amount of signal processing performed by coded Doppler demultiplexer 212d (the details will be described below).
For example, in radar transmitter 100d, modulated signal emitter 102 periodically emits sawtooth-shaped modulated signals for VCO control. Here, the transmission period is represented by Tr.
Transmission frequency controller 112d controls, based on orthogonal code element index OC_INDEX output from coder 106d in phase rotation amount setter 104d, center frequency fc of a frequency-modulated signal (chirp signal) output from VCO 103d for each transmission period Tr.
For example, transmission frequency controller 112d sets center frequency fc of a frequency-modulated signal to be output from VCO 103d to f0 in transmission period Tr corresponding to OC_INDEX=1. Further, transmission frequency controller 112d sets center frequency fc of a frequency-modulated signal to be output from VCO 103d to f0+Δf in transmission period Tr corresponding to OC_INDEX=2. Likewise, transmission frequency controller 112d sets center frequencies fc of the frequency-modulated signals to be output from VCO 103d to (f0+2Δf) to (f0+(Loc−1)Δf) in transmission periods Tr where OC_INDEX=3 to Loc, respectively.
With this control, transmission frequency controller 112d can perform control such that, for example, center frequencies fc of chirp signals are periodically changed over every code transmission period (Loc×Tr).
VCO 103d outputs frequency-modulated signals to phase rotators 107 and radar receiver 200d (for example, mixer section 204), based on the output of transmission frequency controller 112d and the output of modulated signal emitter 102.
In
Further, VCO 103d outputs a chirp signal having center frequency fc(Loc+1) equal to f0 during the (Loc+1)-th transmission period Tr #(Loc+1). Likewise, VCO 103d outputs chirp signals having center frequencies fc(Loc+2) to fc(2Loc) equal to (f0+Δf) to (f0+(Loc−1)Δf) during the (Loc+2)-th transmission period Tr #(Loc+2) to the (2Loc)-th transmission period Tr #(2Loc), respectively.
Likewise, VCO 103d outputs a chirp signal having center frequency fc(m) equal to f0+mod(m−1, Loc)Δf during the m-th transmission period Tr #m.
That is, in
In this manner, in
The respective chirp signals may be, for example, chirp signals having the same frequency-modulation bandwidth Bw in time width TA of a range gate. In the example illustrated in
Other operation of radar transmitter 100d illustrated in
Next, an example of the operation of radar receiver 200d of radar apparatus 10d will be described.
In radar receiver 200d, the processing performed by antenna channel processors 201 on signals received by receive antennas 202, and the operation of subsequent CFAR section 211 are similar to the operation in Embodiment 1. In radar receiver 200d, furthermore, a direction estimation process performed by direction estimator 213d using the output of coded Doppler demultiplexer 212d is also similar to the operation in Embodiment 1.
In radar receiver 200d, for example, the operation of coded Doppler demultiplexer 212d and a conversion process for Doppler velocity information of the target, which is performed by direction estimator 213d, is different from that in Embodiment 1.
The conversion of distance information R(fb) based on the beat frequency index is similar to that in Embodiment 1, and direction estimator 213d may output, for example, based on equation 31, distance information R(fb) using beat frequency index (or distance index) fb.
The following describes an example of the operation of coded Doppler demultiplexer 212d, which is different from that in Embodiment 1.
For example, when radar transmission signals for which center frequencies fc of chirp signals are changed by f0, f0+Δf, f0+(Loc−1)Δf for every transmission period Tr within each code transmission period (Loc×Tr) are used, radar reflected waves obtained in response to the transmission of chirp signals whose center frequencies fc are f0, f0+Δf, and f0+(Loc−1)Δf are input to the first, second, . . . , and Loc-th Doppler analyzers 210 as reception signals.
Thus, the center frequencies of the respective radar reflected waves to be input to Loc Doppler analyzers 210 are the same. Accordingly, Loc Doppler analyzers 210 perform Doppler analysis on reception signals that are radar reflected waves obtained when chirp signals having the same center frequency fc are transmitted.
In contract, center frequencies fc of the chirp signals differ in Loc Doppler analyzers 210. Accordingly, to correct phase rotation caused by the time difference of Doppler analysis among Loc Doppler analyzers 210 to separate code-multiplexed signals, coded Doppler demultiplexer 212d uses following equation 81-12 instead of equation 35 described in <(1) code separation process> in Embodiment 1.
Equation 81-12 is different from equation 35 in the additional term for correcting the phase rotation caused by the change in center frequency:
Here, R(fb_cfar) denotes a distance estimation value calculated from beat frequency index fb_cfar in accordance with equation 31.
Equation 81-12 is derived from the following. For example, when the center frequency of a transmission chirp signal for the first Doppler analyzer 210 is used as a reference, the center frequencies of transmission chirp signals for the second to Loc-th Doppler analyzers 210 are different by Δf, . . . , and (Loc−1)Δf, respectively. Accordingly, phase rotation amounts in reflected-wave arrival time (2R(fb_cfar)/Co) from R(fb_cfar) are different.
That is, when the output of the first Doppler analyzer 210 is used as a phase reference, the phase rotation amount for the noc-th Doppler analyzer 210 is given by 2π(noc−1)Δf×(2R(fb_cfar)/Co). To compensate for the phase rotation, in equation 81-12, the term
is derived. Here, noc=1, . . . , Loc.
In this manner, when center frequencies fc of chirp signals are periodically changed over a plurality of chirp transmission periods within code transmission period (Loc×Tr) (or transmission periods, the number of which is equal to a divisor of Loc), the period for switching Doppler analyzers 210 for each code element matches the timing at which center frequencies fc of chirp signals are changed in a round. For example, the period in which center frequencies fc of chirp signals are periodically changed over a plurality of chirp transmission periods matches the period for switching Doppler analyzers 210 for each code element, thereby allowing coded Doppler demultiplexer 212d to uniquely decide a phase correction value for each of Loc Doppler analyzers 210. This enables coded Doppler demultiplexer 212d to easily perform phase correction in (1) code separation process. Therefore, coded Doppler demultiplexer 212d can easily perform the code separation process.
Next, an example of the operation of direction estimator 213d that is different from that in Embodiment 1 (for example, a conversion process for Doppler velocity information of the target) will be described.
For example, direction estimator 213d may calculate Doppler velocity information vd(fTARGET) of the target detected based on Doppler frequency fTARGET of the target determined by coded Doppler demultiplexer 212d in accordance with following equation 81-13 and output Doppler velocity information vd(fTARGET) as a position measurement result.
As described above, in Variation 10, in radar apparatus 10d, center frequencies fc of chirp signals are periodically changed over a plurality of transmission periods. Even in this case, as in Embodiment 1, radar apparatus 10d (for example, MIMO radar) can extend the effective Doppler frequency bandwidth to 1/(Tr) and can extend the Doppler frequency (relative velocity) detection range with no ambiguity. In addition, radar apparatus 10d can reduce mutual interference between code-multiplexed signals to about the noise level. According to Variation 10, therefore, radar apparatus 10d can improve target-object sensing accuracy over a wider Doppler frequency range.
Variation 10 is applicable not only to Embodiment 1 but also to Variations 1 to 7 of Embodiment 1, and respective similar effects can be achieved. For example, in any of Variations 1 to 7 of Embodiment 1, radar transmission signals including chirp signals whose center frequencies fc are changed in the manner illustrated in
According to Variation 10, furthermore, the timing at which center frequencies fc of chirp signals are changed in a round matches code transmission period (Loc×Tr). Accordingly, for example, the timing at which center frequencies fc of chirp signals are changed in a round can be made to match the period for switching Doppler analyzers 210 for each code element (in other words, the transmission period of a code used for code multiplexing). This enables coded Doppler demultiplexer 212d to easily perform phase correction.
While Variation 10 has described the case where the timing at which center frequencies fc of chirp signals are changed in a round matches code transmission period (Loc×Tr), the timing at which center frequencies fc of chirp signals are changed in a round may be equal to (divisor of (Loc×Tr)). When the timing at which center frequencies fc of chirp signals are changed in a round is a period given by (divisor E of Loc×Tr), that is, ε transmission periods (ε×Tr), coded Doppler demultiplexer 212d may use equation 81-14 instead of equation 81-12, where E denotes a divisor of Loc and satisfies ε>1.
Further, Variation 10 has described the case where radar transmission signals are used for which center frequencies fc of chirp signals are changed by an integer multiple of Δf, such as f0, f0+Δf, f0+(Loc−1)Δf for each transmission period Tr, in code transmission period (Loc×Tr). However, the change width of center frequencies fc is not limited to a frequency of an integer multiple of Δf, and center frequencies fc may be set to vary by any frequency.
For example, radar transmission signals may be used for which center frequencies fc of chirp signals are changed by f0, f0+Δf1, f0+Δf2, . . . , and f0+ΔfLoc-1 for every transmission period Tr in code transmission period (Loc×Tr). Here, Δf1, Δf2, . . . , and ΔfLoc-1 are variable frequency values of center frequencies fc of chirp signals in respective transmission periods Tr in code transmission period (Loc×Tr). In this case, coded Doppler demultiplexer 212d may use following equation 81-15 instead of equation 81-12 in the code separation process.
Embodiment 1 has described the case where both Doppler multiplexing transmission and coding are used. In contrast, this embodiment describes a case where both Doppler multiplexing transmission and time multiplexing are used.
In this embodiment, for example, when Doppler-multiplexed transmission signals are transmitted in a multiplexed manner also using time multiplexing (hereinafter referred to as “time-Doppler multiplexing”), the numbers of time multiplexes for the respective Doppler-multiplexed transmission signals are set to be non-uniform. This enables radar apparatus to identify transmit antennas and determine the presence or absence of Doppler aliasing based on reception signals of time-Doppler multiplexed signals. As a result, the effective Doppler frequency bandwidth can be extended to 1/(Tr), and the Doppler frequency (relative velocity) detection range with no ambiguity can be extended.
[Configuration of Radar Apparatus]
The following describes, as an example, a configuration of a radar scheme (also referred to as, for example, chirp pulse transmission (fast chirp modulation)) that uses frequency-modulated pulse waves such as chirp pulses. Note that the modulation scheme is not limited to that for frequency modulation. For example, an exemplary embodiment of the present disclosure is also applicable to a radar scheme that uses a pulse compression radar configured to transmit a pulse train after performing phase modulation or amplitude modulation.
Radar apparatus 20 includes radar transmitter (transmitting branch) 300 and radar receiver (receiving branch) 400.
[Configuration of Radar Transmitter 300]
Radar transmitter 300 includes radar transmission signal generator 101, phase rotation amount setter 301, phase rotators 304, transmission controllers 305, and transmit antennas 306.
Phase rotation amount setter 301 sets phase rotation amounts (in other words, phase rotation amounts used for Doppler multiplexing) for phase rotators 304. Phase rotation amount setter 301 includes, for example, Doppler shift setter 302 and time multiplexer 303.
For example, Doppler shift setter 302 sets phase rotation amounts ϕndm corresponding to Doppler shift amounts DOPndm to be applied to radar transmission signals (for example, chirp signals) and outputs phase rotation amounts ϕndm to time multiplexer 303. Here, ndm=1, . . . , NDM. NDM denotes the number of Doppler multiplexes (in other words, different numbers of Doppler shifts).
In Radar apparatus 20, since time multiplexing performed by time multiplexer 303 is also used, the number of Doppler multiplexes NDM may be set to be smaller than the number Nt of transmit antennas 306 used for multiplexing transmission. The number of Doppler multiplexes NDM is greater than or equal to, for example, 2.
Doppler shift amounts DOP1, DOP2, . . . , and DOPDM are assigned different phase rotation amounts by, for example, dividing a phase rotation range greater than or equal to 0 and less than 2π. For example, phase rotation amount ϕndm for applying Doppler shift amount DOPndm is assigned, as given by equation 1 (the angle is expressed in radian).
The assignment of phase rotation amounts for applying Doppler shift amounts DOP1, DOP2, . . . , and DOPDM is not limited to that in this assignment method. For example, the assignment of phase rotation amounts given by equation 1 may be shifted. For example, phase rotation amounts may be assigned such that ϕndm=2π(ndm)/NDM. Alternatively, phase rotation amounts ϕ1, ϕ2, . . . , and ϕDM may be randomly assigned to Doppler shift amounts DOP1, DOP2, . . . , and DOPDM using a phase-rotation-amount assignment table.
Time multiplexer 303 sets, for phase rotation amounts ϕ1, . . . , and ϕNDM for applying NDM Doppler shift amounts output from Doppler shift setter 302, “time-multiplexed Doppler phase rotation amounts” based on the number of time multiplexes NTD, and outputs time-multiplexed Doppler phase rotation amounts” to phase rotators 304.
In time multiplexer 303, the number of time multiplexes (hereinafter referred to as the number of time-Doppler multiplexes) for time-multiplexing a Doppler multiplexed signal using the ndm-th Doppler shift amount DOPndm output from Doppler shift setter 302 is represented by “NDOP_TD(ndm)”. Here, ndm=1, . . . , NDM.
The following describes an example of the operation of time multiplexer 303.
Time multiplexer 303 sets the number of time-Doppler multiplexes NDOP_TD(ndm) so that, for example, the sum of the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), . . . , and NDOP_TD(NDM) to be used to time-multiplex the respective Doppler multiplexed signals is equal to the number Nt of transmit antennas 306 used for multiplexing transmission. In other words, time multiplexer 303 sets the number of time-Doppler multiplexes NDOP_TD(ndm) so as to satisfy following equation 82. This enables radar apparatus 20 to perform multiplexing transmission in the Doppler frequency domain and the time domain (hereinafter referred to as time-Doppler multiplexing transmission) using Nt transmit antennas 306.
Here, time multiplexer 303 sets, for example, the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), . . . , and NDOP_TD(NDM) so as to include different numbers of time-Doppler multiplexes in the range greater than or equal to 1 and less than or equal to NTD. For example, time multiplexer 303 sets the numbers of time-Doppler multiplexes such that not all of the numbers of time-Doppler multiplexes are set to the number of time multiplexes NTD, but at least one of the numbers of time-Doppler multiplexes is set to 1 (in other words, no multiplexing). In other words, time multiplexer 303 sets the numbers of time-Doppler multiplexes for Doppler multiplexed signals to be non-uniform. With this setting, for example, radar apparatus 20 can individually separate and receive signals transmitted from plural transmit antennas 306 in a time-Doppler multiplexed manner through reception processing described below.
Time multiplexer 303 sets, in the m-th transmission period Tr, time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) given by following equation 83 for phase rotation amount ϕndm for applying the ndm-th Doppler shift amount DOPndm and outputs time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) to phase rotators 304.
Here, the subscript “ndop_td(ndm)” represents an index less than or equal to the number of time-Doppler multiplexes NDOP_TD(ndm) for phase rotation amount ϕndm for applying Doppler shift amount DOPndm. For example, ndop_td(ndm)=1, . . . , NDOP_TD(ndm).
For example, as given by equation 83, time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) provides a constant phase rotation amount for applying Doppler shift amount DOPndm for the duration of NTD transmission periods, the number of which is equal to the number of time multiplexes used for time multiplexing (for example, NTD×Tr).
Further, time multiplexer 303 outputs, for each transmission period (Tr), time multiplexing index TD_INDEX to radar receiver 400 (output switching section 401 described below). Further, time multiplexer 303 outputs, for each transmission period Tr, time multiplexing index TD_INDEX and the number of time-Doppler multiplexes NDOP_TD(ndm) to transmission controllers 305.
TD_INDEX represents a time multiplexing index that indicates a transmission period in the duration of NTD transmission periods (NTD×Tr), the number of which is equal to the number of time multiplexes, in other words, indicates a transmission duration or transmission timing for time multiplexing. TD_INDEX cyclically varies in the range of 1 to NTD for each transmission period (Tr), as given by following equation 84.
TD_INDEX=mod(m−1,NTD)+1 (Equation 84)
Here, mod(x, y) denotes a modulo operator and is a function that outputs the remainder of x divided by y. Further, m=1, . . . , Nc. Nc denotes the number of transmission periods used for radar position determination (the number of transmissions of radar transmission signals). The number Nc of transmissions of radar transmission signals is set to an integer multiple (Ntdslot times) of NTD. For example, Nc=NTD×Ntdslot.
Next, an example method for setting the numbers of time-Doppler multiplexes NDOP_TD(ndm) for Doppler multiplexed signals to be non-uniform using time multiplexer 303 will be described.
For example, time multiplexer 303 sets the number of time multiplexes NTD satisfying the condition below. For example, the number of time multiplexes NTD and the number of Doppler multiplexes NDM satisfy the following relationship for the number Nt of transmit antennas 306 used for multiplexing transmission.
(Number of time multiplexes NTD)×(number of Doppler multiplexes NDM)>number Nt of transmit antennas used for multiplexing transmission
For example, among the numbers of time multiplexes NTD and the numbers of Doppler multiplexes NDM satisfying the above-described condition, the use of a combination having a smaller value of the product (NTD×NDM) is desirable in terms of both characteristics and complexity of circuit configuration. Note that among the numbers of time multiplexes NTD and the numbers of Doppler multiplexes NDM satisfying the above-described condition, a combination having a smaller value of the product (NTD×NDM) is not restrictive, and any other combination may be applied.
For example, in a case where Nt=3, the combination of NDM=2, NTD=2 is desirable.
In this case, in the assignment of Doppler shift amounts DOP1 and DOP2 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1) and NDOP_TD(2) are set in two combinations, namely, the combination of NDOP_TD(1)=2 and NDOP_TD(2)=1, and the combination of NDOP_TD(1)=1 and NDOP_TD(2)=2.
Further, for example, in a case where Nt=4, the combination of NDM=3 and NTD=2 or the combination of NDM=2 and NTD=3 is desirable.
For example, in a case where Nt=4, NDM=3, and NTD=2, in the assignment of Doppler shift amounts DOP1, DOP2, and DOP3 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), and NDOP_TD(3) are set in three combinations, namely, the combination of NDOP_TD(1)=2, NDOP_TD(2)=1, and NDOP_TD(3)=1, the combination of NDOP_TD(1)=1, NDOP_TD(2)=2, and NDOP_TD(3)=1, and the combination of NDOP_TD(1)=1, NDOP_TD(2)=1, and NDOP_TD(3)=2.
For example, in a case where Nt=4, NDM=2, and NTD=3, in the assignment of Doppler shift amounts DOP1 and DOP2 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1) and NDOP_TD(2) are set in two combinations, namely, the combination of NDOP_TD(1)=3 and NDOP_TD(2)=1, and the combination of NDOP_TD(1)=1 and NDOP_TD(2)=3.
Further, for example, in a case where Nt=5, the combination of NDM=3 and NTD=2 is desirable.
For example, in a case where Nt=5, NDM=3, and NTD=2, in the assignment of Doppler shift amounts DOP1, DOP2, and DOP3 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), and NDOP_TD(3) are set in three combinations, namely, the combination of NDOP_TD(1)=2, NDOP_TD(2)=2, and NDOP_TD(3)=1, the combination of NDOP_TD(1)=2, NDOP_TD(2)=1, and NDOP_TD(3)=2, and the combination of NDOP_TD(1)=1, NDOP_TD(2)=2, and NDOP_TD(3)=2.
Further, for example, in a case where Nt=6 or 7, the combination of NDM=4 and NTD=2 is desirable.
For example, in a case where Nt=6, NDM=4, and NTD=2, in the assignment of Doppler shift amounts DOP1, DOP2, DOP3, and DOP4 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), NDOP_TD(3), and NDOP_TD(4) are set in following six combinations.
For example, in a case where Nt=7, NDM=4, and NTD=2, in the assignment of Doppler shift amounts DOP1, DOP2, DOP3, and DOP4 and time multiplexes, the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), NDOP_TD(3), and NDOP_TD(4) are set in following four combinations.
Likewise, for example, in a case where Nt=8 or 9, the combination of NDM=5 and NTD=2 is desirable. Further, for example, in a case where Nt=10, the combination of NDM=6 and NTD=2 is desirable. The number Nt of transmit antennas 306 is not limited to that in the example described above, and an exemplary embodiment of the present disclosure is also applicable to Nt=11 or more.
Next, an example of how time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) is set will be described.
For example, a description will be given of a case where in time multiplexer 303, the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, and NTD=2. In this case, for example, if the numbers of time-Doppler multiplexes are set such that NDOP_TD(1)=1 and NDOP_TD (2)=2, time multiplexer 303 sets time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 85 to 87 and outputs time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) to phase rotators 304. Here, m=1, . . . , Nc.
[102]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,ϕ1,ϕ1,2ϕ1,2ϕ1,3ϕ1,3ϕ1, . . . } (Equation 85)
[103]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,ϕ2,ϕ2,2ϕ2,2ϕ2,3ϕ2,3ϕ2, . . . } (Equation 86)
[104]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }={0,0,ϕ2,ϕ2,2ϕ2,2ϕ2,3ϕ2,3ϕ2, . . . } (Equation 87)
In equations 86 and 87, time-multiplexed Doppler phase rotation amounts ψ1, 2(m) and ψ2, 2(m) are the same phase rotation amount, and the corresponding signals are subjected to, for example, an operation of temporally shifting transmission timings (time multiplexing) in transmission controllers 305 described below.
Here, as an example, a description will be given of a case where the phase rotation amount for applying Doppler shift amount DOPndm is given by ϕndm=2π(ndm−1)/NDM in equation 1, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, and phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π, are used. That is, intervals of phase rotation amounts ϕndm for applying Doppler shift amounts DOPndm are equal. In this case, time multiplexer 303 sets coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) given by following equations 88 to 90 and outputs coded Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) to phase rotators 304. Here, m=1, . . . , Nc. Here, a modulo operation for 2π is performed, and results are expressed in radians ranging from 0 or more to less than 2π (the same applies to the following description).
[105]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,0,0,0,0,0,0, . . . } (Equation 88)
[106]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,π,π,0,0,π,π, . . . } (Equation 89)
[107]
{ψ2,2(1),ψ2,2(2),ψ2,2(3),ψ2,2(4),ψ2,2(5),ψ2,2(6),ψ2,2(7),ψ2,2(8), . . . }={0,0,π,π,0,0,π,π, . . . } (Equation 90)
As given by equations 88 to 90, when a phase rotation amount is set to ϕndm=2π(ndm−1)/NDM, into which 2π is equally divided, time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) are changed in transmission periods given by NDM×NTD=2×2=4.
As given by equations 88 to 90, furthermore, the number of phases used for a phase rotation amount for applying a Doppler shift amount (for example, two, namely, 0 and π) is equal to the number of Doppler shift amounts used for multiplexing transmission (in other words, the number of Doppler multiplexes) NDM=2.
Further, for example, a description will be given of a case where in time multiplexer 303, the number of transmit antennas used for multiplexing transmission Nt=6, the number of Doppler multiplexes NOM=4, and NTD=2. In this case, for example, if the numbers of time-Doppler multiplexes are set such that NDOP_TD (1)=1, NDOP_TD (2)=1, NDOP_TD (3)=2, and NDOP_TD (4)=2, time multiplexer 303 sets time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 91 to 96 and outputs time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 304. Here, m=1, . . . , Nc.
[108]
{ψ1,1(1),ψ1,1(2),ψ1,1(3),ψ1,1(4),ψ1,1(5),ψ1,1(6),ψ1,1(7),ψ1,1(8), . . . }={0,0,ϕ1,ϕ1,2ϕ1,2ϕ1,3ϕ1,3ϕ1, . . . } (Equation 91)
[109]
{ψ1,2(1),ψ1,2(2),ψ1,2(3),ψ1,2(4),ψ1,2(5),ψ1,2(6),ψ1,2(7),ψ1,2(8), . . . }={0,0,ϕ2,ϕ2,2ϕ2,2ϕ2,3ϕ2,3ϕ2, . . . } (Equation 92)=
[110]
{ψ1,3(1),ψ1,3(2),ψ1,3(3),ψ1,3(4),ψ1,3(5),ψ1,3(6),ψ1,3(7),ψ1,3(8), . . . }={0,0,ϕ3,ϕ3,2ϕ3,2ϕ3,3ϕ3,3ϕ3, . . . } (Equation 93)
[111]
{ψ2,3(1),ψ2,3(2),ψ2,3(3),ψ2,3(4),ψ2,3(5),ψ2,3(6),ψ2,3(7),ψ2,3(8), . . . }={0,0,ϕ3,ϕ3,2ϕ3,2ϕ3,3ϕ3,3ϕ3, . . . } (Equation 94)
[112]
{ψ1,4(1),ψ1,4(2),ψ1,4(3),ψ1,4(4),ψ1,4(5),ψ1,4(6),ψ1,4(7),ψ1,4(8), . . . }={0,0,ϕ4,ϕ4,2ϕ4,2ϕ4,3ϕ4,3ϕ4, . . . } (Equation 95)
[113]
{ψ2,4(1),ψ2,4(2),ψ2,4(3),ψ2,4(4),ψ2,4(5),ψ2,4(6),ψ2,4(7),ψ2,4(8), . . . }={0,0,ϕ4,ϕ4,2ϕ4,2ϕ4,3ϕ4,3ϕ4, . . . } (Equation 96)
In equations 93 and 94, time-multiplexed Doppler phase rotation amounts ψ1, 3(m) and ψ2, 3(m) are the same phase rotation amount, and the corresponding signals are subjected to an operation of temporally shifting transmission timings (time multiplexing) in transmission controllers 305 described below. Also in equations 95 and 96, time-multiplexed Doppler phase rotation amounts ψ1, 4(m) and ψ2, 4(m) are the same phase rotation amount, and the corresponding signals are subjected to an operation of temporally shifting transmission timings (time multiplexing) in transmission controllers 305 described below.
Here, as an example, the phase rotation amount for applying Doppler shift amount DOPndm is given by ϕndm=2π(ndm−1)/NDM, and phase rotation amount ϕ1 for applying Doppler shift amount DOP1, which is equal to 0, phase rotation amount ϕ2 for applying Doppler shift amount DOP2, which is equal to π/2, phase rotation amount ϕ3 for applying Doppler shift amount DOP3, which is equal to π, and phase rotation amount ϕ4 for applying Doppler shift amount DOP4, which is equal to 3π/2, are used. In this case, time multiplexer 303 sets time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) given by following equations 97 to 102 and outputs time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) to phase rotators 304. Here, m=1, . . . , Nc.
As given by equations 97 to 102, when a phase rotation amount is set to ϕndm=2π(ndm−1)/NDM, into which 2π is equally divided, time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), ψ1, 3(m), ψ2, 3(m), ψ1, 4(m), and ψ2, 4(m) are changed in transmission periods given by NDM×NTD=4×2=8.
As given by equations 97 to 102, furthermore, the number of phases used for a phase rotation amount for applying a Doppler shift amount (for example, four, namely, 0, π/2, π, and 3π/2) is equal to the number of Doppler shift amounts used for multiplexing transmission (in other words, the number of Doppler multiplexes) NDM=4.
While the description has been given of, as an example, the setting of phase rotation amounts in a case where the number Nt of transmit antennas 306 is equal to 3 and the number of Doppler multiplexes NDM is equal to 2 and in a case where the number Nt of transmit antennas 306 is equal to 6 and the number of Doppler multiplexes NDM is equal to 4, the number Nt of transmit antennas 306 and the number of Doppler multiplexes NDM are not limited to the values described above. For example, the number of phases used for a phase rotation amount for applying a Doppler shift amount may be equal to the number NDM of Doppler shift amounts used for multiplexing transmission.
The foregoing description has been given of a method for setting phase rotation amounts using time multiplexer 303.
In
The sum of the numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), . . . , and NDOP_TD(NDM) is set to be equal to the number Nt of transmit antennas 306, and Nt time-multiplexed Doppler phase rotation amounts are respectively input to Nt phase rotators 304.
Each of Nt phase rotators 304 applies time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) input thereto to a chirp signal output from radar transmission signal generator 101 for each transmission period Tr. In the following, phase rotator 304 that applies time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) is represented by “phase rotator PROT #[ndop_td(ndm), ndm]”.
The outputs of Nt phase rotators 304 are output to Nt transmission controllers 305. For example, the output of phase rotator PROT #[ndop_td(ndm), ndm] is output to transmission controller TXCTRL #[ndop_td(ndm), ndm]. Here, ndm=1, . . . , NDM, and ndop_td(ndm)=1, . . . , NDOP_TD(ndm).
In
In the following, whether transmission controller TXCTRL #[ndop_td(ndm), ndm] is in transmission on or transmission off state for each time multiplexing index TD_INDEX is displayed using “transmission controller response STATE_TXCTRL #[ndop_td(ndm), ndm]=[state(1), state(2), . . . , state(NTD)]”.
An element of STATE_TXCTRL #[ndop_td(ndm), ndm], namely, state(TD_INDEX), indicates a transmission on or transmission off state for TD_INDEX=1, . . . , NTD. For example, state(TD_INDEX) may indicate 1 for transmission on, and state(TD_INDEX) may indicate 0 for transmission off.
For example, when the number of time-Doppler multiplexes NDOP_TD(ndm)=1 for a Doppler multiplexed signal using Doppler shift amount DOPndm, ndop_td(ndm)=1, and the number of time multiplexes is 1. Thus, transmission controller TXCTRL #[1, ndm] outputs all transmission signals for time multiplexing indices TD_INDEX=1, . . . , NTD (transmission on). Accordingly, when the number of time-Doppler multiplexes NDOP_TD(ndm)=1 for a Doppler multiplexed signal using Doppler shift amount DOPndm, the transmission controller response is represented by STATE_TXCTRL #[1, ndm]=[1, 1, . . . , 1].
In contrast, for example, when the number of time-Doppler multiplexes NDOP_TD(ndm)=NTD for a Doppler multiplexed signal using Doppler shift amount DOPndm, ndop_td(ndm)=1, . . . , NTD, and the number of time multiplexes is NTD. Thus, transmission controller TXCTRL #[ndop_td(ndm), ndm] outputs a transmission signal when time multiplexing index TD_INDEX=ndop_td(ndm) (transmission on). Accordingly, when the number of time-Doppler multiplexes NDOP_TD(ndm)=NDOP_TD(ndm)=NTD for a Doppler multiplexed signal using Doppler shift amount DOPndm, the transmission controller response STATE_TXCTRL #[ndop_td(ndm), ndm] is 1 (transmission on) when TD_INDEX=ndop_td(ndm), and is 0 (transmission off) when TD_INDEX=ndop_td(ndm).
For example, when the number of time-Doppler multiplexes NDOP_TD(1)=4 (=NTD) for a Doppler multiplexed signal using Doppler shift amount DOP1, the transmission controller response for transmission controller TXCTRL #[1, 1] is represented by STATE_TXCTRL #[1, 1]=[1, 0, 0, 0], the transmission controller response for transmission controller TXCTRL #[2, 1] is represented by STATE_TXCTRL #[2, 1]=[0, 1, 0, 0], the transmission controller response for transmission controller TXCTRL #[3, 1] is represented by STATE_TXCTRL #[3, 1]=[0, 0, 1, 0], and the transmission controller response for transmission controller TXCTRL #[4, 1] is represented by STATE_TXCTRL #[4, 1]=[0, 0, 0, 1].
The outputs of transmission controllers 305 (referred to as, for example, time-Doppler multiplexed signals) are amplified to defined transmission power and are then radiated into a space from Nt transmit antennas 306 in a transmit array antenna section. In other words, phase rotation amounts corresponding to Doppler shift amounts are applied to radar transmission signals, which are then time-multiplexed and transmitted from plural transmit antennas 306.
In the following, transmit antenna 306 that radiates the output of transmission controller TXCTRL #[ndop_td(ndm), ndm] into a space is represented by “transmit antenna Tx #[ndop_td(ndm), ndm]”. Here, ndm=1, . . . , NDM, and ndop_code(ndm)=1, . . . , NDOP_CODE(ndm).
For example, in a case where the number of transmit antennas used for multiplexing transmission Nt=3, the number of Doppler multiplexes NDM=2, NTD=2, and numbers of time-Doppler multiplexes are set such that NDOP_TD(1)=1 and NDOP_TD(2)=2, time-multiplexed Doppler phase rotation amounts ψ1, 1(m), ψ1, 2(m), and ψ2, 2(m) are output from time multiplexer 303 to phase rotators 304 for every transmission period.
For example, phase rotator PROT #[1, 1] applies, for each transmission period, phase rotation amount ψ1, 1(m) to a chirp signal generated in radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[1, 1] is output in a manner given by following expression 103 based on transmission controller response STATE_TXCTRL #[1, 1]=[1,1] for transmission controller TXCTRL #[1, 1], and is transmitted from transmit antenna Tx #[1, 1]. Here, cp(t) denotes a chirp signal for each transmission period.
exp[jψ1,1(1)]cp(t),exp[jψ1,1(2)]cp(t),exp[jψ1,1(3)]cp(t), . . . ,exp[jψ1,1(Nc)]cp(t) (Expression 103)
Likewise, phase rotator PROT #[1, 2] applies, for each transmission period, phase rotation amount ψ1, 2(m) to a chirp signal generated in radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[1, 2] is output in a manner given by following expression 104 based on transmission controller response STATE_TXCTRL #[1, 2]=[1,0] for transmission controller TXCTRL #[1, 2], and is transmitted from transmit antenna Tx #[1, 2].
exp[jψ1,2(1)]cp(t),exp[jψ1,2(2)]cp(t),exp[jψ1,2(3)]cp(t), . . . ,exp[jψ1,2(Nc)]cp(t) (Expression 104)
Likewise, phase rotator PROT #[2, 2] applies, for each transmission period, phase rotation amount ψ2, 2(m) to a chirp signal generated in radar transmission signal generator 101 for each transmission period. The output of phase rotator PROT #[2, 2] is output in a manner given by following expression 105 based on transmission controller response STATE_TXCTRL #[2, 2]=[0,1] for transmission controller TXCTRL #[2, 2], and is transmitted from transmit antenna Tx #[2, 2].
exp[jψ2,2(1)]cp(t),exp[jψ2,2(2)]cp(t),exp[jψ2,2(3)]cp(t), . . . ,exp[jψ2,2(Nc)]cp(t) (Expression 105)
The foregoing description has been given of an example of how time-multiplexed Doppler phase rotation amount ψndop_td(ndm), ndm(m) is set.
In this embodiment, accordingly, plural transmit antennas 306 are associated with combinations (in other words, assignment) of Doppler shift amounts DOPndm and transmission durations for time multiplexing (for example, TD_INDEX) such that in each of the combinations, at least one of Doppler shift amount DOPndm or the transmission duration for time multiplexing (for example, TD_INDEX) is different. In this embodiment, furthermore, the number of time multiplexes (in other words, the number of time-Doppler multiplexes NDOP_TD(ndm)) corresponding to each Doppler shift amount DOPndm in combinations of Doppler shift amounts DOPndm and transmission durations is different.
For example, in this embodiment, Nt transmit antennas 306 include at least plural transmit antennas 306 from which transmission signals that are time-multiplexed are transmitted, and at least one transmit antenna 306 from which a transmission signal that is not time-multiplexed is transmitted. In other words, radar transmission signals transmitted from radar transmitter 300 include at least a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes NDOP_TD(ndm) is set to the number of time multiplexes NTD, and a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes NDOP_TD(ndm) is set to be smaller than the number of codes NTD.
[Configuration of Radar Receiver 400]
In
The z-th signal processor 206 includes NTD Doppler analyzers 402-1 to 402-NTD. For example, data is input to the ntd-th Doppler analyzer 402 in the z-th signal processor 206 for every NTD transmission periods (NTD×Tr) by using output switching section 401. Accordingly, the ntd-th Doppler analyzer 402 performs Doppler analysis for each distance index fb using data of Ntdslot transmission periods among Nc transmission periods. Here, ntd denotes the index for time multiplexing, and ntd=1, . . . , NTD.
For example, when Ntdslot is a power of 2, FFT processing is applicable in Doppler analysis. In this case, the FFT size is Ntdslot, and a maximum Doppler frequency without causing aliasing derived from the sampling theorem is ±1/(2NTD×Tr). Further, the Doppler frequency interval for Doppler frequency index fs is 1/(Ntdslot×NTD×Tr), and the range of Doppler frequency index fs is given by fs=−Ntdslot/2, . . . , 0, . . . , Ntdslot/2−1.
The following description will be given of a case where Ntdslot is a power of 2, as an example. When Ntdslot is not a power of 2, for example, data with zero padding can be used to perform FFT processing with a number of data sizes (FFT sizes) equal to a power of 2. In the FFT processing, Doppler analyzer 402 may perform multiplication by a window function coefficient such as of the Han window or the Hamming window. The application of a window function can suppress sidelobes around the beat frequency peak.
For example, output VFTzntd(fb, fs) of Doppler analyzers 402 in the z-th signal processor 206 is given by following equation 106, where j is the imaginary unit and z=1 to Na.
In
For example, CFAR section 403 performs two-dimensional CFAR processing with the distance axis and the Doppler frequency axis (corresponding to the relative velocity) or CFAR processing that is a combination of one-dimensional CFAR processing operations by power addition of outputs VFTzntd(fb, fs) of Doppler analyzers 402 in the first to Na-th signal processors 206, for example, as given by following equation 107. As the two-dimensional CFAR processing or the CFAR processing that is a combination of one-dimensional CFAR processing operations, for example, processing disclosed in NPL 2 may be applied.
CFAR section 403 adaptively sets a threshold and outputs distance index fb_cfar that provides received power greater than the threshold, Doppler frequency index fs_cfar, and received power information PowerFT(fb_cfar, fs_cfar) to time-Doppler demultiplexer 404.
For example, when phase rotation amount ϕndm for applying Doppler shift amount DOPndm is determined using equation 1, intervals ΔFD of Doppler shift amounts in the Doppler frequency domain, which are output from Doppler analyzers 402, are equal, where ΔFD=Ntdslot/NDM. Accordingly, in the outputs of Doppler analyzers 402, peaks are detected at intervals of ΔFD for Doppler-shift multiplexed signals in the Doppler frequency domain. When phase rotation amount ϕndm is determined using equation 1, ΔFD may not sometimes be an integer depending on Ntdslot and NDM. In this case, following equation 108 may be used to achieve ΔFD having an integer value. The following describes a reception processing operation using ΔFD having an integer value.
Further, CFAR section 403 may divide the respective outputs of Doppler analyzers 402 into ranges at Doppler shift amount intervals ΔFD and perform CFAR processing (for example, Doppler domain compression CFAR processing) after, as given by following equation 109, providing power addition (for example, Doppler domain compression) for the respective ranges while matching peak positions of Doppler-shift multiplexed signals. Here, fs_comp=−Ntdslot/2, . . . , −Ntdslot/2+ΔFD−1.
This can compress the Doppler frequency range for the CFAR processing to 1/NDM to reduce the amount of CFAR processing and can simplify the circuit configuration. In addition, CFAR section 403 enables power addition for NDM Doppler-shift multiplexed signals, resulting in SNR being improved by about (NDM)1/2. As a result, the radar sensing performance of radar apparatus 20 can be improved.
As described above, for example, CFAR section 403, which uses Doppler domain compression CFAR processing, adaptively sets a threshold and outputs distance index fb_cfar that provides received power greater than the threshold, Doppler frequency index fs_comp_cfar, and received power information PowerFT(fb_cfar, fs_comp_cfar+(nfd−1)×ΔFD) for Doppler frequency indices (fs_comp_cfar+(nfd−1)×ΔFD) of NDM Doppler multiplexed signals, where nfd=1, . . . , NDM, to time-Doppler demultiplexer 404.
Note that phase rotation amount ϕndm for applying Doppler shift amount DOPndm is not limited to that given in equation 1. CFAR section 403 can apply Doppler domain compression CFAR processing, for example, if Doppler-shift multiplexed signals have phase rotation amounts ϕndm such that peaks are detected at constant intervals in the Doppler frequency domain output from Doppler analyzers 402.
Next, an example of the operation of time-Doppler demultiplexer 404 illustrated in
Time-Doppler demultiplexer 404 separates signals transmitted in a time-Doppler multiplexed manner, using the output of Doppler analyzer 402, based on the outputs of CFAR section 403, namely, distance index fb_cfar, Doppler frequency index fs_comp_cfar, and received-power information PowerFT(fb_cfar, fs_comp_cfar±(nfd−1)×ΔFD) for Doppler frequency indices (fs_comp_cfar±(nfd−1)×ΔFD) of NDM Doppler multiplexed signals, where nfd=1, . . . , NDM, and performs discrimination of transmit antennas 306 and discrimination of Doppler frequencies (or Doppler velocities or relative velocities).
As described above, for example, time multiplexer 303 in phase rotation amount setter 301 does not set all of NDM numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), . . . , and NDOP_TD(NDM) to NTD, but sets at least one of the numbers of Doppler multiplexes to 1. For example, time-Doppler demultiplexer 404 detects (1) a Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1, and performs (2) discrimination of transmit antennas 306 and discrimination of Doppler frequencies of the target based on the detected Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1.
The following describes the processes (1) and (2) described above, which are performed by time-Doppler demultiplexer 404.
<(1) Process of Detecting Time-Doppler Multiplexed Signal for which Number of Time-Doppler Multiplexes is Set to 1>
There are NDM candidate correspondences between NDM time-Doppler multiplexed signals and the outputs of the respective Doppler analyzers 402 for Doppler frequency indices (fs_comp_cfar), (fs_comp_cfar+ΔFD), (fs_comp_cfar+2ΔFD), . . . , and (fs_comp_cfar+(NDM−1)ΔFD) of NDM time-Doppler multiplexed signals for distance index fb_cfar output from CFAR section 403.
For example, if Doppler shift amount DOPndm set in Doppler shift setter 302 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, there are NDM candidate correspondences with cyclically shifted elements, as follows, with consideration given to Doppler aliasing. Here, patterns of the candidate correspondences are numbered DopCase=1 to NDM.
For example, DopCase=1 indicates a correspondence among Doppler shift amounts in the initial state (when the relative velocity to the target is zero). For example, more aliasing components are included as the relative velocity of the target increases in a direction in which the distance to the target decreases, and the resulting correspondences are associated with DopCase=2, 3, . . . , NDM. In other words, DopCase=NDM, NDM-1, . . . , 2 is associated as the relative velocity of the target increases in a direction in which the distance to the target increases.
Here, a table indicating the position of DOPndm counting from the beginning of each DopCase (the position (or order) of DOPndm in DopCase) can be created in advance, based on the Doppler shift amounts set in Doppler shift setter 302. In the following, DOPposi(DOPndm, DopCase) denotes the operator that outputs the position of DOPndm counting from the beginning of each DopCase. For example, in the above-described example of DopCase, DOPposi(DOP1, 1)=1, DOPposi(DOP1, 2)=2, DOPposi(DOP1, NDM)=NDM, DOPposi(DOP2, 1)=2, DOPposi(DOP2, 2)=3, and DOPposi(DOP2, NDM)=1.
Time-Doppler demultiplexer 404 performs a detection process on, for example, the outputs of Doppler analyzers 402 in the z-th signal processor 206 indicated by Doppler frequency indices (fs_comp_cfar+(nfd−1)×ΔFD) of NDM time-Doppler multiplexed signals for distance indices fb_cfar output from CFAR section 403 to detect Doppler multiplexed signals for which the numbers of time-Doppler multiplexes are set to 1.
Here, time-Doppler multiplexed signals for which the numbers of time-Doppler multiplexes are set to 1 are transmitted using the same transmit antenna 306, and time-Doppler multiplexed signals for which the numbers of time-Doppler multiplexes are set to NTD are transmitted using different transmit antennas 306.
Time-Doppler demultiplexer 404 performs coherent addition processing on, for example, the outputs of Doppler analyzers 402 for each time multiplexing index TD_INDEX using DopCase as a parameter, with consideration given to the presence or absence of aliasing. Time-Doppler demultiplexer 404 utilizes, for example, DopCase for which a coherent addition value that realizes maximum power is obtained, as the most likely pattern of Doppler shift amounts to identify transmit antennas.
The detection process may be performed on all of nfd=1, . . . , NDM for all the candidates of DopCase=1, . . . , NDM. Note that the position of, in each DopCase, a Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 is set by Doppler shift setter 302 and is known to radar apparatus 20. Accordingly, time-Doppler demultiplexer 404 performs the following detection process to reduce the operation amount of the separation process.
For example, when a time-Doppler multiplexed signal for which the number of time Doppler multiplexes is set to 1 is Doppler shift amount DOPndm1, time-Doppler demultiplexer 404 performs coherent addition processing on the outputs of Doppler analyzers 402 using candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1 within each DopCase, with consideration given to the presence or absence of Doppler aliasing. Then, time-Doppler demultiplexer 404 detects DopCase that realizes maximum power as follows, based on the result of coherent addition processing with consideration given to the presence or absence of Doppler aliasing.
(a) Coherent Addition Processing when Outputs of Doppler Analyzers 402 Contain No Doppler Aliasing
Coherent addition signal ADDCOHz(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) given by following equation 110 represents a coherent addition value for the output of Doppler analyzer 402 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) in the z-th signal processor 206.
In equation 110, in the exp term, since the sampling times for the outputs of Doppler analyzers 402 for each time multiplexing index are shifted, phase correction is performed in accordance with Doppler frequency index (fs_comp_cfar±(DopCase−1)×ΔFD-DOP1). In equation 110, furthermore, DopCase=1, . . . , NDM. Further, Doppler shift amount DOPndm set in Doppler shift setter 302 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, and DOP1 falls within the range of fs_comp=−Ncode/2, . . . , −Ncode/2+ΔFD−1 in the initial state (when the relative velocity to the target is zero). Accordingly, time-Doppler demultiplexer 404 calculates an amount of phase correction using, for example, DOP1 as a reference in equation 110.
When a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 is Doppler shift amount “DOPndm1”, for DopCase=1, . . . , NDM, coherent addition signals for candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1 are obtained using equation 110. Accordingly, time-Doppler demultiplexer 404 obtains the outputs of a total of NDM coherent addition signals. For example, time-Doppler demultiplexer 404 calculates coherent addition signals for all the receive antennas z=1, . . . , Na in accordance with equation 110 and calculates time separation signal power sum Pow_ADDCOH(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD), as given by following equation 111.
(b) Coherent Addition Processing when Outputs of Doppler Analyzers 402 Contain Doppler Aliasing
If the outputs of Doppler analyzers 402 contain Doppler aliasing, time-Doppler demultiplexer 404 performs coherent addition processing using phase correction (the exp term) used in the coherent addition processing as phase correction with consideration given to Doppler aliasing.
For example, in a case where NTD=2, coherent addition signal ADDALIASz(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) given by following equation 112 represents a coherent addition value with consideration given to Doppler aliasing for the output of Doppler analyzer 402 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) in the z-th signal processor 206.
In equation 112, in the exp term, since the sampling times for the outputs of Doppler analyzers 402 for each time multiplexing index are shifted, phase correction is performed in accordance with Doppler frequency index (fs_comp_cfar±(DopCase−1)×ΔFD-DOP1) and phase correction (exp{−j2π(ntd−1)/NTD}) to address aliasing is performed. In equation 112, furthermore, DopCase=1, . . . , NDM. Further, Doppler shift amount DOPndm set in Doppler shift setter 302 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, and DOP1 falls within the range of fs_comp=−Ncode/2, . . . , −Ncode/2+ΔFD−1 in the initial state (when the relative velocity to the target is zero). Accordingly, time-Doppler demultiplexer 404 calculates an amount of phase correction using DOP1 as a reference in equation 112.
When a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 is Doppler shift amount DOPndm1, for DopCase=1, . . . , NDM, coherent addition signals with consideration given to Doppler aliasing for candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1 are obtained using Equation 112. Accordingly, time-Doppler demultiplexer 404 obtains the outputs of a total of NDM coherent addition signals. For example, time-Doppler demultiplexer 404 calculates coherent addition signals for all the receive antennas z=1, Na using equation 112 and calculates coherent addition power sum Pow_ADDALIAS(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) with consideration given to Doppler aliasing, as given by following equation 113.
For example, in a case where NTD=3, coherent addition signal ADDALIASz(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) given by following equations 114 and 115 represents a coherent addition value with consideration given to Doppler aliasing for the output of Doppler analyzer 402 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) in the z-th signal processor 206.
In equation 114, in the exp term, since the sampling times for the outputs of Doppler analyzers 402 for each time multiplexing index are shifted, phase correction is performed in accordance with Doppler frequency index (fs_comp_cfar±(DopCase−1)×ΔFD−DOP1) and phase correction (exp{−j2πSign(fest)(ntd−1)/NTD}) to address aliasing is performed. In equation 114, furthermore, DopCase=1, . . . , NDM. Further, Doppler shift amount DOPndm set in Doppler shift setter 302 is represented by DOP1<DOP2< . . . <DOPDM-1<DOPDM, and DOP1 falls within the range of fs_comp=−Ncode/2, . . . , −Ncode/2+ΔFD−1 in the initial state (when the relative velocity to the target is zero). Accordingly, time-Doppler demultiplexer 404 calculates an amount of phase correction using DOP1 as a reference in equation 114.
When a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 is Doppler shift amount DOPndm1, for DopCase=1, . . . , NDM, coherent addition signals with consideration given to Doppler aliasing for candidate DOPposi(DOPndm1, DopCase) including Doppler shift amount DOPndm1 are obtained using Equation 114. Accordingly, time-Doppler demultiplexer 404 obtains the outputs of a total of NDM coherent addition signals. For example, time-Doppler demultiplexer 404 calculates coherent addition signals for all the receive antennas z=1, . . . , Na using equation 114 and calculates coherent addition power sum Pow_ADDALIAS(fb-efar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) with consideration given to Doppler aliasing, as given by following equation 116.
In a case where NTD=3, coherent addition values with consideration also given to multiplexed aliasing determination are calculated.
Also in a case where NTD>3, coherent addition value Pow_ADDCOH(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) and coherent addition value Pow_ADDALIAS(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) with consideration given to Doppler aliasing are calculated.
As described above, time-Doppler demultiplexer 404 performs coherent addition processing on the outputs of Doppler analyzers 402, with consideration given to the presence or absence of Doppler aliasing.
Time-Doppler demultiplexer 404 detects coherent addition values that provide maximum power for DopCase=1, . . . , NDM, based on coherent addition value Pow_ADDCOH(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase)−1)×ΔFD) and coherent addition value Pow_ADDALIAS(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase)−1)×ΔFD) with consideration given to Doppler aliasing. In the following, the index number of DopCase that provides a coherent addition value that realizes maximum power is represented by “DopCase_max”.
For example, when coherent addition value Pow_ADDCOH(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) indicates maximum power, time-Doppler demultiplexer 404 determines “absence of aliasing”.
For example, when coherent addition value Pow_ADDALIAS(fb_cfar, fs_comp_cfar±(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) with consideration given to Doppler aliasing indicates maximum power, time-Doppler demultiplexer 404 determines “presence of aliasing”.
A reduction in the SNR of a reception signal that is a reflected wave may make it difficult to discriminate detection of the maximum coherent addition power level due to the influence of noise. To address this difficulty, time-Doppler demultiplexer 404 may introduce a determination condition and perform a process of, for example, adopting a determination result (in other words, detection result) when the determination condition is satisfied and removing (in other words, not adopting) a determination result when the determination condition is not satisfied. This can reduce the probability of erroneous detection of the noise component and the like. For example, time-Doppler demultiplexer 404 may adopt, as a determination result, detection value PMAX of received power having maximum coherent addition value indicating the presence or absence of Doppler aliasing when satisfying PMAX>LEVDETECT×PowerFT(fb_cfar, fs_comp_cfar) LEVDETECT is a determination threshold. LEVDETECT is a real number satisfying 0<LEVDE0TECT<1.
<(2) Process for Discrimination of Transmit Antennas 306 and Discrimination of Doppler Frequencies of Target>
For example, time-Doppler demultiplexer 404 detects Doppler frequency fTARGET of the target in the range of −1/(2×Tr)≤fTARGET<1/(2×Tr).
For example, time-Doppler demultiplexer 404 discriminates transmit antennas 306 and Doppler frequencies of the target, based on the determination result of DopCase (for example, DopCase_max) and the determination result of the presence or absence of aliasing.
[Case (a): Case without Aliasing]
Doppler Frequency Determination for Target:
Time-Doppler demultiplexer 404 determines that the Doppler frequency index of the target is (fs_comp_cfar+(DopCase_max−1)×ΔFD−DOP1). For example, the Doppler frequency interval for Doppler frequency index fs_comp_cfar is given by 1/(Ntdslot×NTD×Tr). Therefore, time-Doppler demultiplexer 404 determines that Doppler frequency fTARGET of the target is (fs_comp_cfar (DopCase_max−1)×ΔFD−DOP1)/(Ntdslot×NTD×Tr).
Transmit Antenna Determination:
Time-Doppler demultiplexer 404 determines, for transmit antenna 306 with the number of time-Doppler multiplexes being 1, that, for example, coherent addition signal ADDCOHz(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) in the z-th signal processor 206 is a reception signal for transmit antenna Tx #[1, ndm1].
Further, time-Doppler demultiplexer 404 determines, for transmit antenna 306 with the number of time-Doppler multiplexes being NTD, that the output of Doppler analyzer 402 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm, DopCase_max)−1)×ΔFD) (see, for example, following equations 117 and 118) in the z-th signal processor 206 is a reception signal for transmit antenna Tx #[ntd, ndm]. In equation 117, the exp term is a term for correcting the Doppler phase caused by time multiplexing. Further, ndm=1, . . . , NDM, but excluding ndm1. Further, ntd=1, . . . , NTD.
[Case (b): Case with Abasing]
Doppler Frequency Determination for Target:
Time-Doppler demultiplexer 404 determines that the Doppler frequency index of the target is fs_comp_cfar (DopCase_max−1)×ΔFD)−DOP1−Ncode×Sign(fs_comp_cfar (DopCase_max−1)×ΔFD−DOP1). For example, the Doppler frequency interval for Doppler frequency index fs_comp_cfar is given by 1/(Ntdslot×NTD×Tr). Therefore, time-Doppler demultiplexer 404 determines that Doppler frequency fTARGET of the target is (fs_comp_cfar+(DopCase_max−1)×ΔFD−DOP1−Ntdslot×Sign(fs_comp_cfar+(DopCase_max−1)×ΔFD−DOP1))/(Ntdslot×NTD×Tr). Sign(x) is a sign function and is a function for real numbers x, providing as output 1 when x>0, 0 when x=0, and −1 when x<0.
In this manner, when the outputs of Doppler analyzers 402 contain Doppler aliasing, time-Doppler demultiplexer 404 determines a Doppler frequency of the target, with consideration given to the aliasing component (for example, Ntdslot×Sign(fs_comp_cfar+(DopCase_max−1)×ΔFD−DOP1)).
Transmit Antenna Determination:
Time-Doppler demultiplexer 404 determines, for transmit antenna 306 with the number of time-Doppler multiplexes being 1, that, for example, coherent addition signal ADDALIASz(fb_cfar, fs_comp_cfar+(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) with consideration given to aliasing for distance index fb_cfar and Doppler frequency index (fs_comp_cfar+(DOPposi(DOPndm1, DopCase_max)−1)×ΔFD) in the z-th signal processor 206 is a reception signal for transmit antenna Tx #[1, ndm1].
Further, time-Doppler demultiplexer 404 determines, for transmit antenna 306 with the number of time-Doppler multiplexes being NTD, that, for example, the output of Doppler analyzer 402 for distance index fb_cfar and Doppler frequency index (fs_comp_cfar±(DOPposi(DOPndm, DopCase_max)−1)×ΔFD) (see, for example, following equations 119 and 120) in the z-th signal processor 206 is a reception signal for transmit antenna Tx #[ntd, ndm]. In equation 119, the exp term is a term for correcting the Doppler phase caused by time multiplexing. Further, ndm=1, . . . , NDM, but excluding ndm1. Further, ntd=1, . . . , NTD.
As described above, radar apparatus 20 does not set all of NDM numbers of time-Doppler multiplexes NDOP_TD(1), NDOP_TD(2), . . . , and NDOP_TD(NDM) to NTD, but sets at least one of the numbers of time-Doppler multiplexes to 1 (in other words, no time multiplexing). With this setting, for example, a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to NTD is transmitted at a set transmission timing (for example, timing of transmission on), whereas a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 is transmitted for each transmission period, for example. Therefore, a reception signal of a time-Doppler multiplexed signal for which the number of time-Doppler multiplexes is set to 1 has a maximum value for a coherent addition value with consideration given to the presence or absence of Doppler aliasing. Time-Doppler demultiplexer 404 can utilize this characteristic to discriminate transmit antennas 306 and to determine a Doppler frequency of the target including the presence or absence of aliasing.
Accordingly, in this embodiment, the range over which a Doppler frequency is detectable without ambiguity can be extended to the range greater than or equal to −1/(2Tr) and less than 1/(2Tr). For example, when a single antenna is used for transmission, the range over which a Doppler frequency is detectable without ambiguity is the range greater than or equal to −1/(2Tr) and less than 1/(2Tr). In this embodiment, therefore, even when plural transmit antennas 306 are used, radar apparatus 20 can detect a Doppler frequency without ambiguity in a manner similar to that when a single antenna is used for transmission.
The foregoing description has been given of an example of the operation of time-Doppler demultiplexer 404.
In
For example, direction estimator 405 generates virtual receive array correlation vector h(fb_cfar, fs_comp_cfar), based on the Doppler frequency determination result for the target and the transmit antenna determination result, and performs a direction estimation process.
As described above, in this embodiment, radar apparatus 20 applies phase rotation amounts corresponding to Doppler shift amounts to radar transmission signals and time-multiplexes radar transmission signals (for example, Doppler multiplexed signals) to transmit radar transmission signals (for example, time-Doppler multiplexed signals) from plural transmit antennas 108 in a multiplexed manner. In this embodiment, plural transmit antennas 108 are associated with combinations of Doppler shift amounts (DOPndm) and transmission durations for time multiplexing (for example, transmission timings) such that in each of the combinations, at least one of Doppler shift amount (DOPndm) or the transmission duration for time multiplexing (for example, transmission timing) is different. In this embodiment, furthermore, the number of time multiplexes corresponding to each Doppler shift amount in combinations of Doppler shift amounts and transmission durations for time multiplexing is different. In other words, the numbers of time multiplexes for the respective Doppler-multiplexed transmission signals are set to be non-uniform.
Radar apparatus 20 can determine, based on, for example, a reception signal of each time-Doppler multiplexed signal (for example, coherent addition signal), transmit antenna 306 associated with the time-Doppler multiplexed signal (in other words, the combination of Doppler shift amount and transmission duration) and the presence or absence of Doppler aliasing (for example, DopCase and the like). This enables radar apparatus 20 to appropriately determine a Doppler frequency of the target even in the presence of Doppler aliasing.
According to this embodiment, therefore, radar apparatus 20 can extend the effective Doppler frequency bandwidth to 1/(Tr) and can extend the Doppler frequency (relative velocity) detection range with no ambiguity. Accordingly, radar apparatus 20 can improve target-object sensing accuracy over a wider Doppler frequency range.
In this embodiment, furthermore, time-Doppler multiplexing, which is performed using both Doppler multiplexing and time multiplexing, can reduce the number of Doppler multiplexes compared with the use of only Doppler multiplexing in multiplexing transmission. This can increase the intervals of phase rotation amounts for applying Doppler shifts, thereby, for example, relieving the accuracy requirements (phase modulation accuracy) for the phase shifters and achieving the cost reduction effect of an RF section, including reduction of the man-hours required for the adjustment of the phase shifters.
In this embodiment, furthermore, since time-Doppler multiplexing is performed using both Doppler multiplexing and time multiplexing, radar apparatus 20 performs, for each time multiplexing index, Fourier frequency analysis (FFT processing) to detect the Doppler frequency (detect the relative velocity). Accordingly, for example, compared with Fourier frequency analysis (FFT processing) to detect the Doppler frequency (detect the relative velocity) using only Doppler multiplexing in multiplexing transmission, the FFT size is (1/number of time multiplexes NTD) and the number of FFTs is increased by (the number of time multiplexes NTD) times. According to this embodiment, therefore, the operation reduction effect of FFT processing can be achieved, and the effect of simplification of the circuit configuration and cost reduction can also be achieved.
The foregoing description has been given of an exemplary embodiment of the present disclosure.
For example, Embodiment 1 has described the case where in an example in which the number of code multiplexes is set to be smaller than NCM, the number of coded Doppler multiplexes is 1 (in other words, no code multiplexing). However, in an exemplary embodiment of the present disclosure, for example, when the number of coded Doppler multiplexes is set to be smaller than NCM, the number of coded Doppler multiplexes may be set to be greater than or equal to 2.
In a radar apparatus according to an exemplary embodiment of the present disclosure, a radar transmitter and a radar receiver may be individually arranged in physically separate locations. In a radar receiver according to an exemplary embodiment of the present disclosure, a direction estimator and any other component may be individually arranged in physically separate locations.
Further, in an exemplary embodiment of the present disclosure, the values used (such as the number Nt of transmit antennas, the number Na of receive antennas, the number of Doppler multiplexes NDM, the number of codes NCM, and the number of time multiplexes NTD) are examples, and these values are not restrictive.
A radar apparatus according to an exemplary embodiment of the present disclosure includes, for example, a central processing unit (CPU), a storage medium such as a read only memory (ROM) that stores a control program, and a work memory such as a random access memory (RAM), which are not illustrated. In this case, the functions of the sections described above are implemented by the CPU executing the control program. However, the hardware configuration of the radar apparatus is not limited to that in this example. For example, the functional sections of the radar apparatus may be implemented as an integrated circuit (IC). Each functional section may be formed as an individual chip, or some or all of them may be formed into a single chip. In the foregoing description, the term “section” used to indicate each constituent element may be interchangeably referred to as any other term such as “circuit (circuitry)”, “device”, “unit”, or “module”.
Various embodiments have been described with reference to the drawings hereinabove. Obviously, the present disclosure is not limited to these examples. Obviously, a person skilled in the art would arrive variations and modification examples within a scope described in claims, and it is understood that these variations and modifications are within the technical scope of the present disclosure. Each constituent element of the above-mentioned embodiments may be combined optionally without departing from the spirit of the disclosure.
The above embodiments have been described with an example of a configuration using hardware, but the present disclosure can be realized by software in cooperation with hardware.
Each functional block used in the description of each embodiment described above is typically realized by an LSI, which is an integrated circuit. The integrated circuit controls each functional block used in the description of the above embodiments and may include an input terminal and an output terminal. The LSI may be individually formed as chips, or one chip may be formed so as to include a part or all of the functional blocks. The LSI herein may be referred to as an IC, a system LSI, a super LSI, or an ultra LSI depending on a difference in the degree of integration.
However, the technique of implementing an integrated circuit is not limited to the LSI and may be realized by using a dedicated circuit, a general-purpose processor, or a special-purpose processor. In addition, a Field Programmable Gate Array (FPGA) that can be programmed after the manufacture of the LSI or a reconfigurable processor in which the connections and the settings of circuit cells disposed inside the LSI can be reconfigured may be used.
If future integrated circuit technology replaces LSIs as a result of the advancement of semiconductor technology or other derivative technology, the functional blocks could be integrated using the future integrated circuit technology. Biotechnology can also be applied.
A radar apparatus according to an exemplary embodiment of the present disclosure includes a plurality of transmit antennas that transmit a plurality of transmission signals using a multiplexing transmission, and a transmission circuit that applies phase rotation amounts corresponding to combinations of Doppler shift amounts and code sequences to the plurality of transmission signals, wherein each of the combinations of the Doppler shift amounts and the code sequences has at least one different from other combination, and wherein the number of multiplexes of the code sequence corresponding to at least one of the Doppler shift amounts in the combinations is different from the number of multiplexing of code sequences corresponding to the remaining Doppler shift amount.
In an exemplary embodiment of the present disclosure, the number of phases used for phase rotation amounts for applying the Doppler shift amounts is smaller than the number of the plurality of transmit antennas.
In an exemplary embodiment of the present disclosure, intervals of phase rotation amounts for applying the Doppler shift amounts are equal.
In an exemplary embodiment of the present disclosure, the number of phases used for phase rotation amounts for applying the Doppler shift amounts is equal to the number of the Doppler shift amounts that are used for the multiplexing transmission.
In an exemplary embodiment of the present disclosure, the radar apparatus further includes a plurality of receive antennas that receive reflected wave signals, the reflected wave signals being the plurality of transmission signals reflected from a target, and reception circuit that determines a transmit antenna corresponding to each of the reflected wave signals among the plurality of transmit antennas and determines whether the reflected wave signal contains aliasing in a Doppler frequency domain, based on a signal subjected to code separation for the reflected wave signal using a corresponding one of the code sequences.
In an exemplary embodiment of the present disclosure, the transmit antennas have a sub-array configuration.
In an exemplary embodiment of the present disclosure, the transmission circuit sets intervals of the Doppler shift amounts to vary for each of frames in which the plurality of transmission signals are transmitted.
In an exemplary embodiment of the present disclosure, the transmission circuit sets transmission periods of the plurality of transmission signals to vary for each of frames in which the plurality of transmission signals are transmitted.
In an exemplary embodiment of the present disclosure, the transmission circuit multiplies each of the plurality of transmission signals by a pseudo-random code sequence.
In an exemplary embodiment of the present disclosure, the transmission circuit sets an association between each of the plurality of transmit antennas and one of the combinations of the Doppler shift amounts and the code sequences to vary for each of frames in which the plurality of transmission signals are transmitted.
In an exemplary embodiment of the present disclosure, the transmission circuit sets an association between each of the plurality of transmit antennas and one of the Doppler shift amounts to vary for each of transmission periods of the plurality of transmission signals.
In an exemplary embodiment of the present disclosure, a radar apparatus includes a plurality of transmit antennas that transmit transmission signals, and a transmission circuit that applies phase rotation amounts to the transmission signals to perform multiplexing transmission to transmit the transmission signals from the plurality of transmit antennas in a multiplexed manner, wherein the number of phases used for the phase rotation amounts is smaller than the number of the plurality of transmit antennas.
In an exemplary embodiment of the present disclosure, the number of the plurality of transmit antennas is any of 3, 6, and 7.
In an exemplary embodiment of the present disclosure, a radar apparatus includes a plurality of transmit antennas that transmit transmission signals, and a transmission circuit that applies phase rotation amounts corresponding to Doppler shift amounts to the transmission signals and time-multiplexes the transmission signals to perform multiplexing transmission to transmit the transmission signals from the plurality of transmit antennas in a multiplexed manner, wherein the plurality of transmit antennas are associated with combinations of the Doppler shift amounts and transmission durations for the time multiplexing such that in each of the combinations, at least one of the Doppler shift amount or the transmission duration for the time multiplexing is different, and the number of multiplexes for the time multiplexing corresponding to each of the Doppler shift amounts in the combinations is different.
In an exemplary embodiment of the present disclosure, the number of phases used for phase rotation amounts for applying the Doppler shift amounts is equal to the number of the Doppler shift amounts that are used for the multiplexing transmission.
In an exemplary embodiment of the present disclosure, intervals of phase rotation amounts for applying the Doppler shift amounts are equal.
A radar apparatus according to an exemplary embodiment of the present disclosure includes a plurality of transmit antennas that transmit transmission signals, and a transmission circuit that applies phase rotation amounts corresponding to Doppler shift amounts and code sequences to the transmission signals to perform multiplexing transmission to transmit the transmission signals from the plurality of transmit antennas in a multiplexed manner, wherein each of the transmission signals is a chirp signal, and a center frequency of the chirp signal is changed for each of transmission periods of the transmission signals or each of transmission periods of the code sequences.
A radar apparatus according to an exemplary embodiment of the present disclosure includes a plurality of transmit antennas that transmit transmission signals, and a transmission circuit that applies phase rotation amounts corresponding to Doppler shift amounts and code sequences to the transmission signals to perform multiplexing transmission to transmit the transmission signals from the plurality of transmit antennas in a multiplexed manner, wherein the transmission signals are chirp signals, and center frequencies of the chirp signals are changed in a round in a period that is a divisor multiple of a code length of each of the code sequences relative to each of transmission periods of the transmission signals.
While various embodiments have been described herein above, it is to be appreciated that various changes in form and detail may be made without departing from the sprit and scope of the invention(s) presently or hereafter claimed.
This application is entitled and claims the benefit of Japanese Patent Application No. 2019-110522, filed on Jun. 13, 2019 and Japanese Patent Application No. 2019-221249, filed on Dec. 6, 2019, the disclosure of which including the specification, drawings and abstract is incorporated herein by reference in its entirety.
The present disclosure is suitable as a radar apparatus for wide-angle range sensing.
Number | Date | Country | Kind |
---|---|---|---|
2019-110522 | Jun 2019 | JP | national |
2019-221249 | Dec 2019 | JP | national |
This is a continuation of U.S. patent application Ser. No. 16/898,856, filed Jun. 11, 2020, which claims the benefit of priority of Japanese Pat. Appl. No. 2019-110522, filed Jun. 13, 2019, and of Japanese Pat. Appl. No. 2019-221249, filed Dec. 6, 2019. The entire disclosure of each of the above-identified documents is incorporated herein by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
20080303711 | Matsuoka | Dec 2008 | A1 |
20110074620 | Wintermantel | Mar 2011 | A1 |
20130135140 | Kishigami | May 2013 | A1 |
20170254879 | Tokieda, I | Sep 2017 | A1 |
20190170873 | Reid | Jun 2019 | A1 |
20190357204 | Matsumura | Nov 2019 | A1 |
20200292663 | Bai | Sep 2020 | A1 |
20210377095 | Nayebi et al. | Dec 2021 | A1 |
20220052715 | Ma et al. | Feb 2022 | A1 |
20220107402 | Kishigami | Apr 2022 | A1 |
Number | Date | Country |
---|---|---|
2828186 | Aug 2012 | CA |
2008-304417 | Dec 2008 | JP |
2011-526371 | Oct 2011 | JP |
2014-119344 | Jun 2014 | JP |
2017-177539 | Oct 2017 | JP |
2016031882 | Mar 2016 | WO |
Entry |
---|
Cadzow, “Direction-of-Arrival Estimation Using Signal Subspace Modeling”, IEEE Transactions on Aerospace and Electronic Systems, vol. 28, No. 1, Jan. 1992, pp. 64-79. |
Kronauge et al., “Fast Two-Dimensional CFAR Procedure”, IEEE Transactions on Aerospace and Electronic Systems, vol. 49, No. 3, Jul. 2013, pp. 1817-1823. |
Li et al., “MIMO Radar with Colocated Antennas”, IEEE Signal Processing Magazine, vol. 24, Issue: 5, Sep. 2007, pp. 106-114. |
Winkler, “Novel Waveform Generation Principle for short-range FMCW-Radars”, In Proc. German Microw. Conf., 2009, pp. 1-4. |
English language translation of Office Action from Japan Patent Office (JPO) in Japanese Patent Appl. No. 2019-221249, dated Mar. 7, 2023. |
Number | Date | Country | |
---|---|---|---|
20220381903 A1 | Dec 2022 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 16898856 | Jun 2020 | US |
Child | 17870383 | US |