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

Data rate measurement method with selectable range in linear optical sampling

Open Access Open Access

Abstract

Linear optical sampling (LOS) is one of the most powerful techniques for high-speed signal monitoring. To measure the data-rate of signal under test (SUT) in optical sampling, multi-frequency sampling (MFS) was proposed. However, the measurable data-rate range of the existing method based on MFS is limited, which makes it very difficult to measure the data-rate of high-speed signals. To solve the above problem, a range selectable data-rate measurement method based on MFS in LOS is proposed in this paper. Through this method, the measurable data-rate range can be selected to match the data-rate range of SUT and the data-rate of SUT can be measured precisely, independently of the modulation format. What’s more, the sampling order can be judged using the discriminant in the proposed method, which is key for plotting eye diagrams with correct time information. We experimentally measure the baud-rates of PDM-QPSK signal from 800 MBaud to 40.8 GBaud in different ranges and judge the sampling orders. The relative error of measured baud-rate is less than 0.17% while the error vector magnitude (EVM) is less than 0.38. Compared with the existing method, under the same sampling cost, our proposed method realizes the selectivity of the measurable data-rate range and the judgment of sampling order, greatly extends the measurable data-rate range of SUT. Hence, the data-rate measurement method with selectable range has great potential for high-speed signal data-rate monitoring.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

With the development of information technology, the data transmission rate of optical communication system is increasing rapidly. The single channel data rate of optical transmission has reached 800 Gbit/s [1] or even higher [2,3,4], but monitoring high-speed signals is facing challenges. The traditional electrical-sampling based optical signal analyzer needs to use ultra-high speed photoelectric detectors and analog-to-digital conversion to complete high-speed signal monitoring, so this way has high requirements for hardware, high manufacturing costs, and is difficult to realize due to electronic bottleneck. Linear optical sampling (LOS) is a powerful method for high-speed optical signal monitoring [5,6]. It uses an optical pulse with low repetition frequency to sample the high-speed signal and reconstruct eye diagram and constellation diagram through software synchronization [7]. LOS realizes the linear optical process by coherently mixing the short optical pulse of the mode-locked fiber laser (MFL) with the optical signal. After linear mixing, LOS performs photoelectric conversion and sampling with low-speed balanced photodetector (BPD). Therefore, LOS obtains high-speed signal characteristics through undersampling scheme, which uses low-speed equipment to reduce hardware overhead. Many scholars have made contributions to improving the representation ability of LOS system, improving the accuracy of eye diagram reconstruction, and reducing the computational complexity [8,9,10,11]. With the wide application of wavelength division multiplexing (WDM), the fiber-optics-frequency-comb (FOFC) is proposed to replace the narrowband MFL as an optical sampling source to extend the range of operation wavelength for the LOS-based optical modulation analyzer (OMA). Since the use of BPDs introduces time bias and amplitude bias for LOS system, Liao proposed a scheme to mitigate the enhanced sampling error arising in the non-ideal response of a balanced photodetector [12]. But this procedure still introduces noise due to the constraints of the unideal response of BPDs. To overcome such issue, Zhe Yu presented a bias balance detection (BBD) scheme and demonstrated a 32 GBaud PDM-QPSK experiment for a FOFC based LOS system [13]. Meanwhile, MFL with a fixed repetition frequency has a theoretical upper limit of the measurement range [19] and the low repetition frequency of MLF increases the sampling period, which increases the constructing time of the eye diagram. To increase the repetition frequency, our team proposed a dual-pulse mixing (DPM) based LOS system using fiber delay lines with the multiplied optical pulse, which extends the measurement range and reduces error vector magnitude (EVM) bias by 9.1% [14]. All of the above studies have improved the performance of LOS in terms of measurement range, accuracy and reliability, respectively. However, in the above studies, the entire process of LOS is transparent to the data-rate, which does not provide information on the length and order of time for constructing eye diagram. In order to comprehensively monitor the optical signal, the ability of optical signal monitoring system to measure different data-rates of signal under test (SUT) is preferable.

To realize the data-rate measurement in optical sampling, multi-frequency sampling (MFS) was proposed. Xiangyu Wu proposed a bit-rate measurement method for optical performance monitoring, by alternatively tuning the repetition frequency of the optical pulse and differentially software-synchronized processing algorithm [15]. However, only several SUT with known discrete points in the low-speed range can be correctly measured like 10 Gbit/s and 40 Gbit/s in their experiment due to the existence of blind region for measurement [16]. Focus on this problem, Lin Zuo proposed a method to avoid the blind region. Through three times of optical sampling with slightly different repetition frequencies, this method calculates twice and takes the maximum result as the correct one. Almost all the bit-rate values in the whole low-speed range can be correctly measured, which achieved the monitoring of an unknown bit-rate in the low-speed range truly [17,18]. However, they found that after avoiding the blind region, the real measurable data-rate range was only a quarter of what was originally proposed in literature [15]. This makes it more difficult to monitor the data-rate of higher speed signal. There has been no major progress in linear optical sampling for the data-rate measurement of SUT in recent years. What’s more, in conventional LOS, the recovered eye diagram waveform may be the mirror image of the real waveform because sampling order cannot be judged, resulting in the error of information on time order in eye diagram. Until now, as far as we know, no method that can judge the sampling order in LOS has been proposed.

In this work, we propose a range selectable data-rate measurement method based on MFS in LOS. We derived the formula and found that the measurable range of data-rate can be greatly extended by adding a frequency selection factor. Firstly, select a measurable data-rate range, then use three laser pulses with different repetition frequencies to sample SUT for three times. Using the data obtained from each sampling, the required parameters can be obtained in software synchronization after applying a serious of DSP algorithms. Then the data-rate of SUT can be converted and the sampling order can be judged using our proposed formula, which is the key to plot the eye diagram with correct time information. We experimentally measure the baud-rates of QPSK signal from 0.8 to 40.8 GBaud and judge the sampling orders, when the EVM is less than 0.38, the relative error of the measured baud-rate is less than 0.17%. The rest of this paper is organized as follows. The theoretical mode and principle of range selectable data-rate measurement method is introduced in Section 2. The methods of avoiding blind regions and judging the sampling order are discussed in Section 3. The experiment for a triple-frequency LOS system is demonstrated and the performance of range selectable data-rate measurement method is analyzed in Section 4. Finally, the conclusion is presented.

