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

Simplified model for nonlinear noise calculation in coherent optical OFDM systems

Open Access Open Access

Abstract

A simplified closed form expression for the noise power due to four-wave mixing in coherent OFDM systems is derived. The proposed model is in very good agreement with the exact model. The derived analytical expressions can be used in performance evaluation of systems employing CO-OFDM with any number of subcarriers and/or as an integral part of physical layer aware routing algorithms.

© 2014 Optical Society of America

1. Introduction

Coherent optical frequency division multiplexing (CO-OFDM) demonstrated the potential of increasing spectral efficiency, making possible optical transmission with per channel data rates beyond 400 Gb/s for future high-speed networks [1, 2]. In such systems, no dispersion compensation mechanisms are implemented whilst the subcarriers are very densely packed in order to maximize spectral efficiency. As a result, there is significant FWM noise [35] which in densely spaced CO-OFDM systems dominates over SPM and XPM.

In this paper, we provide a closed-form approximate expression for the noise power of FWM in OFDM systems the exact solution of which is given in [6]. For a large number of subcarriers the computational time of the model in [6] is significant, making it cumbersome in performance calculations whilst it is useless in Routing and Spectral Assignment (RSA) algorithms that take into account physical layer performance [7]. In prior art, approximations have been reported [3, 4]. However, in [3] there is a lower limit in the number of symbols/s or in the number of subcarriers for a given subcarrier spacing for which the model can be used whilst in [4] fiber dispersion and subcarrier spacing are not taken into account. The main advantages of our analytical model are: 1) there is no lower limit in the number of symbols/s used 2) the expression for the nonlinear noise holds for an arbitrary subcarrier spacing, 3) the closed form expression derived in this work allows to generalise results for the behavior of FWM and ease its introduction to planning/design tools.

2. Simplified FWM model derivation

The system under investigation is the same as in [3] which consists of a single polarization system with many wavelength channels, each of them OFDM modulated. This system can be evaluated as a huge single band OFDM channel provided that the bandwidth of each wavelength channel is much larger than the guard band between them. From [6], and following the same approach as [3] which is generally valid for the central subcarrier, the total noise power due to the FWM in a subcarrier i (Nsub/2iNsub/2), is:

PNLi=Dijk218γ2Leff2eaLPij=Nsub/2Nsub/2k=Nsub/2Nsub/2PjPk11+(2πλ2acΔf2(ik)(jk)D)2sin2{Ns2πλ2cΔf2(ik)(jk)DL/2}sin2{2πλ2cΔf2(ik)(jk)DL/2}(1+4eaLsin2{2πλ2cΔf2(ik)(jk)DL/2}(1eaL)2)
with Nsub the number of subcarriers, Px the power of subcarrier x,Dijk the degeneracy factor (equal to 6 for non-degenerate products and 3 for degenerate products), γ=2πn2/(λAeff) the nonlinear coefficient of the fiber, Aeff the effective core area of the fiber, D the local dispersion parameter, Leff=(1eaL)/a the effective length of the fiber, L the span length, Δf the subcarrier spacing and Ns the number of spans. Equation (1) for the worst case subcarrier frequency (which is the central one with i = 0) and Px=Px, is simplified to:
PNL={Dijk218γ2Leff2P3Ns2Nsub2,Nsub2Nsπλ2Δf2DL4a1c<1Dijk218γ2Leff2P3Nsa1cλ2Δf2DL(1Log[4a1cNsub2Nsπλ2Δf2DL]),1Nsub2Nsπλ2Δf2DL4a1cπNsub4Dijk218γ2Leff2P3Ns2(Nsuba1cNsλ2Δf2DLLog[4Nsubπ]),Nsub2Nsπλ2Δf2DL4a1cπNsub4
Equation (2) is obtained as follows. The phased array term due to the interference between the FWM products in Nsspans is given as
sin2{Nsπλ2Δf2(ik)(jk)DL/c}sin2{πλ2Δf2(ik)(jk)DL/c}
The constructive interferences depend on the term
πλ2Δf2(ik)(jk)DL/c
such as Eq. (3) equals Ns2 when Eq. (4) equals to 0,±π,±2π.... In this case principal maxima occur, producing strong contributions to the total FWM power. Moreover, for a typical OFDM system the majority of the produced FWM terms is around the first maximum. An example is shown in Fig. 1 where the statistical distribution of all of the produced FWM terms for 256 subcarriers is plotted against Eq. (4) for Δf=100MHz and L=100km. As it is obvious the main FWM power is produced by subcarriers that constructively interfere, with their interference located around the first maximum. In parallel, the term
11+(2πλ2acΔf2(ik)(jk)D)2
shows us how strong the interference between the subcarriers in a single span is and gets its maximum value which is 1 around the central maximum while around the kth principal maximum its value is much lower and equals 11+(2πk/aL)2,k=1,2,.... For example, in a system with L=100km and a=0.2dΒ/km the value of the second principal maximum located in ±π is 0.35 and around the third principal maximum located in ±2π is 0.118. This shows that the interference effect is much weaker away from the central principal maximum so it can be ignored. Figure 2 shows the FWM efficiency which is the product between the FWM coefficient for a single span calculated in Eq. (5) and the FWM coefficient between Ns spans calculated in Eq. (3). The number of spans is 10. As a result, by taking into account the two facts that 1) the occurrence of the FWM terms is much higher around the central principal maximum and 2) the constructive interference is much stronger around the central principal maximum we can omit all other principal maxima by using the Taylor approximation for
sin2{πλ2Δf2(ik)(jk)DL/c}(πλ2Δf2(ik)(jk)DL/c)2
for Eq. (4) around 0. As a result Eq. (3) can be simplified to
Ns2sinc2{Nsπλ2Δf2(ik)(jk)DL/c}
Furthermore, the term
4eaLsin2{πλ2Δf2(ik)(jk)DL/c}(1eaL)2
can be omitted for spans longer than 30 km.

 figure: Fig. 1

