Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Real-time multi-target ranging based on chaotic polarization laser radars in the drive-response VCSELs

Open Access Open Access

Abstract

According to the principle of complete chaos synchronization and the theory of Hilbert phase transformation, we propose a novel real-time multi-target ranging scheme by using chaotic polarization laser radar in the drive-response vertical-cavity surface-emitting lasers (VCSELs). In the scheme, to ensure each polarization component (PC) of the master VCSEL (MVCSEL) to be synchronized steadily with that of the slave VCSEL, the output x-PC and y-PC from the MVCSEL in the drive system and those in the response system are modulated by the linear electro-optic effect simultaneously. Under this condition, by simulating the influences of some key parameters of the system on the synchronization quality and the relative errors of the two-target ranging, related operating parameters can be optimized. The x-PC and the y-PC, as two chaotic radar sources, are used to implement the real-time ranging for two targets. It is found that the measured distances of the two targets at arbitrary position exhibit strong real-time stability and only slight jitter. Their resolutions are up to millimeters, and their relative errors are very small and less than 2.7%.

© 2017 Optical Society of America

1. Introduction

Laser radar is generalized widely to apply in target detection, imaging, ranging, speed measurement and other fields [1–6]. The laser radars (lidars) for target ranging are mainly emitted from pulse laser and modulated continuous wave (CW) laser. In the short-pulse ranging techniques, the resolution is usually in the range of several meters, mainly determined by the pulse bandwidth [7, 8]. However, short pulse technology has many shortcomings, such as high manufacturing cost, bulky, heavy, low efficiency, and difficult to operate and maintain [9, 10]. In comparison, the modulated CW laser has many advantages such as compact size, low cost, high efficiency, easy operation and so on. However, the ranging resolution of the laser is in the range of tens of meters since that it is limited by the code rate and modulation speed [11-12]. To achieve a better resolution, an expensive external modulator has to be used for high-speed modulation. The chaotic radar ranging has many advantages over short-pulse and CW laser ranging, such as low probability of intercept, high range and velocity resolution, strong anti-interference ability, easy generation and low cost [13–24]. And it plays an important role in modern radar research.

The chaotic lidar (CLIDAR) generated by using the nonlinear dynamics of semiconductor laser has a higher resolution and is easy to generate and amplify, benefiting from the broad bandwidth of the optical chaos [13, 14]. By using the CLIDAR, the resolution of centimeter magnitude can be easily achieved. The CLIDAR for target ranging has an obvious advantage over CW lidar, such as that it does not need the generator for high-rate code and the modulator with high speed. In 2001, based on the semiconductor laser with optical feedback, Myneni K et al. firstly proposed a CLIDAR ranging scheme [4], which expands a new application of optical chaos in measuring range and velocity of target. In the next few years, Jia-Ming Liu and his associates considered a chaotic laser as a light source, and proposed a novel chaotic radar system [19]. In the system, using the correlation between the signal wave fed back from the measured target and a delayed reference wave, they realized the ranging for target, and achieved stably resolution in centimeter. In the same year, they experimentally achieved the target ranging, and the resolution of the ranging is of 9 cm [20]. In 2008, considering two plane mirrors as two-target, Bingjie Wang and his associates demonstrated the reliability of real-time ranging of the two targets by using CLIDAR, and implemented a resolution of 9 cm between two targets [21]. In 2015, they experimentally explored the target detection and ranging technology in loss media [22]. Moreover, Pengfei Du et al. experimentally realized the remote target ranging by using CLIDAR [23]. So far, as we know, most of CLIDAR ranging schemes are realized by using the correlation between the signal waveform and the reference waveform [13–24]. In these schemes, the ranging resolution can reach several centimeters, and be further improved by using the correlation principle if laser chaotic bandwidth can be further enhanced. Recently, some schemes for laser chaotic bandwidth enhancement have been experimentally and theoretically implemented [25–32]. For example, the chaos bandwidth of 14GHz is experimentally generated by optical heterodyning of two external-cavity laser diodes [27-28]. Optical chaos extending up to 26.5 GHz bandwidth has been achieved by injecting the chaos light from a DFB laser with conventional feedback into a fiber ring resonator [29]. The chaos bandwidth can be extended up to 20GHz by using optical injection [30], and further enlarged up to 32.3 GHz using dual wavelength injection [31]. Replacing conventional optical feedback with phase-conjugate feedback can improve the chaos bandwidth. In the range of achievable phase-conjugate mirror reflectivities, the bandwidth can be further enlarged up to about 12GHz, and its increase reaches 27% [32]. However, how much range resolution can be improved by enhancing the chaotic bandwidth needs to be further explored. On the other hand, due to the interference of the spontaneous emission noise and channel noise, the quality of the correlation can be weakened, which seriously affect the accuracy of the target ranging. Based on the principle of complete chaos synchronization (CCS), these problems can be solved when the chaotic radar in a system of the drive-response lasers is applied in target ranging. The reason is that the resolution of target ranging mainly depends on the CCS quality of the drive-response lasers. The distance of the target at arbitrary position can be measured under steady CCS. It is noted that the CCS quality of the drive-response lasers has very strong robustness to noise [33-34]. The complete synchronization can be used to play the role of the noise filter. Therefore, the noise has little effect on target ranging. However, in the system of the drive-response VCSELs with the spin flip mode (SFM) [33–35], since that the random anticorrelated hopping between the two orthogonal polarization components (PC) emitted by each VCSEL is induced by the feedback [36] or the injection current [37], the simultaneously chaotic synchronization for each polarization component (PC) is extremely unstable and very poor, making the resolution to be reduced seriously, and the real-time stability to be poor in target ranging. In our previous work [38-39], we proposed a novel scheme to control the CCS quality in the drive-response VCSELs by electro-optic (EO) modulation, where the CCS quality of each PC appears periodic variation and can be stably controlled simultaneously when the applied electric field is fixed at a certain value. Motivated by these, according to the principle of CCS, the theory of Hilbert phase transform and linear electro-optical (EO) effect, we present a novel scheme of real-time multi-target ranging by using chaotic polarization lidar (CPLIDAR) in the system of the drive-response VCSELs, where two PCs can be achieved stability and well CCS by using linear EO effect in periodically poled LiNbO3 (PPLN) crystal. Under the completely synchronized condition, we further explore the resolution, the real-time stability, and the relative errors of two target ranging.

2. Theory and model

Figure 1 presents the scheme diagram of multi-target ranging by using CPLIDAR in the drive-response VCSELs. Here, the upper part is the drive system; the lower part is the response system; the master VCSEL (MCSEL) is drive laser; the slave VCSEL (SVCSEL) is response laser. In the drive system, with the bias current and the feedback strength fixed at certain value, the MVCSEL emits two chaotic PCs [see Eqs. (1) - (3)]. Since that their polarized direction along the direction of the x-axis and the y-axis, they are, respectively, defined as the x-PCM and the y-PCM, where the subscript M means the MVCSEL. The feedback cavity for the y-PCM is composed of the mirror 1 (M1), M2 and M6 (FC1); the feedback cavity for the x-PCM is formed by the M3-M5 (FC2). The microwave signal m1 (m2) is modulated to the x-PCM (y-PCM) by the MOD1 (MOD2), where the MOD is EO modulator. The neutral density filter 1 (NDF1) and the NDF2 are used to adjust the optical feedback strength. The transmitting x-polarization optical antenna (TXPOA) is used to transmit the chaotic x-PCM with m1 (the probe signal 1). The chaotic y-PC with m2 (the probe signal 2) is transmitted by the transmitting y-polarization optical antenna (TYPOA). In the response system, the chaotic x-PCM’ and y-PCM’ are from the receiving x-polarization optical antenna (RXPOA) and the receiving y-polarization optical antenna (RYPOA), respectively. The optical amplifier 1 (OA1) and the OA2 are used to compensate the loss of the probe signal 1and the probe signal 2 in the antennas and the channels. Under certain conditions, the SVCSEL can emit two PCs, which are defined as the x-PCS and the y-PCS [see Eqs. (4) - (6)]. Here, the subscript S denotes the SVCSEL.

 figure: Fig. 1