2. Theoretical mode and principle of range selectable data-rate measurement method

LOS system generally adopts a MFL to generate a pulse with a repetition frequency of about one hundred megahertz to realize under-sampling for SUT [5]. Assuming that the data-rate of SUT is B and the repetition frequency of sampling pulse is f, the frequency reduction coefficient D can be calculated by the following formula [15]:

$$D = \frac{B}{f} = M + X,$$
where $M\textrm{ = }\lfloor{{B / f}} \rfloor$ is the integer part of the frequency reduction coefficient, which is recorded as the frequency reduction factor, and X is the decimal part of the frequency reduction coefficient, recorded as the frequency reduction remainder. Define the equivalent step length of sampling as:
$$dt = \frac{1}{f} - \frac{M}{B}.$$

The pulse samples the signal equivalently on the envelope with dt as the time step length [4]. Assuming that the signal period is T, the equivalent sampling points in one signal period are:

$$n = \frac{T}{{dt}}.$$

For an N-point-long sampled data stream, the number of equivalently scanned bit slots S can be obtained by the following formula:

$$S = \frac{N}{n} = \frac{{Ndt}}{T}.$$

After one sampling in LOS, without knowing B, the number of equivalently scanned bit slots S can be accurately obtained in the chirp-Z transform based software synchronization method in DSP [10]. After S is obtained, the N-point-long sampled data stream can be divided and superimposed to reconstruct the eye diagram [6].

In MFS based LOS, a mode-locked fiber laser with adjustable repetition frequency is used to generate two pulses with different repetition frequencies, recorded respectively as ${f_1},{f_2}$. For the convenience of discussion, ${f_1} > {f_2}$ is taken. The SUT with data-rate B is sampled by these two pulses, and (5) can be obtained:

$$\frac{{{S_1}}}{{B{N_1}}} = \frac{1}{{{f_1}}} - \frac{{{M_1}}}{B}$$
$$\frac{{{S_2}}}{{B{N_2}}} = \frac{1}{{{f_2}}} - \frac{{{M_2}}}{B}.$$

The subscripts 1 and 2 of the above variables represent that the variable is obtained after the SUT was sampled by pulse with repetition frequency ${f_1}$ and ${f_2}$. Because ${f_1} > {f_2}$, (6) can be obtained from Eq. (1):

$${M_2} = {M_1} + P,$$
where P is the frequency selection factor, which is a non-negative integer. Then (7) is obtained:
$${D_1} = \frac{B}{{{f_1}}} = {M_1} + {X_1}$$
$${D_2} = \frac{B}{{{f_2}}} = {M_1} + P + {X_2}.$$

The sufficient condition for P to be a non-negative integer is:

$${D_2} - {D_1} = \frac{B}{{{f_2}}} - \frac{B}{{{f_1}}} \ge P.$$

Equivalent to:

$$B > \frac{{{f_1}{f_2}P}}{{\Delta f}}.$$

Then the data-rate B can be converted according to (5) and (6):

$$B = \left|{\frac{{{{{S_1}} / {{N_1}}} - {{{S_2}} / {{N_2}}} - P}}{{{1 / {{f_1}}} - {1 / {{f_2}}}}}} \right|.$$

Formula (10) can be obtained only when ${M_2} = {M_1} + P$ is established, the necessary condition is:

$$\left|{\frac{B}{{{f_2}}} - \frac{B}{{{f_1}}} - P} \right|< 1.$$

Combining (9) and (11), we can get:

$$\frac{{f_1}{f_2}P}{\Delta f} < B < \frac{{f_1}{f_2}(P+1)}{\Delta f}.$$

It can be found that when ${f_1},{f_2}$ is determined, the measurable data-rate range can be selected by the frequency selection factor P. So, we can select a suitable P to match the data-rate range of the SUT. The reported method in [17] is actually one case of our proposed method when manually selecting P = 0. The measurable data-rate range of SUT is limited by (12) and the range that does not meet with (12) is called the first blind region for measurement, in which the data-rate cannot be calculated correctly.

3. Blind region avoiding and sampling order judgement

3.1 Blind region avoiding

In the actual simulation, it is found that even if ${M_2} = {M_1} + P$ is met, there will be wrong results caused by the fact that the frequency reduction remainder ${X_1},{X_2}$ are in (0,0.5) and (0.5,1) respectively, or both are in (0.5,1), which are the second and third blind regions. In order to explain the formation for the second and third blind regions, it is necessary to explore the characteristics of fast Fourier transform (FFT) in software synchronization.

In coarse synchronization, due to the conjugate symmetry of FFT, there will be two completely symmetrical peaks in the spectrum of the sampled signal after variance transform, as shown in Fig. 1.

 figure: Fig. 1.

Fig. 1. Spectrum of sampled signal after variance transform.

Download Full Size | PDF

In fact, it is feasible to take the abscissa of any of these two peaks as the number of equivalently scanned bit slots S for eye diagram reconstruction, and eye diagrams with equal Q value can be obtained, as shown in Fig. 2. This is because the down frequency remainder $X = {S / N}$. When ${X_1},{X_2}$ are symmetrical about 0.5, it takes the same number of steps to sample a period of SUT. But the abscissa point of first peak is always taken as the number of equivalently scanned bit slots S in the software synchronization process. This means that when processing two optical sampled data streams with equal length and symmetrical down frequency remainders about 0.5, they will have the same S after software synchronization, rather than the original two S. Then it is impossible to distinguish which is the correct number of equivalently scanned bit slots from software synchronization.

 figure: Fig. 2.

