Apparatus and method for the complete characterization of optical devices including loss, birefringence and dispersion effects

[0001] This is a continuation of U.S. patent application Ser. No. 10/005,819, filed on Dec. 14, 2001, the entire contents of which are incorporated by reference.

[0002] This application claims priority from provisional application Ser. No. 60/255,077 filed Dec. 14, 2000 and now abandoned.

[0003] 1. Field of Application

[0004] The present invention generally relates to the measurement of the optical characteristics of optical waveguide devices and, in particular, relates to devices for characterizing the optical effect of various devices upon an optical system.

[0005] 2. Discussion of Prior Art

[0006] Waveguide devices such as Bragg gratings, interleavers, couplers, isolators, etc. require accurate characterization if they are to be used in deployable optical fiber networks. Generally this characterization has been accomplished using a few simple specifications such as insertion loss, bandwidth, polarization dependent loss, etc. The devices actually have responses that are significantly more complicated that these specifications can describe. Although the simplification of the characterization is necessary to allow for the use of interchangeable parts, and rational specifications for manufacturers to meet, obtaining a full and complete measurement of the device from which any desired specification can be calculated would have benefits to both manufacturers and their customers.

[0007] First of all, full device characterization would permit computer system models that are much more accurate, and would inherently contain aspects of the device that have not been explicitly specified. Second, in cases where manufacturers use different specifications, or consumers require a new set of specifications, devices could be accurately compared, or new specifications developed from the complete measurement function stored in a database. For most devices this measurement function would consist of four complex functions that would occupy less than 200 kB of storage. Finally, if the complete characterization of the device can be achieved with a single instrument, then final testing of devices would be greatly simplified while remaining completely general.

[0008] It is worthwhile to describe some of the fundamental differences between the characteristics of fiber-optic links, and fiber-optic components (excluding fiber). Fiber-optic links are generally very long (>50 km), and have very broad spectral features that vary on the scale of tens of nanometers. Because of the great lengths involved and high launch powers used to over come the loss, nonlinear interactions of light in the fiber are of significant interests. These nonlinear properties are not easily reduced to a simple transfer function matrix. In very long fiber-optic links such a submarine links, optical amplifiers are used to regenerate the signal and overcome loss. These amplifiers often operate in saturation which is another non-linear effect not captured in a simple transfer function matrix. Generally, fiber-optic links are not well characterized by simple transfer functions.

[0009] The relatively short path lengths (<1 m) within fiber optic components prevents any significant non-linear behavior. Because many passive fiber-optic components are filters used to separate channels closely spaced in wavelength, they are designed to have rapid variations in their properties, as a function of wavelength. Optical filters often display variations on the scale of tens of picometers instead of tens of nanometers as seen in fiber optic links. Not only does the transmission amplitude of these functions vary rapidly, but the phase response does as well, and so, in high bit-rate applications, the phase response must be accurately measured as well. Also, due to the fabrication techniques, many filters display some variation as a function of polarization. Again, in high bit-rate systems, this variation must be accurately characterized. It is in these last two measurements, optical phase, and polarization dependence, that the inventions presented here excel.

[0010] The underlying technology used was initially developed by Glombitza and Brinkmeyer, and published in 1993 (“Coherent Frequency-Domain Reflectometry for Characterization of Single-Mode Integrated-Optical Waveguides,” Journal of Lightwave Technology, Vol. 11, No. 8, August 1993). No patents were filed by either author with regard to this technology. Patents of a similar nature are U.S. Pat. No. 5,082,368 “Heterodyne optical time domain reflectometer,” which employed an acousto-optic modulator, and U.S. Pat. No. 4,674,872 “Coherent reflectometer,” which employed an optical phase shifter, and has expired due to nonpayment of maintenance fees. A co-pending U.S. patent application Ser. No. 09/606,120 filed Jun. 16, 2000 by Froggatt and Erdogan for an invention entitled “Single Laser Sweep Full S-Parameter Characterization of Fiber Bragg Gratings,” herein incorporated by reference, expands upon the published and unpatented work by Glombitza and Brinkmeyer to describe how the measurement of the phase and amplitude of the time domain response of a system can be used to determine the frequency domain (spectral) response of the system as well via an inverse Fourier transform. U.S. Pat. No. 5,896,193 for “Apparatus for testing an optical component” also describes an interferometer much like that described by Glombitza and Brinkmeyer although the patent filing on Feb. 14, 1997 post dates the Glombitza and Brinkmeyer article publication date of August 1993 by over three years.

[0011] The Jones Matrix Representation

[0012] An optical fiber component supporting two polarization modes can be fully described using a wavelength dependent Jones matrix. This matrix describes the transfer of energy from the device input to the device output, encompassing the full effects of polarization dependence. The Jones Matrix relates the input electric field to the output electric field by:

[0013] $\begin{array}{cc}\left[\begin{array}{c}{E}_{s,\mathrm{out}}\\ E_{p,\mathrm{out}}\end{array}\right]=\left[\begin{array}{cc}a& b\\ c& d\end{array}\right]\left[\begin{array}{c}{E}_{s,\mathrm{in}}\\ E_{p,\mathrm{in}}\end{array}\right]& \mathrm{Eqn}.1\end{array}$
where E_S and E_P are complex electric field amplitudes of the two orthogonal fields used to represent the total electric field. These orthogonal fields may be any orthogonal state and linear states oriented at 90 degrees are the most common. However, left hand and right hand circular would be an equally valid choice. No particular pair of orthogonal states is needed for the present invention.

[0014] This four element two-by-two matrix is completely general and includes all aspects of the device behavior. If the four complex numbers (a, b, c & d) can be measured as a function of frequency, then the device has been completely characterized. It is therefore an object of this invention to measure these four quantities, and from this basic physical measurement derive any and all of the desired parameters describing the device.

[0015] It is an object of the present invention to supply a controlled light input to a fiber optic device and then compare the output from the device with the input to characterize the effect of the device on such light.

[0016] In the broadest sense, the present invention involves the supplying of light having an optical frequency content over the frequency range of interest where the light has a polarization state which varies linearly with the frequency. This light is provided to the device under test. The light is affected by the device and is output therefrom. The affected output light is combined in a combiner with the input light over the frequency range of interest. The combined light from the combiner is detected and digitized over a plurality of measurements in the frequency range. A microprocessor derives power curves from the digitized combiner output. The microprocessor calculates the Fourier transform of the respective curves and derives from the Fourier transforms four arrays of constants for the Jones matrix characterizing the device.

[0017] In a specific embodiment, a source of variable frequency light which could be a tunable frequency laser is provided where the light has a similar polarization. The light is supplied to a “spinner” which provides output light having a polarization which varies linearly with the changing frequency of the light. The light from the spinner is supplied as an input to the device under test and also to a beam combiner. The resultant light from the device under test is also supplied to the input of the combiner and its output is supplied to a polarization beam splitter. After having the combined light split into two orthogonal polarization states, the two orthogonally polarized beams of light are supplied to respective detectors. The outputs from the detectors is digitized over a plurality of measurements taken at differing frequencies and respective curves power are derived. A programmed microprocessor calculates the Fourier transform for each of the curves and derives four arrays of constants for a Jones matrix characterizing the device under test.