Fig. 1 Scheme diagram of the multi-target ranging by using CPLIDAR based on the drive-response VCSELs. Here, M: mirror; PPLN: periodic poled LiNbo3; IS: isolator; FR: Faraday rotator; HWP: half wave plate; MVCSEL: master VCSEL; SVCSEL: slave VCSEL; um and us: the bias current of MVCSEL and SVCSEL, respectively; E1 and E2: applied electric field; BS: beam splitter; PBS: polarization beam splitter; m1 and m2: modulated microwave signal;m1 andm2: demodulated microwave signals; MOD: modulator; TYPOA: transmitting y-polarization optical antenna; TXPOA: transmitting x-polarization optical antenna; RXPOA: receiving x-polarization optical antenna; RYPOA: receiving y-polarization optical antenna; NDF: neutral density filter; DEL: delayer; F: filter; HT: Hilbert transform; PD: photodiode; CD: distance calculation; DIV: divider; SUB: subtractor; OA: optical amplifier.

Download Full Size | PDF

In the system, how the CPLIDAR signals are transmitted and received in space needs to be further explained as follows. So far, typically, the CPLIDAR signals are transmitted and received by radio antenna. Based on this traditional approach, each chaotic PC should be converted to radio signal by one photodiode in the drive system. To realize the CCS, in the response system, the radio signal reflected from each target should be further converted to the CPLIDAR signal that is almost similar to that from the M-VCSEL, which is very difficultly implemented in technology. Recently, using the surface plasmon theory, several researchers have experimentally and theoretically implemented a variety of optical nanoantennas that could efficiently convert confined optical energy to free-space light, and vice versa [40–43]. Some polarized optical nanoantennas have been further proposed in experiment and theory [42-43]. Therefore, it is possible to solve the above-mentioned problem by using optical antenna. For this purpose, we choose optical antenna to transmit and receive CPLIDAR signal. On the other hand, to avoid the mutual interference between the two probe signals when they are transmitted in space, and received by the RXPOA and RYPOA simultaneously, we consider the TXPOA and the TYPOA to transmit the probe signals 1 and 2, respectively, and use the RXPOA and the RYPOA to receive two reflected probe signals from the T1 and T2, respectively [see Fig. 1]. Since that the polarization of the two probe signals are perpendicular to each other, the mutual interference between them transmitting in space can be suppressed. Moreover, the RXPOA and the RYPOA can isolate the probe signal 2 and the probe signal 1, respectively, owing to the reason that the polarization of the two receiving antennas are perpendicular to each other.

In the following, we consider the probe signal 1 as an example to illustrate the implementation of the target 1 (T1) ranging. The probe signal 1 is split into two beams by the beam splitter 1 (BS1). One of two beams is firstly fed back by the FC2, and then modulated by linear EO effect in the PPLN1 crystal, finally injected into the MVCSEL. The other one is transmitted to the T1 by the TXPOA. After being reflected or scattered by the T1, it is received by the RXPOA. The probe signal 1 from the TXPOA is split into two beams by the BS3 when the channel noise in it is filtered by the filter 1 (F1). One beam of the probe signal 1, firstly undergoes EO modulation in the PPLN2, then injected into the SVCSEL. The chaotic wave emitted by the SVCSEL is decomposed into the x-PCS and the y-PCS by polarization beam splitter (PBS). The x-PCM’ can be completely synchronized with the x-PCS by the control of the applied electric field. The signal m1 can be decoded from the synchronous division between the x-PCM’ and the x-PCS by the divider1 (DIV1), where the decoded signal is named as m1. The phase of m1 and m1 are extracted by the Hilbert phase transform. The distance of the T1 can be calculated based on their phase difference [see Eqs. (21)-(22)]. Similarly, the ranging for the T2 can also be implemented. In the scheme, the propagating delay time of the probe signals 1 and 2 are defined as τc1 and τc2, respectively. We obtain τc2>τc1 when the distance of the T2 is more than that of the T1. However, under this condition, the two PCs cannot be completely synchronized. To achieve their complete synchronizations, we place a delayer between the BS3 and the RXPOA. The probe signal 1 is delayed for time τT by the delayer. According to the theory of the CCS [33] in the drive-response lasers, τc2=τc1+τT is considered here. And τc=τc2 is set for the convenience of discussion.

In 1995, Miguel et al. put forward the spin flip model (SFM) of VCSEL [44], where the laser cavity is designed into cylindrical symmetry. Consequently the output of VCSEL may be two orthogonal PCs (the x-PC and the y-PC) due to the weak material and cavity anisotropies, whose dynamic behaviors are sensitive to the polarization of the laser field. In particular, the VCSEL with optical feedback or optical injection may generate the two orthogonal chaotic PCs [36-37], the chaotic behavior of the x- PC or the y-PC becomes more complex due to the linear phase anisotropy and the amplitude anisotropy. So the chaotic behaviors of each PC of VCSEL output are different from that of the single polarization mode of the edge emitting laser output. According to the SFM [44], the rate equations of the MVCSEL and the SVCSEL need to be modified when they are subject to the injection of the output two PCs from the MVCSEL modulated by linear EO effect [45]. For the MVCSEL, the rate equations are

