Range-angle Decoupled Transmit Beamforming with Frequency Diverse Array
DOI: 10.12000/JR16113
-
Abstract:
It has been shown that Frequency Diverse Arrays (FDA) exhibit a range-angle dependent beam steering feature by employing a uniform frequency increment across the array elements. However, this beam pattern generates maxima at multiple range values, possibly leading to loss of signal-to-interference-plus-noise ratio when the interferences are located at any of the maxima. Herein, we prove that the beam pattern of FDA is range-periodic and propose the basic criteria for the FDA configuration to decouple the range and angle. In an illuminated space, a single-maximum beam pattern can be obtained by configuring the frequency increment between the elements. Specific examples have been discussed herein, and the simulation results verify the proposed theory.
-
1. Introduction
Frequency Diverse Arrays (FDA) radar has recently drawn much attention among the researchers. FDA differs from the traditional phased array by using a small frequency increment across the array elements, which results in a range-angle-dependent beam pattern. FDA radar system is first proposed in Ref. [1]. In FDA radar, a uniform interelement frequency offset is applied across the array elements. FDA radar with uniform small and large frequency offset frequency has been investigated in Refs. [1–7]. Small frequency offset has been exploited to generate range-dependent beampattern, while large frequency offset can get independent echoes from the target.
Unlike the phased array, the range-angle dependency of the FDA beampattern allows the radar system to focus the transmit energy in a desired range-angle space. This unique feature of FDA helps to suppress the range-dependent interferences[8] and increases the received SINR consequently. Especially for the mainlobe interference and clutter, the FDA can achieve a significant improvement in SINR against the phased array because the FDA provides the increased Degrees Of Freedom (DOFs) in range domain. However, the FDA beampattern is shown to be periodic in range and time[2], which goes to maximum at multiple time and range values. With this multiple-maximum beampattern, the resulting SINR will be deteriorated when the interferers are located at any of the maxima. To improve SINR, FDA with Time-Dependent Frequency Offset (TDFO-FDA) was proposed to achieve a time-independent beampattern at the target location[9]. Nevertheless, the proposed beampattern is still periodic in range which will result in the loss of SINR. A nonuniformly spaced linear FDA with linear incremental frequency increment has also been studied in Ref. [10], and a nonrepeating beampattern has been obtained for range-angle imaging of targets. A uniformly spaced linear FDA with Logarithmically (Log-FDA) increasing frequency offset is proposed in Ref. [11]. The proposed strategy provides a nonperiodic beampattern with the single-maximum in space. In Ref. [12], the beampattern of FDA who transmits the pulsed signal has been studied. Lately, few more publications have done some work in decoupling the range-angle dependent beampattern of FDA[13–16]. All these papers only address the properties of the FDA beampattern, and they do not study the common rule for the FDA configuration to form a single-maximum transmit beampattern.
With the pioneer work on FDA radar, we aim to decouple the range and angle in the beampattern and provide a nonperiodic beampattern with the single-maximum in the illuminated range-angle space. In this paper, we propose a basic criteria for the FDA configuration, in which the element spacing and frequency increment are configurable, to form a single-maximum beampattern through mathematical analysis. This single-maximum beampattern, unlike the multiple-maxima beampattern, can help to further suppress range-dependent interferences, causing improved SINR and increased detect ability. The proposed rule for the FDA configuration will be helpful to design the FDA.
The rest of the paper is organized as follows. In Section 2, the basic FDA model has been described and the basic criterion is derived for the FDA configuration to form a single-maximum beampattern through mathematical analysis. Moreover, several specific conditions are introduced. In Section 3, the beampattern has been plotted for the specific conditions discussed in Section 2. Finally, in Section 4 we conclude the paper.
2. Design and Mathematical Analysis of FDA
2.1 System description
Consider an array of M transmit elements, we assume that the waveform radiated from each antenna element is identical with a frequency increment, as shown in Fig. 1. The radiated frequency from the m-th element is
fm=f0+Δfm, m=0,···,M−1 (1) where Dfm is the frequency increment of m-th element with reference to the carrier frequency f0. Specifically, Df0 = 0.
Considering a given far-field point, the phase of the signal transmitted by the m-th element can be represented by
ψm=2πfm(t−rmc) (2) where c and rm are the speed of light and the distance between the m-th element and the observed point, respectively.
The range difference between individual elements is approximated by
rm=r0−dmsinθ0 (3) where θ0 is the desired angle, dm is the spacing between the m-th element and the first element. Specifically, d0 = 0.
So the phase difference between the m-th element and the first element is
Δψm=ψm−ψ0=2π(fm(t−rmc)−f0(t−r0c))=2πΔfmt−2πfmdmsinθ0c+2πΔfmr0c (4) In Eq. (4), the third term is important because it shows that the FDA radiation pattern depends on both the range and the frequency increment. Taking the first element as the reference for the array, the steering vector can be expressed as
a(θ,r,t)=[1,···,exp(−j2π−Δfm(tc+r)+fmdmsinθc)]T (5) where [·]T denotes the transpose operator.
In the pusled-FDA, for
t∈[−Te2,Te2] , Te is the pulsewidth, the maximum value of phase variance[12] during the pulse duration can be derived asξ=maxm2πΔfmTe (6) When the phase variance
ξ is small enough, the beampattern of the pulsed-FDA can be viewed as quasi-static. Actually, in practical radar systems, duty cycle other than pulse duration is often used to describe the characterization of the pulsed-waveform. Then, the phase varianceξ can be further written as|ξ|=maxm2π|Δfm|frdt=2πdtmaxmρm≪1 (7) where dt is the duty cycle, which is usually small. Eq. (7) holds when
maxmρm is also small.So under the condition Eq. (7), the time t can be neglected. So Eq. (5) can be simplified as
a(θ,r)=[1,···,exp(−j2π−Δfmr+fmdmsinθc)]T (8) 2.2 Transmit beampattern analysis
Throughout this paper, we assume a narrow-band system where the propagation delays manifest as phase shifts to the transmitted signals and Eq. (7) is satisfied. To steer the maximum at an expected target location (
θ 0, r0), the complex weights are configured as a(θ 0, r0), so the transmit beampattern can be expressedAF(θ,r)=|aH(θ0,r0)a(θ,r)|=|M−1∑m=0exp(j2πΔfm(r−r0)−fmdm(sinθ−sinθ0)c)| (9)
where [·]H denotes the conjugate transpose operator.
It is easy to see that the beam direction will vary as a function of the range and angle, which means the beampattern is range-angle dependent. Since the beampattern is coupled in the range and angle, the target’s range and angle cannot be estimated directly by the FDA beamformer output. Note that the beampattern is also related to Dfm and dm, so the desired single-maximum beampattern can be obtained by setting the proper Dfm and dm.
In Eq. (9), when m = 0, the exponential term is equal to 1. To obtain the maximum value of the beampattern, the exponential terms should be all equal to 1 for
m=1,2,···,M−1 . So the phase of the exponential term should be the integral multiple of2π , which can be expressed asΔfm(r−r0)−f0dm(sinθ−sinθ0)−Δfmdm(sinθ−sinθ0)c=Lm (10) where Lm is an integer, e.g. Lm = 0, ±1,···,
m=1,2,···,M−1 .The Eq. (10) can be rewritten as
r=Lmc+f0dm(sinθ−sinθ0)Δfm+dm(sinθ−sinθ0)+r0 (11) Since
dm(sinθ−sinθ0)≪r0 , the termdm(sinθ −sinθ0) in Eq. (8) can be neglected. The curves formed by Eq. (11) will be called as range-angle distribution curves throughout the paper. Then Eq. (11) can be approximately expressed asΔfm(r−r0)−f0dm(sinθ−sinθ0)=Lmc, m=1,2,···,M−1 (12) Rewrite Eq. (12) into matrix form as
Ax=b (13) where
A=[Δf1−f0d1Δf2−f0d2⋮⋮ΔfM−1−f0dM−1] ,x=[r−r0sinθ−sinθ0] ,b=[L1cL2c⋮LM−1c] , Lm = 0, ±1, ···,m=1,2,···, M – 1.To decouple the range and angle, the beam-pattern should have the unique maximum point in the range-angle distribution diagram, which means the Eq. (13) has the unique solution (
θ 0, r0). The necessary and sufficient condition of that the Eq. (13) has the unique solution isrank(A)=rank(˜A)=2 (14) where rank (·) is the rank of a matrix,
˜A=(A,b) . Whenrank(A)=rank(˜A)=1 , the Eq. (13) has infinite solutions, corresponding to the conventional FDA condition, which will be discussed in detail later.To satisfy rank (A) = 2,
dm≠PΔfm , P is a constant.To satisfy
rank(A)=rank(˜A) ,Lm=Qdm orLm=SΔfm , Lm is an integer, and assume that Q, S are the minimum non-zero constants to satisfy the equations. Sincedm≠PΔfm must be satisfied, the two equations cannot be hold at the same time but when Lm = 0,∀m=1,2,···, M–1.In the following, under the condition of
dm≠PΔfm,Lm=0,±1 ,···, we make a summary with different parameter configurations:(1) when Lm = 0, the Eq. (13) has the unique solution (
θ 0, r0);(2) when
Lm=Qdm≠0 ,Lm≠SΔfm , the Eq. (10) has the solutions(arcsin(sin(θ0)− kQλ0),r0),λ0=cf0,k=0,±1,±2,··· . When|sin(θ0)−kQλ0|≤1 , the angle grating lobes will occur at anglearcsin(sin(θ0)−kQλ0) in the beampattern. Otherwise, the Eq. (13) has no solution, resulting in no angle grating lobes;(3) when
Lm=SΔfm≠0,Lm≠Qdm , the Eq. (13) has the(θ0,r0+Skc) ,k=0,±1,±2,··· , which means the range grating lobes will occur atr0+Skc in the beampattern.Note that in array theory, when the adjacent element spacing is less than half the wavelength, the angle grating lobes will never appear. If
Qλ0=±1 andθ0=0∘ , where the element spacing isλ 0, then the grating lobes will occur at angle ±90°. ButLm=SΔfm can always be satisfied since Lm is an integer whose range is[−∞,+∞] . The range grating lobes will always occur at ranger0+Skc in the beampattern. The position of the range grating lobe changes with different S. For example, S = 0.01, the distance between the grating lobe and the mainlobe is 3000 km. So ifr0±Sc is out of the illuminated range space[Rmin,Rmax] , the beampattern has a single-maximum point (θ 0, r0) in the illuminated range space, which means the range and angle are decoupled.So we can conclude that the beampattern of the FDA is always range-periodic, the grating lobes will always occur at range
r0+Skc ,k=0,±1,±2,··· . To obtain the single-maximum beampattern in the illuminated range space[Rmin,Rmax] , the designing criteria for the FDA isdm≠PΔfm , and2Sc>Rmax−Rmin ,|sin(θ0)±Qλ0|>1 ,k=0,±1,±2,··· ,Lm=0,±1,··· ,m=1,2,··· ,M – 1, P is a constant, Q, S are the minimum non-zero constants to satisfy the equationsLm= SΔfm≠0 andLm=Qdm≠0 .Once the range and angle is decoupled, the target’s range and angle can be estimated directly by the FDA beamformer output. Also the 2-dementional MUSIC algorithm[16] for estimating the target’s range and angle can be used as well.
2.3 Specific examples
For the conventional FDA,
Δfm=mΔf , dm = md, Df and d are configurable parameters to control the frequency increment and the element spacing. When Lm = nm,n=0,±1,±2,··· ,m=1,2,···,M−1 , we can getrank(A)=rank (˜A)=1 . Under this circumstance, the Eq. (13) has infinite solutions. The solutions form the range-angle dependent curves, which have expressions as:r=f0dΔfsinθ−f0dsinθ0Δf+r0+ncΔf (15) In Eq. (15), the expression is no longer related to m, which means the range-angle curve for different element coincides with each other, as depicted in Fig. 2(a). In the range-angle distri-bution diagram, the curve is periodic in range, and the range difference between the adjacent curves is c/Df. The corresponding beampattern is depicted as Fig. 3(a).
For the Expf-FDA, whose frequency incre-ment is exponentially increased,
Δfm=(bm −1)Δf , dm = md. b is a configurable constant. The range solution rises when b gets larger. The range-angle distribution curves and the beampattern are depicted in Fig. 2(b) and Fig. 3(b), respectively.Likewise, for the Logd-FDA, whose element spacing is logarithmically increased,
Δfm=mΔf ,dm=log(m+1)d. When Lm = nm,n=0,±1, ±2,··· ,m=1,2,···,M−1 , we can getrank(A)= rank(˜A)=2 , so the range grating lobes will occur in the beampattern as depicted in Fig. 3(c), and the range difference between the grating lobe and the mainlobe is c/Df. The corresponding range-angle distribution diagram is depicted in Fig. 2(c).For another kind of FDA, called Expf-Logd-FDA, where the frequency increment is exponentially increased and the element spacing is logarithmically increased,
Δfm=(bm−1)Δf ,dm=log(m+1)d . The range-angle distribution curves and beampattern are depicted in Fig. 2(d) and Fig. 3(d), respectively.3. Simulations, Results, and Discussions
Bampattern expressed in Eq. (9) and the range-angle distribution curves expressed in Eq. (11) were simulated and plotted for different kinds of FDA discussed in Section 2. The results are discussed and compared with the different kinds of FDA. The illuminated range space is (0 km, 800 km]. To generate these plots, the values of the configurable parameters have been taken as listed in Tab. 1. To avoid angle grating lobes, the parameter d is less than half the wavelength.
Table 1. Parameters for simulationsParameter Value Parameter Value Element number M 8 d 0.1 m Reference frequency f0 1 GHz b 1.4 {\Delta}f 1 kHz Desired point
( \theta_0, r0)(0°, 400 km) The range-angle distribution curves are de-picted in Fig. 2. The curves with different color represent the range-angle distribution for different elements except the reference element in the FDA. It is easy to see that the range-angle distribution of the elements are the same in conventional FDA, and the range difference between the adjacent curves is c/Df = 300 km, which is consistent with the analysis in Section 2. For the Expf-FDA and Expf-Logd-FDA,
Δfm =(1.4m−1)Δf , the range between the first grating lobe and the mainlobe is 3e6 km, which is not in the illuminated space. So the curves of the elements in the range-angle distribution diagram have the single intersection point, which means that the beampattern has the single-maximum point in the illuminated range space. But for the Logd-FDA, when Lm = nm, the distance between the grating lobe and the mainlobe is c/Df = 300 km according to the analysis in Section 2. So in the range-angle distribution diagram, the curves have 3 intersection points, the range difference between adjacent intersection points is 300 km.The beampatterns are depicted in Fig. 3. Similar to the range-angle curve in Fig. 2, the beampattern of the conventional FDA is a range-angle-dependent beampattern. For the Expf-FDA and Expf-Logd-FDA, the beampatterns have a single-maximum point in the illuminated space corresponding to the single intersection point in the range-angle distribution diagram. For the Logd-FDA, the beampattern has 3 range grating lobes corresponding to the 3 intersection points in the range-angle distribution diagram. For the latter 3 FDAs, the contour of the beampattern is an ellipse, which is because the range-angle distribution curves are tightly distributed. The direction of the major axis of the ellipse is corresponding to the most tightness distribution direction of the range-angle distribution curves.
In the meanwhile, we analyze the range and angle solutions of different FDAs in Fig. 4. The Expf-FDA and the Exp-Logd-FDA have the same range solution, and the range solution of conventional FDA and Logd-FDA are the same. That is because the bandwidth across the whole array is the same for each pair. The same situation occurs for the angle solution, if the arrays have the same array aperture, they have the same angle solutions.
We choose the target position at (400 km, 20°), and the transmit beampattern is depicted in Fig. 5. Similar to Fig. 3, we can see that the single maximum beam is formed at the target position.
4. Conclusions
In this paper, we propose a basic criteria for the FDA configuration to provide a single-maximum beampattern in the illuminated space. The single-maximum beampattern can be generated by configuring the element spacing and frequency increment of the FDA. Through the analysis, we can find out that the beampattern of FDA is always range-periodic. Results show that a single-maximum beampattern can be generated with the corresponding criteria by choosing the proper frequency increment. As this paper describes the transmitter only, designing an appropriate receiver is our future work.
-
Table 1. Parameters for simulations
Parameter Value Parameter Value Element number M 8 d 0.1 m Reference frequency f0 1 GHz b 1.4 {\Delta}f 1 kHz Desired point
( \theta_0, r0)(0°, 400 km) -
[1] Antonik P, Wicks M C, and Griffiths H D. Range dependent beamforming using element level waveform diversity[C]. International Waveform Diversity Design Conference, Las Vegas, NV, USA, Jan. 2006: 22–27. [2] Secmen M, Demir S, and Hizal A. Frequency diverse array antenna with periodic time modulated pattern in range and angle[C]. IEEE Radar Conference, Boston, MA, USA, Apr. 2007: 427–430. [3] Antonik P, Wicks M C, and Griffiths H D. Multi-mission, multi-mode waveform diversity[C]. IEEE Radar Conference, Verona, NY, USA, Apr. 2006: 580–582. [4] Zhuang L and Liu X Z. Precisely beam steering for frequency diverse arrays based on frequency offset selection[C]. International Radar Conference, Bordeaux, France, Oct. 2009: 1–4. [5] Chen Y G, Li Y T, and Wu Y H. Research on the linear frequency diverse array performance[C]. IEEE International Conference on Signal Processing, Beijing, China, Oct. 2010: 2324–2327. [6] Sammartino P F, Backer C J, Griffiths H D. Frequency diverse MIMO techniques for radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(1): 201–222. DOI: 10.1109/TAES.2013.6404099 [7] Xu J W, Liao G S, and Zhu S Q. Receive beamforming of frequency diverse array radar systems[C]. 31th URSI General Assembly and Scientific Symposium (URSIGASS), Beijing, China, Aug. 2014: 1–5. [8] Wang W Q. Range-angle dependent transmit beampattern synthesis for linear frequency diverse arrays[J]. IEEE Transactions on Antennas Propagation, 2014, 61(8): 4073–4081. [9] Khan W, Qureshi I M. Frequency diverse array radar with time dependent frequency offset[J]. IEEE Antennas and Wireless Propagation Letters, 2014, 13: 758–761. DOI: 10.1109/LAWP.2014.2315215 [10] Wang W Q, So H C, Shao H. Nonuniform frequency diverse array for range-angle imaging of targets[J].IEEE Sensors Journal, 2014, 14(8): 2469–2476. DOI: 10.1109/JSEN.2014.2304720 [11] Khan W, Qureshi I M, Saeed S. Frequency diverse array radarwith logarithmically increasing frequency offset[J]. IEEE Antennas and Wireless Propagation Letters, 2015, 14: 499–502. DOI: 10.1109/LAWP.2014.2368977 [12] Xu Y H, Shi X W, Xu J W. Range-angle-dependent beamforming of pulsed frequency diverse array[J]. IEEE Transactions on Antennas Propagation, 2015, 63(7): 3262–2367. DOI: 10.1109/TAP.2015.2423698 [13] Gao K, Wang W Q. Decoupled frequency diverse array range-angle-dependent beampattern synthesis using non-linearly increasing frequency offsets[J]. IET Microwaves,Antennas&Propagation, 2016, 10(8): 880–884. [14] Gao K, Wang W Q. Transmit beamspace design for multi-carrier frequency diverse array sensor[J]. IEEE Sensors Journal, 2016, 16(14): 5709–5714. DOI: 10.1109/JSEN.2016.2573379 [15] Shao H, Dai J, Xiong J. Dot-shaped range-angle beampattern synthesis for frequency diverse array[J]. IEEE Antennas and Wireless Propagation(published online) , 2016. DOI: 10.1109/LAWP.2016.2527818 [16] Schmidt R O. Multiple emitter location and signal parameter estimation[J]. IEEE Transactions on Antennas Propagation, 1986, 34(1): 276–280. 期刊类型引用(7)
1. 赵研,陶海红,畅鑫. 天线阵列近场FDA技术研究. 雷达科学与技术. 2023(05): 505-510 . 百度学术
2. 王博,陈楚舒,张君鹏,姜杰. 密度锥削阵FDA雷达方向图特性分析. 雷达科学与技术. 2021(04): 455-459+466 . 百度学术
3. 于雷,何峰,董臻,粟毅,张永胜,吴曼青. 一种基于非线性调频信号和空域编码的FDA雷达波形设计方法. 雷达学报. 2021(06): 822-832 . 本站查看
4. 王博,谢军伟,张晶,葛佳昂. 频控阵阵列模型及参数分析. 火力与指挥控制. 2020(05): 25-32 . 百度学术
5. 项靖,舒勤. 基于卡方分布的频率分集阵列解耦波束形成. 信息技术与信息化. 2020(06): 145-149 . 百度学术
6. 吴秉坤,章新华,李鹏. 基于频率多样性的圆阵方位估计方法研究. 舰船科学技术. 2020(17): 166-170 . 百度学术
7. 陈浩,韩宏川,李荣锋,张雪. 频率分集阵列及其研究进展. 现代防御技术. 2019(01): 63-74 . 百度学术
其他类型引用(6)
-