[0037] For simplicity of understanding, similar components will be similarly numbered throughout the several views.

[0038] The inventions described here can readily measure the matrix above, however, the basis vectors used to describe the electric field is left unknown. This topic, and the reason that it can be neglected is described in the section “Unitary Translation and Polarization Characteristics”.

[0039] Unitary Translation and Polarization Characteristics

[0040] Telecommunications systems use optical fiber that supports two polarization modes. Some effort has been expended to develop single-polarization mode optical fiber, but the losses remained too high to permit their use in long-distance transmission systems, and so, two polarization modes are supported in all currently installed fiber-optic links, and will continue to be supported in all links installed in the foreseeable future.

[0041] A “supported” mode is one that will propagate over long distances. The two modes supported in telecom links can be described as “degenerate” meaning that light in each mode propagates down the fiber at nearly the same speed. The difference in the speed of propagation between the modes is referred to as polarization-mode-dispersion. For short lengths of optical fiber (<100 m), this difference in propagation velocity is so small that the delay difference is essentially zero. Over very long transmission distances (>50 km), the very small difference in propagation speed can cause significant differences in the delay between polarization modes. This accumulated difference now forms the barrier to increased data rates in optical fiber. As a result, methods to compensate for polarization mode dispersion are now at the center of many optical fiber communications research programs.

[0042] When the polarization-mode-dispersion in an optical fiber is kept small, the two modes are considered to be degenerate. A consequence of degeneracy is that energy can very easily move between the supported modes in the fiber. This freedom of movement between the supported modes in the optical fiber means that the polarization state in optical fiber is essentially uncontrolled. Mild bends, external stresses, and thermal changes can all change the state of polarization of light in optical fiber.

[0043] Although the input and output polarization states for a segment of fiber can be considered to be independent, the lossless nature of fiber transmission does impose some restrictions on the output state of the light. First of all, the amplitude of the output light must be the same as the input light (which is essentially the definition of lossless). Less obvious, but of great importance is the maintenance of orthogonality, i.e., when the dot product of two fields is zero.

[0044] Orthogonal polarization states do not interfere with one another, and by combining the states in different phases and amplitudes, we can represent all states of polarization. Pairs of orthogonal states, therefore, represent an important physical construct. In a lossless system, such as an optical fiber, orthogonal inputs result in orthogonal outputs. Any section of optical fiber can then be modeled as a “unitary translation.”

[0045] What does the concept of unitary translation mean to optical fiber systems? A good analogy is that of a three dimensional object in space that is rotated. Take for instance, a cat. If I have a cat and rotate it by some amount left to right, and then another amount front to back, potentially every point on the cat changes location, however, all of the important properties about the cat remain the same. Equivalently, if an arbitrary section of optical fiber is placed between an instrument and a device, the measurement of the device is rotated by two unknown angles, but the properties of the device are all measured correctly. And because polarization states are not controlled in fiber-optic communications systems, the end user of the device cannot be certain of the device's “orientation” (with respect to the incident polarization) in the system, and must be assured that the device will work for all orientations.

[0046] The concept of unitary translation can be extended further to measurement error. Any measurement error that can be represented as a multiplication by unitary matrices then becomes irrelevant to the measurement. The ability to fold such measurement errors, or unknowns, into the unknown fiber lead is what allows the present invention to be constructed in a tractable and economical way.

[0047] FIG. 1 is a block diagram of a generic version of an optical vector analyzer in accordance with the present invention. The device under test (DUT) 10 is supplied with light from source 12 (indicated by the dotted line box) where the light has an optical frequency content over the frequency range of interest (for the particular DUT) and the light has a polarization state which varies linearly with the frequency. In one preferred embodiment, the light supplied by the source 12 is provided by an optical source 14 (either a narrow linewidth source such as a tunable laser tunable over the bandwidth of interest or a broadband source of light including the bandwidth of interest) which supplies light to the device entitled “spinner” 16 (discussed in detail below), which provides the light with a polarization state which varies linearly with the frequency. In one embodiment the source 14 is a narrow linewidth tunable source such as a tunable laser, or a broadband source such as an LED or incandescent filament light. Although optical fiber is shown in the various drawings, it is understood that it could be replaced by a free-space collimated beam. Additionally, while couplers (indicated by ellipses) are used, the could be replaced by bulk-optic beam splitters.

[0048] The spinner 16 causes the state of polarization of the light to rotate through a great circle around the Poincaré sphere as function of wavelength. The relationship between the angular change of the polarization state thought the great circle and the frequency of the light must be linear. The period of the variation (i.e. the frequency change needed for 2 pi rotation) must be greater than the spectral resolution of the acquisition system. In the case of a tunable laser source, this is the instantaneous linewidth of the source, and, in the case of an OSA (Optical Spectrum Analyzer, a widely used commercially available device), this is the resolution bandwidth of the OSA.

[0049] The output from source 12 is provided to a reference path 18 and to the input of the DUT 10. The reference path is a method of transmitting some portion of the light from the source 12 to the beam combiner 20. This path must have a Jones Matrix with a determinant that is not close to zero. In other words, this path cannot be a polarizer, and, is preferably a lossless optical transmission such as a fiber, but it need not be perfectly or substantially lossless. The matrix of the reference path must be invertible in a practical sense. The transmission path to the DUT need only satisfy the same criteria as the reference path matrix. The transmission path from the DUT should be substantially lossless, or lossless to the degree of the accuracy required by the measurement.

[0050] The output of the DUT 10 along with the output of the reference path 18 is supplied to the input of a beam combiner 20 (indicated by dotted line box). The combiner adds the light output from the DUT 10 and the reference path 18 for a plurality of frequencies in the frequency range of interest and splits the sum into two orthogonal polarizations. In preferred embodiments, the beam combiner 20 can be a coupling structure 22, such as a bulk-optic beam splitter or a fused tapered coupler, which supplies a combined light signal to a polarizer 24 which breaks the field into two orthogonal polarization states.

[0051] The two orthogonal polarizations from beam combiner 20 are supplied to the microprocessor 26 (indicated by dotted line box). In a preferred embodiment, the microprocessor 26 includes a spectral acquisition device 28 and a mathematical processor 30. The function of the spectral acquisition device 28 is to measure the power as a function of wavelength of the two signals from the polarizer. These signals may be acquired sequentially or in parallel depending upon the application and embodiment.

[0052] The spectral acquisition device 28, if the optical source is a narrow linewidth tunable source such as a laser, will be a pair of detectors from which the power is sampled as a function of the frequency of the tunable source. Also, if the source is a narrow linewidth tunable source, then it may be necessary to send a small part of the light from the source to the spectral acquisition unit so that the wavelength may be accurately measured (as shown by the dotted line 13 connecting source 12 and microprocessor 26 in FIG. 1). If the source is a broadband optical source, then the spectral acquisition device 28 will need to be a spectral measurement system such as an OSA. If the spectral measurement is made by an OSA, then it will not be necessary to send a portion of the light from the source to the OSA.

[0053] The mathematical processing unit 30 includes a digitizer to convert the spectral measurements made in device 26 to digital form if necessary and is programmed to derive curves fitting the spectral outputs, where the curves are optical power versus optical frequency curves. Fourier transforms of each of the curves are calculated by the processor. From the Fourier transforms, four arrays of constants for a Jones matrix characterizing the DUT are calculated. Knowing the four arrays of constants of a Jones Matrix for the DUT allows total knowledge of how and in what fashion any device under test will affect an input signal and thus affect a fiber optic system.

