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

Optimal signal constellation design for ultra-high-speed optical transport in the presence of nonlinear phase noise

Open Access Open Access

Abstract

In this paper, we first describe an optimal signal constellation design algorithm suitable for the coherent optical channels dominated by the linear phase noise. Then, we modify this algorithm to be suitable for the nonlinear phase noise dominated channels. In optimization procedure, the proposed algorithm uses the cumulative log-likelihood function instead of the Euclidian distance. Further, an LDPC coded modulation scheme is proposed to be used in combination with signal constellations obtained by proposed algorithm. Monte Carlo simulations indicate that the LDPC-coded modulation schemes employing the new constellation sets, obtained by our new signal constellation design algorithm, outperform corresponding QAM constellations significantly in terms of transmission distance and have better nonlinearity tolerance.

© 2014 Optical Society of America

1. Introduction

The fiber nonlinearities can be considered as one of the most important factor in coherent long-haul optical transmission system [1]. Due to fiber Kerr nonlinear effects, the refractive index of optical fiber increases with optical intensity to slightly slow down the propagation speed, inducing intensity depending nonlinear phase shift, which includes the nonlinear phase noise (NLPN) because of the noise component in optical intensity introduced by amplifier [2]. In order to improve the performance of the long-haul optical transport systems with coherent detection, the nonlinear phase noise should be considered in order to increase data rates and/or extend the transmission distance of optical transmission systems.

The optimal design of a signal constellation, i.e., placing M constellation points in the complex plane such that symbol error probability (SEP) or a specific criterion is minimized under an average or peak power constraint, is a classical problem in communication theory [3]. However, only few results are known about the constellation design for optical channel with nonlinear effects. A maximum likelihood (ML) detector for phase-shift keying modulated signals and two-stage detector for 16-QAM constellation has been proposed by Lau and Kahn [4]. Beygi has proposed 16-point ring constellation sets for combating the effect of NLPN [5]. Also, design of APSK constellations for coherent optical channels with NLPN is proposed by Hager in [6]. Although these methods can achieve better SEP performance, there is no consideration about an overall algorithm used to design the optimal signal constellation in the presence of phase noise, applicable to arbitrary signal constellation size.

In this paper, a new optimal constellation design algorithm is proposed that is applicable to coherent detection system in the presence of either linear or nonlinear phase noise. As an extension of the optimum signal constellation design (OSCD), which is introduced in [7], the proposed algorithm uses the cumulative log-likelihood function instead of the Euclidian distance, which is applicable to both linear and nonlinear phase noise-dominating scenarios. The optimum source distribution used in the algorithm is generated by maximizing the channel capacity based on well-known Arimoto-Blahut algorithm [8]. In this scheme, the proposed algorithm can interact with channel impairments and log-likelihood function calculation by using the same criterion. Also, an LDPC coded modulation scheme, representing the generalization of [9], suitable for use in combination with constellations obtained by this algorithm is proposed. Monte Carlo simulations indicate that the signal constellations obtained by the proposed algorithm outperforms corresponding QAM counterparts significantly in both linear and nonlinear phase noise dominated scenarios. Notice that in [7] the constellations obtained for ASE noise dominated scenario are evaluated in the channel dominated by ASE noise. In this paper, we design the signal constellations for the phase noise dominated channel, which is a key difference. Simulation results indicate that new constellations we obtained outperform previous constellations for ASE noise dominated scenario and significantly outperform QAM constellations in both linear or nonlinear phase noise dominated channel models.

The paper is organized as follows. In section 2, we introduce the channel model suitable for both linear and nonlinear phase noise dominated channels. In section 3, we present the proposed optimal signal constellation design algorithm for these two channels. In section 4, we introduced the cumulative log-likelihood function and our proposed coded modulation scheme employing the new signal constellation sets obtained by algorithm introduced in section 3. The results of the Monte Carlo simulation are summarized in section 5. Concluding remarks can be found in section 6.

2. Channel model

In this section, we introduce the equivalent coherent optical channel model suitable for both linear and nonlinear phase noise channels. The phase noise dominated channel modeling is a well-investigated topic in optical transmission [10], [11]. Also, the statistical description of the residual error after carrier recovery is introduced in [12], [13]. In this paper, the model we used [6] corresponds to one particular situation, when either the carrier phase estimation (CPE) is imperfect or the nonlinear phase noise is imperfectly compensated for. The proposed signal constellation design algorithm is applicable to the phase noise channel models [10]- [13] as well.