Fig. 2. Eye diagram of 8GBaud QPSK signal. (a) Eye diagram constructed by taking the abscissa of first peak as S. (b) Eye diagram constructed by taking the abscissa of second peak as S. (a) and (b) share the equal Q value, Q = 8.162 dB.

Download Full Size | PDF

In practice, there are four situations as shown in (13) for the distribution of the down frequency remainder after two times of sampling,

$$\begin{aligned} 0 < {X_1} < 0.5, \quad 0.5 \lt {X_2} \lt 1 \\ 0 < {X_2} < 0.5, \quad 0.5 \lt {X_1} \lt 1 \\ 0 < {X_1} < 0.5, \quad 0 \lt {X_2} \lt 0.5\\ 0.5 < {X_1} < 1, \quad 0.5 \lt {X_2} \lt 1. \end{aligned}$$

Only when $0 < {X_1} < 0.5$, $0 < {X_2} < 0.5$, the measured data-rate is correct by taking the abscissa of the first peak as S. In the other three cases, due to the wrong selection of the abscissa, the measured data-rate will be less than or equal to the correct one, which is proved in detail in Appendix A.

To resolve the above issue, a method to avoid the blind region based on three times of sampling is proposed in this paper. Take ${f_1} > {f_2} > {f_3}$, ${f_1} - {f_2} = {f_2} - {f_3}\textrm{ = }\Delta f$, and the data-rate B of SUT should meet the condition in (14) at the same time:

$$\begin{aligned} P < {D_2} - {D_1} & = \frac{B}{{{f_2}}} - \frac{B}{{{f_1}}} < P + 0.25 \\ P < {D_3} - {D_2} & = \frac{B}{{{f_3}}} - \frac{B}{{{f_2}}} < P + 0.25 \\ 2P < {D_3} - {D_1} & = \frac{B}{{{f_3}}} - \frac{B}{{{f_1}}} < 2P + 0.5 \end{aligned}$$
equivalent to:
$$\begin{aligned} \frac{{{f_1}{f_2}P}}{{\Delta f}} < B < \frac{{{f_1}{f_2}(P + 0.25)}}{{\Delta f}} \\ \frac{{{f_2}{f_3}P}}{{\Delta f}} < B < \frac{{{f_2}{f_3}(P + 0.25)}}{{\Delta f}} \\ \frac{{P{f_1}{f_3}}}{{\Delta f}} < B < \frac{{{f_1}{f_3}(P + 0.25)}}{{\Delta f}}. \end{aligned}$$

Then take the intersection of the ranges defined by formulas in (15) as:

$$\frac{{{f_1}{f_2}P}}{{\Delta f}} < B < \frac{{{f_2}{f_3}(P + 0.25)}}{{\Delta f}}.$$

When the data-rate B of SUT meets the condition of (16), ${S_1}$, ${S_2}$ and ${S_3}$ are calculated respectively in software synchronization after three times of sampling. Then use the data obtained from the first and second sampling, and the data obtained from the second and third sampling to calculate the data-rate respectively

$${B_{12}} = \left|{\frac{{{{{S_1}} / {{N_1}}} - {{{S_2}} / {{N_2}}} - P}}{{{1 / {{f_1}}} - {1 / {{f_2}}}}}} \right|$$
$${B_{21}} = \left|{\frac{{{{{S_2}} / {{N_2}}} - {{{S_1}} / {{N_1}}} - P}}{{{1 / {{f_2}}} - {1 / {{f_1}}}}}} \right|.$$

The data-rate calculated according to (17) is recorded as ${B_{12}}$, and the data-rate calculated according to (18) is recorded as ${B_{21}}$. Similarly, ${B_{23}}$ and ${B_{32}}$ can be obtained. These four results will finally be used to determine the correct data-rate and judge the sampling sequence.

To analyze how to avoid the blind region, the six different distributions of frequency reduction coefficient D on the number axis are plotted as shown in Fig. 3.

 figure: Fig. 3.

Fig. 3. Distributions of frequency reduction coefficient D on the number axis. In each distribution, the results of B are different. (a) B12 wrong, B21 correct, B23 wrong, B32 correct. (b) B12 wrong, B21 correct, B23 wrong, B32 wrong. (c) B12 wrong, B21 wrong, B23 correct, B32 wrong. (d) B12 correct, B21 wrong, B23 correct, B32 wrong. (e) B12 correct, B21 wrong, B23 wrong, B32 wrong. (f) B12 wrong, B21 wrong, B23wrong, B32 correct.

Download Full Size | PDF

In the six cases shown in Fig. 3, at least one correct B can be found from ${B_{12}},{B_{21}},{B_{23}}$ and ${B_{32}}$. As it is proved in Appendix A that when in the blind region, the measured result is always less than or equal to the correct one, we only need to take the maximum of ${B_{12}},{B_{21}},{B_{23}}$ and ${B_{32}}$ as the correct data-rate B, that is:

$$B = \max ({B_{12}},{B_{21}},{B_{23}},{B_{32}}).$$

3.2 Sampling order judgement

When $X \in ({0, 0.5} )$, the sampling order is sequential, and the recovered waveform in eye diagram is the same as the real waveform. But when $X \in ({0.5,1} )$, the sampling order is reverse, which will cause the recovered waveform to be the mirror image of the real waveform, resulting in the error of time order, as shown in the following Fig. 4.

 figure: Fig. 4.

Fig. 4. SUT is drawn as a blue solid line, while the two optical sampling pulses are drawn as a red solid line and an orange solid line. Because their down frequency remainders are located at (0,0.5) and (0.5, 1) respectively, their reconstructed eye images are mirror images of each other.

Download Full Size | PDF