ddt(EMx(t)EMy(t))=k(1+ia){[NM(t)1]}(EMx(t)EMy(t))±k(1+ia)inM(t)(EMy(t)EMx(t))(ra+irp)(EMx(t)EMy(t))+kf(EMx(tτ)EMy(tτ))×exp(iω0τ)[1+m1(tτ)1+m2(tτ)]+(βspγeNMζxβspγeNMζy),
dNM(t)dt=re{NM(t)μM+NM(t)(|EMx(t)|2+|EMy(t)|2)+inM(t)[EMy(t)EMx*(t)EMx(t)EMy*(t)]},
dnM(t)dt=rsnM(t)re{nM(t)(|EMx(t)|2+|EMy(t)|2)+iNM(t)[EMy(t)EMx*(t)EMx(t)EMy*(t)]}.
For the SVCSEL, the rate equations can be described by the following equations:
ddt(ESx(t)ESy(t))=k(1+ia){[NS(t)1]}(ESx(t)ESy(t))±k(1+ia)inS(t)(ESy(t)ESx(t))(ra+irp)(ESx(t)ESy(t))+kinj(Epx(tτc)Epy(tτc))×(exp(iω0τc)exp(iω0τc))+(βspγeNζxβspγeNζy),
dNS(t)dt=re{NS(t)μS+NS(t)(|ESx(t)|2+|ESy(t)|2)+inS(t)[ESy(t)ESx*(t)ESx(t)ESy*(t)]},
dnS(t)dt=rsnS(t)re{nS(t)(|ESx(t)|2+|ESy(t)|2)+iNS(t)[ESy(t)ESx*(t)ESx(t)ESy*(t)]},
where the subscripts x and y represent the x-PC and the y-PC, respectively; E is the complex amplitude of light field; N is the total carrier concentration; n is the difference in concentration between carriers with spin-up and spin-down; k=1/(2τp),τp is the carrier photon lifetime; γeis the nonradiative carrier relaxation rate; a is the line width enhancement factor; γsis the spin relaxation rate; γaand γp is the linear dichroism and the linear birefringence, respectively; τ is the round-trip time in the external cavity; τc is the propagation time of light from the MVCSEL to the SVCSEL; kf is the optical feedback strength; kinj is the optical injection strength; ω0 is the center frequency of the MVCSEL and the SVCSEL; L is the length of the PPLN1 or the PPLN2; μM and μS are the normalized bias current of the MVCSEL and the SVCSEL, respectively; the noise strength parameter D is defined as γeNβSP, with βSP being the spontaneous emission factor; ξx and ξy are two independent complex Gaussian white noise events with zero mean and 1 variance, and their time correlation is that ξi*(t)ξj*(t)2δijδ(tt); m1 and m2 both are the modulated microwave signals. In addition, EMx(tτ) and EMY(tτ), respectively, are the complex amplitude of the output x-PCM and the y-PCM from the PPLN1. In the output of the PPLN2 crystal, their complex amplitudes are Epx(tτ) and Epy(tτc), respectively. As shown in Fig. 1, the polarization direction of the x-PCM is along the direction of the o-light in the PPLN1, the y-PC is aligned with the e-light by the FR1 and the HWP1. Under these conditions, considering the delay time probe signals 1 and 2 as the original inputs of the o-light and the e-light, respectively, we obtain
Ux,y(0,tτ)=ω0VSATLνcn1,2EMx,My(tτ)(1+m1,2(tτ)).
Similarly, in the PPLN2, for the propagating delay probe signals 1 and 2, we have
Ux,y(0,tτc)=ω0VSATLνcn1,2EMx,My(t)(1+m1,2(tτc)),
where Ux and Uy are the amplitude of the o-light and the e-light, respectively; ħ is the Planck constant; SA is the effective area of the light spot; V is the volume of the active layer of the VCSEL; υc is the light velocity in a vacuum; TL=2ngvc/Lv is the round trip time in the laser cavity, Lv is the length of the laser cavity, ng is the effective refractive index of the laser active layer; n1 and n2 are the undisturbed refractive indices of the o-light and the e-light, respectively. With the phase mismatch and the weak second-order nonlinear effect, the analytical solutions of the wave-coupling equations of the linear EO effect for the two PCs in the PPLN1 and the PPLN2 are written as [45-46]:
Ux,y(L,t)=ρx,y(L,t)exp(iβ0L)exp[iϕx,y(L,tt0)],
where t0=τ or τc, and
ρx,y(L,tt0)={Ux,y2(0,tt0)cos2(νL)+[γUx,y(0,tt0)d1,3Uy,x(0,tt0)ν]2sin2(νL)}12.
ϕx,y(L,tt0)=tan1[±γUx,y(0,tt0)d1,3Uy,x(0,tt0)vUx,y(0,tt0)tan(vL)].
with
β0=Δkd2d42,
ν=(Δk+d2d4)2+4d1d32,
γ=d4d2Δk2,
d1=k02n1n2reff1E0f1,
d2=k02n1reff2E0f0,
d3=k02n1n2reff1E0f1,
d4=k02n2reff3E0f0,
where L is the length of the PPLN1 or the PPLN2; the wave vector mismatch Δk=kxky+K1, K1=2π/Λ is the first-order reciprocal lattice vector of the PPLN, Λ is the poled period, kx=2πn1vc/ω0 and ky=2πn2vc/ω0 denote the wave vector of the o-light and the e-light at ω0, respectively. k0 is the wave vector of light in vacuum. Here, the K1 is considered to be very close to the wave vector mismatch kxky. Those components that make little contribution to the EO effect are neglected because of the phase mismatch. The effective EO coefficients in Eqs. (13)a) and (13d) are given as follow.
reff1=j,k,l(εijεkk)ajrjklbkcl,
reff2=j,k,l(εijεkk)ajrjklakcl,
reff3=j,k,l(εijεkk)bjrjklbkcl,
where j,k,l=1,2,3 (the same below); εij=njj2 and εkk=nkk2 both are the diagonalizable electric permittivity tensor elements; rijk is the EO tensor elements; a and b are the unit vector of the o-light and the e-light, respectively; c is the unit vector of the applied electric field E; a=(sinφ,cosφ,0) and bb=(cosθcosφ,cosθ,sinφ,sinθ) because that PPLN is a uniaxial crystal, where θ and φ are the polar angle and the azimuth angle, respectively; c = (0, 1, 0) when the applied electric field E0 is along the direction of the z-axis in the crystal. In addition, f0=2R1is the zero-order coefficient of the Fourier series of the crystal structure function; f+1=[1cos(±2πR)+isin(±2πR)]/(±iπ) are the first-order positive and negative Fourier coefficients, respectively; R=i+/(i++i) is the duty radio, i+is the length of the positive domain, i is the length of the negative domains. When the MVCSEL are subject to the injection of the output x-PCM and y-PCM from the PPLN1 crystal, we have
EMx,My(tτ)=SATLνcn1,2ω0VUx,y(L,tτ).
Similarly, while the output x-PCM and y-PCM from the PPLN2 crystal are injected into the SVCSEL, we obtain

Epx,py(tτc)=SATLνcn1,2ω0VUx,y(L,tτc).

3. Results and discussions

In the system of the drive-response VCSELs, the accuracy and the stability of target ranging mainly depends on the CCS quality between the x-PCM’ (y-PCM’) from -the MVCSEL output and the x-PCS (y-PCS) from the SVCSEL output. To describe the CCS quality, the correlation functions of the two PCs are introduced, and defined as follows.

ρx,y=[IMx,My(tΔt)IMx,My(tΔt)][ISx,Sy(tΔt)ISx,Sy(tΔt)]{[IMx,My(tΔt)IMx,My(tΔt)]2[[ISx,Sy(tΔt)ISx,Sy(tΔt)]2]}1/2,
Where IMx,My(t)=|EMx,My(t)|2 and ISx,Sy(t)=|ESX,SY(t)|2; the symbol < > denotes the time average (the same below); Δt=τcτ. The coefficient ρ ranges from 0 to 1. With the bigger ρ, the two PCs have higher synchronization quality. The case that ρ = 1 denotes that the two chaotic PCs can be achieved complete synchronization when [33].
ESx,Sy(t)=EMx,My(tΔt).
Therefore, the two demodulation signals are
m1,2'(t-Δt)=(1+m1,2(tΔt))EMx,My(tΔt)ESx,Sy(t)1.
We consider two modulated signals as sine waves, namely,
m1,2(t)=Asin(ωt+ϕ0),
where ω is the signal angular frequency, A is the amplitude. Under CCS [see Eq. (17)], the two demodulation signals are derived from Eqs. (17) - (18) as follows.
m1,2'(tΔt)=Asin[ω(tΔt)+ϕ1,2(t)].
Here, ϕ1(t) and ϕ2(t) are the phase of m1' and m2'; ϕ = ϕ1 = ϕ2 since that τc2=τc1+τT. The analytic signal of m1'(t) and m2'(t) are firstly calculated by the following equation
ψ1,2(tΔt)=m1,2'(tΔt)+jm˜1,2'(tΔt)=Aejϕ(t),
where m˜1'(t) and m˜2'(t) are the Hilbert transform of m1'(t) and m2'(t). From Eq. (21), the phase of m1' and m2'are obtained as follows.
ϕ(t)=arctanm˜1,2'(tΔt)m1,2'(tΔt).
From Eq. (22), the delay time Δt is written as
Δt=Δϕ(t)ω,
where Δϕ(t)=ϕ(t)ϕ0. Combined with Δt=τcτ and Eq. (23), the distance of the T1 and the T2 are derived as

d1=(Δt+ττT)c2,
d2=(Δt+τ)c2.

