Methods for Estimating Modal Bandwidth Spectral Dependence

Abstract

Methods for estimating the Effective Modal Bandwidth (EMB) of laser optimized Multimode Fiber (MMF) at a specified wavelength, λS, based on the measured EMB at a first reference measurement wavelength, λM. In these methods the Differential Mode Delay (DMD) of a MMF is measured and the Effective Modal Bandwidth (EMB) is computed at a first measurement wavelength. By extracting signal features such as centroids, peak power, pulse widths, and skews, as described in this disclosure, the EMB can be estimated at a second specified wavelength with different degrees of accuracy. The first method estimates the EMB at the second specified wavelength based on measurements at the reference wavelength. The second method predicts if the EMB at the second specified wavelength is equal or greater than a specified bandwidth limit.

Description

BACKGROUND OF THE INVENTION

The present invention relates in general to the field of optical fibers and more specifically, to multimode fibers (MMF) designed for operation at multiple wavelengths. The present invention also relates to the field of modeling, designing, production, sorting and testing of MMFs. More specifically it relates to the estimation of the MMF EMB at multiple wavelengths.

The invention is also related to modal and chromatic dispersion compensation in Vertical Cavity Surface Emitting Laser (VCSEL) based MMF channels [1]. The methods described here can provide an estimation of the skew in radial DMD pulse waveforms (tilt) at different wavelengths which is critical in the field of modal-chromatic dispersion compensation.

The need for higher bandwidth has been mainly driven by the increasing demand for high-speed backbone data aggregation fueled by video transmission, server applications, virtualization, and other emerging data services. Cost, power consumption, and reliability advantages have favored the predominance of short and intermediate reach optical channels employing transmitters utilizing VCSELs operating at 850 nm over MMF. MMF is currently utilized in more than 85% of datacenter installations, and has a larger core diameter than single-mode fiber (SMF), which reduces connection losses, relaxes alignment tolerances, and reduces connectorization cost.

Recently, new modulation technologies for VCSEL-MMF channels such as PAM-4, and Short Wavelength Division Multiplexing (SWDM) [SWDM alliance], has been proposed in order to increase the data rates. Standards organizations, including the Institute of Electrical and Electronics Engineers (IEEE) working group 802.3cd and the T11 Technical Committee within the International Committee for Information Technology Standards (INCITS) PI-7, are already working on new applications for PAM-4 for optical serial rates over 50 Gb/s per wavelength.

The SWDM concept is similar to the Coarse Wavelength Division Multiplexing (CWDM), already used for SMF channels operating in the 1310 nm spectral region. SWDM requires the specification of the minimum EMB at the wavelengths limits of the operating spectrum (e.g. 850 nm and 953 nm).

The EMB is computed from DMD pulse measurements. The DMD test method, specified within standards organizations [2], describes a procedure for launching a temporally short and spectrally narrow pulse (reference pulse) from a SMF into the core of a MMF at several radial offsets [5]. After propagating through the MMF under test, the pulses are received by a fast photodetector which captures all the MMF core power. The EMB is estimated by the Fourier domain deconvolution of the input pulse from a weighted sum of the received signals for each radial offset launch. The set of weight values utilized in the computation belong to a set of ten representative VCSELs described elsewhere [2]. Due to the test complexity, it is time consuming and the equipment required to perform the test is expensive; EMB test requirements for multiple wavelengths will significantly increase testing time and consequently, increase fiber cost. A method to estimate the EMB from measurements at a single wavelength would therefore, reduce testing time and cost. The challenges to achieve such a method are described below.

FIG. 1, shows a simulation of EMB vs wavelength 100 for a MMF fiber compliant to the OM4 standard. In this figure, we show the EMB 105 has a peak value at λ_p 120. The labels 115 and 125 show the measured and predicted wavelengths, λ_M and λ_S, respectively. The range 110 shows the spectral window in which the fiber can maintain an EMB higher than a specified value, i.e. 4700 MHz km for OM4.