Fig. 1 Statistical distribution of produced FWM terms for 256 subcarriers, Δf = 100MHz and L = 100 km.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 FWM efficiency against Eq. (4).

Download Full Size | PDF

Since in the previous step all principal maxima except the central have been eliminated we can also eliminate the subsidiary maxima located in positions where Eq. (4) equals±3π2Ns,±5π2Ns.... This elimination is justified by the fact that the strength of a FWM term located in the first and the second subsidiary maximum is less than 0.055 and 0.03 respectively of the strength of the central principal maximum for a system with more than 5 spans. So we can approximate the sinc(.) function with a Gaussian function of the form (for all combinations)

j=Nsub/2Nsub/2k=Nsub/2Nsub/2Ns2(e(Nsπλ2Δf2(ik)(jk)DL/c)2a12)2
Variable a1 takes values around 3 for good matching. By expanding Eq. (9) into Taylor series, we get
j=Nsub/2Nsub/2k=Nsub/2Nsub/2Ns2(a12a12+(Nsπλ2Δf2(ik)(jk)DL/c)2)2
The next step is the conversion from summation to integration of Eq. (10) and the mathematical solution of the double integral. This conversion can be justified that as we can approximate integration with the rectangle rule, we can also approximate a sum using integration under similar conditions. In our approximation we assume i=0, so we calculate the nonlinear noise for the central subcarrier which is the worst case subcarrier. We set a new variable n=jk and Eq. (10) becomes
k=Nsub/2Nsub/2n=Nsub/2kNsub/2kNs2(a12a12+(Nsπλ2Δf2knDL/c)2)2
We now convert the double summation to integration and by dividing with (Nsπλ2Δf2kDL/c)2we write Eq. (11) in a more convenient form:
Nsub/2Nsub/2Nsub/2kNsub/2kNs2(a1cNsπλ2Δf2kDL)4(1(a1cNsπλ2Δf2kDL)2+n2)2dndk
By using the reduction formula
1(x2+m2)kdx=x2m2(k1)(x2+m2)k1+2k32m2(k1)1(x2+m2)k1dx
And knowing that
1(x2+m2)dx=1mArcTan[xm]
Equation (12) is integrated into
Nsub/2Nsub/2Ns2(a122(Nsub/2k(Nsπλ2Δf2kDL/c)2(Nsub/2k)2+a12Nsub/2k(Nsπλ2Δf2kDL/c)2(Nsub/2k)2+a12)+a12(Nsπλ2Δf2kDL/c)(ArcTan(Nsπλ2Δf2kDL/ca1(Nsub/2k))ArcTan(Nsπλ2Δf2kDL/ca1(Nsub/2k))))dk
We can further approximate Eq. (15) by taking the fact that ArcTan(x)=ArcTan(x) and by removing k from the terms Nsub/2k and Nsub/2k with maximum absolute error of 0.63 dB within the range 1-20 spans, 100-800 MHz with step 100MHz and 16-1024 subcarriers. So we have
Nsub/2Nsub/2Ns2(a122(Nsub(Nsπλ2Δf2kDL/c)2(Nsub/2)2+a12)+a1(Nsπλ2Δf2kDL/c)ArcTan(Nsπλ2Δf2kDL/ca1(Nsub/2)))dk
The conversion from discrete summation of Eq. (11) to integration leading to Eq. (16) is not valid for the full range of values of k and system parameters. That is because after integration the peak values of the function inside the integral of Eq. (16) against k are not calculated as in discrete summation, resulting in significantly lower noise especially in higher number of subcarriers and large subcarrier spacing. This can be shown in Fig. 3 where inner sum of Eq. (11) and the function inside the integral of Eq. (16) are plotted against various values of k for Ns=10spans, D=17ps/(nmkm), Nsub=64, Δf=800MHz and L=100km. For higher values of k the peak values of discrete summation of Eq. (11) converge to Ns2.To deal with this we calculate the limit when the following equation is lower than 1:

 figure: Fig. 3