Now we analyze how to judge the sampling order. As it can be seen from Fig. 3 that the correct B can be obtained only when the both ${X_1}$ and ${X_2}$ are in $({0,0.5} )$ or $({0.5,1} )$, and at most one in ${B_{ij}}$ and ${B_{ji}}$ is correct. According to this feature, a method to judge the sampling order can be provided. When ${M_2} = {M_1}\textrm{ + }P$ and ${X_1}, {X_2} \in ({0.5,1} )$ are met at the same time, ${X_2} > {X_1}$ can be obtained according to ${B / {{f_2}}} = {M_1} + P + {X_2}$, ${B / {{f_1}}} = {M_1} + {X_1}$ and ${B / {{f_2}}} - {B / {{f_1}}} > P$. Since in software synchronization we take the abscissa of the first peak point as S, the data-rate converted from formula (18) is correct.

Remove the absolute value symbol in (18) to obtain ${B_{Jud}}$, which is called sampling order discriminant:

$${B_{Jud}} = \frac{{{{{S_2}} / {{N_2}}} - {{{S_1}} / {{N_1}}} - P}}{{{1 / {{f_2}}} - {1 / {{f_1}}}}}.$$

When ${X_1}, {X_2} \in ({0.5,1} )$ and P > 0, ${B_{Jud}}$ is negative.

When ${X_1}, {X_2} \in ({0, 0.5} )$, the correct measured data-rate is ${B_{12}}$, and remove the absolute value symbol in (17) to obtain:

$${B_{Jud}} = \frac{{{{{S_1}} / {{N_1}}} - {{{S_2}} / {{N_2}}} - P}}{{{1 / {{f_1}}} - {1 / {{f_2}}}}}.$$

When ${X_1}, {X_2} \in ({0, 0.5} )$, P > 0, ${B_{Jud}}$ is positive.

Then it can be judged that when ${B_{Jud}}$ is positive, the sampling order is sequential, and correspondingly when ${B_{Jud}}$ is negative, the sampling order is reverse. This method is effective only when P > 0 because the result of (10) after exchanging subscripts 1 and 2 are equal when P = 0.

As we known, the frequency selection factor P should be selected manually according to the prior knowledge of the data-rate range of SUT. For example, when f1 = 98.53 MHz, f2= 97.33 MHz, f3 = 96.13 MHz is set, the measurable data-rate range is shown in Table 1. When the data-rate of SUT is more than 7.991 GBaud and less than 9.746 GBaud, P = 1 should be selected during calculating the data-rate. But it should be noted that the width of the range determined by (16) gradually decreases with the increase of P. When the width of the range is equal to 0, there is:

$$\frac{{{f_1}{f_2}{P_{\max }}}}{{\Delta f}}\textrm{ = }\frac{{{f_2}{f_3}({P_{\max }} + 0.25)}}{{\Delta f}}.$$

Tables Icon

Table 1. measurable range corresponding to P

Simplify (22) to obtain:

$${P_{\max }}\textrm{ = }\frac{{{f_3}}}{{8\Delta f}},$$
which shows that there is an upper limit of P,
$$P \le {P_{\max }}.$$

Practically, ${P_{\max }}$ must not be selected to ensure that the width of the measurable range matches the data-rate range of SUT.

In the actual processing of data streams, in addition to data-rate measurement, some DSP algorithms for conventional coherent optical communication are also applied, such as normalization orthogonalization, polarization division demultiplexing, frequency offset estimation, carrier phase recovery, etc. The flow chart of all DSP algorithms applied in data-rate measurement is shown in Fig. 5.

 figure: Fig. 5.

Fig. 5. Schematic of the DSP for LOS system with range selectable data-rate measurement.

Download Full Size | PDF

The range selectable data-rate measurement method is actually a calculation of the transmission frequency of the time-domain waveform, which means that it is transparent to the modulation format. Thus, it is also feasible to measure the data-rate of signals of different modulation formats such as 16QAM, 64QAM, etc.

4. Experimental analysis

We conducted an experiment to verify the data-rate measurement method with selectable range described in Section 2. A custom MFL pump source, with operating frequency of 95.00 MHz to 100.00 MHz and adjustable accuracy of 0.01 MHz (gain spectrum width 20 nm), is used as a seed source. We set the pulse repetition frequency as f1 = 98.53 MHz, f2 = 97.33 MHz, f3 = 96.13 MHz. A fiber laser (EXFO IQS-636) is used to output a continuous wavelength laser with a center wavelength of 1550 nm and a linewidth of 100kHz as an optical carrier. The optical carrier is modulated by IQ Mach Zehnder modulator (IQ MZM), which is driven by the PRBS signal trace from the arbitrary waveform generator (Keysight M8192). The signal generated by the arbitrary waveform generator has been shaped by Raised Cosine filter with the roll-off coefficient α= 0.38. In the LOS system, we use a 2 × 8-wavelength hybrid (Kylia, COH24-X) to realize the mixing of SUT and optical pulse. The data-rates of SUT are set from 800 MBaud to 40.8 GBaud as in Table 2, and each SUT is sampled by pulses with different repetition frequencies for three times. The signal trace after the wavelength hybrid is shown in Fig. 6. After the photoelectric conversion in four BPDs (Thorlabs, BPD470C), a 4-channel data acquisition card (NI5160) is applied to collect the time waveform at a sampling rate of 1.25 G samples/s. After collecting the data, we use the DSP process in Fig. 5 to measure the data-rate and characterize SUT. When processing, the P is selected from 0 to 5 corresponding to the measurable data-rate range shown in Table 1. The baud-rate setting and the corresponding P are shown in Table 2. Since the baud-rate of a QPSK signal is equal to the bit-rate of channel I or channel Q, channel I of SUT is taken for data-rate measurement in experiment. For each SUT, it is feasible to use the data collected from any of the three samples for eye diagram reconstruction. In the experiment we use the data collected by the sampling with the repetition frequency of 96.13 MHz.

 figure: Fig. 6.