Since that the accuracy and the real-time stability of the target ranging depends on the CCS quality of the two PCs, how to choose the reasonable ranges of the crucial parameters is key to making the synchronization of the two PCs to be high quality and strong stability. Here, the key parameters include the bias current u, the applied electric field E, the propagating delay τc, the injection strength kinj and the feedback strength kf. For this purpose, we numerically solve Eqs. (1) - (3) and (4) - (6) by using the four-order Runge-Kutta method. The numerical values of the parameters in calculation are given in Tab. 1. We firstly calculated the maps of the correlation coefficients (ρx and ρy) evolution in the parameter space of u and τc [see Figs. 2(a) and 2(b)], as well as in that of E and τc [see Figs. 2(c) and 2(d)], where kf = kinj = 2 ns−1. In Figs. 2 (a) and 2(b), E is fixed at 0.36kV/mm. u = 1.5 is given in Figs. 2(c) and 2(d). One sees from Figs. 2(a) and 2(b) that, with the increase of u, ρx and ρy appear fast oscillational change between 0 and 1 when u ranges from 1 to 1.2, denoting that the chaotic synchronization between the x-PCM’ and the x-PCS (CS1) and that between the y-PCM’ and the y-PCS (CS2) are both poor. However, with the increase of u from 1.2 to 2, ρx and ρy are identically equal to 1, indicating that the CS1 and CS2 can be steady and good. For example, as shown in Figs. 3(a) and 3(b), ρx=ρy=1 when u = 1.5. At this time, the two PCs are achieved CCS; but for u = 1.18, ρx oscillates between 0.06 and 0.98, and ρy ranges from 0.33 to 0.98; if u is fixed at 1.15, the variation range of ρx is from 0.87 to 1, and ρy jitters between 0.87 and 1. These discussions indicate that the CS1 and CS2 display extremely unstable when u is less than 1.2. The reason is that the symmetry between the rate equations of the drive laser and those of the response one is destroyed because of the simultaneous existence of the two PCs [33, 36-37]. From Figs. 2(c) and 2(d), it is found that, for an arbitrary τc, ρx and ρy vary periodically with E. For example, when E changes in a certain range, such as 0.1~0.2kV/mm, 0.4~0.5kV/mm and 0.7~0.8kV/mm, ρx and ρy oscillate violently between 0 and 1. Under these conditions, the CS1 and CS2 are also jittered violently and very poor. On the other hand, while E varies in some limited ranges, such as 0.2~0.4kV/mm, 0.5 ~0.7kV/mm and 0.8 ~1.0 kV/mm, ρx and ρyidentically equal to 1, indicating that the two PCs can be achieved CCS. Figures 3(c) and 3(d) further describe the real-time change of ρx and ρy under different E. It is concluded from Figs. 2 and 3 that the two PCs can be obtained stable and well synchronized quality when the bias current u and the applied electric field E are controlled in the appropriate range. This can be attributed to the case that the imbalance in the symmetry between the rate equations of the drive laser and those of the response one is offset because of the exchange of energy between the two PCs by EO modulation.

Tables Icon

Table 1. Numerical values for the calculation of the target ranging

 figure: Fig. 2

Fig. 2 Maps of the correlation coefficients ρxand ρy evolution in the parameter space of u and τc, as well as in that of E and τc.

Download Full Size | PDF

 figure: Fig. 3

Fig. 3 Dependence of the correlation coefficients ρx and ρy on the propagating delay time τcunder different parameters such as E and u, where kf=kinj=2ns1 and E = 0.36kV/ mm in (a) and (b); u = 1.5 in (c) and (d).

Download Full Size | PDF

When τc are fixed at 5ns and 7ns, respectively, we give the chaotic behaviors of the output x-PC and y-PC from the MVCSEL and those from the SVCSEL, as shown in Figs. 4 and 5, respectively. Here, E = 0.36kV/mm, u = 1.5. Their time series and spectrums are displayed in the left column and the right column of these two figures, respectively. Their powers are, respectively, defined as follows:

PMx,y=10log10(|EMx,My|2|EMx|2+|EMy|2),PSx,y=10log10(|ESx,Sy|2|ESx|2+|ESy|2).
As shown in Figs. 4 and 5, for the MVCSEL and the SVCSEL, the time traces and the power spectrum distributions of the output two PCs are both in chaotic state. The output x-PCS and y-PCS have, respectively, the same chaotic dynamic behaviors as the delay output x-PCM (the x-PCM’) and the output y-PCM with delay (the y-PCM’), which indicates that each PC of the MVCSEL is completely synchronized with that of the SVCSEL by selecting appropriate E and μ. Moreover, the chaos bandwidth of the x-PCM’ and the y- PCM’ are about 13GHz, which are the same as and that of the x-PCS and the y- PCS. These chaos bandwidths can be further enhanced by increasing feedback strength and injection strength, which can be beneficial to improve the modulation rate of m1 and m2. However, according to the Eqs. (17) - (25), the range resolution mainly depends on the CCS quality of each PC, and is independent of its chaotic bandwidth in the drive-response synchronization system.

 figure: Fig. 4

Fig. 4 Chaotic behavior of the output x-PC and y-PC from the M-VCSEL and those from the SVCSEL with τc=5ns, E = 0.36kV/mm, u = 1.5. The left column: time series; the right column: power spectrum distribution.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 Chaotic behavior of the output x-PC and y-PC from the M-VCSEL and those from the VCSEL when τc=7ns, E = 0.36kV/mm, u = 1.5. The left column: time series; The right column: power spectrum distribution.

Download Full Size | PDF

In the following discussions for chaotic radar ranging, to ensure the two PCs to be in a stable state of complete synchronization, we take E and μ as 0.36kV/mm and 1.5, respectively. Under CCS, for the verification of the reliability and accuracy of chaotic radar ranging for the targets T1 and T2, two sets of the actual distances of the T1 and the T2 are considered as follows: LT1=0.45m and LT2=0.75m; LT1=0.75m and LT2=1.05m. From these distances, we obtain two sets of the propagating delay times as follows: τc1=3ns and τc2=5ns; τc1=5ns and τc2=7ns. The delay time τTin the delayer is considered as 2ns to ensure that the two PCs both have the same propagating delay. In addition, two sinusoidal signals have the same angular frequency, i.e., ω1=ω2=ω, where ω=2π/4(τcτ).

For the above-mentioned two sets of different target distances, Fig. 6 presents the plots of ϕ0 and ϕ, as well as d1 and d2 versus time, where ϕ0 is the initial phase of m1 and m2; ϕ is the phase of m1’ and m2’; d1 and d2 are the measured distance of the T1 and T2, respectively. From Figs. 6(a)-6(d), it is found that ϕ0 and ϕ both show a linear increase in periodicity when the T1 and T2 locate in different distances, respectively. The difference (Δϕ) between ϕ0 and ϕ almost remains unchanged at any time. By using Δϕ from Figs. 6(a) - 6(d) and Eqs. (23) - (25), we further calculate the time traces of d1 and d2 for two sets of different target positions, as shown in Figs. 6(e) and 6(f). One sees from Fig. 6(e) that d1 is slightly jittered between 0.4392 ~0.4572m when LT1=0.45m; IfLT1=0.75m, it varies in a small range of from 0.7488 to 0.7524m. For LT2=0.75m, d2 fluctuates in small range such as 0.7392~0.7527m; if LT2=1.05m its variation range is from 1.0492 to 1.0513m. These results indicate that the two chaotic PCs have very high real-time stability in target distance measurement under CCS, and the ranging resolution reaches millimeter magnitude.

 figure: Fig. 6

Fig. 6 For the targets T1 and T2 at different positions, time traces of the phases ϕ0 and ϕ, as well as the measured distances d1 and d2 when the output two PCs are in complete chaotic synchronization state. Here, (a): ϕ0 and ϕ for LT1 = 0.45m; (b) ϕ0 and ϕ for LT1 = 0.75m; (c): ϕ0 and ϕ for LT2 = 0.75m; (d): ϕ0 and ϕ for LT2 = 1.05m; (e): d1 for LT1 = 0.45m (the blue solid line), and d1 for LT1 = 0.75m (the red solid line); (f): d2 for LT2 = 0.75m (the red dotted line), and d2 for LT2 = 1.05m (the blue dotted line).