Fig. 3 The red curve shows the inner sum of Eq. (11) and the green curve shows the function inside the integral of Eq. (16), both plotted for various values of k.

Download Full Size | PDF

(a122(Nsub(Nsπλ2Δf2kDL/c)2(Nsub/2)2+a12)+a1(Nsπλ2Δf2kDL/c)ArcTan(Nsπλ2Δf2kDL/ca1(Nsub/2)))

Firstly, we convert the ArcTan(.) function of Eq. (17) to complex function since ArcTan(x)=i2(Log(1ix)Log(1+ix)) and we get

a122(Nsub(Nsπλ2Δf2kDL/c)2(Nsub/2)2+a12)+ia12(Nsπλ2Δf2kDL/c)(Log(1iNsπλ2Δf2kDLNsub2a1c)Log(1+iNsπλ2Δf2kDLNsub2a1c))
Next we split the real from the imaginary part of the complex numbers and since iLog(1±ix)=12iLog((1ix)2)Arg(1±ix) we have
a122(Nsub(Nsπλ2Δf2kDL/c)2(Nsub/2)2+a12)+a12(Nsπλ2Δf2kDL/c)(Arg(1iNsπλ2Δf2kDLNsub2a1c)Arg(1+iNsπλ2Δf2kDLNsub2a1c))
Arg(1ix)Arg(1+ix) converges to π for large x which is the range when Eq. (17) approaches 1, so we can write
a122(Nsub(Nsπλ2Δf2kDL/c)2(Nsub/2)2+a12)+a1π2(Nsπλ2Δf2kDL/c)
The first term of Eq. (20) is much smaller than the second and as a result a1π2(Nsπλ2Δf2kDL/c)=1 so k=|a1c2Nsλ2Δf2DL|. Equation (16) can be rewritten as
Nsub/2Nsub/2Ns2(2a12c2Nsub(Nsπλ2Δf2DL)2(1(2a1cNsubNsπλ2Δf2DL)2+k2)+a1c(Nsπλ2Δf2kDL)ArcTan(Nsπλ2Δf2kDLNsub2a1c))dk,Nsub2<a1c2Nsλ2Δf2DLanda1c2Nsλ2Δf2DLa1c2Nsλ2Δf2DLNs2(2a12c2Nsub(Nsπλ2Δf2DL)2(1(2a1cNsubNsπλ2Δf2DL)2+k2)+a1c(Nsπλ2Δf2kDL)ArcTan(Nsπλ2Δf2kDLNsub2a1c))dk+2Ns2(Nsub2a1c2Nsλ2Δf2DL),Nsub2a1c2Nsλ2Δf2DL
The next step is the integration of Eq. (21). This integral can be solved analytically using Eq. (14) and the following equation:
x1x11xArcTan(xy)dx=i(PolyLog(2,iyx1)PolyLog(2,iyx1))
where function PolyLog(2,iyx1) denotes the polylogarithm function with s = 2: Lis(z)=k=1zkks and y a product of variables.

As a result Eq. (21) is solved into