In principle, based on MMF theory, when all the physical parameters of the fiber are known (i.e. dimensions, refractive profile, dopant types and content), the EMB at λ_S can be predicted from the EMB value at λ_M. In practice however, variations in the refractive index design and dopant content during the preform fabrication process produce changes in 100 which prevent the estimation of the EMB at λ_S. FIG. 2 shows simulated MMFs with identical EMB at λ_M 200, but different EMB spectral dependence. Peaks 205, 210, 205 are different and uncorrelated with 200. Moreover, since the spectral windows 220, 225 and 230 are different, an estimation of the EMB at λ_S 240 is not possible.

Shown in FIG. 3. is the EMB at 850 nm and 953 nm for a large number of simulated fibers, represented using rectangle markers 300 with random variations in their refractive index core. The horizontal and vertical axes of this figure represents the EMB at λ_M=850 nm vs. EMB λs=953 nm respectively. A subset of these fibers that meet the TIA-492AAAD OM4 EMB specification are represented by diamonds markers 305. This figure shows the lack of correlation among EMBs at 850 nm and 953 nm. For example, in 310, a measured fiber with EMB=6000 MHz km at 850 nm can have any value from 1500 to 3000 MHz km at 953 nm. Conversely, 315 shows that a MMF with EMB=2000 MHz km at 953 nm can have any value from 200 to 15000 MHz km at 850 nm. This simulation, which was extended for a large range of wavelengths from 800 nm to 1100 nm, clearly shows that there is no direct relationship between the fiber's EMB at a specified wavelength, λ_S, and the EMB at a measured wavelength, λ_M, when λs≠λ_M.

A method that enables the prediction of the EMB at an arbitrary wavelength based on measurements at another wavelength is desirable to reduce testing time and cost of a MMF.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows the EMB for an OM4 standard compliant multimode fiber.

FIG. 2 shows the EMB for several OM4 standard compliant multimode fibers with identical bandwidths at 850 nm.

FIG. 3 shows EMB at 953 nm vs. EMB at 850 nm for 1000 multimode fibers.

FIG. 4 shows a traditional MMF refractive index profile.

FIG. 5 shows an MMF refractive index profile with a trench in the cladding.

FIG. 6 shows a flow diagram for a training method.

FIG. 7 shows a flow diagram for an estimation method.

FIG. 8(a) shows DMD pulses at 850 nm and 953 nm for three fibers of population A.

FIG. 8(b) shows DMD pulses at 850 nm and 953 nm for three fibers of population B.

FIG. 9 shows measured EMB at 953 nm vs. 850 nm for A and B population.

FIG. 10(a) shows centroid features for A and B populations for radial offsets 1-12.

FIG. 10(b) shows centroid features for A and B populations for radial offsets 13-24.

FIG. 11(a) shows peak position features for A and B populations for radial offsets 1-12.

FIG. 11(b) shows peak position features for A and B populations for radial offsets 12-24.

FIG. 12 shows G_C and G_P functions as a function of radial offset for A and B populations.

FIG. 13(a) shows measured and estimated DMD pulses at 953 nm for A population.

FIG. 13(b) shows measured and estimated DMD pulses at 953 nm for B population.

FIG. 14 shows estimated EMB vs. measured EMB after training.

FIG. 15 shows the results for the classification of fibers.

SUMMARY OF THE INVENTION

Methods for estimating the Effective Modal Bandwidth (EMB) of laser optimized Multimode Fiber (MMF) at a specified wavelength, λ_S, based on the measured EMB at a first reference measurement wavelength, λ_M. In these methods the Differential Mode Delay (DMD) of a MMF is measured and the Effective Modal Bandwidth (EMB) is computed at a first measurement wavelength. By extracting signal features such as centroids, peak power, pulse widths, and skews, as described in this disclosure, the EMB can be estimated at a second specified wavelength with different degrees of accuracy. The first method estimates the EMB at the second specified wavelength based on measurements at the reference wavelength. The second method predicts if the EMB at the second specified wavelength is equal or greater than a specified bandwidth limit.

DETAILED DESCRIPTION OF THE INVENTION

The present invention discloses novel methods to estimate the EMB of a MMF at a desired wavelength, from measurements performed at another wavelength. The first method, Method 1, can be used to predict the EMB at an arbitrary wavelength, λ_S, based on an EMB measurement at a different wavelength, λ_M. The second method can be used to evaluate if the EMB at an arbitrary wavelength, λ_S, is equal of greater than a minimum specified threshold. Each method provides different degree of complexity and accuracy.