Fig. 6. Experimental setup for LOS system with range selectable data-rate measurement.

Download Full Size | PDF

Tables Icon

Table 2. Set baud-rate and corresponding P

The measurement results for each SUT are compared with the data-rate set in Table 2, and the accuracy performance of the data-rate measurement method is tested by relative error.

$$\textrm{Relative}_\textrm{error} = {\frac{{|{B^{\prime} - B} |}}{B}^{}}.$$

A comparison of the measured baud-rate and the real baud-rate is shown in Fig. 7(a). The data-rate set according to Table 2 is drawn as a red solid line in Fig. 7(a) while the measured data-rate is drawn as a blue hollow circle. In LOS, the quality of optical signal is characterized by EVM value. In order to test the performance of our proposed method under different signal quality, the EVM values of SUT at each data-rate are also applied for comparison.

 figure: Fig. 7.

Fig. 7. Experimental data analysis. (a) Real baud-rate B and measured baud-rate B. (b) EVM and relative error.

Download Full Size | PDF

In Fig. 7(b), the relative error of the measured data-rate is plotted as a red line, while the EVM value is plotted as a blue line. It can be seen that the quality of SUT is best at about 8 GBaud, and the EVM value gradually increases with the increase of the baud-rate. However, the relative error is not sensitive to the deterioration of signal quality, which shows that the proposed method has better applicability. Figure 7(b) shows that when the EVM is less than 0.38, the relative error of the measured baud-rate is less than 0.17%. This proves that our proposed method has a high accuracy even when the signal quality is very poor, which shows the reliability of the method. The reconstructed eye diagram with measured time information and constellation diagram of 8 GBaud QPSK signal using the repetition frequency of 96.13 MHz is shown in Fig. 8.

 figure: Fig. 8.

Fig. 8. Eye diagram with measured time information and constellation diagram of 8 GBaud QPSK signal, Q = 8.1620 dB, EVM = 0.1503.

Download Full Size | PDF

5. Conclusion

This paper proposes a range selectable data-rate measurement method based on MFS in LOS. In our proposed method, the measurement range can be selected according to the data-rate range of the SUT. In LOS system, a MFL is used to generate pulses with different repetition frequencies to sample the SUT for three times, and then the data acquired by the ADC is used to calculate the data-rate of the SUT in our proposed method. The data-rate of SUT can be measured without blind region in the selected range, and the sampling order can be judged, which is key for plotting eye diagrams with correct time information. Compared with the existing method, under the same sampling cost, our proposed method realizes the selectivity of the measurable data-rate range and the judgment of the sampling order, greatly extend the measurable data-rate range of SUT. We experimentally measure the data-rates of QPSK signal as high as 40.8 GBaud in a triple-frequency LOS system and judge the sampling order. When EVM is lower than 0.38, the relative error of measured baud-rate is less than 0.17%. Thus, the range selectable data-rate measurement method has great potential for future high-speed signal monitoring.

Appendix A

Since ${f_1}$ and ${f_2}$ are equal in status in multi-frequency sampling, for the convenience of discussion, ${f_1} > {f_2}$ is taken, the calculation formula of the measured data-rate ${B^{\prime}}$ is written as follows:

$${B^{\prime}} = \frac{{{f_1}{f_2}|{X_1} - {X_2} - P|}}{{{f_1} - {f_2}}}.$$

In the first case, when ${M_2} = {M_1}\textrm{ + }P$ is satisfied but not $({X_1} - 0.5)({X_2} - 0.5) > 0$, because of ${f_1} > {f_2}$, ${X_2} > {X_1}$ can be obtained according to ${B / {{f_2}}} = {M_1} + P + {X_2}$, ${B / {{f_1}}} = {M_1} + {X_1}$ and ${B / {{f_2}}} - {B / {{f_1}}} > P$, there must be ${X_2} \in ({0.5,1} )$ and ${X_1} \in ({0,0.5} )$. Since the abscissa of the first peak point is taken as the number of equivalently scanned bit slots S in software synchronization, ${X_2}$ should be replaced by $1 - {X_2}$ in (26), and ${B^{\prime}}$ is:

$${B^{\prime}} = \frac{{{f_1}{f_2}|{X_1} - 1 + {X_2} - P|}}{{{f_1} - {f_2}}}.$$

Using $\frac{B}{{{f_1}}} = {M_1} + {X_1}$, $\frac{B}{{{f_2}}} = {M_1} + P + {X_2}$ to eliminate ${X_2},{X_1}$ in (27):

$${B^{\prime}} = \frac{{|B{f_2} - {M_1}{f_1}{f_2} - {f_1}{f_2} + B{f_1} - {M_2}{f_1}{f_2} - P{f_1}{f_2}|}}{{{f_1} - {f_2}}}.$$

In order to prove the measured data-rate ${B^{\prime}}$ is less than or equal to the real data-rate B, it is only necessary to prove:

$$- B \le \frac{{B{f_2} - {M_1}{f_1}{f_2} - {f_1}{f_2} + B{f_1} - {M_2}{f_1}{f_2} - P{f_1}{f_2}}}{{{f_1} - {f_2}}} \le B.$$

Considering ${f_1} - {f_2} > 0$, (29) can be written as:

$$- B({f_1} - {f_2}) \le B{f_2} - {M_1}{f_1}{f_2} - {f_1}{f_2} + B{f_1} - {M_2}{f_1}{f_2} - P{f_1}{f_2} \le B({f_1} - {f_2}).$$

The proof of (30) is equivalent to the proof of (31) and (32):

$$2B - {M_1}{f_1} - {f_1} - {M_2}{f_1} - P{f_1} \le 0$$
$$0 \le - {M_1}{f_2} - {f_2} + 2B - {M_2}{f_2} - P{f_2}.$$

