Devices based on the transmission and processing of optical signals are becoming increasingly common. Computer and communication networks that utilize optical fibers for the transmission of data are now commonplace. Such networks rely on optical fibers and other elements such as light amplifiers, multiplexers, demultiplexers, dispersion compensators, etc. to carry and process the light signals.
The passage of an optical signal through devices such as optical fibers results in a delay that can vary with the state of the polarization of the light. This dependence of the propagation characteristics of light waves on the polarization state is often referred to as polarization mode dispersion (PMD). PMD can result in problems in the transmission of data at very high data rates because portions of the light signal having different polarization states will arrive at slightly different times, and hence, a light pulse used to transmit data will suffer a broadening effect that can lead to interference between successive pulses in a pulse train of data.
A knowledge of the PMD of an optical fiber and optical components is important for a number of reasons. For example, a system designer must know the magnitude of the PMD before setting the maximum data transmission rate over the fiber. Accordingly, methods for measuring PMD have been sought. In one proposed method, four independent measurements of the group delay of an optical fiber or other device under test are made with light of four different polarization states. Since the measured group delay drifts with temperature and is perturbed by vibrations, this technique suffers from environmental sensitivity, especially in long fibers. Further, the PMD is obtained by taking the difference of the large measured group delays. Hence, a small error in the measured group delays leads to a much larger error in the PMD.
The present invention includes a polarization mode dispersion analyzer having a light source, a sensor, an output phase signal analyzer, and a controller. The light source generates a probe light signal that is intensity modulated and also polarization modulated, the light source being adapted to apply the light signal to a device under test. The sensor generates an output phase signal related to the phase of the intensity modulation of an output optical signal leaving the device under test. The output phase signal analyzer measures an amplitude and phase of at least one frequency component of the output phase signal at a frequency related to the polarization modulation. The controller generates a signal indicative of a differential group delay of the device under test utilizing the measured amplitude and phase. In one aspect of the invention, the controller also measures a group delay associated with the device under test. In another aspect of the invention, the probe light signal includes a light signal in which all three Stokes vector polarization components are periodic functions of time that are characterized by one or more modulation frequencies and wherein the at least one frequency component is one of the modulation frequencies.
The manner in which the present invention provides its advantages can be more easily understood with reference to
In one embodiment of the present invention, the probe light source signal is generated by a combination of light source 21, intensity modulator 26, and polarization modulator 22. Light source 21 generates a light signal having a fixed polarization state. The intensity of the optical signal from light source 21 is modulated by intensity modulator 26 at an intensity modulation frequency. The polarization of the light leaving intensity modulator 26 is modulated by polarization modulator 22 prior to the light source signal being applied to a device under test 25. The manner in which the polarization is modulated will be discussed in detail below.
The light leaving device under test 25 is input to a sensor 23 that measures the power in the light signal that leaves device under test 25 as a function of time. The sensor 23 must be capable of receiving optical signals at the intensity modulation frequency. In one embodiment, sensor 23 includes a photodetector and a synchronous detection circuit that generates in-phase and quadrature signals from the output of the photodetector. The synchronous detection circuit may operate at the intensity modulation frequency. Alternatively, the synchronous detection circuit may include an electrical mixer that mixes the electrical frequency down from the intensity modulation frequency to a lower intermediate frequency (IF). Then, the synchronous detection, implemented in hardware or software, may operate at the IF. The in-phase (I) and quadrature (Q) output of the synchronous detector provide an estimate of power (R) and phase (θ) of the received intensity modulated signal. The estimated phase θ(t) is further processed as described below. Any other arrangement that generates the phase θ(t) of the intensity modulated optical signal received by the photodetector at the intensity modulation frequency could be utilized.
The output phase signal θ(t) can be described in terms of polarization modulation frequencies. The output phase signal, θ(t), of sensor 23 is analyzed by a phase-shift analyzer 30 that measures the amplitude and phase of the frequency components contained in θ(t) at a set of frequencies that are determined by polarization modulator 22. Phase-shift analyzer 30 can include a vector spectrum analyzer, a lock-in amplifier, or any other form of synchronous detection, implemented in hardware or in software, that allows simultaneous measurement of amplitude and phase of the frequency components of θ(t) at the set of frequencies in question.
The operation of these components is under the control of controller 24 that also performs the computations needed to provide a measurement of the PMD of device under test 25. Controller 24 can be any data processing system that is also capable of generating the required potentials for the polarization modulator, reading the information from phase shift analyzer 30, and performing the computations described below. General-purpose signal generation and data processing systems or special purpose hardware can be utilized to construct such a controller.
As noted above, intensity modulator 26 modulates the light signal from light source 21 at an intensity modulation frequency fm. The intensity modulation frequency must be high enough to measure phase shifts related to the group delay or the differential group delay (or PMD) of the device under test. The intensity modulation frequency is typically between few MHz and few GHz. In the absence of polarization modulation by polarization modulator 22, θ is a constant. If device under test 25 has detectable PMD and the polarization state of the intensity modulated light signal is modulated appropriately by polarization modulator 22, then the observed phase will be modulated at frequencies determined by the polarization modulation. That is, θ will become a function of time, θ(t). The manner in which the PMD is derived from the observed modulation θ(t) will be explained in detail below after a discussion of the manner in which the polarization of the light signal applied to the device under test is modulated.
The operation of the present invention can be more easily understood in terms of the Stokes vector, which describes the state of polarization of a light signal. The Stokes vector has 4 components, S0-S3. The first component, S0, is the intensity of the light signal and the remaining three components describe the state of polarization of the light signal. The polarization state of the light signal is represented as a vector in a three dimensional space in which the unit vectors along the three axes can be viewed as representing the fraction of the light with various types of polarization. The S1 axis measures the content of linear polarization with the positive values corresponding to horizontally polarized light, and the negative values corresponding to vertically polarized light. The S2 axis measures the content of linear polarization at 45 degrees to the horizontal (or vertical) with the positive values corresponding to +45 degree polarized light, and the negative values corresponding to −45 degree polarized light. Finally, the S3 axis measures the content of circular polarization, the positive values representing right-hand circularly (right circular) polarized light and the negative values representing left-hand circularly (left circular) polarized light. The normalized Stokes vector has all its components normalized with respect to its first element. Thus, the normalized intensity is equal to one. Refer now to
The Stokes vector parameters can be related to the electric field of the light waves. Refer now to
S
0
=|E
x|2+|Ey|2
S
1
=|E
x|2−|Ey|2
S
2=2Re(E*xEy)
S
3=2Im(E*xEy) (1)
It should be noted that the components of the electric field vector are complex values. It can be shown that θ(t) is related to a set of elementary parameters p0, p1, p2, and p3 by the following relationship.
where, after appropriate scaling, p0 represents the group delay; p1 represents the 0° component of the differential group delay; p2, represents the 45° component of the differential group delay, and p3 represents the circular component of the differential group delay. Here, q(t), u(t), and v(t) are normalized Stokes vector components that are modulated, and hence, are functions of time. The parameters, pj, can be related to the elementary matrices of the Jones calculus for describing polarization states. The reader is referred to Bogdan Szafraniec and Douglas Baney, “Elementary Matrix-Based Vector Optical Network Analysis” Journal of Lightwave Technology, Vol. 4, April 2007 for a detailed description of elementary matrices and elementary parameters.
Referring again to
Refer now to
The manner in which the potentials are chosen will be explained in more detail below. For the present discussion, it is sufficient to note that a first periodic waveform is applied between electrodes 53 and 52, and a second periodic waveform is applied between electrodes 53 and 51. Electrode 53 is a reference (ground) electrode. In general, the waveforms have the same period. Over each period of the waveforms, the Stokes vector of the output light traverses a predetermined path (trajectory) on the surface of the Poincare sphere. The path is chosen such that all of the polarization dependent components of the Stokes vector are modulated with sufficient amplitude to measure the corresponding elementary parameter, pj, related to each Stokes vector component during each cycle of the applied waveforms. The path may also be chosen to have its center of gravity in the center of the sphere.
Refer now to
It should be noted, however, that while q(t), u(t) and v(t) are periodic, q(t), u(t) and v(t) cannot each be pure tones simultaneously. For the components to be pure tones, there must be three frequencies, wq, wu, and wv, for which
q(t)=cos(ωqt)
u(t)=cos(ωut+Du)
v(t)=cos(ωvt+Dv)
q
2(t)+u2(t)+v2(t)=1, (3)
where Du and Dv are fixed phase shifts. It can be shown that this system of equations has no solutions with these constraints.
While a solution in which each of the components is a single tone cannot be found, a solution that only depends on three tones is possible. For example,
q(t)=cos(2ωt)
u(t)=(sin(ωt)+sin(3ωt))/2
v(t)=(−cos(ωt)+cos(3ωt)/2 (4)
The above equations satisfy the constraint q(t)2+u(t)2+v(t)2=1 and describe a trajectory that produces only three tones in the signal from sensor 23.
A more detailed discussion of the considerations that go into choosing a trajectory on the Poincare sphere is provided below. For the purposes of the present discussion, it will be assumed that a trajectory for the Stokes vector on the Poincare sphere has been chosen.
Given a trajectory on the Poincare sphere, the potentials that must be applied to the electrodes in the polarization modulator must be determined. To provide these potentials for a known constant input polarization state, a calibration table that maps the voltages on the two electrodes onto polarizations on the Poincare sphere is constructed. For the purposes of this discussion, it will be assumed that the polarization modulator is a modulator such as that shown in
This process can be more easily understood by referring to
Consider a normalized form of the Stokes vector (1, q(t), u(t), v(t)). Since the Stokes vector is modulated using a periodic modulation function, each of its components is also a periodic function. Hence, each component can be represented as a Fourier series. If the Stokes vector executes a closed loop on the Poincare sphere at an angular frequency, ω, then each component can be expanded in a Fourier series with ω as the fundamental frequency. The number of significant harmonics in the series depends on the details of the trajectory chosen on the Poincare sphere. For example, the trajectory described by Eq. (4) has only 3 significant harmonics. In the more general case, mathematically, the components of the polarization dependent components of the Stokes vector can be written in the following form
q(t)=C1+A1,1sin(ωt+φ1,1)+A1,2sin(2ωt+φ1,2)+A1,3sin(3ωt+φ1,3)
u(t)=C2+A2,1sin(ωt+φ2,1)+A2,2sin(2ωt+φ2,2)+A2,3sin(3ωt+φ2,3)
v(t)=C3+A3,1sin(ωt+φ3,1)+A3,2sin(2ωt+φ3,2)+A3,3sin(3ωt+φ3,3) (5)
The constants Ci, φi,j, and Ai,j, where i=1 to 3 and j=1 to N, can be measured experimentally by measuring the output of the polarization modulator for each Stokes' component using an appropriate polarization filter that removes two of the three components. The measurement can be performed using a vector spectrum analyzer, a lock-in amplifier, or any other form of synchronous detection that allows simultaneous measurement of amplitude and phase of individual harmonics. As will become clear from the following discussion, the number of harmonics that are significant, N, must be at least 3.
The above-described polarization modulation patterns all involve expanding the polarization dependent components of the Stokes vector in a harmonic series. That is, each component is expanded in terms of a number of component frequencies in which the component frequencies are integer multiples of some fundamental frequencies. However, as will be discussed in detail below, there are cases in which the polarization dependent components of the Stokes vector can be expanded in a series in which the frequencies are not integer multiples of a predetermined frequency. Hence, in the general case, it will be assumed that
q(t)=C1+A1,1sin(ω1t+φ1,1)+A1,2sin(ω2t+φ1,2)+A1,3sin(ω3t+φ1,3)
u(t)=C2+A2,1sin(ω1t+φ2,1)+A2,2sin(ω2t+φ2,2)+A2,3sin(ω3t+φ2,3
v(t)=C3+A3,1sin(ω1t+φ3,1)+A3,2sin(ω2t+φ3,2)+A3,3sin(ω3t+φ3,3) (6)
As will become clear from the following discussion, there must be at least three frequencies ωj. In the case of a harmonic expansion, ωj=j*ω, where ω is the fundamental frequency.
Referring again to Eq. (2), the measured phase, θ(t), can be re-written in the form
θ(t)=p0+p1q(t)+p2u(t)+p3v(t) (7)
Substituting the expansion of Eq. (6) for q(t), u(t), and v(t)
where N is at least 3. It is convenient to introduce a complex notation that combines constants Ai,j and φi,j into a single complex constant zi,j=Ai,j exp(jφi,j), where j=√{square root over (−1)} is an imaginary number. In the complex notation equation (8) becomes
Hence, now complex θ(t) is modulated at each of the frequencies wj. As noted above, phase-shift analyzer 30 shown in
θ0=p0+p1C1+p2C2+p3C3 (9a)
and
θ1=p1Z1,1+p2Z2,1+p3Z3,1
θ2=p1Z1,2+p2Z2,2+p3Z3,2
θ3=p1Z1,3+p2Z2,3+p3Z3,3 (9b)
It is important to note that the equation (9b) can be rewritten in a matrix form θ=Z·p, where the matrix Z, contains elements zi,j, and θ and p are vectors. In equation (9a), θ0 is the DC component of θ(t). The constants Cj are the average components of the Stokes vector over the modulation path on the Poincare sphere. They describe the unmodulated part of the light. To simplify the following discussion, it will be assumed that the polarization modulation paths on the Poincare sphere are chosen such that the average of each polarization dependent component of the Stokes vector is 0. i.e., Cj=0. This choice corresponds to the depolarized light. Then, from equation (9a), θ0=p1, i.e., the group delay is found from the DC component of the demodulated θ(t). To be precise, the group τg is
where fIM is the frequency of the intensity modulation.
The system of equations shown in Eq. 9(b) will have a solution if the determinant of the matrix Z is non-zero. In this case, the values of the elementary parameters pj can be determined. It is important to note that the reference phase utilized in the synchronous detection of the components of θ(t) must be properly set. In general, the reference phase is set such that the computed values of pj are real numbers. The differential group delay is obtained from these parameters as follows:
DGD=√{square root over (p12+p22+p32)}/(2πfIM) (10)
where fIM is the intensity modulation frequency utilized by intensity modulator 26.
The above-described embodiments utilized only three of the frequencies even in those cases in which the Stokes vector component expansion includes additional harmonics or other frequencies. However, embodiments in which more of the components are utilized to provide an over determined system in which noise is further reduced could be constructed. In addition, if the determinant of Z is 0 for some choice of the 3 frequencies and there are additional frequencies, a matrix constructed from other frequencies may have a non-zero determinant.
In the above-described embodiments the designer determines the trajectory on the Poincare sphere and generates the modulation signals that are applied to the polarization modulator from a calibration model for that polarization modulator. The coefficients of the matrix Z are then measured experimentally. If the determinant of Z is zero, or too small to allow for an accurate solution of the system of equations, a new trajectory on the Poincare sphere is chosen and the process repeated.
Alternatively a known trajectory as that described by the Eq. (4), discussed above, can be used. The trajectory produces only three frequencies. The corresponding Z matrix is:
where j=√{square root over (−1)}. The determinant of the above matrix is equal to j/2. Refer now to
The choice of trajectory from among those that generate matrices that have non-zero determinants can be guided by some general principles that are listed below. Trajectories that generate fewer frequencies for all Stokes vector components are preferred. Only three frequencies are needed to solve for the corresponding coefficients that determine the DGD. The additional harmonics divert energy that would have gone into the harmonics that are being used; hence, trajectories that generate a significant number of additional harmonics are likely to lead to lower signal-to-noise ratios.
The number of harmonics that are generated by any given trajectory may depend on the number of harmonics in the corresponding drive signals that are applied to the electrodes in the polarization modulator. Also, complicated voltage waveforms are more difficult to synthesize, and hence, can lead to more complex driving circuitry for the modulator.
There is also a limit on the voltages that can be generated by the controller and applied to the polarization modulator. Hence, a trajectory on the Poincare sphere must be traversable using voltages that are within some predetermined range of voltages that are determined by the polarization modulator and the controller.
Refer now to
The modulation of the various Stokes vector components are shown in
The above-described trajectories on the Poincare sphere are closed loops, and hence, the modulation frequencies are harmonics of the frequency with which the closed loop is traversed. For the purposes of the present discussion, a path will be defined as being closed if it begins and ends at the same point on the Poincare sphere. This will always be the case when the Stokes vector is a periodic function. In some cases, it may be advantageous to use modulation frequencies that are unrelated frequencies instead of harmonics. For example, such unrelated frequencies could reduce some errors caused by harmonics produced by non-linearities of the power sensor. Trajectories in which the Stokes vector components are modulated in a periodic manner without requiring the trajectory to be closed are possible. An example of such a trajectory is given by
q(t)=cos(2ω1t)
u(t)=(sin(2ω1t−ω2t)+sin(2ω1t+ω2t))/2
v(t)=(cos(2ω1t−ω2t)+cos(2ω1t+ω2t))/2 (14)
with ω1=eω/2 and ω2=ω. Here e is the irrational number, 2.71828 . . . The controller detects modulation at (e−1)ω, eω, and (e+1)ω, where ω is chosen to provide detection at frequencies that are within the range of the analyzer contained within the controller. It should be noted that while the Stokes vector components are described by periodic functions, the trajectory defined by Eq. (14) is not periodic. The trajectory is endless and eventually samples the entire Poincare sphere surface without repeating itself.
The above-described embodiments of the present invention utilize a light source that has a fixed polarization state. The light source can be a tunable laser light source that allows characterization of components over wavelength. The fixed polarization state of the laser source is maintained within the intensity modulator and modulated by the polarization modulator. Highly monochromatic tunable laser sources are very attractive in embodiments in which the device under test includes an optical fiber, optical fiber components, or other devices having fiber interfaces or small dimensions. However, embodiments based on other light sources such as LEDs can also be constructed. If the light source does not provide light with a constant fixed polarization, a polarization filter can be introduced between the light source and the polarization modulator or as part of the input port of the polarization modulator.
The polarization modulators discussed above can modulate the polarization of the fixed polarization light source at frequencies of up to 100 kHz and in some cases 1 MHz. Hence, the DGD measurements can be made in a time that is short compared to the time periods over which the optical properties of the device under test change. In addition, the present invention does not rely on taking differences of large group delays to measure the DGD. Accordingly, the present invention provides a substantial improvement over the prior art methods discussed above.
The above-described embodiments of the present invention output the value of DGD defined in Eq. (10). However, embodiments that provide the values of the individual quantities pj, or some subset thereof, can also be provided. In this regard, it should be noted that the pj, for j=1 to 3 are related to the differential group delays associated with specific polarization states.
Various modifications to the present invention will become apparent to those skilled in the art from the foregoing description and accompanying drawings. Accordingly, the present invention is to be limited solely by the scope of the following claims.