1. Field of the Invention
The present invention relates to a radar device such as a frequency modulated continuous wave (FMCW) radar device used for anticollision of a moving object such as an automobile, preceding-vehicle following travel while keeping a constant inter-vehicle range, or the like, for detecting a relative velocity or a range with respect to a target existing outside the moving object through transmission and reception of a radar wave.
2. Description of the Related Art
As illustrated in
where B is a frequency displacement width of the transmission signal S1, f0 is a center frequency of the transmission signal S1, T is a period of time required for modulation in one period, and C is a light speed.
As described above, the conventional radar device is capable of detecting the range to the target and a range change rate through association (hereinafter, referred to as “pairing”) of the beat frequencies in the up-chirp state and the down-chirp state with each other. The beat frequencies obtained in the up-chirp state and the down-chirp state, respectively, are offset to each other even if the beat frequencies are obtained for the same target. In particular, under the environment in which there exist a plurality of targets, that is, a plurality of beat frequencies, there is a need to determine which beat frequency in the up-chirp state corresponds to the beat frequency in the down-chirp state, and the determination is extremely difficult.
As a countermeasure against the above-mentioned problem, there has been proposed the following FMCW radar device (for example, see Japanese Patent No. 2778864). The conventional FMCW radar device is capable of coping with the environment in which the plurality of targets exist, in such a manner that, in pairing of the beat frequencies obtained in the up-chirp state and the down-chirp state, pairing of the beat frequencies obtained in the up-chirp state and the down-chirp state is performed so that the beat frequencies obtained in each sweep period are arranged in ascending order, and the arrangement is saved.
However, the above-mentioned conventional technology requires pairing of the beat frequencies obtained in the up-chirp state and the down-chirp state with an aim to obtain the range and the velocity with respect to the external target. For that reason, when the frequency in one of those states is not obtained, there appears a target (undetectable target) which cannot be detected because the frequency pair cannot be selected though the target actually exists, or a target (false target) which would not originally exist. This causes reliability of measurement results to deteriorate.
To solve the above-mentioned drawbacks, there has been proposed the following method (for example, see Japanese Patent No. 4186744). In the proposed method, when a given target is first observed, the beat frequencies obtained in the up-chirp state and the down-chirp state are paired with each other to obtain the range to the target and the range change rate. In second and subsequent observations, the range and the range change rate which have been first obtained, and the beat frequency in the up-chirp state or the down-chirp state are directly employed to calculate the range to the target and the range change rate.
However, the conventional technology suffers from the following problems. The conventional radar device disclosed in Japanese Patent No. 4186744 does not require pairing of the beat frequencies, but updates the range and the range change rate without consideration of input of an angle. Therefore, when the angle of the target is changed, the probability that a precision in tracking the target position deteriorates is high. Further, a method of removing an erroneous pair in pairing of the beat frequencies for use in a normal radar device is not taken.
The present invention has been made to solve the above-mentioned problems, and therefore aims at providing a radar device which is capable of improving a precision in tracking a target even if the number of peaks is different in beat frequency between an up-chirp state and a down-chirp state, and of removing an erroneous pair, by using a tracking filter that receives the beat frequencies together with a tracking filter that receives an observation value after pairing.
In particular, according to the conventional technologies, there has been proposed the method of updating the range and the range change rate with a direct input of the beat frequencies in the up-chirp state and the down-chirp state. On the other hand, the present invention proposes a method of estimating the target position accurately to remove the erroneous pair.
A radar device according to the present invention includes: a receiver that receives, as a reception signal, a signal obtained by reflecting a transmission signal periodically increased or decreased in frequency with a constant modulation width by a target; and a beat frequency detector, which is configured to: mix the reception signal and the transmission signal together to generate a beat signal; obtain a first beat frequency distribution according to the beat signal being in an up-chirp state in which the frequency of the transmission signal increases, to thereby specify a first frequency peak of the first beat frequency distribution; and to obtain a second beat frequency distribution according to the beat signal being in a down-chirp state in which the frequency of the transmission signal decreases, to thereby specify a second frequency peak of the second beat frequency distribution. Further, the radar device of the present invention includes: a beat frequency pair selector that produces a pair observation value of the first frequency peak of the first beat frequency distribution and the second frequency peak of the second beat frequency distribution to calculate a range and a Doppler velocity with respect to the target; and an angle measurement processor that calculates an angle of the target based on the pair observation value. Still further, the radar device of the present invention includes: a pair observation value tracking filter that updates a position and a velocity of a track according to the pair observation value including the range, the Doppler velocity, and the angle by means of an existing track; and a beat frequency tracking filter that updates the position and the velocity of the track according to one of the first frequency peak and the second frequency peak by means of the existing track. Yet further, the radar device of the present invention includes: an integration/selection unit that one of integrates the track of the pair observation value tracking filter and the track of the beat frequency tracking filter together, and selects one of the track of the pair observation value tracking filter and the track of the beat frequency tracking filter, as a system track; a system track memory that stores the system track therein; and an abnormal value determination unit that determines, when the pair observation value from the angle measurement processor is not identical with the system track stored in the system track memory, the pair observation value as an abnormal value, and starts no tracking with respect to the abnormal value.
According to the radar device of the present invention, the tracking filter that receives the beat frequencies and the tracking filter that receives the observation value after pairing are used together, thereby enabling the precision in tracking the target to be improved even when the number of peaks in beat frequency is different between the up-chirp state and the down-chirp state, and the erroneous pair to be removed.
In the accompanying drawings:
Hereinafter, a description is given of a radar device according to preferred embodiments of the present invention with reference to the accompanying drawings.
A radar device according to a first embodiment of the present invention is described with reference to
Referring to
The signal processor 1 includes a receiver 11, an A/D converter 12, a beat frequency detector 13, a beat frequency pair selector 14, and an angle measurement processor 15.
The beat frequency tracking filter 2 includes an observation value output determination unit 21, a track input determination unit 22, an up-chirp tracking filter 23, a down-chirp tracking filter 24, and a sub track memory 25.
Referring to
Referring to
Next, an operation of the radar device according to the first embodiment of the present invention is described with reference to the drawings.
The receiver 11 included in the signal processor 1 receives, as a reception signal, a signal obtained by reflecting, by a target, a transmission signal having a frequency periodically increased or decreased with a given modulation width. The A/D converter 12 converts a beat signal with an intermediate frequency, which is output from the receiver 11, into a digital signal.
The beat frequency detector 13 performs frequency analysis through a fast Fourier transform (FFT). That is, the beat frequency detector 13 mixes the reception signal and the transmission signal together to generate a beat signal. The beat frequency detector 13 then obtains a first beat frequency distribution according to the beat signal being in an up-chirp state in which a frequency of the transmission signal increases, and specifies a first frequency peak of the first beat frequency distribution. The beat frequency detector 13 further obtains a second beat frequency distribution according to the beat signal being in a down-chirp state in which the frequency of the transmission signal decreases, and specifies a second frequency peak of the second beat frequency distribution.
A frequency U(t)i of the beat signal in the up-chirp state and a frequency D(t)j of the beat signal in the down-chirp state are extracted. Herein, i and j are representative of the number of peaks obtained by the beat frequency detector 13. It is presumed that to each of the beat frequencies is allocated a time (hereinafter, referred to as “beat frequency observation time”) at which the beat frequency is obtained. As illustrated in
The beat frequency pair selector 14 produces the pair observation value of the first frequency peak of the first beat frequency distribution and the second frequency peak of the second beat frequency distribution to calculate a range and a Doppler velocity with respect to the target.
The angle measurement processor 15 calculates an angle of the target on the basis of the pair observation value from the beat frequency pair selector 14.
The beat frequency obtained by the beat frequency detector 13 is output to the beat frequency tracking filter 2. The beat frequency tracking filter 2 updates, with the use of an existing track, a position and velocity of the track on the basis of the first frequency peak of the first beat frequency distribution or the second frequency peak of the second beat frequency distribution.
The observation value output determination unit 21 included in the beat frequency tracking filter 2 determines whether the input beat frequency has been obtained in the up-chirp state or the down-chirp state, outputs an up-beat frequency to the up-chirp tracking filter 23, and outputs a down-beat frequency to the down-chirp tracking filter 24.
As illustrated in
First, the prediction unit 231 calculates a predicted state vector Xk|k−1 represented by the following expressions (2) and (3) with the use of a updated state vector Xk−1|k−1 represented by the following expression (1) for the existing track (system track or sub track) output from the system track memory 5 (or the sub track memory 25). The prediction unit 231 calculates a state prediction covariance matrix Pk|k−1 represented by the following expression (4) with the use of a updated state covariance matrix Pk−1|k−1 and a drive noise covariance matrix Qk−1. A state transition matrix Φk−1 used for calculation is calculated through the following expression (5) with the use of a difference Δtk−1(=t(k)−t(k−1)) from a time t(k) allocated to the beat frequency newly supplied by the observation value output determination unit 21, which is represented by the following expression (6).
The beat frequency converter 232 converts the track into a beat frequency predicted value fuk|k−1 and a state prediction variance Puk|k−1, as represented by the following expressions (7) to (11).
The correlation unit 233 performs a correlation process of the beat frequency predicted value and the beat frequency. The correlation unit 233 first determines whether or not the observation value fuo of the beat frequency in the up-chirp state at the time tk when the observation value output determination unit 21 executes the output satisfies an inequality represented by the following expression (12). In the expression (12), du is a determination threshold value, and Su is a residual variance of the target which is defined by the following expression (13). In the expression (13), Ak is an measurement error variance.
When there is no beat frequency observation value satisfying the expression (12) at all, the correlation unit 233 does not output the beat frequency to the smoothing unit 234, and the smoothing unit 234 executes a process (hereinafter referred to as “memory track process”) of replacing the updated state value with a predicted value as represented by the following expressions (14) and (15).
xk|k=xk|k−1 (14)
Pk|k=Pk|k−1 (15)
On the other hand, when there is a beat frequency observation value satisfying the expression (12) (when there is a correlation), the correlation unit 233 outputs the beat frequency observation value to the smoothing unit 234, and updates the updated state value of the target. Further, when a plurality of beat frequency observation values satisfy the expression (12), the correlation unit 233 updates the updated state value of the target with the use of a normal correlation algorithm such as nearest neighbor (NN). For example, the updating expression for the updated state value is represented by the following expression (16). In the expression (16), a gain matrix Kuk is calculated by using a theoretical formula such as an extended Kalman filter generally known. FlagEKF has 0 as an initial value, and adds 1 only when there is a correlation, as represented by the following expression (17).
x
k|k
=x
k|k−1
+K
k
u
{f
o
u
−f
k|k−1
u} (16)
FlagEKF=FlagEKF+1 (17)
The smoothing unit 234 outputs the updated state vector, the updated state variance, and the time to the sub track memory 25, and regards the updated state vector and the updated state variance as a sub track. The sub track memory 25 outputs the sub track to the track input determination unit 22.
Subsequently, the track input determination unit 22 outputs the sub track to the down-chirp tracking filter 24. The configuration and operation of the down-chirp tracking filter 24 are identical with those of the up-chirp tracking filter 23 as described above. With a superscript u of a state vector in the expressions (7) to (11) being changed to d, the down-chirp tracking filter 24 operates in the same manner. Therefore, detailed description thereof is omitted.
Assuming that the signal processing time is different as illustrated in
The beat frequency tracking filter 2 updates the sub track with the aid of the beat frequency obtained till the subsequent tracking period. Even when the order of the up-chirp and the down-chirp is replaced with each other in the signal processor 1, it is possible to perform the tracking process according to the order of the input beat frequency. Described above is the processing of the beat frequency tracking filter 2.
Subsequently, processing of the pair observation value tracking filter 3 is described. The pair observation value tracking filter 3, as illustrated in
An operation of the prediction unit 31 is identical with that of the prediction unit 231 of the up-chirp tracking filter 23, and therefore description thereof is omitted. The correlation unit 32 performs the correlation process of the predicted value output from the prediction unit 31 and the pair observation value output from the angle measurement processor 15. The correlation unit 32 first determines whether or not a pair observation value zk at the time tk satisfies an inequality of the following expression (18). In the expression (18), dLKF is a determination threshold value, and S is a residual covariance of the target which is defined in the following expression (19). In the expression (19), Ak is an measurement error covariance.
(zk−Hkxk|k−1)TSk|k−1−1(zk−Hkxk|k−1)≦dLKF (18)
S
k|k−1
=H
k
P
k|k−1
H
k
T
+A
k (19)
H=I4×4 (20)
zk=[Rok sin θokRok cos θok{dot over (R)}ok sin θok{dot over (R)}ok cos θok]T (21)
When it is determined that there is no correlation, the correlation unit 32 does not output the pair observation value to the smoothing unit 33, and the smoothing unit 33 executes the memory track process represented by the expressions (14) and (15). On the other hand, when it is determined that there is a correlation, the correlation unit 32 outputs the pair observation value to the smoothing unit 33 to update the updated state value of the target. The smoothing expression is represented by the following expression (22). A gain matrix Kk is calculated by using a theoretical formula such as a Kalman filter generally known. FlagLKF has 0 as an initial value, and adds 1 only when there is a correlation, as represented by the following expression (23).
x
k|k
=x
k|k−1
+K
k
{z
k
−H
k
x
k|k−1} (22)
FlagLKF=FlagLKF+1 (23)
The pair observation value tracking filter 3 outputs the updated state vector, the updated state covariance matrix, and information on the paired up-beat frequency and down-beat frequency allocated to the pair observation value correlated to the track to the integration/selection unit 4.
Next, an operation of the integration/selection unit 4 is described with reference to
The integration/selection unit 4 first performs checking of the EKF track and the LKF track (Steps 401 to 403). In this time, when both of those tracks are memory tracks, the integration/selection unit 4 regards the system track as a memory track, and registers the EKF track as the system track in the system track memory 5 (Step 411).
xk|kSYS=xk|kEKF (24)
Pk|kSYS=Pk|kEKF (25)
When the LKF track exists, and the EKF track is a memory track, the integration/selection unit 4 assumes that the LKF track is updated by the observation value of the erroneous pair, and deletes the LKF track (Step 408). Further, the integration/selection unit 4 registers the EKF track in the system track memory 5 as the system track (Step 409). When the LKF track is a memory track, and the EKF track exists, the integration/selection unit 4 registers the EKF track in the system track memory 5 as the system track (Step 410).
When both of the EKF track and the LKF track exist, the integration/selection unit 4 evaluates the adequacy of the track. First, the integration/selection unit 4 performs the determination of correlation between the tracks through the following expression (26) (Step 404). In this case, x and P output from the respective tracks correspond to xk|k and Pk|k, respectively. In the expression (26), dTRK represents the determination threshold value.
(xLKF−xEKF)T(PLKF+PEKF)−1(xLKF−xEKF)≦dTRK (26)
When the expression (26) is not satisfied in the above-mentioned correlation determination, the integration/selection unit 4 assumes that the tracks are of no correlation, and registers the EKF track in the system track memory 5 as the system track (Step 406). Then, the integration/selection unit 4 deletes the LKF track (Step 407).
When the expression (26) is satisfied, the integration/selection unit 4 assumes that the tracks are of correlation, and integrates the tracks together (Step 405). As a method involving integrating the tracks together, the integration/selection unit 4 integrates the tracks together through a covariance intersection technique, taking a color property into consideration, for example, as represented by the following expressions (27) and (28). The integration/selection unit 4 registers the integrated track in the system track memory 5 as the system track. In this case, a state vector xSYSk|k and the updated state covariance matrix PSYSk|k of the system track are given by the following expressions (27) and (28), respectively. In the expressions (27) and (28), ω represents a parameter. Further, the integration/selection unit 4 also registers, in the system track memory 5, the information on the pair of up-beat frequency and down-beat frequency which are allocated to the pair observation value correlated with the track of the pair observation value tracking filter 3.
x
k|k
SYS
=P
k|k
SYS
[ω{P
LKF}−1xLKF+(1−ω){PEKF}−1xEKF] (27)
P
k|k
SYS
=[ω{P
LKF}−1+(1−ω){PEKF}−1]−1 (28)
Further, there may be applied a least square integration method as represented by the following expressions (29) and (30) though the color property is not considered. The method has no need to set the parameter ω.
x
k|k
SYS
=P
k|k
SYS
[{P
LKF}−1xLKF+{PEKF}−1xEKF] (29)
P
k|k
SYS
=[{P
LKF}−1+{PEKF}−1]−1 (30)
For the purpose of reducing a calculation load, there is a weighting and integrating method using a trace of the updated state covariance matrix of the track as represented by the following expression (31) (it is assumed that the trace of a matrix A is tr(A)). The method has no need to calculate an inverse matrix though the updated state covariance matrix of the system track cannot be calculated.
Further, with an aim to suppress the calculation load, there are a weighting and integrating method using a predetermined parameter a(0≦a≦1) as represented by the following expression (32), and a weighting and integrating method depending on the update status of the track, as represented by the following expression (33). In this case, nLKF in the expression (33) represents the number of times by which the tracks are input to the integration/selection unit 4 without the LKF track becoming the memory track during N samplings in the past from the present. In the case of the EKF track, nLKF is replaced with nEKF. In particular, when a simple tracking filter such as an α−β filter is employed as the tracking filter, the method is effective with no need to update the updated state covariance.
When the integration/selection unit 4 continuously selects the EKF track because of no correlation between the tracks in Step 404 of
Further, when the LKF track is the memory track in Step 403 of
The abnormal value determination unit 6 regards the pair observation value input from the angle measurement processor being not identical with the information on the up-beat frequency and the down-beat frequency which are input from the system track memory 5, as an abnormal value, and does not start tracking. That is, when the pair observation value from the angle measurement processor 15 is not identical with the system track, the abnormal value determination unit 6 regards the pair observation value as the abnormal value, and does not start tracking with respect to the abnormal value. In addition, the abnormal value determination unit 6 may output, to the beat frequency pair selector 14 within the signal processor 1, the up- and down-frequency pair allocated to the pair observation value used for the system track so as to be paired preferentially.
As described above, according to the first embodiment, in the radar device with the FMCW radar, the tracking precision can be improved as compared with the conventional radar device with the FMCW radar to thereby enable a reduction in occurrence of the erroneous track, by means of the integration/selection unit 4 for integrating the tracks together or selecting one of the tracks, and the abnormal value determination unit 6 for determining the erroneous pair based on the pair correlated with the track. In this case, there are used the beat frequency tracking filter 2 directly receiving the beat frequencies obtained in the up-chirp state and the down-chirp state to update the position and the velocity of the target, and the pair observation value tracking filter 3 receiving the range, the Doppler velocity, and the angle of the target, which are obtained from the pair of the up-chirp state and the down-chirp state to update the position and the velocity of the target, together.
Further, according to the first embodiment, in the radar device with the FMCW radar, the abnormal value determination unit 6 outputs, to the signal processor 1, the pair correlated with the track as a candidate that is preferentially paired, thereby enabling the occurrence of the erroneous pair to be reduced.
Further, according to the first embodiment, in the radar device with the FMCW radar, the tracking precision can be enhanced by the integration/selection unit 4 selecting the track small in residual error when the track of the beat frequency tracking filter 2 is continuously selected even if the tracks of the pair observation value tracking filter 3 and the beat frequency tracking filter 2 are not correlated with each other.
Further, according to the first embodiment, in the radar device with the FMCW radar, when the track of the beat frequency tracking filter 2 is continuously selected, and there is no observation value correlated by the pair observation value tracking filter 3, the integration/selection unit 4 integrates a predicted track of the track of the pair observation value tracking filter 3 and the track of the beat frequency tracking filter 2 together, thereby enabling the target tracking precision to be enhanced.
Further, according to the first embodiment, in the radar device with the FMCW radar, the integration/selection unit 4 executes weighting and integration with the use of a predetermined parameter when integrating the tracks together, thereby enabling the calculation load to be reduced.
Further, according to the first embodiment, in the radar device with the FMCW radar, the integration/selection unit 4 executes weighting and integration through the covariance intersection technique when integrating the tracks together, taking the color property into consideration, thereby enabling the tracking precision to be ensured.
Further, according to the first embodiment, in the radar device with the FMCW radar, the integration/selection unit 4 executes weighting and integration through the least square integration method when integrating the tracks together, thereby enabling the tracking precision to be ensured without setting the parameter.
Further, according to the first embodiment, in the radar device with the FMCW radar, the integration/selection unit 4 executes weighting and integration with the use of trace of the updated state covariance matrix of each track when integrating the tracks together, thereby enabling the calculation load to be reduced.
Further, according to the first embodiment, in the radar device with the FMCW radar, the integration/selection unit 4 executes weighting and integration with the use of the number of updating times of the tracking for each track when integrating the tracks together, thereby enabling the calculation load to be reduced.
A radar device according to a second embodiment of the present invention is described with reference to
Referring to
Referring to
Next, an operation of the radar device according to the second embodiment of the present invention is described with reference to the drawings.
The radar device according to the second embodiment is configured to input the angle of the pair observation value to the beat frequency tracking filter 2A to execute the tracking process.
The angle tracking filter 26, as illustrated in
The angle converter 262 converts the track into an angle predicted value θk|k−1 and a state prediction variance Pθk|k−1, as represented by the following expressions (34) to (36).
The correlation unit 263 performs a process of correlating the predicted value output from the angle converter 262 with the pair observation value output from the observation value output determination unit 21. First, the correlation unit 263 determines whether or not the pair observation value zk at the time tk satisfies an inequality of the following expression (37). In the expression (37), dθ is a determination threshold value, and Sθ is a residual covariance of the target which is defined in the following expression (38).
When it is determined that there is no correlation, the correlation unit 263 does not output the pair observation value to the smoothing unit 264, and the smoothing unit 264 executes the memory track process represented by the above-mentioned expressions (14) and (15).
On the other hand, when it is determined that there is a correlation, the correlation unit 263 outputs the angle observation value to the smoothing unit 264, and updates the updated state value of the target.
x
k|k
=x
k|k−1
+K
k
θ{θo−θk|k−1} (39)
FlagEKF=FlagEKF+1 (40)
The smoothing unit 264 outputs the updated state vector, the updated state variance, and the time to the sub track memory 25, and regards the updated state vector and the updated state variance as a sub track. The sub track memory 25 outputs the sub track to the track input determination unit 22. Other processing is identical with that in the above-mentioned first embodiment, and therefore description thereof is omitted.
As described above, according to the second embodiment, in the radar device with the FMCW radar, the beat frequency tracking filter 2A directly receives the angle of the pair observation value in addition to the up-chirp and the down-chirp to update the position and the velocity of the target, thereby enabling the tracking precision to be improved.
A radar device according to a third embodiment of the present invention is described with reference to
Referring to
Referring to
Next, an operation of the radar device according to the third embodiment of the present invention is described with reference to the drawings.
The radar device according to the third embodiment is configured to input, when angle measurement is executed every beat frequency to obtain an angle, the angle to the beat frequency tracking filter 2B to perform the tracking process.
As illustrated in
The beat frequency angle converter 232B calculates the beat frequency predicted value fuk|k−1 and the state prediction variance puk|k−1 thereof, and the angle predicted value θk|k−1 and the state prediction variance pθk|k−1 thereof. Formula for calculation is omitted.
The correlation unit 233 performs a process of correlating the predicted value output from the beat frequency angle converter 232B with the beat frequency and the angle observation value. First, the correlation unit 233 determines whether or not the beat frequency observation value fuo and the angle observation value θo at the time tk satisfy an inequality of the following expression (41). In the expression (41), d is a determination threshold value, Su is a residual covariance defined in the above-mentioned expression (13), and Sθ is a residual covariance defined in the above-mentioned expression (38).
When it is determined that there is no correlation, the correlation unit 233 does not output the beat frequency and the angle observation value to the smoothing unit 234, and the smoothing unit 234 executes the memory track process represented by the above-mentioned expressions (14) and (15).
On the other hand, when it is determined that there is a correlation, the correlation unit 233 outputs the beat frequency and the angle observation value to the smoothing unit 234, and updates the updated state value of the target.
x
k|k
=x
k|k−1
+K
k
u,θ
{z
ok
−z
k|k−1} (42)
zok=[fouθo]T (43)
zk|k−1=[fk|k−1uθk|k−1]T (44)
FlagEKF=FlagEKF+1 (45)
Next, the smoothing unit 234 outputs the updated state vector, the updated state variance, and the time to the sub track memory 25, and regards the updated state vector and the updated state variance as a sub track. The sub track memory 25 outputs the sub track to the track input determination unit 22. Other processing is identical with that in the above-mentioned first embodiment, and therefore description thereof is omitted.
As described above, according to the third embodiment, in the radar device with the FMCW radar, the beat frequency tracking filter 2B directly receives, in addition to the up-chirp and the down-chirp, the angle allocated thereto to update the position and the velocity of the target, thereby enabling the tracking precision to be improved.
Number | Date | Country | Kind |
---|---|---|---|
2009-104154 | Apr 2009 | JP | national |