By introducing phase code into multicarrier orthogonal frequency division multiplex signal, the multicarrier phase coded (MCPC) radar signal possesses a good spectrum utilization rate and can achieve a good combination of narrowband and wideband processing. Radar pulse train signal not only reserves the high range resolution of monopulse signal, but also has the same velocity resolution performance as continuous wave signal does. In this study, we use the chaotic biphase code generated by Chebyshev mapping to conduct a phase modulation on MCPC pulse train so as to design two different types of multicarrier chaotic phase coded pulse train signal. The ambiguity functions of the two pulse train signals are compared with that of P4 code MCPC pulse train. In addition, we analyze the influences of subcarrier number, phase-modulated bit number, and period number on the pulse train’s autocorrelation performance. The low probability of intercept (LPI) performance of the two signals is also discussed. Simulation results show that the designed pulse train signals have a thumbtack ambiguity function, a periodic autocorrelation side lobe lower than P4 code MCPC pulse train, and excellent LPI performance, as well as the feature of waveform diversity.
1. Introduction
The MCPC signal based on OFDM technique [1] is a new wideband radar signal which has attracted much attention recently. It has a flexible structure, a thumbtack ambiguity function. This signal not only has the advantages of narrowband radar but also can synthesize the wideband signal with orthogonal narrowband signals so as to achieve multichannel separation and quick processing. Starting from OFDM signal, Levanon and Mozeson [2] analyzed the structure, spectrum, ambiguity function, autocorrelation, and power spectrum of monopulse, continuous wave, and pulse train of MCPC signal. Besides, the methods for decreasing the peak-to-mean envelope power ratio (PMEPR) of these signals are elaborated. Studies were performed on the waveform design, pulse compression, target detection, and imaging of MCPC signal in [3]. Reference [4] introduced a method of least squares for the allocation of a proper phase to obtain desired ambiguity function. In literature [5], the features of MCPC signal and frequency stepped signal were combined to propose the multicarrier phase coded frequency stepped radar signal, which has a high range resolution yet has the defect of range-Doppler coupling.
In order to extract Doppler information, monopulse needs to be accumulated to obtain the pulse train signal. The pulse train signal reserves the high range resolution of monopulse signal while maintaining the velocity resolution performance of continuous wave signal [6]. However, the traditional pulse train signal has a poor performance in measuring remote and high speed targets because of range or velocity ambiguity, and the antijamming performance is also low. MCPC signal has a flexible structure, and different pulse train signals can be obtained using different phase coding methods. Studying the MCPC pulse train signal with complex modulation can improve the detection performance and LPI performance of radar.
Like the noise signal, the chaotic signal exhibits continuous power spectrum, initial-value sensitivity, ergodicity, and aperiodicity, and therefore chaotic signal possesses good range and velocity resolutions, thumbtack ambiguity function, excellent antijamming ability, and electromagnetic compatibility. Thus, chaotic signal avoids the problem of range ambiguity. Moreover, it overcomes many shortcomings of the noise signal, such as the difficulties of generation, replication, and application. These properties of chaotic signal conform to the low probability of intercept (LPI) and antijamming requirements for modern radar.
In this study, the multicarrier chaotic phase coded pulse train signal is designed using the chaotic biphase code for phase modulation on MCPC pulse train by different modes. The ambiguity function as well as the autocorrelation performance of the signal is simulated, and the LPI performance is also analyzed.
2. Structure of MCPC Pulse Train Signal
The complex envelope of the monopulse MCPC signal is given by
(1)f(t)=∑n=1N∑m=1Mωnan,ms[t-(m-1)tb]hhhhhhhhh×exp(j2π(n-1)Δft),
where N is the subcarrier number and M is the phase-modulated bits. Consider an,m=ejφn,m, where φn,m is the mth phase element of the nth subcarrier. One has s(t)≡1 for 0≤t<tb and zero elsewhere. The frequency difference between two adjacent subcarriers Δf is set equal to the inverse of the bit duration tb yielding orthogonal frequency division multiplexing. ωn=|ωn|ejθn (∑n=0N-1|ωn|2=1) is the complex weight associated with the nth subcarrier, and |ωn| is the frequency weighted amplitude. θn is the frequency weighted phase, also known as the initial phase. The structure of the MCPC signal is illustrated schematically in Figure 1.
Structure of MCPC signal.
The schematic diagram shows that the MCPC signal consists of N sequences transmitted simultaneously on N subcarriers. Each sequence contains M phase-modulated bits.
The general expression for the MCPC pulse train complex envelope fTr(t) is
(2)fTr(t)=∑p=1P∑n=1N∑m=1Mwnan,m,ps[t-(m-1)tb-(p-1)Tr]hhhhhhhhhhh×exp[j2π(n-1)(t-(p-1)Tr)Δf],
where P is the period number and Tr is pulse repetition interval. Consider an,m,p=ejφn,m,p, where φn,m,p is the mth phase element of the nth sequence in pth pulse. The rest of the parameters are the same as those of monopulse MCPC signal mentioned above.
3. Design of Multicarrier Chaotic Phase Coded Radar Pulse Train Signal
Chaos is considered as a phenomenon of random behavior in certainty nonlinear dynamical system [7]. Chaotic mapping can generate plenty of signals [8] whose rules are difficult to grasp by the interference side. Considering the balance, correlation, and the flatness of power spectrum in band, this study uses Chebyshev mapping [9] to generate the biphase coded sequence.
Chebyshev mapping is defined as
(3)xn+1=cos(karccos(xn)),
where x0∈(-1,1) and Chebyshev mapping order k=4.
The chaotic signal generated by chaotic maps is decimal. This decimal chaotic signal can be converted to binary sequence using the threshold comparison approach. Suppose Cn is the binary code after quantization, xn is the original decimal chaotic sequence, and x-=limN→∞(1/N)∑n=1Nxn is the mean of the sequence, which is the threshold. The binary code Cn is obtained by the following expression:
(4)Cn={1,xn≥x--1,xn<x-n=0,1,2,….
Because there are many parameters in MCPC pulse train signal, the phase code of the subcarrier in each pulse can be changed to obtain MCPC pulse train signals with different properties. Based on the differences of the phase codes of each pulse and subcarrier in MCPC pulse train signals, the classification of MCPC pulse train signal is made as shown in Table 1.
Classification of MCPC pulse train signal.
IIS MCPC
INS MCPC
NIS MCPC
NNS MCPC
Phase code of each pulse
Identical
Identical
Nonidentical
Nonidentical
Phase code of each subcarrier
Identical Sequence
Nonidentical Sequence
Identical Sequence
Nonidentical Sequence
Considering the features of the chaotic sequence mentioned above, we use the chaotic sequence to perform phase code modulation on the MCPC pulse train so as to generate chaotic modulation NNS MCPC pulse train signal. Because different initial values can result in chaotic sequences with different properties, all the chaotic sequences mentioned in the following refer to the chaotic biphase sequence with maximum main-to-sidelobe ratio as selected by the principle of autocorrelation maximization [10]. We design the following two chaotic modulated NNS MCPC pulse train signals with P periods, N subcarriers, and M bits through different chaotic modulation methods.
(1) The Chebyshev mapping generates L=P×N×M chaotic biphase codes, which are expressed in set form by {x1,x2,x3,…xL}. Suppose the phase coded set of MCPC pulse train signal is φn,m,p, where p=1,2,…P, n=1,2,…N, and m=1,2,…M. The phase of each subcarrier in MCPC pulse train is modulated by
(5)φn,m,p=xMN(p-1)+M(n-1)+m.
Therefore, the first type of chaotic modulation NNS MCPC pulse train is obtained and denoted by Chaos_NNS MCPC_I. Figure 2 shows the structure of Chaos_NNS MCPC_I pulse train.
Structure of Chaos_NNS MCPC_I pulse train.
(2) Each time a pulse is transmitted, the Chebyshev mapping generates L=N×M chaotic biphase codes with different initial values. The codes are expressed in set form by {x1,x2,x3,…xL}. Suppose the phase coded set of a single MCPC pulse is φn,m, where n=1,2,…N and m=1,2,…M. The phase of each subcarrier in the present transmitted pulse is modulated based on
(6)φn,m=xM(n-1)+m.
Therefore, the second type of chaotic modulation NNS MCPC pulse train is obtained and denoted by Chaos_NNS MCPC_II. The structure of Chaos_NNS MCPC_II pulse train is shown as Figure 3.
Structure of Chaos_NNS MCPC_II pulse train.
4. Derivation of Ambiguity Function
The ambiguity function is an important tool for analyzing radar signals. The ambiguity χ(τ,v) of a signal f(t) is given by
(7)χ(τ,υ)=∫-∞+∞f(t)f*(t+τ)exp(j2πvt)dt,
where τ is the time shift variable and v is the frequency shift. By substituting (2) into (7), the ambiguity function of NNS MCPC pulse train is obtained:
(8)χ(τ,v)=∫-∞+∞fTr(t)fTr*(t+τ)exp(j2πvt)dt=∫-∞+∞∑p=1P∑n=1N∑m=1Mwnan,m,phhhhhhhhhhhhhh×s[t-(m-1)tb-(p-1)Tr]hhhhhhhhhhhhhh×exp[j2π(n-1)(t-(p-1)Tr)Δf]×∑p1=1P∑n1=1N∑m1=1Mwn1*an1,m1,p1*hhhhhhhhhhhhhh×s*[t+τ-(m1-1)tb-(p1-1)Tr]hhhhhhhhhhhhhh×exp[-j2π(n1-1)hhhhhhhhhhhhhhhhhhhh×(t+τ-(p-1)Tr)Δf]hhhhhhhhhhhhhh×exp(j2πvt)dt=∑p=1P∑n=1N∑m=1M∑p1=1P∑n1=1N∑m1=1Mwnwn1*an,m,pan1,m1,p1*hhhhhhhhhhhhhhhhhhhh×exp{j2πTrΔfhhhhhhhhhhhhhhhhhhhhhhhhhh×[(p1-1)(n1-1)hhhhhhhhhhhhhhhhhhhhhhhhhhhhh-(p-1)(n-1)]}hhhhhhhhhhhhhhhhhhhh×exp[-j2π(n1-1)τΔf]×∫-∞+∞exp{j2π[(n-n1)Δf+v]t}hhhhhhhhh×s[t-(m-1)tb-(p-1)Tr]hhhhhhhhh×s*[t+τ-(m1-1)tb-(p1-1)Tr]dt.
Let (n-n1)Δf+v=Fd, t-(m-1)tb-(p-1)Tr=t′, and the integral in (8) can be simplified as
(9)∫-∞+∞exp{j2π[(n-n1)Δf+v]t}hhhh×s[t-(m-1)tb-(p-1)Tr]hhhhh×s*[t+τ-(m1-1)tb-(p1-1)Tr]dt=∫-∞+∞exp{j2πFd[t′+(m-1)tb+(p-1)Tr]}hhhhhh×s(t′)s*[t′+τ+(m-m1)tb+(p-p1)Tr]dt′=exp{j2πFd[(m-1)tb+(p-1)Tr]}×∫-∞+∞exp(j2πFdt)×s(t)hhhhhhhhh×s*[t+τ+(m-m1)tb+(p-p1)Tr]dt.
Since s(t)≡1, for 0≤t<tb, and zero, elsewhere, (8) is nonzero only when |τ+(m-m1)tb+(p-p1)Tr|<tb. Let τ+(m-m1)tb+(p-p1)Tr=ts; then (9) can be abbreviated as
(10)exp{j2πFd[(m-1)tb+(p-1)Tr]}×∫-∞+∞exp(j2πFdt)×s(t)s*(t+ts)dt.
If ts>0, the integral of (10) becomes
(11)∫0tb-tsexp(j2πFdt)dt=ej2πFdtj2πFd|0tb-ts=ejπFd(tb-ts)sin[πFd(tb-ts)]πFd;
if ts<0, the integral of (10) becomes
(12)∫-tstbexp(j2πFdt)dt=ej2πFdtj2πFd|-tstb=ejπFd(tb-ts)sin[πFd(tb+ts)]πFd.
According to (11) and (12), the result for the integral in (10) can be expressed as
(13)∫-∞+∞exp(j2πFdt)×s(t)s*(t+ts)dt=ejπFd(tb-ts)sin[πFd(tb-|ts|)]πFd,|τ+(m-m1)tb+(p-p1)Tr|<tb.
So the ambiguity function of NNS MCPC pulse train is
(14)χ(τ,v)=∑p=1P∑n=1N∑m=1M∑p1=1P∑n1=1N∑m1=1Mwnwn1*an,m,pan1,m1,p1*hhhhhhhhhhhhhhhhhhhh×exp{j2πTrΔf[(p1-1)(n1-1)hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh-(p-1)(n-1)]}hhhhhhhhhhhhhhhhhhhh×exp[-j2π(n1-1)τΔf]hhhhhhhhhhhhhhhhhhh×exp{j2πFd[(m-1)tbhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh+(p-1)Tr]}hhhhhhhhhhhhhhhhhhh×ejπFd(tb-ts)sin[πFd(tb-|ts|)]πFd,|τ+(m-m1)tb+(p-p1)Tr|<tb=χauto(τ,v)+χcross(τ,(n-n1)Δf+v),|τ+(m-m1)tb+(p-p1)Tr|<tb,
where
(15)χauto(τ,v)=∑p=1p=p1P∑n=1n=n1N∑m=1m=m1M|wn|2|an,m,p|2hhhhhhhhhhh×exp[-j2π(n-1)τΔf]hhhhhhhhhhh×exp{j2πv[(m-1)tb+(p-1)Tr]}hhhhhhhhhhh×ejπv(tb-ts)sin[πv(tb-|ts|)]πv,|τ|<tb,
is the auto-ambiguity function of Chaos_NNS MCPC pulse train and
(16)χcross(τ,(n-n1)Δf+v)=∑p=1p≠p1P∑n=1n≠n1N∑m=1m≠m1M∑p1=1P∑n1=1N∑m1=1Mwnwn1*an,m,pan1,m1,p1*hhhhhhhhhhhhhhhhhhhhh×exp{j2πTrΔfhhhhhhhhhhhhhhhhhhhhhhhhhhh×[(p1-1)(n1-1)hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh-(p-1)(n-1)]}hhhhhhhhhhhhhhhhhhhh×exp[-j2π(n1-1)τΔf]hhhhhhhhhhhhhhhhhhhh×exp{j2πFd[(m-1)tbhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh+(p-1)Tr]}hhhhhhhhhhhhhhhhh×ejπFd(tb-ts)sin[πFd(tb-|ts|)]πFd,hhhhhhhhhhhhhhh|τ+(m-m1)tb+(p-p1)Tr|<tb,
is the cross-ambiguity function of Chaos_NNS MCPC pulse train.
For Chaos_NNS MCPC_I pulse train, its ambiguity function is (14), where
(17)an,m,p=ejφn,m,p=ejxMN(p-1)+M(n-1)+m=ejcos(karccos(xMN(p-1)+M(n-1)+m-1)).
For Chaos_NNS MCPC_II pulse train, its ambiguity function is (14), where
(18)an,m,p=ejφn,m,p=ejxM(n-1)+m=ejcos(karccos(xM(n-1)+m-1)).
Furthermore, we note that different p has different initial values in Chaos_NNS MCPC_II pulse train.
5. Experimental Results5.1. Ambiguity Function
The simulation compares the ambiguity of the INS MCPC and NNS MCPC pulse train signals based on P4 code with that of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II proposed in this study. Suppose P=5, N=7, M=13, tb=M×10-6 s, and the duty ratio is 33%. The phase code of INS MCPC pulse train is the P4 code cyclic shift. Figure 4 shows the ambiguity function of INS MCPC based on P4 code.
Ambiguity function of INS MCPC based on P4 code.
A typical example of NNS MCPC pulse train signal is the COCS (consecutive ordered cyclic shift) P4 code pulse train [2], in which the modulation sequence of all subcarriers of each pulse and all pulses of the same subcarrier are P4 code cyclic shift. Figure 5 shows the ambiguity function of COCS P4 MCPC pulse train.
Ambiguity function of COCS P4 MCPC pulse train.
Since the phase codes in NNS MCPC pulse train signal are different, the periodic side lobe ratio of this signal at integer multiple of Tr is obviously lower than that obtained for INS MCPC pulse train signal.
Figure 6 shows the ambiguity function of the Chaos_NNS MCPC_I proposed in this paper with initial values 0.35.
Ambiguity function of Chaos_NNS MCPC_I.
Figure 7 shows the Chaos_NNS MCPC_II proposed in this paper. The initial values are [0.75,-0.52,-0.74,0.15,0.63].
Ambiguity function of Chaos_NNS MCPC_II.
It can be seen from Figures 6 and 7 that the two signals proposed in this paper possess a thumbtack ambiguity function since randomness is introduced into NNS MCPC pulse train signal in each subcarrier and pulse. Thus, it has high range and velocity resolution. The main lobe near the origin is narrow and the periodic side lobe at integer multiple of Tr is even lower than that obtained with P4 COCS MCPC pulse train signal. That means the correlation in each period attenuates due to randomness. Although Chaos_NNS MCPC_I and Chaos_NNS MCPC_II pulse train have about the same autocorrelation performance, Chaos_NNS MCPC_II increases the complexity and flexibility of signal design because of different chaotic initial values in each pulse, resulting in higher waveform diversity.
5.2. Analysis on Autocorrelation Performance
There are many parameters in multicarrier chaotic phase code pulse train. Changing the number of subcarriers, bits, and periods will affect the autocorrelation performance of the signal. The simulation result is as follows.
(1) Influence of Subcarrier Number on Autocorrelation. While the numbers of periods and bits are kept unchanged (P=5, M=13), the number of subcarriers is increased to three times (N=21) the original number (the sampling rate is also changed to three times the original value). The corresponding autocorrelation performance of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II is shown in Figure 8.
Influence of subcarrier number on autocorrelation.
(2) Influence of Bit Number on Autocorrelation. While the numbers of periods and subcarriers are kept unchanged (P=5, N=7), the numbers of bits is increased to three times (M=39) the original number. The corresponding autocorrelation performance of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II is shown in Figure 9.
Influence of bit number on autocorrelation.
(3) Influence of Period Number on Autocorrelation. While the numbers of subcarriers and bits are kept unchanged (N=7, M=13), the numbers of periods is increased to three times (P=15) the original value. The corresponding autocorrelation performance of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II is shown in Figure 10.
Influence of period number on autocorrelation.
Figures 8 and 9 show that the length of chaotic sequence increases whenever N or M grows. Thus, the randomness increases and the autocorrelation performance improves further. In addition, the autocorrelation of Chaos_NNS MCPC_II is more improved than that of Chaos_NNS MCPC_I. This is because, under the condition of a short chaotic sequence, increasing the length more obviously improves the randomness. It can be seen from Figure 10 that the chaotic sequence becomes longer as the number of periods increases, and the autocorrelation performance improves. However, due to the increase of period number, optimal initial values are running out of options, which make the autocorrelation performance of chaotic sequence become worse. Therefore, the improvement of autocorrelation performance of Chaos_NNS MCPC_II signal is not apparent. We can improve the signal autocorrelation performance under large period number by changing the type of chaotic mapping.
5.3. Analysis on LPI Performance
The intercept factor of radar [11] is written as
(19)α=A1TB,
where A is a parameter related to radar and reconnaissance receiver and can be regarded as a constant. TB denotes the time-bandwidth product. TB is an important factor that affects the LPI performance of radar. The larger, the better LPI performance is. B=N/tb and T=Mtb are the MCPC pulse train signal’s bandwidth and time width, respectively. So the time-bandwidth product is
(20)TB=N×M,
and intercept factor is
(21)α=A1N×M.
The LPI performance of MCPC pulse train can be improved by raising the numbers of subcarriers and bits. On the other hand, the power spectra of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II are shown in Figure 11. The figure suggests that the two signals have a flat power spectrum, making it difficult to detect by reconnaissance receivers using the pulse accumulation method [12], which raise the difficulty of being intercepted and identified by enemies. Moreover, because of the randomness of chaos and the complex modulation mode of the signal, the LPI performance is further improved.
Power spectra of Chaos_NNS MCPC_I and Chaos_NNS MCPC_II.
6. Conclusions
Radar waveform design is one of the main measures to raise the LPI performance of radar and is also an important aspect of radar system design. In this study, the chaotic sequence is introduced into the MCPC pulse train so as to design the Chaos_NNS MCPC_I and Chaos_NNS MCPC_II pulse train signal based on chaotic biphase modulation. We compare the ambiguity function of the signal proposed in this paper with that of MCPC pulse train modulated by P4 code and analyze the influences of the subcarrier number, bit number, and period number on the autocorrelation performance of the two signals. Finally, we discuss the LPI performance of the two signals. The multicarrier chaotic phase coded pulse train signal designed in this paper possesses thumbtack ambiguity function and low periodic autocorrelation side lobe. Also, it has a good LPI performance and the characteristic of waveform diversity, which enhances the detection ability and viability of radar.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Aeronautical Science Foundation of Aviation Industry Corporation of China under Grant no. 20112053018.
LevanonN.Multifrequency complementary phase-coded radar signal200014762762842-s2.0-003447548910.1049/ip-rsn:20000734LevanonN.MozesonE.2004New York, NY, USAJohn Wiley & SonsBinD.2011Changsha, ChinaNational University of Defense Technology(Chinese)SebtM. A.SheikhiA.NayebiM. M.Orthogonal frequency-division multiplexing radar signal design with optimised ambiguity function and low peak-to-average power ratio2009321221322-s2.0-6034911972910.1049/iet-rsn:20080106HuoK.JiangW.-D.LiX.MaoJ.-J.A new OFDM Phase-Coded Stepped-Frequency radar signal and its characteristic20113336776832-s2.0-7995550122710.3724/SP.J.1146.2010.00532SkolnikM. I.20033rdMcGrawHillBoccalettiS.GrebogiC.LaiY.-C.ManciniH.MazaD.The control of chaos: theory and applications200032931031972-s2.0-0034178680CarrollT. L.Adaptive chaotic maps for identification of complex targets2008242562622-s2.0-5004909569710.1049/iet-rsn:20070128KohdaT.TsunedaA.SakaeT.Chaotic binary sequences by Chebyshev maps and their correlation propertiesProceedings of the 2nd IEEE International Symposium on Spread Spectrum Techniques and Applications (ISSTA '92)19926366LindnerJ.Binary sequences up to length 40 with best possible autocorrelation function197511215072-s2.0-0016562524SchleherD. C.Low probability of intercept radarProceedings of the International Radar Conference1985346349BabuP.StoicaP.Spectral analysis of nonuniformly sampled data—a review20102023593782-s2.0-7634909603210.1016/j.dsp.2009.06.019