These methods can be used for the design and manufacturing processes of MMF that have a core and a cladding where the index of refraction of the cladding is less than that of the core. The core has a gradient index of refraction which varies from a peak value at the center of the core to a minimum value at the cladding interface following a predominant alpha-profile function to minimize modal dispersion [JLT 2012]. Refractive index profiles for two types of MMF are shown in FIGS. 4 and 5. In FIG. 4 a traditional MMF refractive index profile is shown. The profile 400 does not present any abrupt discontinuity inside the core or inside the cladding. The propagating mode groups of this fiber are shown in 410. In FIG. 5 the refractive index profile 500 abruptly changes in the cladding due to the refractive index trench 520 introduced to provide lower bending loss. Labels 510 and 515 shows some of the propagating and leaking mode groups respectively.

Waveguide theory for alpha-profile fibers has been well developed [ref]. The theory can enable the modeling of fiber DMD behavior over a broad range of wavelengths, when the profiles and dopants concentrations are known. In practice however, due to manufacturing variations the designed “optimum” refractive index profile is distorted deterministically and randomly. Very small alterations in 400 or 500, basically change the way the mode groups 410, 510 interact with the variations in refractive index, which destroys or reduces the correlations among DMDs at different wavelengths as it was showed in FIG. 3.

Method 1

This method, can be used to predict the EMB at an arbitrary wavelength, λ_S, based on an EMB measurement at a different wavelength, λ_M. The method was developed based on the inventors' realization that in order to increase the correlation among EMB measurements at λ_M, and a second wavelength, λ_S, a new approach that fully utilizes the information provided by the measured DMD waveforms is required. The method proposed here uses the DMD pulse waveform information at λ_M, such as centroids, peak position, width, shapes, energy per radial offset, and skews, to predict the EMB at a second wavelength. Statistical and signal processing techniques disclosed here, allow us to extract and utilize those parameters to distort the DMD pulse waveforms acquired at λ_M, to predict the DMD pulse waveforms at λ_S. This method which requires a training of the algorithm, enables the prediction of EMBs at different wavelengths from one measurement. FIGS. 6 and 7 show the block diagrams for the training and estimation processes respectively. For illustrative purposes, we use an example to describe both methods.

Training for Method 1

In 600, the populations of TIA-492AAAD standards compliant OM4 fibers from two suppliers (A and B), which use different manufacturing processes are selected. It is understood that the population used here is only an example and is not restricted to any specific number of fiber suppliers. In 602, we select a subset of 24 fibers from manufacturer A and 12 from manufacturer B for training. In 604, the DMD of all fibers are measured at the first measurement wavelength, λ_M=850 nm, and the second specified wavelength, λ_S, which in this example is taken to be 953 nm. These measurements are stored in the array y(r,t,λ) for analysis. FIGS. 8(a) and 8(b) show the DMD radial pulses for three MMF from each population at 850 nm (dark trace) and 953 nm (lighter trace). FIGS. 8(a) and 8(b) show that most of the fibers have similar DMD pulses at low radial offset for both wavelengths. For population A, the DMD pulse shapes are very different at larger radial offset for the two wavelengths.

The EMBs computed from the measured DMD pulses for the A and B populations at 850 nm and 953 nm are shown in FIG. 9. These measurements agree with simulation results showed in FIG. 3, which indicates that EMBs at different wavelengths are uncorrelated.

In step 606 of FIG. 6, the main features of the DMD pulses at each wavelength are extracted. This process captures the main characteristics required to describe the DMD pulses at each radial offset and each wavelength for post-processing and analysis. As an illustrative example, here we extract the centroid, mean power, peak power value and position, and the root mean square (RMS) width. The centroid feature is computed using,

$\begin{array}{cc}{C}_{r,\lambda }=\frac{\sum _k{t}_{k}y\left(r,t_k,\lambda \right)}{\sum _{k}y\left(r,t_k,\lambda \right)}& \left(1\right)\end{array}$

where r is the radial offset index that relates the position of the single-mode launch fiber to the MMF core center axis during the DMD measurement, t is the discrete length normalized temporal, k is the time index. The variable t and k are related to the number of temporal samples simulated or acquired from the oscilloscope during DMD measurements at a given wavelength. The mean power is computed by,