[0054] Details of the “spinner” 16 are shown in FIG. 2 and comprises an fiber-optic network constructed from Polarization Maintaining (PM) fiber 32, a splitter such as a PM coupler 34, a polarization beam combiner 36. In this embodiment the optical source 14 is a tunable laser which is desirable because the power vs. wavelength will be initially measured in 10 to 50 femtometer steps. The light from the laser is coupled to the polarization maintaining fiber 32 such that all of the light is in a single polarization mode. The light is then evenly split, using the polarization maintaining (PM) coupler 34. Half of the light travels through a long path Lp, and half of the light travels through a short path Sp. Light from the two paths is then recombined using the polarization spinner 36.

[0055] The spinner 16 puts the light from one path into the x-polarization, and light from the other path into the orthogonal y-polarization. Because the path lengths of the two components are unequal, the relative phase between the x-polarized light and the y-polarized light is strongly wavelength dependent. As a result the polarization state of the light exiting the polarization beam combiner 36 varies from linear to elliptical to right hand circular to elliptical, etc. as shown in the 9 graphs of FIG. 3. Graph 1 illustrates the evolution of the polarization state of the light exiting the polarization beam combiner 36. The curves show the path that the electric field of the light traverses at optical frequencies. Note that the states include both circular polarization states, and two orthogonal linear states.

[0056] If the polarization state is viewed on the Poincaré sphere as shown in FIG. 4, it will trace a great circle through the poles of the sphere. As the light passes through the fiber following the beam splitter, the polarization states will be changed. The changes in polarization caused by the single-mode fiber will be lossless, and as a consequence, the states will always trace a great circle. As shown in FIG. 4, the Poincaré sphere representation of the states through which the light exiting the polarization-beam-combiner traverses as the laser is tuned. The great circle traversed by the light has poles that are defined by the states of polarization of the light in each of the “spinner” arms. The resultant light output on the single mode (SM) fiber from combiner 36 of the spinner 16 has a polarization state which varies linearly with the frequency of the light input to the spinner.

[0057] In order to monitor the source of light, in a preferred embodiment some of the light exiting the polarization beam combiner 36 of the spinner 16 is picked off by the single-mode coupler 38 as shown in FIG. 5. This picked off light is again split by a second single-mode coupler 39 (for example a 90/10 couple) and provides two light samples. One light sample (90% of the original beam) is used to trigger the acquisition of power at the final signal detectors in the spectral acquisition block 28 with the trigger detector 50. The trigger signal is generated by passing the light through a long Michelson interferometer 40 comprised of coupler 42, the trigger length of fiber Lt, and the two Faraday Rotator Mirrors (FRM's) 44. FRM's 44 are used in the interferometer to prevent polarization fading. In addition, the other light sample (10% of the original beam) is passed through a reference gas cell 46 and gas cell detector 48 to calibrate the laser wavelength to a National Institutes of Standards and Technology (NIST) traceable standard every sweep.

[0058] FIG. 5 illustrates the elements making up the spectral acquisition block 28 as shown by the dotted line. In FIG. 1, where the light sample is provided to the spectral acquisition block 28 along the dotted line 13, such line is shown in FIG. 5 as the line connecting couplers 38 & 39. In FIG. 5, the Michelson interferometer 40 is connected to the trigger detector 50 which produces sampling signals at very precise frequency intervals. The FRM's 44 prevent fringe fading for fiber positions. The Gas Cell 46 serves as a NIST traceable standard for both the frequency increment of the Michelson interferometer trigger, and the absolute wavelength of the measurement. As shown, the light that is not picked off is split by a 50/50 coupler 41 (also known as a 3 dB coupler) and passed both through the DUT 10, and the reference path 18.

[0059] FIG. 6 is a specific form of the generic invention shown in FIG. 1. The light from the two paths is then recombined in the coupler 22 in a preferred embodiment using a single-mode coupler 43. The recombined light is then passed to the polarizer 24 which includes an all-fiber polarization controller 52. This controller 52 forms an arbitrary lossless matrix that transforms one polarization state to another, without loss. The polarization controller could be located at any position between the spinner and the beam splitter 54 (except for along the DUT path). This transformation constitutes a “mapping” function. Adjusting the controller changes this matrix, but maintains its lossless nature. Mapping transforms the great-circle circumscribed on the Poincaré-sphere by the incident light to any other great circle. After passing through the polarization controller 52, the recombined light is then incident on a polarization-beam-splitter 54. At the beam splitter 54, the light is split into its x- and y-polarization components. These components are then incident on the s and p detectors 56 and 58, corresponding to the x and y components, respectively.

[0060] While the detectors 56 & 58 are shown outside of the dotted line encompassing the polarizer 24, they could be included within the polarizer and their output signal provided to the spectral acquisition block 28 as shown in FIG. 5. Alternatively, they could be included in the spectral acquisition block 28. The location is not important—just the fact that the output of the detectors 56 & 58, along with the outputs of the trigger detector 50 and the gas cell detector 48, are provided to the microprocessor 26.

[0061] The polarization controller 52 can be adjusted to one of two possible alignments. The adjustment should be made with no light coming through the measurement path. In one polarization controller alignment—referred to as the “nulled alignment”—the two linear polarization states are mapped to the two linear polarization states of the polarization-beam-splitter 54 when they pass through the reference path. In this alignment, the relative phase of the two states of polarization at the polarization-beam-coupler 43 does not affect the power splitting ratio at the polarization-beam-splitter 54 for reference path light, and as a result, no fringes are observed when the laser is tuned.

[0062] In the other polarization controller alignment—referred to as the “maximum contrast alignment”—the great circle circumscribed by the reference light is made to intersect with the two linear polarization states of the beam splitter 54. Whenever the reference polarization state passes through one of the linear states of the beam splitter, all of the reference power is incident on one detector, and none is incident on the other detector. The reference power then oscillates between the s and p detectors 56 & 58, respectively, going to zero on one detector when the power is maximized on the other detector. This output in a preferred embodiment is provided to the spectral acquisition unit 28.

[0063] If we look at the Fourier transform of the fringes, we get a time domain response, with the time domain response of each Jones matrix element arriving at a different time. By properly selecting the complex valued points associated with each matrix element, we can then use an inverse Fourier transform to compute the frequency domain (spectral) Jones matrix elements. While each term is multiplied by a phase error, the result of this phase error is equivalent to an arbitrary section of lossless optical fiber added to the device, and thus has no effect on the measured properties of the device.

[0064] FIGS. 7A and 7B are illustrations of the amplitude versus frequency of the Fourier Transform of the interference fringes measured by the network in FIG. 6 for the “nulled” alignment. FIGS. 8A and 8B are illustrations of the amplitude of the Fourier Transform of the interference fringes measured by the network in FIG. 6 for the “maximized” alignment.

[0065] Theoretical Analysis

[0066] We can begin by representing the DUT or “device” as a simple diagonal 2×2 matrix with complex, frequency-dependent entries,

[0067] $\begin{array}{cc}\stackrel{\stackrel{\_}{\_}}{M_D}=\left[\begin{array}{cc}{m}_1\left(\omega \right)& 0\\ 0& m_2\left(\omega \right)\end{array}\right]& \mathrm{Eqn}.2\end{array}$

[0068] In order to measure the complete response of the DUT, it is necessary that two orthogonal polarizations be incident upon the device. These states cannot be incident simultaneously, because the linear superposition of these states will form a new, single polarization. In order to perform a rapid measurement that will completely characterize the device, it is desirable to use a polarization state that changes rapidly with frequency, such that the device response does not appreciably change for one complete cycle through polarization states. An optical device to produce this rapidly varying polarization using a PM splitter and a polarization spinner 36, has previously been disclosed in FIG. 2 as a “spinner” and would cost on the order of about $3000.

[0069] In the “spinner” of FIG. 2, the fiber-optic network produces a polarization state that is strongly dependent on wavelength. The use of Polarization Maintaining (PM) fiber and a polarization beam coupler 43 guarantees orthogonality of the recombined states, and results in a well-defined set of states through which the polarization evolves as the frequency of the light is changed.

[0070] The normalized field emerging from this device is then given by,

[0071] $\begin{array}{cc}{E}_0=\frac{1}{\sqrt{2}}\left[\begin{array}{c}1\\ e^{-i\omega {\tau }_p}\end{array}\right]& \mathrm{Eqn}.3\end{array}$

[0072] In order to meet the criteria above that the polarization state vary rapidly enough that the response of the DUT does not change appreciably through a single cycle, the delay, τ_p, in the network above must be greater than the impulse response length of the device, τ_h. The fiber network that transfers this field from the 50/50 beam splitter 41 to the DUT will have largely uncontrollable polarization properties. It is desirable to impose two restrictions on this network. First, it have no Polarization Dependent Loss (PDL), and second, the matrix describing the network be constant as a function of frequency over the range of frequencies that is of interest. Of these two requirements, the zero PDL criteria is the most difficult to satisfy. The matrix for such a network is,

[0073] $\begin{array}{cc}\stackrel{\stackrel{\_}{\_}}{M_T}=\left[\begin{array}{cc}{e}^{i\varphi }\cos \theta & -e^{-i\psi }\sin \theta \\ e^{i\psi }\sin \theta & e^{-i\varphi }\cos \theta \end{array}\right]e^{-i\omega {\tau }_T}& \mathrm{Eqn}.4\end{array}$
There will be an arbitrary lossless matrix between the polarization spinner 36, and the device, and then a second arbitrary lossless matrix between the device and the polarizing beam splitter 54, and so the total response matrix is given by,
{double overscore (M_R)}={double overscore (M_T2)}·{double overscore (M_D)}·{double overscore (M_TI)}  Eqn. 5
which yields a new matrix in which each entry may be nonzero,

[0074] $\begin{array}{cc}\stackrel{\stackrel{\_}{\_}}{M_R}=\left[\begin{array}{cc}a& b\\ c& d\end{array}\right]& \mathrm{Eqn}.6\end{array}$

[0075] FIG. 9 is an illustration of the optical path and some associated lossless translation matrices (rot [1], rot [2], and rot [3]). The electric field also undergoes a translation through a lossless matrix from the spinner 16 to the polarizing beam splitter 54. This reference field is then given by,

[0076] $\begin{array}{cc}{\stackrel{\rightharpoonup }{E}}_{\mathrm{ref}}=\left[\begin{array}{cc}{e}^{i\varphi }\cos \theta & -e^{-i\psi }\sin \theta \\ e^{i\psi }\sin \theta & e^{-i\varphi }\cos \theta \end{array}\right]\left[\begin{array}{c}1\\ e^{-i\omega {\psi }_p}\end{array}\right]\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}& \mathrm{Eqn}.7\end{array}$
The measured field is given by,