Download Full Size | PDF

In addition, the relative errors RE1 and RE2 are introduced to describe the accuracy of the target ranging, and defined as follows:

RE1,2=|<d1,2>LT1,2|LT1,2×100%,
where the symbol “| |” indicates absolute value. We quantize the relative errors, making them easier for observation. State 1 indicates that RE1 and RE2 are no more than 5%; state 2: 5%≤RE1,2≤ 20%; state 3: 20%≤RE1,2≤ 50%; State 4: 50% < RE1,2≤80%; state 5: RE1,2>80%.

Figure 7 gives the maps of RE1 and RE2 evolutions with τc in the parameter space of μ and τc . [see Figs. 7(a) and 7(b)], as well as in the parameter space of E and τc [see Figs. 7(c) and 7(d)]. As shown in Figs. 7(a) and 7(b), if E is fixed at 0.36kV/mm, since that the two PCs have been achieved CCS [see Figs. 2(c)-2(d)], RE1 and RE2 are almost no more than 5% when μ≥1.5 and τc ≥4ns,. With μ fixed at 1.5, RE1 and RE2 vary periodically with E if τc ≥4ns. For example, RE1 and RE2 are more than 20% while E varies in several limited ranges, such as 0.38-0.5, and 0.68-0.8, due to the deterioration of the CS1 and CS2 [see Figs. 2(c)-2(d)]; they are almost no more than 5% if E exists in some regions, such as 0.18-0.38, 0.48-0.68, and 0.78-0.98, owing to the good quality of the CS1 and CS2 [see Figs. 2(c)-(d)]. These mean that the small relative errors can be obtained by controlling E.

 figure: Fig. 7

Fig. 7 Maps of the relative errors RE1 and RE2 evolutions with the propagating time τc in the parameter space of μ and τc [(a) and (b)], where E = 0.36kV/mm. The maps of their evolutions in the parameter space of E and τc [(c) and (d)], where μ=1.5 ; (e)-(h): the partial enlargement of (a)-(d) in turn.

Download Full Size | PDF

Figure 8 further present the dependences of the relative errors, the signal-to-noise ratios (SNRs), and the correlation coefficients on the propagating time under μ=1.5 and E=0kV/mm. Here, the SNR1 and SNR2 are defined as

SNR1,2=10log10(PS1,2/Pn1,2),
where Ps1 and Ps2 are, respectively, expressed as

 figure: Fig. 8

Fig. 8 Dependences of the relative errors (RE1 and RE2), the signal-to noise ratios (SNR1 and SNR2), and the correlation coefficients (ρx and ρy) on the propagating time τc, where μ=1.5, E=0kV/mm.

Download Full Size | PDF

Ps1,2=[EMx,My(tΔt)/ESx,Sy(t)]2.

Pn1 and Pn2 are the corresponding spontaneous radiated noise power, and written as

Pn1,2=(Dξx,y)2.
It is found from Fig. 8 that if E=0kV/mm, ρx and ρy show irregular oscillation between 0 and 1, which indicates that the CS1 and CS2 are very unstable and seriously deteriorated. Moreover, the SNR of two PCs appear extremely unstable and violent oscillating [see Figs. 8(c) and 8(d)]. For example, SNR1 oscillates between −33.175dB and 5.182 dB, SNR2 varies irregularly between −30.955dB and −26.711dB, denoting that the SNR of the x-PC has serious fluctuation, while the fluctuation range for the SNR of the y-PC is relatively small since that the x-PC has large synchronization error and the y-PC possesses smaller one. Under these conditions, the variation range of RE1 is from 0 to 126.954%, and that of RE2 is between 0 and 100.555%. However, ρx and ρy identically equal to 1 when E = 0.36kV/mm [see Figs. 9 (a) and 9(b)], indicating that the two PCs are achieved stable CCS. Under this condition, their SNRs are low and only slight jitter [see Figs. 9(c) and 9(d)]. For example, the SNR1 is slightly jittered between −32.963dB and −33.043dB; the SNR2 varies slightly between −30.883dB and −30.949dB. As a result, the fluctuation range of RE1 and RE2 are greatly reduced. Their variation ranges are 0~2.696% and 0~2.661%, respectively [see Figs. 9(e) and 9(f)].

 figure: Fig. 9

Fig. 9 Dependences of the relative errors (RE1 and RE2), the signal-to noise ratios (SNR1 and SNR2), and the correlation coefficients (ρx and ρy) on the propagating time τc, where E = 0.36kV/mm and μ=1.5.

Download Full Size | PDF

From [33–35], it is found that the CCS quality of each PC depends not only on the parameters discussed above, but also on other key parameters of the system such as γP, γs, kf, king and τ. Therefore, the variation of these four different parameters can seriously influence the relative errors (RE1 and RE2). Figure 10 displays the evolutions of RE1 and RE2 as well as ρx and ρy with these parameters under E = 0.36kV/mm, τc = 7ns, μ=1.5 and kf = kinj. As shown from Figs. 10(a1) and 10(a2), there appear very low relative errors when γP is small. But large relative errors will occurs if γP is large. For example, if γP varies from 1ns−1 to 2.837ns−1, RE1 and RE2 are both small and become 1.787%, due to the reason that ρx and ρy are identically equal to 1. With the increase of γP from 2.837ns−1 to 6.327ns−1, RE1 fluctuates between 0.1217% and 97.43%, due to the fluctuation range of ρx between 0.679 and 1. RE2 irregularly varies between 0.1615% and 109.05%, owing to the variation range of ρy from 0.7569 to 1. If γP increases from 6.327ns−1 to 10 ns−1, RE1 and RE2 both exhibit violently fluctuation since that ρx and ρy fluctuates over a wide range. It is found from Figs. 10(b1) and 10(b2) that the relative errors have less values and vary in small range when increases from 50ns−1 to 100ns−1. For example, RE1 varies in a small range of from 1.177% to γs 1.944%, and RE2 changes from 0.6059% to 3.195% when γs increases from 50ns−1 to 100ns−1, owing to the reason that ρx and ρy are almost close to 1. These indicate that the relative errors have strong robustness when γs is large. From Figs. 10(c1) and 10(c2), one sees that if kf varies from 1ns−1 to 1.5ns−1, RE1 fluctuates between 0.4505% and 63.54%, and RE2 oscillates between 0.4193% and 53.34% because that ρx and ρy change between 0.3636 and 1. When kf increases from 1.5ns−1 to 2.796 ns−1, RE1 and RE2 almost equal to 1.787% since that ρx and ρy both are of 1. With kf further increasing from 1.5ns−1 to 3ns−1, RE1 and RE2 begin to fluctuate again. These show that the large relative errors occur in the region of the smaller values of kf or that of the larger values of kf. As further seen from Figs. 10(d1) and 10(d2), RE1 and RE2 are small and less than 2.7% when ρx and ρy equal to 1 in the parameter space of τ. However, although ρx and ρy are slightly reduced, but also RE1 and RE2 have very large values if τ is no more than 4ns. For example, with ρx being 0.9958 and 0.9988, RE1 is 28.34% and 21.11% when τ is fixed at 1.102ns and 3.99ns, respectively. But if τ is fixed at 5.796ns, RE1 is small and of 0.3281% when ρx equals to 0.9985. This denotes that the relative errors for these two targets are very sensitive and weak robust to the CCS if short delay occurs in the master laser feedback loop.

 figure: Fig. 10

Fig. 10 Evolutions of the correlation coefficients (ρx and ρy), and the relative errors (RE1 and RE2) with the four different key parameters when E = 0.36kV/mm, τc = 7ns, μ=1.5 and kf = kinj. Here, (a1): ρx and ρyαγP ; (a2): RE1 and RE2αγP; (b1): ρx and ρyαγs ; (b2): RE1 and RE2 αγs; (c1): ρx and ρyα kf; (c2): RE1 and RE2 αkf; (d1): ρx and ρyατ ; (d2): RE1 and RE2 ατ.