2.1 Coherent optical channel dominated by the linear phase noise

The equivalent channel model for coherent detection, upon compensation of linear and nonlinear impairments and CPE, can be represented as:

rk=s(ak,θk)+zk,    rk=[rk(1)rk(i)rk(N)]T
s(ak,θk)=ejθk[ak(1)ak(i)ak(N)]T,  zk=[zk(1)zk(i)zk(N)]T
where rk(i) is the ith component of observation vector at kth symbol interval, ak(i) is the ith coordinate of transmitted symbol (at kth symbol interval), and zk is the corresponding noise vector (with Gaussian-like distribution of components). The θk denotes the residual phase error (at kth time instance) due to laser phase noise, nonlinear phase noise, and imperfect CPE. N denotes the number of components in the observation vector and j denotes the imaginary unit. In polarization-division multiplexing (PDM), Eqs. (1) and (2) correspond to either x- or y- polarization state. In 4-D signaling, the components above represent projections along in-phase and quadrature basis functions corresponding to x- and y- polarization. The model above is applicable to few-mode fiber (FMF) applications as well. For instance, to describe the laser phase noise and imperfect CPE, the Wiener phase noise model can be used:
θk=(θk1+Δθk) mod 2π
where Δθk is zero-mean Gaussian process of variance δΔθ=2πΔfTs. The Ts denotes the symbol duration and Δf denotes either linewidth or frequency offset. The phase slips can also be modeled by Markov-like process of certain memory. The probability density function (PDF) of the phase increment in Eq. (3) is given by
pΔΘ(Δθk)=n=p(0,δΔθ2,Δθkn2π)
where p(0,δΔθ2,Δθkn2π) denotes the Gaussian PDF of zero-mean, variance δΔθ2, and argument Δθkn2π. The resulting noise process is Gaussian-like with the power spectral density of N0 so that the corresponding conditional probability density function is given by:
pR(r|ak,θk)=e||rks(ak,θk)||2N0/(πN0) 
For non-Gaussian channels, we will need to use the method of histograms to estimate the conditional probability density function pR(r|ak,θk).

2.2 Coherent optical channel dominated by the nonlinear phase noise

In this case, our channel model is nonlinear phase noise model with discrete amplification for the finite number of fiber spans. When the optical signal is periodically amplified by EDFAs, the nonlinear phase noise is unavoidably added to the optical signal and accumulated as the number spans increases. For convenience, we consider the discrete memoryless channel model, which is introduced in [1] and can be described as follows:

Y=(X+Z)ejΦNL
where XX is the channel input, Z is the total additive noise, and Y is the channel observation. In each fiber span, the overall nonlinear phase shift ΦNL is given by
ΦNL=0LγP(z)dz=γLeffP 
where P is the launch power and γ is the nonlinear Kerr-parameter. For a fiber span length of L with attenuation coefficient of α, the power evolution is described as P(z)=PeαL and the effective length is defined as
Leff=1eαLα
For a system with NA fiber spans, the overall nonlinear phase noise is given by:
ΦNL= γLeff{|E0+n1|2+|E0+n1+n2|2++|E0+n1++nNA|2}
where E0 is the baseband representation of the transmitted electric field, nk is independent identically distributed zero-mean circular Gaussian random complex variable with variance δ02. The total additive noise at the end of all fiber segments has the variance δ2E[Z2]=2NAδ02 and can be calculated as [14]
δ2=2nsphναΔνNA L 
where all parameters are summarized in the Table 1 [4].

Tables Icon

Table 1. Parameters of the System under Study

In this channel model, the variance of the phase noise is dependent on the channel input and the channel is specified by the number of spans, transmission length, and the launch power. In Fig. 1, we show the received constellation diagrams for 16-QAM constellation with K=25, L=2000Km for different powers P.

 figure: Fig. 1

Fig. 1 Constellation diagrams of received signal with 16-QAM constellation, L = 2000km, for different launch powers.

Download Full Size | PDF

3. Optimal signal constellation design for phase noise dominated channels