[0077] $\begin{array}{cc}{\stackrel{\rightharpoonup }{E}}_{\mathrm{meas}}=\left[\begin{array}{cc}a& b\\ c& d\end{array}\right]\left[\begin{array}{c}1\\ e^{-i\omega {\tau }_p}\end{array}\right]\frac{e^{-i\omega {\tau }_D}}{4}& \mathrm{Eqn}.8\end{array}$
and the field at the beam splitter 54 is given by

[0078] $\begin{array}{cc}\begin{array}{c}{\stackrel{\rightharpoonup }{E}}_{\mathrm{split}}=\left[\begin{array}{cc}{e}^{i\varphi }\cos \theta & -e^{-i\psi }\sin \theta \\ e^{i\psi }\sin \theta & e^{-i\varphi }\cos \theta \end{array}\right][\begin{array}{c}1\\ e^{-i\omega {\tau }_p}\end{array}]\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\\ \left[\begin{array}{cc}a& b\\ c& d\end{array}\right][\begin{array}{c}1\\ e^{-i\omega {\tau }_p}\end{array}]\frac{e^{-i\omega {\tau }_D}}{4}\end{array}& \text{Eqn. 9}\end{array}$
Multiplying the vector and matrix in the first term gives,

[0079] $\begin{array}{cc}\begin{array}{c}{\stackrel{\rightharpoonup }{E}}_{\mathrm{split}}=\left[\begin{array}{c}{e}^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \\ e^{i\psi }\sin \theta +e^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta \end{array}\right]\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\\ \left[\begin{array}{cc}a& b\\ c& d\end{array}\right]\left[\begin{array}{c}1\\ e^{-i\omega {\tau }_p}\end{array}\right]\frac{e^{-i\omega {\tau }_D}}{4}\end{array}& \mathrm{Eqn}.10\end{array}$
and the multiplying the vector and matrix in the second term gives,

[0080] $\begin{array}{cc}\begin{array}{c}{\stackrel{\rightharpoonup }{E}}_{\mathrm{split}}=\left[\begin{array}{c}{e}^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \\ e^{i\psi }\sin \theta +e^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta \end{array}\right]\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\\ \left[\begin{array}{c}a+b{e}^{-i\omega {\tau }_p}\\ c+d{e}^{-\mathrm{i\omega }{\tau }_p}\end{array}\right]\frac{e^{-i\omega {\tau }_D}}{4}\end{array}& \mathrm{Eqn}.11\end{array}$
The fields are then split at the beam splitter 54, and sent to separate detectors 56, 58 for power measurement; these fields are given by,

[0081] $\begin{array}{cc}{E}_s=\left(e^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \right)\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\left(a+b{e}^{-i\omega {\tau }_p}\right)\frac{e^{-i\omega {\tau }_D}}{4}& \mathrm{Eqn}.12\end{array}$
and

[0082] $\begin{array}{cc}{E}_p=\left(e^{i\psi }\sin \theta +e^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta \right)\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\left(c+d{e}^{-i\omega {\tau }_p}\right)\frac{e^{-i\omega {\tau }_D}}{4}& \mathrm{Eqn}.13\end{array}$
A signal proportional to the power incident on the s-detector 56 is given by multiplying the field by its complex conjugate,