$\begin{array}{cc}{\mathrm{Ymean}}_{r,\lambda }=\frac{\sum _{k}y\left(r,t_k,\lambda \right)}{\sum _k{t}_k}& \left(2\right)\end{array}$

The peak power is computed using,

$\begin{array}{cc}{\mathrm{Ypeak}}_{r,\lambda }={\max }_t\left(y\left(r,t_k,\lambda \right)\right)& \left(3\right)\end{array}$

where max_t(.) is a function that finds the maximum of the DMD pulses for each radial offset and for each wavelength. The peak position is computed using,

$\begin{array}{cc}{P}_{r,\lambda }=\mathrm{find\_peak}\left(y\left(r,t,\lambda \right)\right)& \left(4\right)\end{array}$

where, find peak is a function that finds the maximum value of the DMD pulses for each radial offset and for each wavelength. The RMS width of the pulse for each radial offset is computed,

$\begin{array}{cc}{W}_{r,\lambda }=\sqrt{\frac{\sum _k{\left(t_k-C_{r,\lambda }\right)}^{2}y\left(r,t,\lambda \right)}{\sum _{k}y\left(r,t_k,\lambda \right)}-{\left(T_{REF}\right)}^2}& \left(5\right)\end{array}$

where T_rEF is the RMS width of the reference pulse used for the measurement.

The features extracted from DMD measurements at λ_M, are used to predict features at λ_S, based on the model described in equations (6-8).

$\begin{array}{cc}{C}_{r,{\lambda }_s}=\left(1+I_C\left(r\right)\right)C_{r,{\lambda }_M}+F_C\left({\lambda }_M,{\lambda }_S\right)G_C\left(r\right)& \left(6\right)\end{array}$

where C_r,λS, and C_r,λM represent the centroids per radial offset at λ_M and λ_S, I_C(.,.),F_C(.,.),G_C(.) is the set of polynomial functions that describe the relationship between centroids at those wavelengths.

$\begin{array}{cc}{P}_{r,{\lambda }_s}=\left(1+I_P\left(r\right)\right)P_{r,{\lambda }_M}+F_P\left({\lambda }_M,{\lambda }_S\right)G_P\left(r\right)& \left(7\right)\end{array}$

where P_r,λS, and C_r,λM represent the centroids per radial offset at λ_M and λ_S., I_P(.,.), F_P(.,.), G_P(.) is the set of polynomial functions that describe the relationship between peak positions at those wavelengths.

$\begin{array}{cc}{W}_{r,{\lambda }_s}=\left(1+I_W\left(r\right)\right)W_{r,{\lambda }_M}+F_W\left({\lambda }_M,{\lambda }_S\right)G_W\left(r\right)& \left(8\right)\end{array}$

where W_r,λS, and W_r,λM represent the centroids per radial offset at λ_M and λ_S., I_W(.,.), F_W(.,.), G_W(.) is the set of polynomial functions that describe the relationship between widths at those wavelengths.

The F(.,.) functions are solely dependent on the measured and targeted wavelength. These functions accommodate for chromatic effects in the refractive index and material. The G(.) functions are solely dependent on radial offsets and accommodate for relationships between the group velocity of DMD pulses at different radial offset in the fiber core. The I(.) functions, dependent on the radial offset, accommodates for mode transition due to the change of wavelengths.

In step 608, the features extracted from the measured DMD pulses at the two wavelengths are used to find the coefficients of the polynomial functions described above (6-8). Standard curve fitting techniques are applied as described in [3]. For the samples used in this example, FIGS. 10 (a) and (b), show the centroid features for 850 nm and 953 nm for radial offsets from 1 to 24 microns for the two fiber populations A (red) and B (blue). FIGS. 11 (a) and (b), show the peak positions for 850 nm and 953 nm for radial offsets from 1 to 24 microns for the two fiber populations. For these samples, F(850,953) was 16 ps/μm/km for population A and 13.3 ps/μm/km for population B. The functions G_C(.,.) and G_P(.,.) for a cubic polynomial curve fitting is shown in FIGS. 12 (a) and 12 (b) for fiber populations A and B respectively. Similarly, curves for the other features described above (1-5) are obtained.