In the presence of phase noise, we can use an algorithm similar to OSCD algorithm [7], but now changing the optimization criterion. Instead of minimizing the mean square error, we define the cumulative log-likelihood function and get the optimal signal constellation that maximizes this function. Namely, the Euclidean distance receiver is optimum only for the AWGN channel. For the linear or nonlinear phase error dominated channel, the Arimoto-Blahut algorithm can be applied to find the optimal source distribution. The resulting OSCD constellation optimized from the training sequences generated from the optimal source distribution will make the input symbols to approximately have similar optimal source distribution, as the channel model is Gaussian plus phase noise. Our proposed algorithm, named here as a log-likelihood ratio-based OSCD (LLR-OSCD) algorithm, can be formulated as follows.

  • 0) Initialization: Choose the signal constellation that will be used for initialization of size M.
  • 1) Generate the training sequence from the optimum source distribution, denoted as {xj; j=0,  , n1}.
  • 2) Group the samples from this sequence into M clusters. The membership to the cluster is determined based on LLR of sample point and candidate signal constellation points from previous iteration. Each sample point is assigned to the cluster with the largest LLR. Given the m th subset (cluster) with N candidate constellation points, denoted as A^m={yi ; i=1,  , N}, find the log-likelihood (LL) function of partition P(A^m)={Si ; i=1,  , N}, as follows
    LLm=LL({A^m,P(A^m)})=n1k=0n1maxyA^mLL(xk,y)

    The function LL(xj,y) represents the cumulative log-likelihood function defined as

    LL(xk,y)=1NSi=1NS{xk1Re[(y1+y2j)ej×PNi]}2+{xk2Im[(y1+y2j)ej×PNi]}22δ2

    where NS denotes the number of phase noise samples and the corresponding phase noise sample is denoted as PNi. xk1 and xk2 denote the first and second coordinates of the point xk. Similarly, y1 and y2 denote the coordinates of the received point y. The equation above is written by having in mind a two-dimensional (2D) signal constellation design, but it can straightforwardly be generalized for arbitrary dimensionality.

  • 3) Determine the new signal constellation points as center of mass for each cluster.

    Repeat the steps 2)-3) until convergence is achieved.

The PNi can be generated by target channel model introduced in Section 2 for both linear and nonlinear phase noise case. That is to say, if PNi is generated by linear phase noise model, then the resulting constellation is obtained for linear phase noise dominated scenario. Meanwhile, if PNi is generated by the nonlinear phase noise model, the result is obtained for nonlinear phase noise dominated scenario. The number of noise sample is chosen as 50 and the length of training sequence is 104 for 8-ary constellation and 105 for 16-ary constellation. We use QAM constellation as a starting constellation in initialization step.

As an illustration, in Fig. 2 we provide the result of the proposed algorithm. The proposed constellations, which can be named as LLR-OSCD constellations, for linear phase noise model are shown in Fig. 2(a), by setting the frequency offset × symbol duration product to 103. Meanwhile, the constellation, which is called NL-OSCD, for nonlinear phase noise is also shown in Fig. 2(b). The result is generated at launch power of –6 dB for the transmission length of 2000km.

 figure: Fig. 2

Fig. 2 The optimized 2D 16-ary signal constellations for: (a) linear phase noise model, (b) nonlinear phase noise model (after 2000 km of SMF and launch power of – 6dBm).

Download Full Size | PDF

4. Proposed LDPC coded modulation scheme

For convenience, we can use the likelihood function defined as

L(ak,θk)=pR(r|ak,θk)pR(r|ak=0) 
If the sequence of L=T/Ts statistically independent symbols, a=[a1aL]T, is transmitted, the corresponding likelihood function will be:
l(a,θ)=l=1LL(ak,θk)
To avoid numerical overflow problems, the log-likelihood function should be used instead, so we have that
l(a,θ)=logL(a,θ)
Clearly, the maximum likelihood (ML) approach leads to exponential increase in complexity as sequence length L increases. Here we propose a different strategy, which is inspired by Monte Carlo integration method. Namely, to calculate the log-likelihood function we will need to perform the following L -dimensional numerical integration:
l(a)=log{exp[l(a,θ)]pΘ(θ)dθ}
Instead of numerical integration, we propose to use Monte Carlo integration. Namely, by using the Monte Carlo integration, the log-likelihood function l(a) can be estimated as:
l(a)=logEθ{exp[l(a,θ)]}
where the expectation averaging Eθ is performed for different phase noise realizations. This method is particularly simple for memoryless phase noise processes, Wiener phase noise process and cyclic slip phase noise process described as a Markov process of reasonable memory. It can be shown that complexity of this method is O{(m2+S)Nr}, where m is the channel memory, S is the sequence length and Nr is the number of phase noise realizations. Compared to the ML method, whose complexity is O{ML}, where M is the signal constellation size, complexity is significantly lower for long sequences to be detected. This method requires the knowledge of Markov phase noise process, which is quite easy to characterize by training. In particular, for the Wiener phase noise process and memoryless phase noise process, the Gaussian noise generator is only needed. Also, we can also use the nonlinear phase noise model introduced in Section 2 to generate the phase noise process.