Download Full Size | PDF

From the above results, it is concluded that the relative errors of the ranging for the two targets are seriously dependent on the quality of the CS1 and CS2. When the CS1 and CS2 exhibit unstable and are deteriorated, the ranging accuracies for the two targets are very low and their stability are very poor. It is noted that the two PCs can be achieved steady CCS by using EO modulation and selecting the appropriate system parameters. Under this condition, the ranging accuracy of the target at any position is greatly improved, and its stability is further strengthened. Furthermore, two chaotic PCs both have high and stable SNR.

4. Conclusions

We propose a novel real-time multi-target ranging scheme by using CPLIDAR in the drive-response VCSELs, based on the EO modulation theory, complete chaotic synchronization theory and Hilbert transform principle. Here, we discuss the effects of some crucial parameters, such as the bias current, the applied electric field, and the propagating delay, the feedback or injection strength, the spin relaxation, the birefringence, the delay in the master feedback loop, on the chaotic synchronization quality and stability. By using EO modulation and selecting the system parameters reasonably, the two chaotic polarization components with arbitrary propagating delay can be achieved steady complete synchronization. Under this condition, the real-time stability for the target ranging is fairly strong and only slight jitter. The resolution of the target ranging reaches millimeter magnitude. It is important to note that the relative error of the target at any position is very small and less than 2.7%. These results have potential application in the real-time multi-target ranging with high resolution and anti-noise jamming.

Funding

National Natural Science Foundation of China (NSFC) (Grant No. 61475120); Innovative Projects in Guangdong Colleges and Universities, China (Grant No. 2015KTSCX146); and the Youth Foundation of Wuyi university (No. 2014zk08).

References and links

1. Y. N. Ji and M. J. Zhang, “Ultra-wideband microwave-photonic chaotic radar for remote ranging,” International Conference on Optical Communication. 121(1), 1–3 (2013). [CrossRef]  

2. M. U. Piracha, D. Nguyen, D. Mandridis, T. Yilmaz, I. Ozdur, S. Ozharar, and P. J. Delfyett, “Range resolved lidar for long distance ranging with sub-millimeter resolution,” Opt. Express 18(7), 7184–7189 (2010). [CrossRef]   [PubMed]  

3. M. C. Amann, T. Bosch, M. Lescure, R. Myllylä, and M. Rioux, “Laser ranging: a critical review of usual techniques for distance measurement,” Opt. Eng. 40(1), 10–19 (2001). [CrossRef]  

4. K. Myneni, T. A. Barr, B. R. Reed, S. D. Pethel, and N. J. Corron, “High precision ranging using a chaotic laser pulse train,” Appl. Phys. Lett. 78(11), 1496–1498 (2001). [CrossRef]  

5. F. Pellen, P. Olivard, Y. Guern, J. Cariou, and J. Lotrian, “Radio frequency modulation on an optical carrier for target detection enhancement in sea-water,” J. Phys. D Appl. Phys. 34(7), 1122–1130 (2001). [CrossRef]  

6. P. Weibring, T. Johansson, H. Edner, S. Svanberg, B. Sundnér, V. Raimondi, G. Cecchi, and L. Pantani, “Fluorescence lidar imaging of historical monuments,” Appl. Opt. 40(33), 6111–6120 (2001). [CrossRef]   [PubMed]  

7. G. N. Pearson, P. J. Roberts, J. R. Eacock, and M. Harris, “Analysis of the performance of a coherent pulsed fiber lidar for aerosol backscatter applications,” Appl. Opt. 41(30), 6442–6450 (2002). [CrossRef]   [PubMed]  

8. A. Rybaltowski and A. Taflove, “Signal-to-noise ratio in direct-detection mid-infrared random-modulation continuous-wave lidar in the presence of colored additive noise,” Opt. Express 9(8), 389–399 (2001). [CrossRef]   [PubMed]  

9. D. M. He and G. G. L. Seet, “Underwater lidar imaging (UWLI) system in short turbid water tank,” Proc. SPIE 4416, 272–275 (2001). [CrossRef]  

10. M. P. Bristow, “Suppression of afterpulsing in photomultipliers by gating the photocathode,” Appl. Opt. 41(24), 4975–4987 (2002). [CrossRef]   [PubMed]  

11. N. Takeuchi, N. Sugimoto, H. Baba, and K. Sakurai, “Random modulation CW lidar,” Appl. Opt. 22(9), 1382–1386 (1983). [CrossRef]   [PubMed]  

12. N. Takeuchi, H. Baba, K. Sakurai, and T. Ueno, “Diode-laser random-modulation CW lidar,” Appl. Opt. 25(1), 63–67 (1986). [CrossRef]   [PubMed]  

13. F. Y. Lin and J. M. Liu, “Nonlinear dynamics of a semiconductor laser with delayed negative optoelectronic feedback,” IEEE J. Quantum Electron. 39(4), 562–568 (2003). [CrossRef]  

14. J. S. Lawrence and D. M. Kane, “Injection locking suppression of coherence collapse in a diode laser with optical feedback,” Opt. Commun. 167(1-6), 273–282 (1999). [CrossRef]  

15. T. Deng, G. Q. Xia, J. J. Chen, X. Tang, X. D. Lin, X. Yang, S. W. Hua, and Z. M. Wu, “Experimental investigation on nonlinear dynamics of 1550nm VCSEL simultaneously subject to orthogonal optical injection and negative optoelectronic feedback,” Laser Phys. 27(4), 045402 (2017). [CrossRef]  

16. J. J. Chen, Z. M. Wu, L. Fan, X. Tang, X. D. Lin, T. Deng, and G. Q. Xia, “Polarization bistability in a 1550nm vertical-cavity surface-emitting laser subject to variable polarization optical injection,” IEEE Photonics J. 9(2), 1502309 (2017). [CrossRef]  

17. J. J. Chen, Z. M. Wu, T. Deng, X. Tang, X. Yang, and G. Q. Xia, “Current-and feedback-induced state bistability in a 1550nm-VCSEL with negative optoelectronic feedback,” IEEE Photonics J. 9(1), 1500310 (2017). [CrossRef]  

18. F. Alonge, M. Branciforte, and F. Motta, “A novel method of distance measurement based on pulse position modulation and synchronization of chaotic signals using ultrasonic radar systems,” IEEE Trans. Instrum. Meas. 58(2), 318–329 (2009). [CrossRef]  

19. F. Y. Lin and J. M. Liu, “Chaotic lidar,” IEEE J. Sel. Top. Quantum Electron. 10(5), 991–997 (2004). [CrossRef]  

20. F. Y. Lin and J. M. Liu, “Chaotic radar using nonlinear laser dynamics,” IEEE J. Quantum Electron. 40(6), 815–820 (2004). [CrossRef]  

21. B. J. Wang, Y. C. Wang, L. Q. Kong, and A. B. Wang, “Multi-target real-time ranging with chaotic laser radar,” Chin. Opt. Lett. 6(11), 868–870 (2008). [CrossRef]  

22. B. J. Wang, H. Xu, P. Yang, L. Liu, and J. X. Li, “Target detection and ranging through lossy media using chaotic radar,” Entropy (Basel) 17(4), 2082–2093 (2015). [CrossRef]  

23. P. F. Du, D. X. Geng, W. Wang, and M. L. Gong, “Laser detection of remote targets applying chaotic pulse position modulation,” Opt. Eng. 54(11), 114102 (2015). [CrossRef]  

24. X. Ai, R. Nock, J. G. Rarity, and N. Dahnoun, “High-resolution random-modulation cw lidar,” Appl. Opt. 50(22), 4478–4488 (2011). [CrossRef]   [PubMed]  