[0083] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=[\left(e^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \right)\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\\ \left(a+b{e}^{-i\omega {\tau }_p}\right)\frac{e^{-i\omega {\tau }_D}}{4}]·[(e^{-i\varphi }\cos \theta -\\ e^{i\left(\psi +\omega {\tau }_p\right)}\sin \theta )\frac{e^{i\omega {\tau }_{\mathrm{ref}}}}{2}+\left(a^{*}+b^{*}{e}^{i\omega {\tau }_p}\right)\frac{e^{i\omega {\tau }_D}}{4}]\end{array}& \mathrm{Eqn}.14\end{array}$
Expanding this gives,

[0084] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left(1-e^{i\left(\psi +\varphi +\omega {\tau }_p\right)}\cos \theta \sin \theta -e^{-i\left(\psi +\varphi +\omega {\tau }_p\right)}\cos \theta \sin \theta \right)+\\ \frac{1}{16}\left({a}^2+{b}^2+{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}+a^{*}b{e}^{-i\omega {\tau }_p}\right)+\\ \frac{1}{8}\left(e^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \right)\left(a^{*}+b^{*}{e}^{i\omega {\tau }_p}\right)e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}+\\ \frac{1}{8}\left(e^{-i\varphi }\cos \theta -e^{i\left(\psi +\omega {\tau }_p\right)}\sin \theta \right)\left(a+b{e}^{-i\omega {\tau }_p}\right)e^{-i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}\end{array}& \mathrm{Eqn}.15\end{array}$
Reducing the expression to only its real parts,