Our proposed LDPC coded modulation scheme used in the simulations below is shown in the Fig. 3. The b independent data are first encoded by an (n,k) LDPC encoder and written in row-wise fashion into b×n block interleaver. Then either NL-OSCD or LLR-OSCD can be used for modulation and 2D (I/Q) modulator performs the electrical-to-optical conversion. The polarization division multiplexing is used. On receiver side, after polarization diversity receiver, which provides the projection along both polarizations and in-phase and quadrature channels, we perform compensation of linear and nonlinear impairments as well as CPE. We then calculate the log-likelihood function based on Eq. (17), by using the Monte Carlo integration approach introduced above.

 figure: Fig. 3

Fig. 3 The LDPC coded modulation scheme with Monte Carlo integration to evaluate LLRs. PBS/C: polarization beam splitter/combiner, LPF: low-pass filter, ADC: analog-to-digital converter, APP: a posteriori probability.

Download Full Size | PDF

5. Numerical results

The LDPC coded modulation scheme intorduced in Section 4 is used in simulation. The symbol rate is set to Rs=32.25GS/s and the quasi-cyclic LDPC (16935, 13550) code of rate 0.8 is used. All the results are obtained for 3 outer iterations and 20 inner (LDPC decoder) iterations.

5.1 Simulation for linear phase noise channel

The BER performance for LDPC-coded polarization-division multiplexed OSCD and LLR-OSCD constellation sets are summarized in the Fig. 4. We have shown the performance different LLR-OSCD constellation sets optimized with different initial constellations (put in brackets) and different number of phase noise samples, denotes as NS. The constellation size is 16 and the frequency offset × symbol duration product is set to 103. The green line denotes the BER performance of LLR-OSCD that uses the IPQ (Iterative Polar Quantization) as initial constellation, while the black line is optimized from QAM. It is evident that the 16-ary LLR-OSCD optimized from QAM performs the best for 500 phase noise samples. The proposed LLR-OSCD with (NS = 500) outperforms conventional QAM by 1.1 dB and OSCD by 0.2 dB.

 figure: Fig. 4

Fig. 4 BER performance for proposed LLR-OSCDs for 16-ary 2D constellation.

Download Full Size | PDF

5.2 Simulation for the nonlinear phase noise channel

In simulations related to the nonlinear phase noise channel, we first fix the number of spans and the span length (the span length is 80km for all simulations), while vary the launch power in order to determine the optimum launch power (minimizing the BER), as illustrated in Fig. 5. The Fig. 5 clearly indicates the existence of the optimal launch power for the total transmission distance of 4000 km for both 8-ary NL-OSCD and 8-QAM. It can also be noticed that the nonlinear tolerance of 8-NL-OSCD is better. NL-OSCD (1, 2000) curve denotes that this constellation has been obtained by employing the algorithm described in Sec. 3 with the launch power of 1 dBm and for total transmission distance of 2000 km. We can see that NL-OSCDs with different parameters have similar optimal power around 1 dBm. Then we fix the span length and launch power in order to see the performance of different constellations against the total transmission length. The total transmission distance that can be achieved by employing NL-OSCDs and LDPC coding is around 6700 km, as shown in Fig. 6.

 figure: Fig. 5

Fig. 5 BER vs. launch power for uncoded 8-ary NL-OSCD and 8-QAM.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 BER vs transmission length plot for LDPC-coded 8-ary NL-OSCD and 8-star-QAM.

Download Full Size | PDF

In Fig. 6, we use LLR-OSCD to denote the constellation designed for linear phase noise dominated scenario. It is obvious to see that the NL-OSCD outperforms 8-ary star-QAM for almost 500 km. Notice that LLR-OSCD performs a little worse than the QAM, which indicates that LLR-OSCD constellation obtained for linear phase noise channel model is not applicable to the nonlinear phase noise channels. We then perform the similar study for 16-ary NL-OSCD and the results are summarized in Fig. 7 and Fig. 8.

 figure: Fig. 7

Fig. 7 BER vs. launch power for uncoded 16-ary NL-OSCD and 16-QAM.

Download Full Size | PDF

 figure: Fig. 8

Fig. 8 BER vs transmission length plot for LDPC-coded 16-ary NL-OSCD and 16-QAM.