25. R. M. Nguimdo, M. Khoder, J. Danckaert, G. Van der Sande, and G. Verschaffelt, “Fast phase response and chaos bandwidth enhancement in semiconductor lasers subject to optical feedback and injection,” Opt. Lett. 39(20), 5945–5948 (2014). [CrossRef]   [PubMed]  

26. A. Quirce, A. Vally, and H. Thienpont, “Enhancement of chaos bandwidth in VCSELs induced by simultaneous orthogonal optical injection and optical feedback,” IEEE J. Quantum Electron. 52(10), 1–9 (2016). [CrossRef]  

27. A. B. Wang, B. J. Wang, L. L. Yuancai, and K. Alan Shore, “Optical heterodyne generation of high-dimensional and broadband white chaos,” IEEE J. Sel. Top. Quantum Electron. 21(6), 531–540 (2015). [CrossRef]  

28. C. H. Cheng, Y. C. Chen, and F. Y. Lin, “Chaos time delay signature suppression and bandwidth enhancement by electrical heterodyning,” Opt. Express 23(3), 2308–2319 (2015). [CrossRef]   [PubMed]  

29. A. Wang, Y. Wang, Y. Yang, M. Zhang, H. Xu, and B. Wang, “Generation of flat-spectrum wideband chaos by fiber ring resonator,” Appl. Phys. Lett. 102(3), 031112 (2013). [CrossRef]  

30. A. Uchida, T. Heil, Y. Liu, P. Davis, and T. Aida, “High-frequency broad-band signal generation using a semiconductor laser with a chaotic optical injection,” IEEE J. Quantum Electron. 39(11), 1462–1467 (2003). [CrossRef]  

31. M. Zhang, T. Liu, P. Li, A. Wang, J. Zhang, and Y. Wang, “Generation of broadband chaotic laser using dual-wavelength optically injected fabry-pérot laser diode with optical feedback,” Photo. Technol. Lett. IEEE 23(24), 1872–1874 (2011). [CrossRef]  

32. É. Mercier, D. Wolfersberger, and M. Sciamanna, “High-frequency chaotic dynamics enabled by optical phase-conjugation,” Sci. Rep. 6(1), 18988 (2016). [CrossRef]   [PubMed]  

33. I. Gatare, M. Sciamanna, A. Locquet, and K. Panajotov, “Influence of polarization mode competition on the synchronization of two unidirectionally coupled vertical-cavity surface-emitting lasers,” Opt. Lett. 32(12), 1629–1631 (2007). [CrossRef]   [PubMed]  

34. D. Z. Zhong and Z. M. Wu, “Complete chaotic synchronization mechanism of polarization mode of VCSEL with anisotropic optical feedback,” Opt. Commun. 282(8), 1631–1639 (2009). [CrossRef]  

35. Z. D. Zhong, T. Deng, and G. L. Zheng, “Manipulation of the complete chaos synchronization in dual-channel encryption system based on polarization division multiplexing,” Wuli Xuebao 63(7), 070504 (2014).

36. M. Sciamanna, K. Panajotov, H. Thienpont, I. Veretennicoff, P. Mégret, and M. Blondel, “Optical feedback induces polarization mode hopping in vertical-cavity surface-emitting lasers,” Opt. Lett. 28(17), 1543–1545 (2003). [CrossRef]   [PubMed]  

37. M. Virte, K. Panajotov, H. Thienpont, and M. Sciamanna, “Deterministic polarization chaos from a laser diode,” Nat. Photonics 7(1), 60–65 (2012). [CrossRef]  

38. D. Zhong, Y. Ji, and W. Luo, “Controllable optoelectric composite logic gates based on the polarization switching in an optically injected VCSEL,” Opt. Express 23(23), 29823–29833 (2015). [CrossRef]   [PubMed]  

39. D. Z. Zhong, W. Luo, and G. L. Xu, “Controllable all-optical stochastic logic gates and their delay storages based on the cascaded VCSELs with optical-injection,” Chin. Phys. B 25(9), 094202 (2016). [CrossRef]  

40. L. Yousefi, “Highly directive hybrid plasmonic leaky wave optical nanoantenna,” Prog. Electromagnetics Res. 50, 85–90 (2014). [CrossRef]  

41. M. A. Panahi, L. Yousefi, and M. Shahabadi, “Highly directive hybrid plasmonic leaky-wave optical antenna with controlled side-lobe level,” J. Lightwave Technol. 33(23), 4791–4798 (2015). [CrossRef]  

42. G. Rui, R. L. Nelson, and Q. Zhan, “Circularly polarized unidirectional emission via a coupled plasmonic spiral antenna,” Opt. Lett. 36(23), 4533–4535 (2011). [CrossRef]   [PubMed]  

43. Y. B. Xie, Z. Y. Liu, Q. J. Yang, G. H. Sun, X. J. Zhang, and Y. Y. Zhu, “Controlling the state of polarization via optical nanoantenna feeding with surface plasmon polaritons,” Appl. Phys. Lett. 108(13), 131102 (2016). [CrossRef]  

44. Q. Feng, J. V. Moloney, and J. V. Moloney, “Light-polarization dynamics in surface-emitting semiconductor lasers,” Phys. Rev. A 52(2), 1728–1739 (1995). [CrossRef]   [PubMed]  

45. W. L. She and W. K. Lee, “Wave coupling theory of linear electrooptic effect,” Opt. Commun. 195(1), 303–311 (2001). [CrossRef]  

46. G. L. Zheng, H. C. Wang, and W. L. She, “Non-collinear quasi-phase-matched linear electro-optic effect in periodically poled LiNbO3 and its applications,” J. Opt. A, Pure Appl. Opt. 10(1), 015102 (2008). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1
Fig. 1 Scheme diagram of the multi-target ranging by using CPLIDAR based on the drive-response VCSELs. Here, M: mirror; PPLN: periodic poled LiNbo3; IS: isolator; FR: Faraday rotator; HWP: half wave plate; MVCSEL: master VCSEL; SVCSEL: slave VCSEL; um and us: the bias current of MVCSEL and SVCSEL, respectively; E1 and E2: applied electric field; BS: beam splitter; PBS: polarization beam splitter; m1 and m2: modulated microwave signal; m 1 and m 2 : demodulated microwave signals; MOD: modulator; TYPOA: transmitting y-polarization optical antenna; TXPOA: transmitting x-polarization optical antenna; RXPOA: receiving x-polarization optical antenna; RYPOA: receiving y-polarization optical antenna; NDF: neutral density filter; DEL: delayer; F: filter; HT: Hilbert transform; PD: photodiode; CD: distance calculation; DIV: divider; SUB: subtractor; OA: optical amplifier.
Fig. 2
Fig. 2 Maps of the correlation coefficients ρ x and ρ y evolution in the parameter space of u and τ c , as well as in that of E and τ c .
Fig. 3
Fig. 3 Dependence of the correlation coefficients ρx and ρy on the propagating delay time τ c under different parameters such as E and u, where k f = k inj =2n s 1 and E = 0.36kV/ mm in (a) and (b); u = 1.5 in (c) and (d).
Fig. 4
Fig. 4 Chaotic behavior of the output x-PC and y-PC from the M-VCSEL and those from the SVCSEL with τ c =5ns, E = 0.36kV/mm, u = 1.5. The left column: time series; the right column: power spectrum distribution.
Fig. 5
Fig. 5 Chaotic behavior of the output x-PC and y-PC from the M-VCSEL and those from the VCSEL when τ c =7ns, E = 0.36kV/mm, u = 1.5. The left column: time series; The right column: power spectrum distribution.
Fig. 6
Fig. 6 For the targets T1 and T2 at different positions, time traces of the phases ϕ0 and ϕ, as well as the measured distances d1 and d2 when the output two PCs are in complete chaotic synchronization state. Here, (a): ϕ0 and ϕ for L T 1 = 0.45m; (b) ϕ0 and ϕ for L T 1 = 0.75m; (c): ϕ0 and ϕ for L T 2 = 0.75m; (d): ϕ0 and ϕ for L T 2 = 1.05m; (e): d1 for L T 1 = 0.45m (the blue solid line), and d1 for L T 1 = 0.75m (the red solid line); (f): d2 for L T 2 = 0.75m (the red dotted line), and d2 for L T 2 = 1.05m (the blue dotted line).
Fig. 7
Fig. 7 Maps of the relative errors RE1 and RE2 evolutions with the propagating time τ c in the parameter space of μ and τ c [(a) and (b)], where E = 0.36kV/mm. The maps of their evolutions in the parameter space of E and τ c [(c) and (d)], where μ=1.5 ; (e)-(h): the partial enlargement of (a)-(d) in turn.
Fig. 8
Fig. 8 Dependences of the relative errors (RE1 and RE2), the signal-to noise ratios (SNR1 and SNR2), and the correlation coefficients ( ρ x and ρ y ) on the propagating time τ c , where μ=1.5, E=0kV/mm.
Fig. 9
Fig. 9 Dependences of the relative errors (RE1 and RE2), the signal-to noise ratios (SNR1 and SNR2), and the correlation coefficients ( ρ x and ρ y ) on the propagating time τ c , where E = 0.36kV/mm and μ=1.5.
Fig. 10
Fig. 10 Evolutions of the correlation coefficients ( ρ x and ρ y ), and the relative errors (RE1 and RE2) with the four different key parameters when E = 0.36kV/mm, τc = 7ns, μ=1.5 and kf = kinj. Here, (a1): ρ x and ρ y α γ P ; (a2): RE1 and RE2 α γ P ; (b1): ρ x and ρ y α γ s ; (b2): RE1 and RE2 α γ s ; (c1): ρ x and ρ y α kf; (c2): RE1 and RE2 αkf; (d1): ρ x and ρ y ατ ; (d2): RE1 and RE2 ατ.