Ns2(a1cNsπλ2Δf2DL(ArcTan(Nsπλ2Δf2DLNsub24a1c)ArcTan(Nsπλ2Δf2DLNsub24a1c))+ia1cNsπλ2Δf2DL(PolyLog(2,iNsπλ2Δf2DLNsub24a1c)PolyLog(2,iNsπλ2Δf2DLNsub24a1c)),Nsub2<a1c2Nsλ2Δf2DLandNs2(a1cNsπλ2Δf2DL(ArcTan(πNsub4)ArcTan(πNsub4))+ia1cNsπλ2Δf2DL(PolyLog(2,iπNsub4)PolyLog(2,iπNsub4)))+2Ns2(Nsub2a1c2Nsλ2Δf2DL),Nsub2a1c2Nsλ2Δf2DL
Knowing that ArcTan(x)=ArcTan(x) Eq. (23) is simplified to:
Ns2(2a1cNsπλ2Δf2DLArcTan(Nsπλ2Δf2DLNsub24a1c)+ia1cNsπλ2Δf2DL(PolyLog(2,iNsπλ2Δf2DLNsub24a1c)PolyLog(2,iNsπλ2Δf2DLNsub24a1c))),Nsub2<a1c2Nsλ2Δf2DLandNs2(2a1cNsπλ2Δf2DLArcTan(πNsub4)+ia1cNsπλ2Δf2DL(PolyLog(2,iπNsub4)PolyLog(2,iπNsub4)))+2Ns2(Nsub2a1c2Nsλ2Δf2DL),Nsub2a1c2Nsλ2Δf2DL
ArcTan(x) for x>1 converges to π2 and for x1 converges to x while i(PolyLog(2,ix)PolyLog(2,ix)) converges to πLog(1x) for x>1 and to x for x1 . Since πNsub4>1 for all cases we can write
Ns2Nsub2,Nsπλ2Δf2DLNsub24a1c<1Nsa1cλ2Δf2DL(1Log(4a1cNsub2Nsπλ2Δf2DL)),Nsπλ2Δf2DLNsub24a1c1Ns2(Nsuba1cNsλ2Δf2DLLog(4Nsubπ)),NsubNsλ2Δf2DLa1c1
As a result, the nonlinear noise for the worst case subcarrier or the central subcarrier with i = 0 is simplified to Eq. (2). The performance of the system is dominated by the central subcarriers giving us an upper bound for the expected BER. The approach of the central subcarrier is considered acceptable since only the few outer subcarriers have significantly improved performance.

3. Results and discussions

The accuracy between the approximated model of Eq. (2) and the model in [3] against the model of Inoue in Eq. (1) can be shown in Fig. 4 where the absolute approximation error in dB is plotted against various values of subcarrier spacing and number of subcarriers. This figure depicts the maximum absolute error within the range of 1-20 spans with length of 100 km andD=17ps/(nmkm). We can see from Fig. 4 that this error between 1 and 20 spans is less than 1.75 dB for total bandwidth B=ΔfNsub=100GHz and less than 1.1 dB for B=50GHz while exceeds 2 dB in higher system bandwidths. This happens because as the subcarrier spacing increases, the number of the produced FWM terms located in the second and third principal maxima is higher so they cannot be ignored. On the other hand, the accuracy of our model is much better than the well known model in [3] the inaccuracy of which within the range of 1-20 spans exceeds 2 dB for large subcarrier spacing while showing good agreement for subcarrier spacing up to 200 MHz in all cases. This inaccuracy is expected and justified by the fact that in model [3] in order to solve the double integral analytically the “dense subcarrier” assumption [3, Eq. (11)-(12)] is made. Practically, this assumption holds when the subcarrier spacing Δfis lower than 200 MHz making the model unsuitable for larger subcarrier spacing. As a result, Eq. (2) is a fast and reliable formula for calculation of the total noise power due to FWM in coherent OFDM systems for a wide range of system parameters.

 figure: Fig. 4

Fig. 4 Comparison of the Absolute Approximation Error (dB) of Eq. (2) and model in [3] against Eq. (1) for various number of subcarriers.

Download Full Size | PDF

Finally, the validity of the approximation of Eq. (2) is tested for systems with fractional spans (with spans longer than 30 km) by averaging the span lengths and maintaining the number of spans:

L=L1+L2+...+LNsNs
In Fig. 5 we can see the comparison between [6, Eq. (10)] and our approximation in a system with 10 spans with variable lengths which in the current example was Ln = 3 spans of 40 km, 3 spans of 80 km and 4 spans of 100 km. As it is evident, a maximum deviation of 1.25 dB is observed which is acceptable. Furthermore, the results indicate that the number of spans is by far more important than the span length for FWM noise generation as expected since non-linear effects occur mainly about over the first twenty kilometres. After that length, the optical power is very low (due to fiber losses) for the non-linear effects to occur.

 figure: Fig. 5

Fig. 5 Comparison of our model against Inoue’s model with fractional spans [6, Eq. (10)] for 10 spans and 100 MHz and 200 MHz subcarrier spacing.

Download Full Size | PDF

4. Conclusion

An analytical approximate formula for calculating the nonlinear noise in CO-OFDM systems has been derived. This formula is valid regardless of the number of subcarriers and has been benchmarked with the results produced from the exact solution against sub-carrier frequency spacing, number of spans and span length for the central (worst) subcarrier in the spectrum. This model gives us the advantage of the fast calculation of the nonlinear noise providing a maximum deviation of 1.75 dB from the exact model for total bandwidth of 100 GHz and less than 1.1 dB for total bandwidth of 50 GHz while showing lower inaccuracy than previous approximated models.

Acknowledgments

We would like to thank Dr. A. Hadjifotiou for valuable discussions and suggestions. This research was funded by the Operational Program “Education and Lifelong Learning” of the Greek National Strategic Reference Framework (NSRF) Research Funding Program: THALES PROTOMI, grant number MIS 377322.

References and links

1. Q. Yang, A. Al Amin, X. Chen, Y. Ma, S. Chen, and W. Shieh, “428-Gb/s single-channel coherent optical OFDM transmission over 960-km SSMF with constellation expansion and LDPC coding,” Opt. Express 18(16), 16883–16889 (2010). [CrossRef]   [PubMed]  

2. B. Zhu, S. Chandrasekhar, X. Liu, and D. W. Peckham, “Transmission performance of a 485-Gb/s CO-OFDM superchannel with PDM-16QAM subcarriers over ULAF and SSMF-based links,” IEEE Photon. Technol. Lett. 23(19), 1400–1402 (2011). [CrossRef]  

3. X. Chen and W. Shieh, “Closed-form expressions for nonlinear transmission performance of densely spaced coherent optical OFDM systems,” Opt. Express 18(18), 19039–19054 (2010). [CrossRef]   [PubMed]  

4. A. J. Lowery, S. Wang, and M. Premaratne, “Calculation of power limit due to fiber nonlinearity in optical OFDM systems,” Opt. Express 15(20), 13282–13287 (2007). [CrossRef]   [PubMed]  

5. X. Zhu and S. Kumar, “Nonlinear phase noise in coherent optical OFDM transmission systems,” Opt. Express 18(7), 7347–7360 (2010). [CrossRef]   [PubMed]  

6. K. Inoue, “Phase-mismatching characteristic of four-wave mixing in fiber lines with multistage optical amplifiers,” Opt. Lett. 17(11), 801–803 (1992). [CrossRef]   [PubMed]  

7. C. T. Politi, V. Anagnostopoulos, and A. Stavdas, “PLI-aware routing in regenerated mesh topology optical networks,” J. Lightwave Technol. 30(12), 1960–1970 (2012). [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 (5)

Fig. 1
Fig. 1 Statistical distribution of produced FWM terms for 256 subcarriers, Δf = 100MHz and L = 100 km.
Fig. 2
Fig. 2 FWM efficiency against Eq. (4).
Fig. 3
Fig. 3 The red curve shows the inner sum of Eq. (11) and the green curve shows the function inside the integral of Eq. (16), both plotted for various values of k.
Fig. 4
Fig. 4 Comparison of the Absolute Approximation Error (dB) of Eq. (2) and model in [3] against Eq. (1) for various number of subcarriers.
Fig. 5
Fig. 5 Comparison of our model against Inoue’s model with fractional spans [6, Eq. (10)] for 10 spans and 100 MHz and 200 MHz subcarrier spacing.

Equations (26)

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

P N L i = D i j k 2 18 γ 2 L e f f 2 e a L P i j = N s u b / 2 N s u b / 2 k = N s u b / 2 N s u b / 2 P j P k 1 1 + ( 2 π λ 2 a c Δ f 2 ( i k ) ( j k ) D ) 2 sin 2 { N s 2 π λ 2 c Δ f 2 ( i k ) ( j k ) D L / 2 } sin 2 { 2 π λ 2 c Δ f 2 ( i k ) ( j k ) D L / 2 } ( 1 + 4 e a L sin 2 { 2 π λ 2 c Δ f 2 ( i k ) ( j k ) D L / 2 } ( 1 e a L ) 2 )
P N L = { D i j k 2 18 γ 2 L e f f 2 P 3 N s 2 N s u b 2 , N s u b 2 N s π λ 2 Δ f 2 D L 4 a 1 c < 1 D i j k 2 18 γ 2 L e f f 2 P 3 N s a 1 c λ 2 Δ f 2 D L ( 1 L o g [ 4 a 1 c N s u b 2 N s π λ 2 Δ f 2 D L ] ) , 1 N s u b 2 N s π λ 2 Δ f 2 D L 4 a 1 c π N s u b 4 D i j k 2 18 γ 2 L e f f 2 P 3 N s 2 ( N s u b a 1 c N s λ 2 Δ f 2 D L L o g [ 4 N s u b π ] ) , N s u b 2 N s π λ 2 Δ f 2 D L 4 a 1 c π N s u b 4
sin 2 { N s π λ 2 Δ f 2 ( i k ) ( j k ) D L / c } sin 2 { π λ 2 Δ f 2 ( i k ) ( j k ) D L / c }
π λ 2 Δ f 2 ( i k ) ( j k ) D L / c
1 1 + ( 2 π λ 2 a c Δ f 2 ( i k ) ( j k ) D ) 2
sin 2 { π λ 2 Δ f 2 ( i k ) ( j k ) D L / c } ( π λ 2 Δ f 2 ( i k ) ( j k ) D L / c ) 2
N s 2 sin c 2 { N s π λ 2 Δ f 2 ( i k ) ( j k ) D L / c }
4 e a L sin 2 { π λ 2 Δ f 2 ( i k ) ( j k ) D L / c } ( 1 e a L ) 2
j = N s u b / 2 N s u b / 2 k = N s u b / 2 N s u b / 2 N s 2 ( e ( N s π λ 2 Δ f 2 ( i k ) ( j k ) D L / c ) 2 a 1 2 ) 2
j = N s u b / 2 N s u b / 2 k = N s u b / 2 N s u b / 2 N s 2 ( a 1 2 a 1 2 + ( N s π λ 2 Δ f 2 ( i k ) ( j k ) D L / c ) 2 ) 2
k = N s u b / 2 N s u b / 2 n = N s u b / 2 k N s u b / 2 k N s 2 ( a 1 2 a 1 2 + ( N s π λ 2 Δ f 2 k n D L / c ) 2 ) 2
N s u b / 2 N s u b / 2 N s u b / 2 k N s u b / 2 k N s 2 ( a 1 c N s π λ 2 Δ f 2 k D L ) 4 ( 1 ( a 1 c N s π λ 2 Δ f 2 k D L ) 2 + n 2 ) 2 d n d k
1 ( x 2 + m 2 ) k d x = x 2 m 2 ( k 1 ) ( x 2 + m 2 ) k 1 + 2 k 3 2 m 2 ( k 1 ) 1 ( x 2 + m 2 ) k 1 d x
1 ( x 2 + m 2 ) d x = 1 m A r c Tan [ x m ]
N s u b / 2 N s u b / 2 N s 2 ( a 1 2 2 ( N s u b / 2 k ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 k ) 2 + a 1 2 N s u b / 2 k ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 k ) 2 + a 1 2 ) + a 1 2 ( N s π λ 2 Δ f 2 k D L / c ) ( A r c Tan ( N s π λ 2 Δ f 2 k D L / c a 1 ( N s u b / 2 k ) ) A r c Tan ( N s π λ 2 Δ f 2 k D L / c a 1 ( N s u b / 2 k ) ) ) ) d k
N s u b / 2 N s u b / 2 N s 2 ( a 1 2 2 ( N s u b ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 ) 2 + a 1 2 ) + a 1 ( N s π λ 2 Δ f 2 k D L / c ) A r c Tan ( N s π λ 2 Δ f 2 k D L / c a 1 ( N s u b / 2 ) ) ) d k
( a 1 2 2 ( N s u b ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 ) 2 + a 1 2 ) + a 1 ( N s π λ 2 Δ f 2 k D L / c ) A r c Tan ( N s π λ 2 Δ f 2 k D L / c a 1 ( N s u b / 2 ) ) )
a 1 2 2 ( N s u b ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 ) 2 + a 1 2 ) + i a 1 2 ( N s π λ 2 Δ f 2 k D L / c ) ( L o g ( 1 i N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) L o g ( 1 + i N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) )
a 1 2 2 ( N s u b ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 ) 2 + a 1 2 ) + a 1 2 ( N s π λ 2 Δ f 2 k D L / c ) ( A r g ( 1 i N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) A r g ( 1 + i N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) )
a 1 2 2 ( N s u b ( N s π λ 2 Δ f 2 k D L / c ) 2 ( N s u b / 2 ) 2 + a 1 2 ) + a 1 π 2 ( N s π λ 2 Δ f 2 k D L / c )
N s u b / 2 N s u b / 2 N s 2 ( 2 a 1 2 c 2 N s u b ( N s π λ 2 Δ f 2 D L ) 2 ( 1 ( 2 a 1 c N s u b N s π λ 2 Δ f 2 D L ) 2 + k 2 ) + a 1 c ( N s π λ 2 Δ f 2 k D L ) A r c Tan ( N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) ) d k , N s u b 2 < a 1 c 2 N s λ 2 Δ f 2 D L a n d a 1 c 2 N s λ 2 Δ f 2 D L a 1 c 2 N s λ 2 Δ f 2 D L N s 2 ( 2 a 1 2 c 2 N s u b ( N s π λ 2 Δ f 2 D L ) 2 ( 1 ( 2 a 1 c N s u b N s π λ 2 Δ f 2 D L ) 2 + k 2 ) + a 1 c ( N s π λ 2 Δ f 2 k D L ) A r c Tan ( N s π λ 2 Δ f 2 k D L N s u b 2 a 1 c ) ) d k + 2 N s 2 ( N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L ) , N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L
x 1 x 1 1 x A r c Tan ( x y ) d x = i ( P o l y L o g ( 2 , i y x 1 ) P o l y L o g ( 2 , i y x 1 ) )
N s 2 ( a 1 c N s π λ 2 Δ f 2 D L ( A r c Tan ( N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) A r c Tan ( N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) ) + i a 1 c N s π λ 2 Δ f 2 D L ( P o l y L o g ( 2 , i N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) P o l y L o g ( 2 , i N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) ) , N s u b 2 < a 1 c 2 N s λ 2 Δ f 2 D L a n d N s 2 ( a 1 c N s π λ 2 Δ f 2 D L ( A r c Tan ( π N s u b 4 ) A r c Tan ( π N s u b 4 ) ) + i a 1 c N s π λ 2 Δ f 2 D L ( P o l y L o g ( 2 , i π N s u b 4 ) P o l y L o g ( 2 , i π N s u b 4 ) ) ) + 2 N s 2 ( N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L ) , N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L
N s 2 ( 2 a 1 c N s π λ 2 Δ f 2 D L A r c Tan ( N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) + i a 1 c N s π λ 2 Δ f 2 D L ( P o l y L o g ( 2 , i N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) P o l y L o g ( 2 , i N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c ) ) ) , N s u b 2 < a 1 c 2 N s λ 2 Δ f 2 D L a n d N s 2 ( 2 a 1 c N s π λ 2 Δ f 2 D L A r c Tan ( π N s u b 4 ) + i a 1 c N s π λ 2 Δ f 2 D L ( P o l y L o g ( 2 , i π N s u b 4 ) P o l y L o g ( 2 , i π N s u b 4 ) ) ) + 2 N s 2 ( N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L ) , N s u b 2 a 1 c 2 N s λ 2 Δ f 2 D L
N s 2 N s u b 2 , N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c < 1 N s a 1 c λ 2 Δ f 2 D L ( 1 L o g ( 4 a 1 c N s u b 2 N s π λ 2 Δ f 2 D L ) ) , N s π λ 2 Δ f 2 D L N s u b 2 4 a 1 c 1 N s 2 ( N s u b a 1 c N s λ 2 Δ f 2 D L L o g ( 4 N s u b π ) ) , N s u b N s λ 2 Δ f 2 D L a 1 c 1
L = L 1 + L 2 + ... + L N s N s
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.