Download Full Size | PDF

The optimal power for 16-ary NL-OSCD is 4 dBm, as shown in Fig. 7, and the total transmission distance exceeds 4400 km, as shown in Fig. 8. The transmission distance can be further increased by employing the turbo equalization principle, which is out of scope of this paper. The Fig. 7 indicates that 16-ary NL-OSCD has better nonlinear tolerance than 16-QAM. On the other hand, Fig. 8 indicates that LDPC-coded 16-ary NL-OSCD outperforms 16-QAM by almost 1000 km and LLR-OSCD by 600 km. The new signal constellations obtained by using the proposed algorithm clearly have better nonlinear tolerance and can significantly extend the total transmission distance of QAM. The main reason for this improvement is due to the fact that during the signal constellation design, the nonlinear nature of the channel is taken into account.

In Fig. 9, we provide an interesting result of our study. If we increase the span length to 100 km, the total transmission length of LDPC-coded 8 ary-NL-OSCD get increased to almost 11000 km with the launch power of −1.5 dB. Clearly, the LDPC-coded 8-ary NL-OSCD outperforms LPDC-coded 16-QM QAM for almost 3000 km. We can conclude that there exists an optimal span length for the system, with which we can achieve better performance. However, in Fig. 9, in addition to constellation optimization, launch power optimization, the span length is optimized as well.

 figure: Fig. 9

Fig. 9 BER vs. transmission length plots for LDPC-coded 16-ary NL-OSCD and 16-QAM when span length is 100 km.

Download Full Size | PDF

6. Conclusion

In this paper, we proposed a new signal constellation designed algorithm, which can be used to design signal constellations suitable for coherent optical channels dominated by either linear or nonlinear phase noise. By using the log-likelihood function as the design criterion in signal constellation design instead of the Euclidean distance, we can take the nonlinear nature of the channel into account. The Monte Carlo simulations demonstrate that for 8-ary and 16-ary cases, LDPC-coded NL-OSCDs outperform the corresponding LDPC-coded QAM constellations for almost 500 km and 1000 km, respectively, for span length of 80 km. The LLR-OSCD outperforms QAM constellation by 1.1 dB for 8-ary, for phase noise dominated channel. We have found that the signal constellation is both launch power and distance dependent, as expected. Further, by optimizing the span length together with the launch power the transmission distance can be further increased. The proposed signal constellation design algorithm can be easily expanded to the multi-dimensional case, which will be the focus of the future work.

Acknowledgments

This work was supported in part by NSF under Grant CCF-0952711.

References and links

1. K. P. Ho, Phase modulated Optical Communication System (Springer, 2005).

2. G. P. Agrawal, Nonlinear Fiber Optics (Academic, 2006).

3. I. Hanzo, W. Webb, and T. Keller, Single- and Multi-carrier Quadrature Amplitude Modulation: Principles and Applications for Personal Communications, WLANs and Broadcasting (Wiley, 2000).

4. A. P. T. Lau and J. M. Kahn, “Signal design and detection in presence of nonlinear phase noise,” J. Lightwave Technol. 25(3), 779–783 (2004).

5. L. Beygi, E. Agrell, and M. Karlsson, “Optimization of 16-point ring constellations in the presence of nonlinear phase noise,” in Optical Fiber Communication Conference/National Fiber Optic Engineers Conference 2011, OSA Technical Digest (CD) (Optical Society of America, 2011), paper OThO4. [CrossRef]  

6. C. Hager, A. Grell, A. Alvarado, and E. Agrell, “Design of APSK constellations for coherent optical channels with nonlinear phase noise,” IEEE Trans. Commun. 61(8), 3362–3373 (2013).

7. I. Djordjevic, T. Liu, L. Xu, and T. Wang, “Optimum signal constellation design for high-speed optical transmission,” in Optical Fiber Communication Conference, OSA Technical Digest (Optical Society of America, 2012), paper OW3H.2. [CrossRef]  

8. T. Cover and J. Thomas, Elements of Information Theory (Wiley, 1991).

9. I. B. Djordjevic and T. Wang, “On the LDPC-coded modulation for ultra-high-speed optical transport in the presence of phase noise,” in Optical Fiber Communication Conference/National Fiber Optic Engineers Conference 2013, OSA Technical Digest (online) (Optical Society of America, 2013), paper OM2B.1. [CrossRef]  