For (31), substitute ${M_2} = {M_1}\textrm{ + }P$ and simplify it to obtain:

$$2B - 2{M_2}{f_1} - {f_1} \le 0.$$

Formula (33) divides $2{f_1}$ by both sides and sort out:

$$\frac{B}{{{f_1}}} - {M_2} \le \frac{1}{2}.$$

Substitute ${B / {{f_1}}} = {M_1} + {X_1}$ into (34) to obtain:

$${M_1} + {X_1} - {M_2} \le \frac{1}{2}.$$

Substitute ${M_2} = {M_1}\textrm{ + }P$ into (35) and simplify it to obtain:

$${X_1} - P \le \frac{1}{2}.$$

According to ${X_1} \in ({0,0.5} )$, when P takes any natural number, (36) is obviously true, so (31) is true.

Then prove (32), substitute ${M_2} = {M_1}\textrm{ + }P$ into (32) and simplify it to obtain (37):

$$2B \ge {f_2} + 2{M_2}{f_2}.$$

Formula (37) divides $2{f_2}$ by both sides to get (38):

$$\frac{B}{{{f_2}}} \ge \frac{1}{2} + {M_2}.$$

Substitute ${B / {{f_2}}} = {M_2} + {X_2}$ into (38) and simplify it to:

$${X_2} \ge \frac{1}{2}.$$

Similarly, according to ${X_2} \in ({0.5,1} )$, easy to know (39) is obviously true, so (32) is true. Then ${B^{\prime}} \le B$ is certified in the first case.

In the second case, when neither ${M_2} = {M_1}\textrm{ + }P$ nor $({X_1} - 0.5)({X_2} - 0.5) > 0$ is satisfied, ${f_1} > {f_2}$, ${B / {{f_2}}} = {M_1} + P + {X_2}$, ${B / {{f_1}}} = {M_1} + {X_1}$ and ${B / {{f_2}}} - {B / {{f_1}}} > P$, so ${X_2} > {X_1}$. It is easy to get ${X_1} \in ({0.5,1} )$, ${X_2} \in ({0,0.5} )$, ${M_2} = {M_1}\textrm{ + }P\textrm{ + }1$. Because the software synchronization takes the abscissa of the first peak point as S, ${X_1}$ should be replaced with $1 - {X_1}$ in (26), and the measured data-rate ${B^{\prime}}$ is:

$${B^{\prime}} = \frac{{{f_1}{f_2}|1 - {X_1} - {X_2} - P|}}{{{f_1} - {f_2}}}.$$

Similarly, use ${B / {{f_1}}} = {M_1} + {X_1}$, ${B / {{f_2}}} = {M_1} + P + {X_2}$ to eliminate ${X_2},{X_1}$ from (40):