Tables (1)

Tables Icon

Table 1 Numerical values for the calculation of the target ranging

Equations (38)

Equations on this page are rendered with MathJax. Learn more.

d dt ( E Mx (t) E My (t) )=k( 1+ia ){[ N M (t)1]}( E Mx (t) E My (t) )±k( 1+ia )i n M (t)( E My (t) E Mx (t) ) ( r a +i r p )( E Mx (t) E My (t) )+ k f ( E Mx (tτ) E My (tτ) ) ×exp(i ω 0 τ)[ 1+ m 1 (tτ) 1+ m 2 (tτ) ]+( β sp γ e N M ζ x β sp γ e N M ζ y ),
d N M (t) dt = r e { N M (t) μ M + N M (t)(| E Mx (t) | 2 +| E My (t) | 2 )+i n M (t) [ E My (t) E Mx * (t) E Mx (t) E My * (t)]},
d n M (t) dt = r s n M (t) r e { n M (t)(| E Mx (t) | 2 +| E My (t) | 2 )+i N M (t) [ E My (t) E Mx * (t) E Mx (t) E My * (t)]}.
d dt ( E Sx (t) E Sy (t) )=k( 1+ia ){[ N S (t)1]}( E Sx (t) E Sy (t) )±k( 1+ia )i n S (t) ( E Sy (t) E Sx (t) )( r a +i r p )( E Sx (t) E Sy (t) )+ k inj ( E px (t τ c ) E py (t τ c ) ) ×( exp(i ω 0 τ c ) exp(i ω 0 τ c ) )+( β sp γ e N ζ x β sp γ e N ζ y ),
d N S (t) dt = r e { N S (t) μ S + N S (t)(| E Sx (t) | 2 +| E Sy (t) | 2 )+i n S (t) [ E Sy (t) E Sx * (t) E Sx (t) E Sy * (t)]},
d n S (t) dt = r s n S (t) r e { n S (t)(| E Sx (t) | 2 +| E Sy (t) | 2 )+i N S (t) [ E Sy (t) E Sx * (t) E Sx (t) E Sy * (t)]},
U x,y (0,tτ)= ω 0 V S A T L ν c n 1,2 E Mx,My (tτ)(1+ m 1,2 (tτ)).
U x,y (0,t τ c )= ω 0 V S A T L ν c n 1,2 E Mx,My (t)(1+ m 1,2 (t τ c )),
U x,y (L,t)= ρ x,y (L,t)exp(i β 0 L)exp[i ϕ x,y (L,t t 0 )],
ρ x,y (L,t t 0 )= { U x,y 2 (0,t t 0 ) cos 2 (νL) + [ γ U x,y (0,t t 0 ) d 1,3 U y,x (0,t t 0 ) ν ] 2 sin 2 (νL) } 1 2 .
ϕ x,y (L,t t 0 )= tan 1 [ ±γ U x,y (0,t t 0 ) d 1,3 U y,x (0,t t 0 ) v U x,y (0,t t 0 ) tan(vL) ].
β 0 = Δk d 2 d 4 2 ,
ν= (Δk+ d 2 d 4 ) 2 +4 d 1 d 3 2 ,
γ= d 4 d 2 Δk 2 ,
d 1 = k 0 2 n 1 n 2 r eff1 E 0 f 1 ,
d 2 = k 0 2 n 1 r eff2 E 0 f 0 ,
d 3 = k 0 2 n 1 n 2 r eff1 E 0 f 1 ,
d 4 = k 0 2 n 2 r eff3 E 0 f 0 ,
r eff1 = j,k,l ( ε ij ε kk ) a j r jkl b k c l ,
r eff2 = j,k,l ( ε ij ε kk ) a j r jkl a k c l ,
r eff3 = j,k,l ( ε ij ε kk ) b j r jkl b k c l ,
E Mx,My (tτ)= S A T L ν c n 1,2 ω 0 V U x,y (L,tτ).
E px,py (t τ c )= S A T L ν c n 1,2 ω 0 V U x,y (L,t τ c ).
ρ x,y = [ I Mx,My (tΔt) I Mx,My (tΔt) ][ I Sx,Sy (tΔt) I Sx,Sy (tΔt) ] { [ I Mx,My (tΔt) I Mx,My (tΔt) ] 2 [ [ I Sx,Sy (tΔt) I Sx,Sy (tΔt) ] 2 ] } 1/2 ,
E Sx,Sy (t)= E Mx,My (tΔt).
m 1,2 ' (t-Δt)= (1+ m 1,2 (tΔt)) E Mx,My (tΔt) E Sx,Sy (t) 1.
m 1,2 (t)=Asin(ωt+ ϕ 0 ),
m 1,2 ' (tΔt)=Asin[ω(tΔt)+ ϕ 1,2 (t)].
ψ 1,2 (tΔt)= m 1,2 ' (tΔt)+j m ˜ 1,2 ' (tΔt)=A e jϕ(t) ,
ϕ(t)=arctan m ˜ 1,2 ' (tΔt) m 1,2 ' (tΔt) .
Δt= Δϕ(t) ω ,
d 1 = (Δt+τ τ T )c 2 ,
d 2 = (Δt+τ)c 2 .
P Mx,y =10 log 10 ( | E Mx,My | 2 | E Mx | 2 + | E My | 2 ) , P Sx,y =10 log 10 ( | E Sx,Sy | 2 | E Sx | 2 + | E Sy | 2 ).
R E 1,2 = | < d 1,2 > L T 1,2 | L T 1,2 ×100%,
SN R 1,2 =10 log 10 ( P S 1,2 / P n 1,2 ),
P s 1,2 = [ E Mx,My (tΔt)/ E Sx,Sy (t)] 2 .
P n 1,2 = (D ξ x,y ) 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.