10. Z. Tao, W. Yan, L. Liu, L. Li, S. Oda, T. Hoshida, and J. C. Rasmussen, “Simple fiber model for determination of XPM effects,” J. Lightwave Technol. 29(7), 974–986 (2011). [CrossRef]  

11. M. Magarini, A. Spalvieri, F. Vacondio, M. Bertolini, M. Pepe, and G. Gavioli, “Empirical modeling and simulation of phase noise in long-haul coherent optical transmission systems,” Opt. Express 19(23), 22455–22461 (2011). [CrossRef]   [PubMed]  

12. G. Colavolpe, A. Barbieri, and G Caire, “Algorithms for iterative decoding in the presence of strong phase noise,” IEEE J. Sel. Area Commun . 23(9) 1748–1757 (2005). [CrossRef]   [PubMed]  

13. M. Magarini, L. Barletta, A. Spalvieri, A. Leven, M. Pepe, and G. Gavioli, “Impact of nonideal phase reference on soft decoding of differentially encoded modulation,” IEEE Photon. Technol. Lett. 24(23), 2179–2182 (2012). [CrossRef]  

14. R.-J. Essiambre, G. Kramer, P. J. Winzer, G. J. Foschini, and B. Goebel, “Capacity limits of optical fiber networks,” J. Lightwave Technol. 28(4), 662–701 (2010). [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 (9)

Fig. 1
Fig. 1 Constellation diagrams of received signal with 16-QAM constellation, L = 2000km, for different launch powers.
Fig. 2
Fig. 2 The optimized 2D 16-ary signal constellations for: (a) linear phase noise model, (b) nonlinear phase noise model (after 2000 km of SMF and launch power of – 6dBm).
Fig. 3
Fig. 3 The LDPC coded modulation scheme with Monte Carlo integration to evaluate LLRs. PBS/C: polarization beam splitter/combiner, LPF: low-pass filter, ADC: analog-to-digital converter, APP: a posteriori probability.
Fig. 4
Fig. 4 BER performance for proposed LLR-OSCDs for 16-ary 2D constellation.
Fig. 5
Fig. 5 BER vs. launch power for uncoded 8-ary NL-OSCD and 8-QAM.
Fig. 6
Fig. 6 BER vs transmission length plot for LDPC-coded 8-ary NL-OSCD and 8-star-QAM.
Fig. 7
Fig. 7 BER vs. launch power for uncoded 16-ary NL-OSCD and 16-QAM.
Fig. 8
Fig. 8 BER vs transmission length plot for LDPC-coded 16-ary NL-OSCD and 16-QAM.
Fig. 9
Fig. 9 BER vs. transmission length plots for LDPC-coded 16-ary NL-OSCD and 16-QAM when span length is 100 km.

Tables (1)

Tables Icon

Table 1 Parameters of the System under Study

Equations (17)

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

r k =s( a k , θ k )+ z k ,     r k = [ r k ( 1 ) r k ( i ) r k ( N ) ] T
s( a k , θ k )= e j θ k [ a k ( 1 ) a k ( i ) a k (N) ] T ,   z k = [ z k ( 1 ) z k ( i ) z k ( N ) ] T
θ k =( θ k1 +Δ θ k ) mod 2π
p ΔΘ ( Δ θ k )= n= p(0, δ Δθ 2 ,Δ θ k n2π)
p R ( r| a k , θ k )= e | | r k s( a k , θ k ) | | 2 N 0 /(π N 0 ) 
Y=( X+Z ) e j Φ NL
Φ NL = 0 L γP( z )dz=γ L eff P 
L eff = 1 e αL α
Φ NL = γ L eff { | E 0 + n 1 | 2 + | E 0 + n 1 + n 2 | 2 ++ | E 0 + n 1 ++ n N A | 2 }
δ 2 =2 n sp hναΔν N A  L 
L L m =LL( { A ^ m ,P( A ^ m ) } )= n 1 k=0 n1 max y A ^ m LL( x k ,y )
LL( x k ,y )= 1 NS i=1 NS { x k1 Re[ ( y 1 + y 2 j ) e j×P N i ] } 2 + { x k2 Im[ ( y 1 + y 2 j ) e j×P N i ] } 2 2 δ 2
L( a k , θ k )= p R ( r| a k , θ k ) p R ( r| a k =0 )  
l( a,θ )= l=1 L L( a k , θ k )
l( a,θ )=logL(a,θ)
l( a )=log{ exp[ l( a,θ ) ] p Θ ( θ )dθ}
l( a )=log E θ { exp[ l( a,θ ) ] }
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.