In 610 the correlations among the features, i.e. the ones shown in FIGS. 1-12 are evaluated. If they are higher than a determined threshold, e.g., 80%, the model is ready to use and the process end in 615. If not, in 612 the signal to noise ratio (SNR) of all DMD measurements are evaluated. If the noise of the measurements is higher than a pre-determine threshold, the measurements need to be repeated. If the SNR is high, but the correlations are low, it is possible that the samples do not represent the fiber population and a new set of samples will be required.

Method 2: Estimation Method

After training, the method for the DMD mapping and estimation, shown in FIG. 7, is ready to use. Here, we use the same example to describe the processes. In 700 the fibers that require EMB estimation are selected. In 702, the model described in (6-8) and the wavelengths (in this case λ_M=850 nm, λ_S,=953 nm) are selected. In 704, the DMD at λ_M is measured. In 706 the features are extracted from the DMD pulse centroids at λ_M using equations (1-5). In 708 the DMD pulses are estimated at λ_S. Next, the model described in equations (6-8) is used to estimate the features C_r,λS,P_r,λS,W_r,λS,Ymean,_λS at λ_S.

The parameter P_r,λS is used to reposition each of the DMD pulses using,

$\begin{array}{cc}{y}_p\left(r,t_k,{\lambda }_S\right)=y(r,t_k-\left(P_{r,{\lambda }_S}-P_{r,{\lambda }_M},{\lambda }_M\right)& \left(9\right)\end{array}$

where the y_P(.,.,.) array represents the estimated DMD pulses after the peak position correction.

The differences between the centroid and peak position are computed at both wavelengths. The variation of these differences are computed as shown,

$\begin{array}{cc}\Delta =\left(C_{r,{\lambda }_S}-P_{r,{\lambda }_S}\right)-\left(C_{r,{\lambda }_M}-P_{r,{\lambda }_M}\right)& \left(10\right)\end{array}$

The parameter Δ is used to estimate the new width and skew of the DMD pulses at λ_S. In the majority of cases, when, λs>λ_M, the DMD pulse width tends to increase. Conversely, when λs<λ_M, the width tends to decrease. The changes in skew and width are corrected using a linear filter as shown,

$\begin{array}{cc}{y}_W\left(r,t_k,{\lambda }_S\right)=\sum _{i=0}^{\mathrm{Ntaps}}{A}_i{y}_p\left(r,t_k-Ki\Delta ,{\lambda }_S\right)& \left(11\right)\end{array}$

where y_W(.,.,.) represents the estimated DMD after equalization, i is the equalizer tap index, Ntaps the number of taps, A_i represents the tap coefficient, K is a scaling factor.

For each fiber, the optimum values of Ntaps, A_i, and K, are found by numerically searching. The constraint conditions or equations for this search are the estimated mean, peak, and the values shown in table I.

TABLE I
∑                          k                                                ⁢                                                                              y                            W                                                    ⁡                                                      (                                                          r                              ,                                                              t                                k                                                            ,                                                              λ                                S                                                                                      )                                                                                                                                                ∑                          k                                                ⁢                                                  t                          k                                                                                      ≤                                          Ymean                                              r                        ,                                                  λ                          S
maxt(yW(r, tk, λS)) ≤ Ypeakr,λS
∑                            k                                                    ⁢                                                                                                                    (                                                                                                      t                                    k                                                                    -                                                                      C                                                                          r                                      ,                                      λ                                                                                                                                      )                                                            2                                                        ⁢                                                                                          y                                W                                                            ⁡                                                              (                                                                  r                                  ,                                  t                                  ,                                                                      λ                                    S                                                                                                  )                                                                                                                                                                                          ∑                            k                                                    ⁢                                                                                    y                              W                                                        ⁡                                                          (                                                              r                                ,                                                                  t                                  k                                                                ,                                                                  λ                                  S                                                                                            )                                                                                                                                            ≤                                          W                                              r                        ,                                                  λ                          S

In 710, the algorithm evaluates if the conditions shown above can be maintained below a pre-determined threshold, e.g., 60% of the estimated constraint' values. If that is not achieved, in 712 the SNR of the DMD measurement is evaluated. Depending on this, the DMD may need to be measured again 704. Otherwise, in 717 it is indicated that the estimation failed. If the conditions compared in 710 are achieved, the algorithm provides the DMD corrected pulses and the estimated EMB is obtained.

FIG. 13, shows the corrected DMD results for populations A and B. In FIG. 14 the estimated and measured bandwidths at λs=953 nm are shown. The correlation for these results is around 80%-90%.

Method 2

This method can be used to predict if the EMB at an arbitrary second wavelength, λ_S, is equal or greater than a specified threshold, EMB_th, based on a DMD measurement at a different wavelength, λ_M. As in the previous case this method utilizes features of the DMD pulse waveforms at λ_M, such as centroids, peak position, width, shapes, energy per radial offset, and skews. The average centroid for positions Rt__start-Rt__end is defined using,

$\begin{array}{cc}{C}_{\mathrm{Top}}\left(R_{\mathrm{T\_start}},R_{\mathrm{T\_end}}\right)=\frac{\sum _{r=R_{\mathrm{T\_end}}}^{R_{\mathrm{T\_start}}}{C}_{r{\lambda }_M}}{R_{\mathrm{T\_end}}-R_{\mathrm{T\_start}}+1}& \left(12\right)\end{array}$

The average centroid for positions RB__start-RB__end is defined using,

$\begin{array}{cc}{C}_{\mathrm{Bottom}}\left(R_{\mathrm{B\_start}},R_{\mathrm{B\_end}}\right)=\frac{\sum _{r=R_{\mathrm{B\_end}}}^{R_{\mathrm{B\_start}}}{C}_{r,{\lambda }_M}}{R_{B\mathrm{\_end}}-R_{B\mathrm{\_start}}+1}& \left(\mathrm{l3}\right)\end{array}$

A function denominated, P-Shift is computed as

$\begin{array}{cc}P-\mathrm{Shift}\left(R_{\mathrm{T\_start}},R_{\mathrm{T\_end}},R_{\mathrm{B\_start}},R_{\mathrm{B\_end}}\right)=C_{\mathrm{Top}}\left(R_{\mathrm{T\_start}},R_{\mathrm{T\_end}}\right)-C_{\mathrm{Bottom}}\left(R_{\mathrm{B\_start}},R_{\mathrm{B\_end}}\right)& \left(14\right)\end{array}$

The slopes using the peak pulse position for two or more radial regions are computed as shown in equation below.

$\begin{array}{cc}P-{\mathrm{Slope\_R}}_k=\frac{1}{L}\frac{\stackrel{\mathrm{R\_end}}{\sum _{r={\mathrm{R\_start}}_k}}\left(P_{r{\lambda }_m}-T_k\right)\left(r-\left({\mathrm{R\_end}}_k+{\mathrm{R\_start}}_k\right)/2\right)}{\sum _{r={\mathrm{R\_start}}_k}^{{\mathrm{R\_end}}_k}{\left(r-\left({\mathrm{R\_end}}_k+{\mathrm{R\_start}}_k\right)/2\right)}^2}& \left(15\right)\end{array}$

where k is the index that represent the selected radial offset regions and

$\begin{array}{cc}{T}_k=\frac{1}{{\mathrm{R\_end}}_k-{\mathrm{R\_start}}_k+1}\stackrel{{\mathrm{R\_end}}_k}{\sum _{r={\mathrm{R\_start}}_k}}{P}_{r{\lambda }_m}& \left(16\right)\end{array}$

The widths for the same k regions that are computed using:

$\begin{array}{cc}M-{\mathrm{Width}}_k=\frac{1}{L}\left[{\max }_{r={\mathrm{R\_start}}_k:R_{en{d}_k}}\left(W\left(r\right)\right)\right]& \left(17\right)\end{array}$

It should be noted for features described in (15-17), the k index can take values from 1 to N_k where N_k<25 r of radial offsets, i.e. 25. In practice, as shown in the algorithm training example described below, low values for Nk, i.e. N_k=2, are enough to provide estimations with low uncertainty.

The training method, which is described below, utilize machine learning techniques to find the radial-offset regions that maximize the difference between parameters such as P_shift, P_slopes and M_widths for two or more population of fibers. One population of fiber will have EMB>EMB_th at λ_S and other populations will not satisfy this constraint. After training the estimation method simply evaluates if the extracted features from MMF under test belong to the regions found during training that satisfy the condition, EMB>EMB_th at λ_S based on the DMD measurements at λ_M.

Training for Method 2

The training process is identical to the one shown in FIG. 6, with exception of steps 606 and 608. We use the same example, starting from 606 of FIG. 6, to illustrate the training method.

In step 606, the main features of the DMD pulses at λ_M, are extracted. Note the differences with the first method which require the computation of the features at each wavelength, λ_M and λ_S. The extracted features are C_r,λM, P_r,λM (centroid, peak and width) using equations (1-5).

In 608, the training is performed. The training is an iterative process that has the goal to maximize a metric or a series of metrics that represents the differences in features of two groups of fibers. One group, Group 1 are composed by the MMFs that have EMB>EMB_th at λ_S and the other group, Group 2 by MMFs that have EMB<EMB_th at λ_S.

Initially, all the MMFs used for training are mapped in a space defined by the P_shift, P_slopes and M_widths. The initial values of the regions utilized in (12-17) which are {R_{B_start},R_{B_end}}, {R_{T_startT}, R_{T_end}}, {R_start_k,R_end_k} are set to random values.

In this example, the utilized metric is a function implemented in C, Python, or Matlab, which computes p-norm distances in the mentioned space, among the MMFs that belong to the groups Group 1 and Group 2.

$\begin{array}{cc}M\left(\text{?}\text{?}\text{?}\text{?},{\mathrm{R\_start}}_1,{\mathrm{R\_end}}_1,\dots ,\text{?}\text{?}\right)=\{(\mathrm{P\_Shift}\mathrm{\_Group}\text{?}\mathrm{P\_Shift}\mathrm{\_Group}\text{?}+\sum _{k=1}^{N_k}{A_{1,k}\left(\mathrm{P\_Slopes}\text{?}-\mathrm{P\_Slopes}\text{?}\right)}^p+,{\hspace{1em}\sum _{k=1}^{N_k}{A_{2,k}\left(\mathrm{M\_Width}\text{?}-\mathrm{M\_Width}\text{?}\right)}^p\}}^{1/p}& \left(18\right)\\ \text{?}\text{indicates text missing or illegible when filed}& \end{array}$

where A_1,k, A_2,k are weight parameters to quantify the relative importance of each features and/or radial offset regions.

In each iteration the coordinate axes are modified by changing the values of {R_{B_start}, R_{B_end}}, {R_{T_start},R_{T_end}},and the set of k parameters {R_start_k,R_end_k}. In addition, the norm parameter p and the weights, can be also optimized in each iteration.

During the optimization process, the values can be changed at random, or in deterministic ways. For example, using the random search algorithms or using gradient methods. The features are recomputed using (12-17) for each new set of regions. The MMFs are mapped in the new space and the utilized metric, i.e. equation (18) is computed. The process continue until the metric is maximized, or until an exhaustive search is produced.

To illustrate how the algorithm improves the metric in each iteration we use a set of 35 MMFs. For sake of simplicity we utilize N_k=2, A_1,1=A_1,2=1, and p=1 and the following simplified version of the metric, (18)

$\begin{array}{cc}M\left(R_{T\mathrm{\_start}},R_{T\mathrm{\_end}},R_{B\mathrm{\_start}},R_{B\mathrm{\_end}},{\mathrm{R\_start}}_1,{\mathrm{R\_end}}_1,{\mathrm{R\_start}}_2,{\mathrm{R\_end}}_2\right)=\left\{{\left(\mathrm{P\_ShiftGroup1}-\mathrm{P\_ShiftGroup2}\right)}^p+\sum _{k=1}^2{A}_{1,k}\left(\mathrm{P\_Slopes}\mathrm{\_Group}{1}_k-\mathrm{P\_Slopes}\mathrm{\_Group}{2}_k\right)\right\}& \left(19\right)\end{array}$

FIG. 15 shows initial mapping of the population for one to 8000 iterations. In the figure. the square markers represent MMF from Group 1 and the circle markers represents MMFs from Group 2. It can be observed that for the initial iterations 1-5000, FIGS. 15 (a), (b), (c), (d) and (e ) it is not possible to differentiate between both populations. After 7000 iterations the algorithm capable of separate MMF from Group 1 and Group 2. The boundaries between the groups, Group 1 and Group 2, in the plane shown in FIG. 15 (i) can be established (see black trace). Based on this classification the optimum radial offset that optimizes the feature extraction from the DMD pulse waveforms was found. The values of the found regions are: 2 to 10 micron radial offsets for the first P-Slope (k=1), 12 to 23 microns for the second P-Slope (k=2). For the P-shift calculation shown in (14), the optimum regions were 2 to 3 microns for C_Top and 18-24 microns for the C_Bottom

The training using the disclosed algorithm demonstrates that the MMFs for Group 1 and Group 2 have distinctive features that can be observed when the optimum set of radial regions to represent them are selected. These results demonstrate a method to predict if EMB>EMB_th at λs based on the DMD measurements at λ_M.

Estimation Method

During training the optimum radial-offset regions to extract the features that optimally represent MMFs that have EMBs>EMB_th at λ_S were found. In the feature-space, see for example FIG. 15(i), the Groups of MMFs that have the desired characteristics can be separated by a line or in general by a polynomial that isolate two regions one for Group 1 and another for Group 2. For the estimation process the features of a MMF are extracted from DMD measurements at λ_M and mapped in the feature-space. If the MMF belongs to the desired regions that produce EMB>EMB_th at λ_S (see FIG. 17 the fiber is accepted. Otherwise the fiber is rejected.

Note that while this invention has been described in terms of several embodiments, these embodiments are non-limiting (regardless of whether they have been labeled as exemplary or not), and there are alterations, permutations, and equivalents, which fall within the scope of this invention. Additionally, the described embodiments should not be interpreted as mutually exclusive, and should instead be understood as potentially combinable if such combinations are permissive. It should also be noted that there are many alternative ways of implementing the methods and apparatuses of the present invention. It is therefore intended that claims that may follow be interpreted as including all such alterations, permutations, and equivalents as fall within the true spirit and scope of the present invention.

Also note that nothing in this disclosure should be considered as limiting and all instances of the invention described herein should be considered exemplary.

Claims

1. A method for generating a mapping model for use in estimating the modal bandwidth of a multimode fiber at a second wavelength (λS), based on a DMD measurement performed at only a first wavelength (λM), comprising: selecting a population of multimode fibers (600);performing a DMD measurement for each multimode fiber of the population of multimode fibers at the first wavelength and the second wavelength (604);extracting at least one signal feature of the DMD measurement from each multimode fiber of the population of multimode fibers at the first wavelength and the second wavelength,the signal feature being at least one of a centroid, mean power, peak power value and position, and root mean square (RMS) width of the DMD measurement at the first wavelength and the second wavelength respectively (608);mapping the DMD measurement of the first and second wavelength against the at least one signal feature to determine coefficients of the mapping model (608).

2. A method for generating a mapping model for use in predicting if the EMB of a multimode fiber (MMF) at an arbitrary wavelength, λS, is equal or greater than a specified bandwidth threshold, EMBth, based on a DMD measurement at a different wavelength, λM, comprising: selecting a population of multimode fibers (600);performing a DMD measurement for each multimode fiber of the population of multimode fibers at the different wavelength, λM, and the arbitrary wavelength, λS (604);extracting at least one signal feature of the DMD measurement from each multimode fiber of the population of multimode fibers at the different wavelength, λM, the signal feature being at least one of a centroid, mean power, peak power value and position, androot mean square (RMS) width of the DMD measurement at the different wavelength, λM, (606);mapping the DMD measurement of the arbitrary wavelength and different wavelength against the at least one signal feature to determine coefficients of the mapping model (608), wherein the coefficients are determined by an iterative process that maximizes a metric that represents the differences in features in a first group of fibers and a second group of fibers, wherein the first group are composed of MMFs that have EMB>EMBth at the arbitrary wavelength λS and second group are composed of MMFs that have EMB<EMBth at the arbitrary wavelength λS.