[0085] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left[1-\sin \left(2\theta \right)\cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]+\\ \frac{1}{16}\left({a}^2+{b}^2+2\mathrm{Re}\left[{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}\right]\right)+\\ \frac{1}{8}\mathrm{Re}[\left(e^{i\varphi }\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta \right)\\ \left(a^{*}+b^{*}{e}^{i\omega {\tau }_p}\right)e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.16\end{array}$
And then expanding the third term,

[0086] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left[1-\sin \left(2\theta \right)\cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]+\\ \frac{1}{16}\left({a}^2+{b}^2+2\mathrm{Re}\left[{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}\right]\right)+\\ \frac{1}{8}\mathrm{Re}[(a^{*}{e}^{i\varphi }\cos \theta +b^{*}{e}^{i\omega {\tau }_p}{e}^{\mathrm{i\varphi }}\cos \theta -e^{-i\left(\psi +\omega {\tau }_p\right)}\\ \sin \theta a^{*}-e^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta b^{*}{e}^{\mathrm{i\omega }{\tau }_p})e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.17\end{array}$
And a little bit or reordering

[0087] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left[1-\sin \left(2\theta \right)\cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]+\frac{1}{16}({a}^2+{b}^2+\\ 2\mathrm{Re}\left[{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}\right])+\frac{1}{4}\mathrm{Re}[(a^{*}{e}^{i\varphi }\cos \theta -b^{*}{e}^{-i\psi }\sin \theta +\\ b^{*}{e}^{i\left(\varphi +\omega {\tau }_p\right)}\cos \theta -a^{*}{e}^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta )e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.18\end{array}$
Now, repeat for the p-polarization

[0088] $\begin{array}{cc}\begin{array}{c}{E}_p·E_p^{*}=[\left(e^{i\psi }\sin \theta +e^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta \right)\frac{e^{-i\omega {\tau }_{\mathrm{ref}}}}{2}+\\ \left(c+d{e}^{-i\omega {\tau }_p}\right)\frac{e^{-i\omega {\tau }_D}}{4}]·[(e^{-i\psi }\sin \theta +\\ e^{i\left(\varphi +\omega {\tau }_p\right)}\cos \theta )\frac{e^{i\omega {\tau }_{\mathrm{ref}}}}{2}+\left(c^{*}+d^{*}{e}^{i\omega {\tau }_p}\right)\frac{e^{i\omega {\tau }_D}}{4}]\end{array}& \mathrm{Eqn}.19\end{array}$
expanding and collecting real parts,

[0089] $\begin{array}{cc}\begin{array}{c}{E}_p·E_p^{*}=\frac{1}{4}\left[1+\sin \left(2\theta \right)\cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]+\frac{1}{16}[{c}^2+\\ {d}^2+2\mathrm{Re}\left[{\mathrm{cd}}^{*}{e}^{i\omega {\tau }_p}\right]]+2\mathrm{Re}[\frac{1}{8}[(e^{i\psi }\sin \theta +\\ e^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta )\left(c^{*}+d^{*}{e}^{i\omega {\tau }_p}\right)e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]]\end{array}& \mathrm{Eqn}.20\end{array}$
expanding the third term again

[0090] $\begin{array}{cc}\begin{array}{c}{E}_p·E_p^{*}=\frac{1}{4}\left[1+\sin \left(2\theta \right)\cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]+\frac{1}{16}[{c}^2+\\ {d}^2+2\mathrm{Re}\left[{\mathrm{cd}}^{*}{e}^{i\omega {\tau }_p}\right]]+\frac{1}{4}\mathrm{Re}[(c^{*}{e}^{i\psi }\sin \theta +\\ d^{*}{e}^{i\omega {\tau }_p}{e}^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta +d^{*}{e}^{i\omega {\tau }_p}{e}^{i\psi }\sin \theta +\\ c^{*}{e}^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta )e^{\mathrm{i\omega }\left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.21\end{array}$
and a final reordering

[0091] $\begin{array}{cc}\begin{array}{c}{E}_p·E_p^{*}=\frac{1}{4}\left[1+\sin \left(2\theta \right)\cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]+\frac{1}{16}[{c}^2+{d}^2+\\ 2\mathrm{Re}\left[c{d}^{*}{e}^{i\omega {\tau }_p}\right]]+\frac{1}{4}\mathrm{Re}[(c^{*}{e}^{i\psi }\sin \theta +d^{*}{e}^{-i\varphi }\cos \theta +\\ d^{*}{e}^{i\left(\psi +\omega {\tau }_p\right)}\sin \theta +c^{*}{e}^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta )e^{i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.22\end{array}$
So, our detected powers for the p and s detectors 58 & 56, respectively are,

[0092] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left[1-\sin \left(2\theta \right)\cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]+\frac{1}{16}({a}^2+{b}^2+\\ 2\mathrm{Re}\left[{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}\right])+\frac{1}{4}\mathrm{Re}[(a{e}^{-i\varphi }\cos \theta -b{e}^{i\psi }\sin \theta +\\ b{e}^{-i\left(\varphi +\omega {\tau }_p\right)}\cos \theta -a{e}^{i\left(\psi +\omega {\tau }_p\right)}\sin \theta )e^{-i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.23\\ \mathrm{and}& \\ \begin{array}{c}{E}_p·E_p^{*}=\frac{1}{4}\left[1+\sin \left(2\theta \right)\cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]+\frac{1}{16}\lfloor {c}^2+{d}^2+\\ 2\mathrm{Re}\left[{\mathrm{cd}}^{*}{e}^{i\omega {\tau }_p}\right]\rfloor +\frac{1}{4}\mathrm{Re}[(c{e}^{-i\psi }\sin \theta +d{e}^{i\varphi }\cos \theta +\\ d{e}^{-i\left(\psi +\omega {\tau }_p\right)}\sin \theta +c{e}^{i\left(\varphi +\omega {\tau }_p\right)}\cos \theta )e^{-i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.24\end{array}$

[0093] Here, we will digress a moment to discuss the parameter, θ, and what might be done about it. FIG. 10 is a schematic of a fiber-optic network for the measurement of the reflection response of a DUT. The path through the polarization controller 52 is the reference path. The polarization controller is adjusted so that the fringes measured at the detectors are maximized. This adjustment must be made with the DUT disconnected, and the laser sweeping. With the fringes maximized, the system is described as being “aligned.” An alternate alignment in which the fringes are nulled is also possible.

[0094] As shown in FIG. 10, we can insert a polarization controller 52 in the reference path, and use it to manipulate the polarization state of the incident field incident on the beam splitter. If we zero out the return signal from the device under test (bend the lead fiber), The fields incident on the detectors are

[0095] $\begin{array}{cc}{E}_s·E_s^{*}=\frac{1}{4}\left[1-\sin \left(2\theta \right)\cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]& \mathrm{Eqn}.25\\ \mathrm{and}& \\ E_p·E_p^{*}=\frac{1}{4}\left[1+\sin \left(2\theta \right)\cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]& \mathrm{Eqn}.26\end{array}$

[0096] By adjusting the polarization controller 52 to achieve full contrast, θ is set to

[0097] $\frac{\pi }{4}+\frac{q\pi }{2}$
where q is an integer and represents an uncertainty in the sign of the trigonometric function. The detected fields are then given by,

[0098] $\begin{array}{cc}\begin{array}{c}{E}_s·E_s^{*}=\frac{1}{4}\left[1\pm \cos \left(\psi +\varphi +\omega {\tau }_p\right)\right]+\frac{1}{16}({a}^2+{b}^2+\\ 2\mathrm{Re}\left[{\mathrm{ab}}^{*}{e}^{i\omega {\tau }_p}\right])+\frac{1}{4\sqrt{2}}\mathrm{Re}[(a{e}^{-i\left(\varphi +k\pi \right)}-b{e}^{i\left(\psi +l\pi \right)}+\\ b{e}^{-i\left(\varphi +k\pi +\omega {\tau }_p\right)}-a{e}^{i\left(\psi +l\pi +\omega {\tau }_p\right)})e^{-i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}\mathrm{and}& \mathrm{Eqn}.27\\ \begin{array}{c}{E}_p·E_p^{*}=\frac{1}{4}\left[1\mp \cos \left(\varphi +\psi +\omega {\tau }_p\right)\right]+\frac{1}{16}\lfloor {c}^2+{d}^2+\\ 2\mathrm{Re}\left[{\mathrm{cd}}^{*}{e}^{i\omega {\tau }_p}\right]\rfloor +\frac{1}{4\sqrt{2}}\mathrm{Re}[(c{e}^{-i\left(\psi -l\pi \right)}+d{e}^{i\left(\varphi -k\pi \right)}+\\ d{e}^{-i\left(\psi -l\pi +\omega {\tau }_p\right)}+c{e}^{i\left(\varphi -k\pi +\omega {\tau }_p\right)})e^{-i\omega \left({\tau }_D-{\tau }_{\mathrm{ref}}\right)}]\end{array}& \mathrm{Eqn}.28\end{array}$
where k and l are arbitrary integers.

[0099] Examining these detected signals it can be seen that the elements of the matrices have been shifted to separate locations in frequency space, with a phase shift. The relative phases between the elements in a Jones matrix are important and must be known to some degree as discussed later. So, by selecting terms based upon their location in frequency space, and, also making a best guess at the polarization spinner delay, τ_p we can extract the following field responses,
E_a=−ae^{i(ψ+lπ+ωΔτp)}  Eqn. 29
E_b=be^{−i(φ+kπ+ωΔτp)}  Eqn. 30
E_ba=ae^−i(φ+kπ)−be^i(ψ+lπ)  Eqn. 31
E_c=ce^{i(φ−kπ+ωΔτp)}  Eqn. 32
E_d=de^{−i(ψ−lπ+ωΔτp)}  Eqn. 33
E_dc=ce^{−i(ψ−lπ)}+de^i(φ−kπ)  Eqn. 34
where we have allowed for an error in our estimate of the spinner delay of Δτ_p. We immediately have expressions for all of the matrix elements, however, an arbitrary phase factor, some of it frequency dependent, is associated with the phase factor.
a=−E_ae^{−i(ψ+lπ+ωΔτp)}  Eqn. 35
b=E_be^{i(φ+kπ+ωΔτp)}  Eqn. 36
c=E_ce^{−i(φ−kπ+ωΔτp)}  Eqn. 37
d=E_de^{i(ψ−lπ+ωΔτp)}  Eqn. 38

[0100] Since we maximized the fringes that oscillate at a frequency of ωτ_p, we can easily determine τ_p for each scan, thus allowing for variations in the spinner path imbalance from scan to scan due to temperature. With the ωτ_p term removed from the exponents, and any static phase error lumped into ψ and φ, we can then construct the matrix,

[0101] $\begin{array}{cc}\left[\begin{array}{cc}a{e}^{-i\left(\varphi +k\pi \right)}& b{e}^{i\left(\psi +l\pi \right)}\\ c{e}^{-i\left(\psi -l\pi \right)}& d{e}^{i\left(\varphi +k\pi \right)}\end{array}\right]=\left[\begin{array}{cc}{e}^{i\frac{\varphi -\psi +\left(k-l\right)\pi }{2}}& 0\\ 0& e^{i\frac{\varphi -\psi +\left(k-l\right)\pi }{2}}\end{array}\right]\left[\begin{array}{cc}a& b\\ c& d\end{array}\right]& \mathrm{Eqn}.39\\ \left[\begin{array}{cc}{e}^{i\frac{\varphi +\psi +\left(k+l\right)\pi }{2}}& 0\\ 0& e^{i\frac{\varphi +\psi +\left(k+l\right)\pi }{2}}\end{array}\right]& \end{array}$

[0102] Note that the measured matrix can be obtained from the original matrix by multiplication with two unitary matrices. This is not a transformation of basis, because the two matrices are not inverses of one another. As a result, diagonalizing the measured matrix does not recover the original matrix. All of the important properties of the matrix, Polarization Mode Dispersion, Polarization Dependent Loss, Spectral Amplitude and Group Delay are all maintained. All that has been altered is the polarization state at which a given minimum or maximum will occur. Since polarization in optical fiber is rarely well controlled, this loss of information regarding the device is not detrimental to the vast majority of applications.

[0103] In the following embodiments, the mathematics describing the device performance is not substantially altered, although the elements used to implement the measurement are altered. These alterations involve the removal of elements from claims and are not obvious.

[0104] No Polarization-Maintaining Fiber Embodiment

[0105] The embodiment described above in FIGS. 1 and 10 contains a PM fiber coupler 22. Proper operation of the device requires that the splitting ratio of the coupler be very well matched. Obtaining PM couplers is difficult and expensive. Obtaining well matched PM couplers is much more difficult and expensive. Polarization-beam-combiners are also difficult and expensive to come by. It would, therefore, be advantageous to eliminate these elements from the invention, and replace them with more easily obtainable devices.

[0106] At the polarization beam combiner 20, the great circle on the Poincaré sphere through which the polarization state rotates is well defined. Once the light has propagated through a few meters of optical fiber, however, the fiber will have transferred the polarization states to a new and completely arbitrary great circle. As a result, the well defined states produced by the PM fiber system serve no useful purpose to the measurement. We are then free to choose any set of orthogonal states to send through the two different length paths.

[0107] As disclosed in U.S. Pat. No. 4,389,090 “Fiber Optic Polarization Controller,” describes a simple apparatus that can be used to map any input polarization state to any other polarization state with little or no loss and little wavelength dependence. The polarization controller consists of three fiber loops that can be tilted independently, and is generally referred to as a set of LeFevre Loops after the inventor.

[0108] If we include a polarization controller of any kind (such as the LeFevre Loops discussed above) in a Mach-Zender Interferometer, we can adjust the controller such that the two states, when recombined at the coupler, are orthogonal to one another. This state of orthogonality is achieved by adjusting the polarization controller until one output of the coupler shows no fringes as the laser is tuned. The lack of interference fringes indicates that the two states of polarization entering the coupler through different paths have no common components, and are, therefore, orthogonal. With the creation of a polarization state that rotates rapidly through a great circle on the Poincaré sphere as a function of wavelength, we have achieved an equivalent function to the PM fiber/polarization-beam-combiner device described above. The remainder of the system then operates exactly as the previously described embodiment.

[0109] FIG. 11 shows an all single-mode (SM) fiber (i.e. not PM fiber) spinner module. By adjusting the polarization controller (LeFevre Loops 60) until no fringes are observed on the nulling detector 62 as the laser is swept, orthogonality of the polarization states in both arms of the “spinner” is guaranteed. The polarization state in one arm is, however, completely arbitrary.

[0110] Calibration of a single Spinner Optical Vector Analyzer

[0111] In order to describe the calibration procedure, it will be necessary to reference specific points of the optical network within one embodiment of the optical vector analyzer. A schematic diagram of the relevant portion of the network appears in FIG. 18 where the couplers are 50/50 3 dB fused-tapered, wavelength flattened couplers. In order to calibrate out the nonideal characterizations of the system, it is necessary to temporarily break one optical path in the spinner embodiment shown in FIG. 18. An optical switch 66 is therefore placed in the c-d path so that this break can be made and the connection re-established for normal operation. Additionally, although the polarization controller is shown in the fiber f-i, it could also be located at other positions such as the k-I fiber.

[0112] Before performing any measurements, and after aligning the optics, two measurement scans are conducted with no DUT attached. The first scan is done with the spinner switch open, and the DC signals on the detectors at points w and x are recorded. These DC signals are equal to the optical power incident on each detector multiplied by the gains of the detectors and their supporting electronics, and are frequency-dependent. We will denote them as
S_w^open(ω)=g_w(ω)P_w^open(ω)  Eqn. 40
and
S_x^open(ω)=g_x(ω)P_x^open(ω)  Eqn. 41
where P, g, and S refer to optical power, gain, and measured signal, respectively, and the subscripts refer to the two detectors. The superscript is a reminder that this information was acquired with the spinner switch open. In order for polarization-dependent loss to be properly calibrated, gw(ω) and gx(ω) must be equal. When the spinner switch is closed and a measurement is made, the optical signal at each detector will contain an AC term and a DC term:
S_w^closed(ω)=g_w(ω)[P_w^closed,dc(ω)+P_w^closed,ac(ω)e^iωτ]  Eqn. 42
S_X^closed(ω)=g_x(ω)[P_x^closed,dc(ω)+P_x^closed,ac(ω)e^iω]  Eqn. 43

[0113] In the above expressions ωis the optical frequency and τ is the difference between the times for an optical signal to propagate from point b to point e and from point c to point d. By taking the Fourier transform of these signals, the AC and DC terms can be separated, and then inverse transforms will bring the information back into the frequency domain. We obtain the following information from this scan:
S_w^closed,dc(ω)=g_w(ω)P_w^closed,dc(ω)  Eqn. 44
S_x^closed,dc(ω)=g_x(ω)P_x^closed,dc(ω)  Eqn. 45
S_w^closed,ac(ω)=g_w(ω)P_w^closed,ac(ω)  Eqn. 46
S_x^closed,ac(ω)=g_x(ω)P_x^closed,ac(ω)  Eqn. 47
where the superscripts now indicate both (a) which peak of the Fourier transform was selected, and (b) that the spinner switch was closed during data acquisition. We now define the following calibration quantities derived from these two scans:

[0114] $\begin{array}{cc}{C}_1\left(\omega \right)=\sqrt{S_w^{\mathrm{open}}}{e}^{i\left(\angle S_w^{\mathrm{closed},\mathrm{ac}}-\angle S_x^{\mathrm{closed},\mathrm{ac}}\right)}& \mathrm{Eqn}.48\end{array}$
C₂(ω)=√{square root over (S_x^open)}  Eqn. 49
C₃(ω)=√{square root over (S_w^closed,dc−S_w^open)}  Eqn. 50
C₄(ω)=√{square root over (S_x^closed,dc−S_x^open)}  Eqn. 51

[0115] Once the above information has been obtained, it will be used to divide all subsequent measurements in the following way. If the raw Optical Vector Analyzer Jones matrix measurement yields the matrix

[0116] $\begin{array}{cc}\lceil \begin{array}{cc}a\left(\omega \right)& b\left(\omega \right)\\ c\left(\omega \right)& d\left(\omega \right)\end{array}\rceil & \mathrm{Eqn}.52\end{array}$
then the matrix that will be used for subsequent calculations will be

[0117] $\begin{array}{cc}\left[\begin{array}{cc}a/C_1& b/C_3\\ c/C_2& d/C_4\end{array}\right]& \mathrm{Eqn}.53\end{array}$

[0118] During a measurement scan, the detector 62 at point h acquires data at the same time as the detectors 58 and 56 at points w and x, respectively. This data has the same form as Eqs. 42 and 43, and the AC component can be extracted in the same way. We will call this result S_h^ac(ω), and multiply the Jones matrix by the phase of S_h^ac(ω) as follows:

[0119] $\begin{array}{cc}\left[\begin{array}{cc}\frac{a}{C_1}{e}^{-i\angle S_h^{\mathrm{ac}}}& \frac{b}{C_3}{e}^{i\angle S_h^{\mathrm{ac}}}\\ \frac{c}{C_2}{e}^{-i\angle S_h^{\mathrm{ac}}}& \frac{d}{C_4}{e}^{i\angle S_h^{\mathrm{ac}}}\end{array}\right]& \mathrm{Eqn}.54\end{array}$

[0120] The next step in calibration is to perform a measurement using a fiber patchcord as the DUT. The resulting matrix, upon being divided by the calibration quantities as in Eqns. 48–51, is stored in memory and is called the reference matrix. It will be denoted as

[0121] $\begin{array}{cc}{\stackrel{\_}{M}}_R=\left[\begin{array}{cc}\frac{a^{\prime }}{C_1}{e}^{-i\angle S_h^{{\mathrm{ac}}^{\prime }}}& \frac{b^{\prime }}{C_3}{e}^{i\angle S_h^{{\mathrm{ac}}^{\prime }}}\\ \frac{c^{\prime }}{C_2}{e}^{-i\angle S_h^{{\mathrm{ac}}^{\prime }}}& \frac{d^{\prime }}{C_4}{e}^{i\angle S_h^{{\mathrm{ac}}^{\prime }}}\end{array}\right]& \mathrm{Eqn}.55\end{array}$
Once the reference matrix has been acquired, the instrument is now considered to be calibrated. All subsequent measured matrices (in the form of Eqn. 54) are multiplied on the right by the inverse of M_R to yield the final measured Jones matrix J for the DUT, where:

[0122] $\begin{array}{cc}\stackrel{\_}{J}=\frac{1}{a^{\prime }{d}^{\prime }-b^{\prime }{c}^{\prime }}\left[\begin{array}{cc}\frac{a}{C_1}{e}^{-i\angle S_h^{\mathrm{ac}}}& \frac{b}{C_3}{e}^{i\angle S_h^{\mathrm{ac}}}\\ \frac{c}{C_2}{e}^{-i\angle S_h^{\mathrm{ac}}}& \frac{d}{C_4}{e}^{i\angle S_h^{\mathrm{ac}}}\end{array}\right]& \text{Eqn. 56}\\ \left[\begin{array}{cc}\frac{d^{\prime }}{C_4}{e}^{i\angle S_h^{{\mathrm{ac}}^{\prime }}}& -\frac{b^{\prime }}{C_3}{e}^{i\angle S_h^{{\mathrm{ac}}^{\prime }}}\\ -\frac{c^{\prime }}{C_2}{e}^{-i\angle S_h^{{\mathrm{ac}}^{\prime }}}& \frac{a^{\prime }}{C_1}{e}^{-i\angle S_h^{{\mathrm{ac}}^{\prime }}}\end{array}\right]& \end{array}$

[0123] All Coupler/Polarization Controller Embodiment

[0124] The polarization beam-splitter 54 remaining in the design in the description above not only adds cost, but places very tight restrictions on the coupler tolerances in the spinner interferometer. A problem with the above embodiment is that a polarization controller 52 and polarizing beam-splitter 54 would be required for each additional receiving channel if multiple detection channels were required.

[0125] These difficulties can be resolved using a double spinner embodiment without polarization control as shown in FIGS. 12 & 13. The light source 14 supplies either a varying frequency light (for example, from a tunable laser source sweeping through the frequency of interest) or a broadband light (including the frequency range of interest) to two spinners 16a and 16b. It is understood that the reference path in the form of the dotted line 13 (shown in FIG. 1) would be needed in the variable frequency light source embodiment, it has been eliminated from FIG. 12 to illustrate a broadband light embodiment with an optical spectrum analyzer (OSA), in this instance a grating spectrometer 68. Other OSA's include a Fourier-Transform spectrometer and a scanning Fabry-Perot spectrometer.

[0126] In the FIGS. 12 & 13 embodiments, the output of each spinner provides light whose polarization varies linearly with frequency, the difference between the two spinners is that the rate of change of polarization with respect to a change in frequency is different. This difference in the rate of change can be defined by a path length difference, or other characteristics of the spinner.

[0127] The outputs of spinners 16a and 16b are applied to the DUT 10, and the reference path 18, respectively. The two outputs are connected to the coupler 22 with the output detected and applied to the microprocessor 26 as before. However, because of the difference in polarization rate of change, there is no need to utilize a polarizer in the processing of the outputs. As with the embodiment of FIG. 1, if the narrow linewidth source, such as a tunable laser, is used, a light sample is supplied to the spectral acquisition block 28.

[0128] One specific embodiment of the double spinner optical vector analyzer is shown in FIG. 13. One Mach-Zender interferometer is used for each spinner. One arm serves as the reference, and the other arm includes the device under test. Each interaction between the four polarizations can be made to have a unique apparent location in the transform of the interferometer fringes. In the FIG. 13 embodiment of the present invention, there is shown a coupler/polarization-controller network for full Jones matrix measurements measured by single device detector 64. Although adjustment of the fringes is not absolutely necessary, in a preferred embodiment, the interference terms occurring at L1, L1+L2, L1−L3 and L1−L3+L2 on detector 62 are adjusted to be non-zero using polarization controller 52.

[0129] Two of the polarization interactions are between different arms of the same spinner interferometer. These interactions are nulled using the polarization controllers in the interferometers (shown as LeFevre Loops 60 to distinguish from the nulling polarization controller 52 although they may all be LeFevre Loops). There are then four interactions remaining, and the four interactions are measurements of the four Jones matrix elements. Each matrix element is scaled by constant. These scaling values can be measured relatively easily, and, since they are determined by coupler splitting-ratios and splice losses, they should be stable over very long periods of time. The ability to calibrate each path loss out of the measurement should make it possible to obtain PDL measurements accurate to better than 0.05 dB.

[0130] FIG. 14 is an illustration of the amplitude versus frequency of the Fourier transform of the two peaks resulting from interference in the smaller interferometers in FIG. 13. By making the differential delay in each interferometer different, due to the increased path length of L1, the interference terms for each is readily separable.

[0131] A further advantage of this measurement method is the ability to use a single detector 64 for the measurement. As a result, gain matching between multiple detectors is not an issue, and relatively low precision (adjustable gain) components can be used. The use of a single device detector 64 also reduces the channel count required for a functioning system from four to three, meaning that a standard 4-channel board will now support a transmission and return-loss measurement.

[0132] FIG. 15 shows the amplitude versus frequency of the Fourier transform of the interference fringes at the device detector 64. Note that the differential delays for each of the small interferometers must be unequal to prevent the overlap of matrix terms, and that the optimal difference is that one delay should be twice the other. The spectral matrix entries are obtained by applying an inverse Fourier transform to each of the time domain responses.

[0133] Including all four elements in a single channel measurement does have the disadvantage of requiring more temporal range, and greater separation between the device, and any possible reflections (primarily connector reflections). In general, if a five nanosecond (5 ns) response is allowed (equivalent to 200 MHz sample spacing), then the fiber leads must be a minimum of 2 meters long. 3 meter leads would permit a 7 ns response time.

[0134] Scalability—Multilple Test Station Architecture

[0135] Because the double-Mach-Zender configuration of FIG. 13 does not have any defined polarization states, and thus no polarization controller associated with the measurement detector, it is easily scalable. Perhaps the first way in which this scaling can be applied would be to measure multiple outputs of a single device such as a Dense Wavelength Division Multiplexer (DWDM), a de-multiplexer, or an optical switch. FIG. 16 is a schematic of a single source, multiple detector, Jones matrix measurement system. Note that each additional detection channels require only a coupler and a detector 64, and no polarizing or polarization control elements.

[0136] A second application of the scalability of the double-Mach-Zender configuration would be to distribute the source light and the measurement light to multiple measurement stations for a plurality of devices under test (DUT 10) in a single facility. In this way the cost of the laser could be spread over multiple measurement stations, resulting in an extremely cost effective solution. FIG. 17 is a schematic of a multiple device/single source Jones Matrix measurement system. Note that each additional detection channel requires only a coupler and a detector, and no polarizing or polarization control elements.

[0137] In view of the various embodiments and discussion of the present invention noted above, many variation and substitutions of elements for those disclosed will be obvious to one of ordinary skill in the art. Accordingly, the invention is limited only by the claims appended hereto.