$${B^{\prime}} = \frac{{|{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - B{f_1} + {M_2}{f_1}{f_2} - P{f_1}{f_2}|}}{{{f_1} - {f_2}}}.$$

To prove ${B^{\prime}} \le B$, just prove:

$$- B \le \frac{{|{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - B{f_1} + {M_2}{f_1}{f_2} - P{f_1}{f_2}|}}{{{f_1} - {f_2}}} \le B.$$

Using simplification steps similar to the first case, (42) can be equivalent to (43) and (44).

$$\frac{B}{{{f_1}}} - {M_1} \le 1$$
$$P + {X_2} \ge 0.$$

Formula (43) is equivalent to ${X_1} \le 1$, which is obviously true from ${X_1} \in ({0.5,1} )$. Formula (44) is established because ${X_2} \in ({0,0.5} )$ and $P \ge 0$. Then ${B^{\prime}} \le B$ in the second case is proved.

In the third case, when both ${M_2} = {M_1}\textrm{ + }P$ and $({X_1} - 0.5)({X_2} - 0.5) > 0$ are satisfied, but ${X_1}, {X_2} \in ({0.5,1} )$, because ${f_1} > {f_2}$, ${B / {{f_2}}} = {M_1} + P + {X_2}$, ${B / {{f_1}}} = {M_1} + {X_1}$ and ${B / {{f_2}}} - {B / {{f_1}}} > P$, there is still ${X_2} > {X_1}$. Since the abscissa of the first peak point is taken as S in software synchronization, ${X_1} - {X_2}$ should be replaced by ${X_2} - {X_1}$ in (26), and the measured data-rate ${B^{\prime}}$ is:

$${B^{\prime}} = \frac{{{f_1}{f_2}|{X_2} - {X_1} - P|}}{{{f_1} - {f_2}}}.$$

Similarly, use ${B / {{f_1}}} = {M_1} + {X_1}$, ${B / {{f_2}}} = {M_1} + P + {X_2}$ to eliminate ${X_2},{X_1}$ in (45):

$${B^{\prime}} = \frac{{|B{f_1} - {M_2}{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - P{f_1}{f_2}|}}{{{f_1} - {f_2}}}.$$

To prove ${B^{\prime}} \le B$, just prove:

$$- B \le \frac{{|B{f_1} - {M_2}{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - P{f_1}{f_2}|}}{{{f_1} - {f_2}}} \le B.$$

Using simplification steps similar to the first case, (47) can be simplified to (48) and (49).

$$B{f_1} - {M_2}{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - P{f_1}{f_2} \le B({f_1} - {f_2})$$
$$B({f_1} - {f_2}) \le B{f_1} - {M_2}{f_1}{f_2} - B{f_2} + {M_1}{f_1}{f_2} - P{f_1}{f_2}.$$

Using simplification steps as those in the first and second cases, (48) can be simplified to:

$$0 \le P.$$

Easy to get that (50) must be established.

Similarly, formula (49) can be reduced to:

$${X_1} \le {X_2}.$$

It can be seen from ${X_1} < {X_2}$ that (51) must be established.

To sum up, ${B^{\prime}} \le B$ is certified in all cases.

Funding

National Key Research and Development Program of China (2018YFB1801705); National Natural Science Foundation of China (61727817, 62021005).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. M. Kong, C. Liu, B. Sang, K. Wang, J. Ding, J. Shi, L. Zhao, W. Zhou, X. Xin, B. Liu, B. Ye, W. Chen, and J. Yu, “Demonstration of 800-gbit/s/carrier tps-64qam wdm transmission over 2,000 km using mimo volterra equalization,” in Optical Fiber Communications Conference and Exhibition (2021), pp. 1–3.

2. F. Pittalà, R. Braun, G. Böcherer, P. Schulte, M. Schaedler, S. Bettelli, S. Calabrò, M. Kuschnerov, A. Gladischq, F. Westphal, C. Xie, R. Chen, Q. Wang, and B. Zheng, “1.71 Tb/s single-channel and 56.51 Tb/s DWDM transmission over 96.5 km field-deployed SSMF,” IEEE Photonics Technol. Lett. 34(3), 157–160 (2022). [CrossRef]  

3. V. Bajaj, F. Buchali, M. Chagnon, S. Wahls, and V. Aref, “Single-channel 1.61 Tb/s optical coherent transmission enabled by neural network-based digital pre-distortion,” in European Conference on Optical Communications (2020), pp. 1–4.

4. F. Wang, H. Yao, J. Wang, T. Mai, X. Xin, and M. Guizani, “Hybrid optical-electrical data center networking: Challenges and solutions for bandwidth resource optimization,” IEEE Commun. Mag. 60(11), 90–96 (2022). [CrossRef]  

5. C. Dorrer, D. Kilper, H. Stuart, G. Raybon, and M. Raymer, “Linear optical sampling,” IEEE Photonics Technol. Lett. 15(12), 1746–1748 (2003). [CrossRef]  

6. C. Dorrer, “Monitoring of optical signals from constellation diagrams measured with linear optical sampling,” J. Lightwave Technol. 24(1), 313–321 (2006). [CrossRef]  

7. M. Westlund, H. Sunnerud, M. Karlsson, and P. Andrekson, “Software-synchronized all-optical sampling for fiber communication systems,” J. Lightwave Technol. 23(3), 1088–1099 (2005). [CrossRef]  

8. H. Sunnerud, M. Skold, M. Westlund, and P. A. Andrekson, “Characterization of complex optical modulation formats at 100 gb/s and beyond by coherent optical sampling,” J. Lightwave Technol. 30(24), 3747–3759 (2012). [CrossRef]  

9. T. Kiatchanog, K. Igarashi, T. Tanemura, D. Wang, K. Katoh, and K. Kikuchi, “Real-time all-optical waveform sampling using a passively mode-locked fiber laser as the sampling pulse source,” in Optical Fiber Communication Conference and the National Fiber Optic Engineers Conference (2006), p. 3.

10. A. Yang, J. Lai, and Y.-N. Sun, “A chirp-z-transform-based software synchronization method for optical performance monitoring,” IEEE Photonics Technol. Lett. 23(22), 1739–1741 (2011). [CrossRef]  

11. J. Xue, Y. Wang, C. Li, J. Liu, L. Han, Q. Zhang, and X. Xin, “Integration Extraction Method Aided Linear Optical Sampling System,” in 20th International Conference on Optical Communications and Networks (2022), pp. 1–3.

12. R. Liao, Z. Wu, S. Fu, S. Zhu, Z. Yu, M. Tang, and D. Liu, “Fiber optics frequency comb enabled linear optical sampling with operation wavelength range extension,” Opt. Lett. 43(3), 439–442 (2018). [CrossRef]  

13. Z. Yu, S. Fu, H. He, Z. Wu, T. Huang, M. Tang, and D. Liu, “Biased balance detection for fiber optical frequency comb based linear optical sampling,” J. Lightwave Technol. 39(11), 3458–3465 (2021). [CrossRef]  

14. F. Wang, R. Gao, Z. Li, Q. Zhang, C. Yu, Q. Tian, Y. Wang, and X. Xin, “Fast linear optical sampling with high repetition-frequency using fiber delay lines,” Opt. Express 30(18), 32895–32907 (2022). [CrossRef]  

15. A. Yang, X. Wu, Y. Qiao, and Y. Sun, “Bit-rate adaptive optical performance monitoring method for fiber communication systems,” Opt. Commun. 284(1), 436–440 (2011). [CrossRef]  

16. A. Yang, X. Liu, J. Lai, L. Zuo, and Y. Sun, “Double frequency optical sampling for high-speed optical fiber communication systems,” in 8th International Conference on Communications and Networking in China (2013), pp. 815–821.

17. L. Zuo, L. Feng, A. Yang, J. Lai, and Y. Sun, “Experimental investigation on bit-rate-adaptive software synchronous optical sampling,” Chin. Opt. Lett. 11(5), 16–19 (2013).

18. L. Zuo, “Research on Key Technologies of Optical Performance Monitoring Based on All Optical Sampling”. Beijing: Beijing Institute of Technology, 2014: 99–118.

19. J. Song, S. Fu, B. Liu, M. Tang, P. Shum, and D. Liu, “Impact of sampling source repetition frequency in linear optical sampling,” IEEE Photonics Technol. Lett. 28(1), 15–18 (2016). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. Spectrum of sampled signal after variance transform.
Fig. 2.
Fig. 2. Eye diagram of 8GBaud QPSK signal. (a) Eye diagram constructed by taking the abscissa of first peak as S. (b) Eye diagram constructed by taking the abscissa of second peak as S. (a) and (b) share the equal Q value, Q = 8.162 dB.
Fig. 3.
Fig. 3. Distributions of frequency reduction coefficient D on the number axis. In each distribution, the results of B are different. (a) B12 wrong, B21 correct, B23 wrong, B32 correct. (b) B12 wrong, B21 correct, B23 wrong, B32 wrong. (c) B12 wrong, B21 wrong, B23 correct, B32 wrong. (d) B12 correct, B21 wrong, B23 correct, B32 wrong. (e) B12 correct, B21 wrong, B23 wrong, B32 wrong. (f) B12 wrong, B21 wrong, B23wrong, B32 correct.
Fig. 4.
Fig. 4. SUT is drawn as a blue solid line, while the two optical sampling pulses are drawn as a red solid line and an orange solid line. Because their down frequency remainders are located at (0,0.5) and (0.5, 1) respectively, their reconstructed eye images are mirror images of each other.
Fig. 5.
Fig. 5. Schematic of the DSP for LOS system with range selectable data-rate measurement.
Fig. 6.
Fig. 6. Experimental setup for LOS system with range selectable data-rate measurement.
Fig. 7.
Fig. 7. Experimental data analysis. (a) Real baud-rate B and measured baud-rate B. (b) EVM and relative error.
Fig. 8.
Fig. 8. Eye diagram with measured time information and constellation diagram of 8 GBaud QPSK signal, Q = 8.1620 dB, EVM = 0.1503.

Tables (2)

Tables Icon

Table 1. measurable range corresponding to P

Tables Icon

Table 2. Set baud-rate and corresponding P

Equations (53)

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

D = B f = M + X ,
d t = 1 f M B .
n = T d t .
S = N n = N d t T .
S 1 B N 1 = 1 f 1 M 1 B
S 2 B N 2 = 1 f 2 M 2 B .
M 2 = M 1 + P ,
D 1 = B f 1 = M 1 + X 1
D 2 = B f 2 = M 1 + P + X 2 .
D 2 D 1 = B f 2 B f 1 P .
B > f 1 f 2 P Δ f .
B = | S 1 / N 1 S 2 / N 2 P 1 / f 1 1 / f 2 | .
| B f 2 B f 1 P | < 1.
f 1 f 2 P Δ f < B < f 1 f 2 ( P + 1 ) Δ f .
0 < X 1 < 0.5 , 0.5 < X 2 < 1 0 < X 2 < 0.5 , 0.5 < X 1 < 1 0 < X 1 < 0.5 , 0 < X 2 < 0.5 0.5 < X 1 < 1 , 0.5 < X 2 < 1.
P < D 2 D 1 = B f 2 B f 1 < P + 0.25 P < D 3 D 2 = B f 3 B f 2 < P + 0.25 2 P < D 3 D 1 = B f 3 B f 1 < 2 P + 0.5
f 1 f 2 P Δ f < B < f 1 f 2 ( P + 0.25 ) Δ f f 2 f 3 P Δ f < B < f 2 f 3 ( P + 0.25 ) Δ f P f 1 f 3 Δ f < B < f 1 f 3 ( P + 0.25 ) Δ f .
f 1 f 2 P Δ f < B < f 2 f 3 ( P + 0.25 ) Δ f .
B 12 = | S 1 / N 1 S 2 / N 2 P 1 / f 1 1 / f 2 |
B 21 = | S 2 / N 2 S 1 / N 1 P 1 / f 2 1 / f 1 | .
B = max ( B 12 , B 21 , B 23 , B 32 ) .
B J u d = S 2 / N 2 S 1 / N 1 P 1 / f 2 1 / f 1 .
B J u d = S 1 / N 1 S 2 / N 2 P 1 / f 1 1 / f 2 .
f 1 f 2 P max Δ f  =  f 2 f 3 ( P max + 0.25 ) Δ f .
P max  =  f 3 8 Δ f ,
P P max .
Relative error = | B B | B .
B = f 1 f 2 | X 1 X 2 P | f 1 f 2 .
B = f 1 f 2 | X 1 1 + X 2 P | f 1 f 2 .
B = | B f 2 M 1 f 1 f 2 f 1 f 2 + B f 1 M 2 f 1 f 2 P f 1 f 2 | f 1 f 2 .
B B f 2 M 1 f 1 f 2 f 1 f 2 + B f 1 M 2 f 1 f 2 P f 1 f 2 f 1 f 2 B .
B ( f 1 f 2 ) B f 2 M 1 f 1 f 2 f 1 f 2 + B f 1 M 2 f 1 f 2 P f 1 f 2 B ( f 1 f 2 ) .
2 B M 1 f 1 f 1 M 2 f 1 P f 1 0
0 M 1 f 2 f 2 + 2 B M 2 f 2 P f 2 .
2 B 2 M 2 f 1 f 1 0.
B f 1 M 2 1 2 .
M 1 + X 1 M 2 1 2 .
X 1 P 1 2 .
2 B f 2 + 2 M 2 f 2 .
B f 2 1 2 + M 2 .
X 2 1 2 .
B = f 1 f 2 | 1 X 1 X 2 P | f 1 f 2 .
B = | f 1 f 2 B f 2 + M 1 f 1 f 2 B f 1 + M 2 f 1 f 2 P f 1 f 2 | f 1 f 2 .
B | f 1 f 2 B f 2 + M 1 f 1 f 2 B f 1 + M 2 f 1 f 2 P f 1 f 2 | f 1 f 2 B .
B f 1 M 1 1
P + X 2 0.
B = f 1 f 2 | X 2 X 1 P | f 1 f 2 .
B = | B f 1 M 2 f 1 f 2 B f 2 + M 1 f 1 f 2 P f 1 f 2 | f 1 f 2 .
B | B f 1 M 2 f 1 f 2 B f 2 + M 1 f 1 f 2 P f 1 f 2 | f 1 f 2 B .
B f 1 M 2 f 1 f 2 B f 2 + M 1 f 1 f 2 P f 1 f 2 B ( f 1 f 2 )
B ( f 1 f 2 ) B f 1 M 2 f 1 f 2 B f 2 + M 1 f 1 f 2 P f 1 f 2 .
0 P .
X 1 X 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.