DEVICE FOR WIDE-BAND SPECTRAL ANALYSIS OF A SIGNAL OF INTEREST

Abstract

A wideband device for the spectral analysis of a signal of interest includes a source designed to generate the signal of interest; an optical splitter element designed to spatially split the signal of interest into a first signal and a second signal; a first frequency-shifting optical cavity comprising a first frequency shifter designed to shift the optical frequency of the first signal by a first frequency f.sub.1 per round trip in the first cavity, the first cavity having a first trip time .sub.1; a second frequency-shifting optical cavity comprising a second frequency shifter designed to shift the optical frequency of the second signal by a second frequency f.sub.2 per round trip in the second cavity, the second cavity having a second trip time .sub.2; the first and the second optical cavity being designed such that a maximum number of round trips of the signal in the first and the second cavity is equal to predetermined N; a detector designed to coherently detect the first signal transmitted by the first cavity and the second signal transmitted by the second cavity and generate a photocurrent proportional to a luminous intensity detected by the detector, an analog low-pass filter designed to filter frequencies of the photocurrent that are lower than min (f.sub.12/; f.sub.2/2) a processor configured to compute a square modulus of the photocurrent filtered by the low-pass filter, from which a temporal representation of frequency information of the signal of interest is determined.

Claims

1. A wideband device (D) for the spectral analysis of a signal of interest comprising: a source (S) designed to generate said signal of interest (Si); an optical splitter element (ES) designed to spatially split said signal of interest into a first signal (V1) and a second signal (V2); a first frequency-shifting optical cavity (DBDF, BDF1) comprising a first frequency shifter (AOM1) designed to shift the optical frequency of the first signal by a first frequency f.sub.1 per round trip in said first cavity, said first cavity having a first trip time .sub.1; a second frequency-shifting optical cavity (DBDF, BDF2) comprising a second frequency shifter (AOM2) designed to shift the optical frequency of the second signal by a second frequency f.sub.2 per round trip in said second cavity, said second cavity having a second trip time .sub.2; the first and the second optical cavity being designed such that a maximum number of round trips of said signal in the first and the second cavity is equal to predetermined N; a detector (PD) designed to coherently detect the first signal (W1) transmitted by the first cavity and the second signal (W2) transmitted by the second cavity and generate a photocurrent (Tr) proportional to a luminous intensity detected by said detector, a low-pass filter (LP) designed to filter frequencies of the photocurrent that are lower than min (f.sub.1/2;f.sub.22/), a processor (UT) configured to compute a square modulus of the photocurrent filtered by said low-pass filter, from which a temporal representation of frequency information of said signal of interest is determined, said frequency information being: a real part of a fractional Fourier transform of said signal of interest, an order of said fractional Fourier transform being set by the value f.sub.1.sub.1f.sub.2.sub.2 when the first cavity and the second cavity are configured to verify the condition f.sub.1.sub.1f.sub.2.sub.2, modulo 1 or a power spectrum of said signal of interest when the first cavity and the second cavity are configured to verify the condition f.sub.1.sub.1=f.sub.2.sub.2 modulo 1.

2. The device as claimed in claim 1, wherein said source comprises a monochromatic continuous-wave laser (CW), an RF source (AM) designed to generate an RF signal s(t), and a modulator (Mod) designed to amplitude-modulate or phase-modulate, using said RF signal s(t), laser radiation generated by said continuous-wave laser, so as to form said signal of interest.

3. The device as claimed in claim 1, wherein the analog low-pass filter (LP) is designed to filter frequencies of said photocurrent that are lower than min [N|f.sub.1f.sub.2|; f.sub.1/2; f.sub.2/2], when the first cavity and the second cavity are configured to verify the condition f.sub.1.sub.1=f.sub.2.sub.2 modulo 1.

4. The device as claimed in claim 1, wherein the first and the second cavity respectively comprise a first and a second amplifier (EDFA, EDFA1, EDFA2) designed to compensate for the losses induced respectively by the first and the second cavity.

5. The device as claimed in claim 1, wherein the first frequency shifter is a first acousto-optic modulator excited by a first local oscillator (OL1) designed to vary said first shift frequency and wherein the second frequency shifter is a second acousto-optic modulator excited by a second local oscillator (OL2) designed to vary said second shift frequency.

6. The device as claimed in claim 1, wherein the first cavity comprises a first controllable delay line (DL1) designed to vary the first trip time .sub.1 and wherein the second cavity comprises a second controllable delay line (DL2) designed to vary the second trip time .sub.2.

7. The device as claimed in claim 1, wherein the first and the second cavity are fiber ring cavities (BDF1, BDF2) comprising respectively a first and a second doped fiber amplifier (EDFA1, EDFA2) and a first and a second optical bandpass filter (BP1, BP2) configured to set said maximum number N of round trips in the first and the second cavity.

8. The device as claimed in claim 7, comprising stabilizing means (ST) for stabilizing the first and the second cavity that are designed to maintain, over time, the coherence of said first signal transmitted by the first cavity with said second signal transmitted by the second cavity.

9. The device as claimed in claim 6, comprising a single ring cavity (BDF), said device furthermore comprising: a first coupler (C1) designed to inject said first signal (V1) into said single cavity in a first direction, a second coupler (C2) designed to inject said second signal (V2) into said single cavity in a second direction, said first cavity corresponding to the single cavity into which the first signal is injected in the first injection direction, said second cavity corresponding to the single cavity into which the second signal is injected in the second injection direction, said single ring cavity comprising: a first circulator (CO1) designed to direct the first signal to a first controllable delay line (DL1) designed to vary the first trip time .sub.1 and comprising said first frequency shifter (AOM1), a second circulator (CO2) designed to direct the second signal to a second controllable delay line (DL2) designed to vary the second trip time .sub.1 and comprising said second frequency shifter (AOM2), a doped fiber amplifier (EDFA), an optical bandpass filter (BP) configured to set said maximum number N of round trips.

10. The device as claimed in claim 1, wherein the first and the second cavity are configured such that 1/|.sub.1.sub.240 GHz.

11. The device as claimed in claim 1, wherein the first and the second cavity are configured such that N is greater than 200.

12. A method for the spectral analysis of a signal of interest (Si) using a first frequency-shifting optical cavity (BDF1) comprising a first frequency shifter (AOM1) having a first trip time .sub.1 and a second frequency-shifting optical cavity (BDF2) comprising a second frequency shifter (AOM2) having a second trip time .sub.2, a maximum number of round trips of said signal of interest in the first and the second cavity being equal to predetermined N, said method comprising the following steps: A. generating said signal of interest; B. spatially splitting said signal of interest into a first signal (V1) and a second signal (V2); C. injecting said first signal into the first frequency-shifting optical cavity (BDF1) and shifting the optical frequency of the first signal by a first frequency f.sub.1 per round trip in said first cavity; D. injecting said second signal into the second frequency-shifting optical cavity (BDF2) and shifting the optical frequency of the second signal by a second frequency f.sub.2 per round trip in said second cavity; E. coherently detecting the first signal transmitted by the first cavity (W1) and the second signal (W2) transmitted by the second cavity and generating a photocurrent (Tr) proportional to a detected luminous intensity, F. filtering frequencies of the photocurrent that are lower than min (f.sub.12/;f.sub.22/), G. computing a square modulus of the filtered photocurrent and determining a temporal representation of frequency information of said signal of interest, said method comprising a step prior to step A, denoted step A0, of: adjusting the first or the second cavity in order to set a difference f.sub.1.sub.1f.sub.2.sub.2 to a desired non-zero value, in order to compute, in step G, a specific order of a real part of a fractional Fourier transform of the signal of interest, or adjusting the first or the second cavity in order to cancel out a difference f.sub.1.sub.1f.sub.2.sub.2, in order to compute, in step G, a power spectrum of said signal of interest.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0057] Other features, details and advantages of the invention will become apparent upon reading the description, given with reference to the appended drawings, which are given by way of example and in which, respectively:

[0058] FIG. 1 shows a schematic view of a spectral analysis device from the prior art,

[0059] FIG. 2 shows a schematic view of a wideband spectral analysis device according to the invention,

[0060] FIG. 3 shows a schematic view of a wideband spectral analysis device according to a first embodiment of the invention,

[0061] FIG. 4 shows a schematic view of a wideband spectral analysis device according to a third embodiment of the invention,

[0062] FIG. 5 shows a schematic view of a wideband spectral analysis device according to a fourth embodiment of the invention,

[0063] FIG. 6 shows a spectral analysis method according to the invention.

[0064] References to the figures, when they are identical, correspond to the same elements.

[0065] The elements are not to scale in the figures unless indicated otherwise.

DETAILED DESCRIPTION

[0066] FIG. 2 is a schematic depiction of a device D for wideband spectral analysis of a signal of interest according to the invention. As will be explained later, by virtue of a first frequency-shifting optical cavity BDF1 and a second frequency-shifting optical cavity BDF2, the device D of the invention makes it possible to determine a temporal representation of frequency information of the signal of interest. According to certain conditions linking these cavities, this temporal representation of the frequency information is the real part of a fractional transform or a power spectrum of the signal of interest.

[0067] Wideband is understood here to mean that the bandwidth of the device is greater than or equal to 20 GHz, preferably greater than or equal to 40 GHz. The parameters controlling the bandwidth of the device will be explained later in the description.

[0068] The device D of the invention comprises a radiation source S designed to generate the signal of interest Si. An optical splitter element ES of the device is designed to spatially split the signal of interest into a first signal V1 and a second signal V2.

[0069] According to one embodiment, the optical path of the signal of interest at the output of the source S is a fiber path and the element ES is a 12 fiber coupler (also called Y coupler). As an alternative, according to another embodiment, the optical path of the signal of interest at the output of the source S is a free-space optical path and the splitter element ES is a beam splitter or a splitter cube.

[0070] The first signal V1 is injected into the first frequency-shifting optical cavity BDF1. This cavity BDF1 has a first trip timer, (round trip time in the cavity) and comprises a first frequency shifter AOM1 designed to shift the optical frequency of the first signal by a first frequency f.sub.1 per round trip in the first cavity.

[0071] Similarly, the second signal V2 is injected into the second frequency-shifting optical cavity BDF2, which has a second trip time .sub.2 and which comprises a second frequency shifter AOM2 designed to shift the optical frequency of the second signal by a second frequency f.sub.2 per round trip in the second cavity.

[0072] The first optical cavity BDF1 and the second optical cavity BDF2 may, indiscriminately and without departing from the scope of the invention, be a linear cavity or a ring cavity, and a free-space cavity or a fiber cavity.

[0073] Essentially, the first and the second optical cavity are designed such that a maximum number of round trips of the signal in the first and the second cavity is equal to a predetermined number N. Thus, calling f.sub.0 the center frequency of the signal of interest, the first and the second cavity each generate a frequency comb comprising respectively the frequencies f.sub.0+nf.sub.1 and f.sub.0+nf.sub.2, with n[1; N]. In concrete terms, the first and the second cavity produce replicas of the signal of interest, which are both time-shifted (by multiples of .sub.1 and .sub.2 respectively) and frequency-shifted (by multiples of f.sub.1 and f.sub.2 respectively).

[0074] According to one embodiment, in order to control this maximum number of round trips N, the first and the second loop respectively comprise a first and a second bandpass filter BP1, BP2 (not shown in FIG. 2 but visible in FIGS. 5 and 6), of bandwidth suitable for transmitting respectively the frequencies f.sub.0+nf.sub.1 and f.sub.0+nf.sub.2, with n[1;N].

[0075] The first frequency shifter AOM1 and the second frequency shifter AOM2 are preferably acousto-optic modulators controlled by a first local oscillator OL1 and a second local oscillator OL2 (not shown in the figures). The excitation frequency generated by the first local oscillator and the second local oscillator makes it possible to vary the first frequency f.sub.1 and the second frequency f.sub.2. As an alternative, the first and the second frequency shifter are single-sideband electro-optic modulators (or SSB MZM, for single-sideband Mach Zehnder modulator).

[0076] The device D furthermore comprises a detector PD designed to coherently detect the first signal W1 transmitted by the first cavity and the second signal W2 transmitted by the second cavity. The detector PD is typically a photodiode or any other photodetector known to those skilled in the art. Preferably, the detector PD is formed by balanced photodiodes in order to detect very small light signal variations.

[0077] This detector PD detects, in real time, a luminous intensity that corresponds to the coherent summation of all replicas of the signal of interest that are both time-shifted and frequency-shifted (the frequency combs) by the first and the second cavity, and then generates a photocurrent Tr proportional to this detected luminous intensity. For the device to work, it is essential that the first signal W1 transmitted by BDF1 is coherent with the signal W2 transmitted by BDF2.

[0078] The photocurrent Tr is then filtered by a low-pass analog filter LP designed to let through frequencies of the photocurrent that are lower than the Nyquist frequencies associated with the first signal W1 and with the second signal W2, that is to say min (f.sub.1/2; f.sub.2/2).

[0079] Finally, the device D comprises a processor UT configured to compute a square modulus of the photocurrent filtered by the filter LP, thus generating a time trace TT. By way of this trace TT, the device D makes it possible to determine a temporal representation of frequency information of the signal of interest.

[0080] In the device D of the invention, the bandwidth is equal to LS.sub.in=1/(.sub.1.sub.2). This characteristic is of great interest because, by minimizing the difference .sub.1.sub.2, it is possible to maximize the bandwidth of the device D.

[0081] The inventors have demonstrated that, depending on the value of the difference f.sub.1.sub.1f.sub.2.sub.2, the trace TT is the real part of the fractional Fourier transform or a power spectrum of the signal of interest. For the sake of clarity, the derivation of the equations explaining this result are presented at the end of this description.

[0082] Thus, according to a first variant of the invention, the first cavity and the second cavity are configured to verify the condition f.sub.1.sub.1=f.sub.2.sub.2 (modulo 1). In this first variant, the photocurrent Tr is a temporal representation of the real part of the Fourier transform of the signal of interest (=/2) and the trace TT is then a temporal representation of the power spectrum of this signal. The device D therefore makes it possible to map the spectrum over time (frequency-to-time mapping) in real time.

[0083] Preferably, in this first variant, the analog low-pass filter LP is designed to filter frequencies of the photocurrent that are lower than a value min [N|f.sub.1f.sub.2|; f.sub.1/2; f.sub.2/2)], before the computing of the square modulus of the photocurrent by the processor UT. Filtering the frequencies lower than N|f.sub.1f.sub.2| makes it possible, if this value is lower than the Nyquist frequencies f.sub.1/2; f.sub.2/2, to process the photocurrent Tr with a processor comprising slower processing electronics (typically 50 Msample/s).

[0084] In this first variant, the frequency resolution of the device is .sub.f=1/(N(.sub.1.sub.2))=LS.sub.in/N. It will thus be understood that, in order to obtain a suitable spectral resolution .sub.f while maintaining a high bandwidth LS.sub.in, it is necessary to maximize the maximum number of return trips in the first and the second cavity BDF1, BDF2. Preferably, in order to obtain a suitable spectral resolution, the first and the second cavity are configured such that a maximum number of round trips N is greater than 200, preferably greater than 500. Typically, for a bandwidth of LS.sub.in=20 GHz, and N=400, the frequency resolution of the device is .sub.f=50 MHz.

[0085] In a second variant of the invention, the first cavity and the second cavity are configured to verify the condition f.sub.1.sub.1f.sub.2.sub.2, modulo 1. This difference should be much less than 1, that is to say that: f.sub.1.sub.1=f.sub.2.sub.2+, ||<<1 (see equations at the end of the description). In this second variant, the trace TT is then a temporal representation of the real part of the fractional Fourier transform (FrFT) of the signal of interest. The order of the fractional transform represented by the device D is set by the difference f.sub.1.sub.1f.sub.2.sub.2. More specifically, when f.sub.1.sub.1=f.sub.2.sub.2 modulo 1 (first variant of the invention), the photocurrent is a temporal representation of the real part of the Fourier transform of the input signal (that is to say order of the FrFT corresponding to /2) while, when f.sub.1.sub.1f.sub.2.sub.2 has a value close to 0 (second variant), the photocurrent represents the real part of the FrFT of the signal of interest. In concrete terms, the order of the FrFT depends on the scale factor (f.sub.2f.sub.1)/(.sub.1.sub.2) (see the equations at the end of the description).

[0086] In this second variant of the invention, it is possible to define a frequency resolution by considering two signals with linear frequency modulation, that are simultaneous and that have the same chirp rate (that is to say the same slope in the time-frequency plane). These two signals therefore differ in terms of starting frequency. By choosing the correct order of the FrFT, the FrFT of these chirped signals consists of two peaks, the spacing between which is proportional to the frequency difference between these two chirped signals. The resolution of the device in the second variant of the invention then corresponds to the smallest value of the frequency difference between the two chirps able to be measured by the device. This definition makes it possible to define a resolution for this system that is a frequency.

[0087] Thus, depending on the value of the difference f.sub.1.sub.1f.sub.2.sub.2, the device D of the invention makes it possible to obtain a real-time temporal representation of two different items of frequency information of the signal of interest: the real part of the FrFT or the power spectrum.

[0088] The temporal resolution of the device (that is to say the time taken by the system to compute and produce the frequency information) is equal to 1/(f.sub.2f.sub.1). It should be noted that the frequency information is obtained in real time, that is to say that the determined frequency information is updated every period 1/(f.sub.2f.sub.1).

[0089] According to one embodiment, the first cavity comprises a first controllable delay line DL1 (not shown) designed to vary the first trip time .sub.1 and/or the second cavity comprises a first controllable delay line DL2 (not shown) designed to vary the second trip time .sub.2. It is thus possible to reduce the difference (.sub.1.sub.2) in order to maximize the bandwidth LS.sub.in of the device D. Preferably, the difference (.sub.1.sub.2) is such that the bandwidth LS.sub.in is greater than or equal to 20 GHz, preferably greater than or equal to 40 GHz. This bandwidth value is not achievable in the device described in the document de Chatellus, Hugues Guillet, Luis Romero Corts, and Jos Azaa. Optical real-time Fourier transformation with kilohertz resolutions. Optica 3.1 (2016): 1-8, where the bandwidth of the device is of the order of 10 MHz as it is set by 1/.sub.1.

[0090] One advantage of the embodiments of the invention comprising delay lines DL1, DL2 and/or frequency shifters AOM1, AOM2 that make it possible to vary f.sub.1 and f.sub.2 is that of making it possible to control the difference f.sub.1.sub.1f.sub.2.sub.2. It is thus possible to switch the operation of the device from the first variant of the invention to the second variant of the invention or vice versa. Moreover, controlling the difference f.sub.1.sub.1f.sub.2.sub.2 makes it possible to choose the order of the FrFT of the signal of interest.

[0091] It should be noted that, in the embodiments of the invention in which the first and the second cavity BDF1, BDF2 respectively comprise a first and second bandpass filter BP1, BP2, the bandwidth LS.sub.in is equal to a spectral width LS.sub.BP of the first and of the second bandpass filter when this spectral width LS.sub.BP is less than 1/(.sub.1.sub.2). In other words, the bandwidth is equal to LS.sub.in=min (1/(.sub.1.sub.2); LS.sub.BP).

[0092] Notably, the device D does not perform any truncation or sampling on the spectrum of the signal of interest that reduces the analyzed part of the spectrum. Thus, unlike some of the abovementioned devices from the prior art, the device D of the invention has a POI of 100% and the frequency information of the signal of interest represented by the trace TT is representative of the entire spectrum of the signal of interest within the bandwidth.

[0093] According to one embodiment, the first and the second cavity respectively comprise a first and a second amplifier EDFA1, EDFA2 (not shown in FIG. 1 but visible in FIG. 4) that are designed to compensate for the losses induced by the first and the second cavity.

[0094] According to a first embodiment of the invention, illustrated in FIG. 3, the source S comprises a monochromatic continuous-wave laser CW, an RF source AM designed to generate an RF signal s(t), and a modulator Mod designed to amplitude-modulate or phase-modulate, using the RF signal s(t), the laser radiation generated by the continuous-wave laser so as to form said signal of interest. In the first embodiment of the invention, the device D makes it possible to determine a temporal representation of frequency information of the RF signal s(t) carried by the continuous-wave laser. Thus, according to the second variant of the invention, the trace TT is then a temporal representation of the FrFT of the RF signal s(t) and, according to the first variant of the invention, the trace TT is then a temporal representation of the convolution product of the power spectrum of the monochromatic laser CW and the power spectrum of the RF signal s(t). When the laser CW has a coherence time long enough for the radiation to remain coherent within the cavities BDF1, BDF2 (typically of the order of around one hundred s), then the trace TT is then a temporal representation of the real part of the spectrum of the RF signal s(t). This first embodiment is of particular interest given the bandwidth LS.sub.in of the device, which may be greater than 40 GHz while at the same time having a POI of 100%, these values not being unachievable using the devices from the prior art carrying out spectral analysis of RF signals.

[0095] In a second embodiment, the source S is a coherent optical source generating a signal of interest Si in the form of light radiation s(t), for example a laser source. In the second embodiment of the invention, the device D makes it possible to determine a temporal representation of frequency information of the signal s(t).

[0096] FIG. 4 schematically illustrates a third embodiment, compatible with the first and the second embodiment, in which the entire optical path of the signal of interest from the source to the detection by the photodiode is fiber-based. The alignment of the device is thus simplified, and the device is less sensitive to impacts and vibrations. In this third embodiment, the signal of interest Si is split, by a Y coupler ES, so as to form the first signal V1 and the second signal V2, both guided in a respective optical fiber. The first and the second cavity BDF1, BDF2 are fiber ring cavities comprising, respectively, a first and a second doped fiber amplifier EDFA1, EDFA2 and a first and a second optical bandpass filter BP1, BP2. The first signal V1 is injected into the first cavity BDF1 by way of a fiber coupler C1 and the second signal V2 is injected into the second cavity BDF2 by way of a fiber coupler C2.

[0097] As explained above, the first and the second doped fiber amplifier EDFA1, EDFA2, for example EDFAs, are designed to compensate for the losses induced in the cavities BDF1, BDF2. The optical bandpass filters BP1, BP2 are configured to set said maximum number N of round trips in the first and the second cavity and to limit the noise resulting from the amplified spontaneous emission of the doped fiber amplifiers.

[0098] Advantageously, in the third embodiment, the device D comprises stabilization means ST (not shown) for stabilizing the first and the second cavity BDF1, BDF2 that are designed to maintain, over time, the coherence between the first signal transmitted by the first cavity and the second signal transmitted by the second cavity. The stabilization device makes it possible to ensure that the radiation transmitted by the two cavities BDF1, BDF2 is mutually coherent. Such a device is well known to those skilled in the art and is for example described in Coherent multi-heterodyne spectroscopy using acousto-optic frequency combs, Opt. Express 26, 13800-13809 (2018).

[0099] FIG. 5 schematically illustrates a fourth embodiment, compatible with the first and the second embodiment, which constitutes an alternative to the third embodiment. This fourth embodiment aims to overcome the problem of maintaining coherence over time between the first signal W1 and the second signal W2 transmitted by the first and the second cavity. For this purpose, the device comprises a single fiber ring cavity BDF configured to form the first and the second cavity BDF1, BDF2 in the injection direction. This configuration with a single counter-propagating ring cavity makes it possible to limit the effects of fiber length fluctuations (vibrations, thermal drifts, etc.), which may be different between the first and the second cavity of the third embodiment, and which lead to a loss of coherence as the round trips progress. It also makes it possible to reduce the number of optical components needed.

[0100] Similarly to the third embodiment, in the device D of the fourth embodiment, the entire optical path of the signal of interest from the source to the detection by the photodiode is fiber-based. The device D comprises a first fiber coupler C1 designed to inject the first signal V1 into the single cavity BDF in a first direction S.sub.1. In addition, the device D comprises a second fiber coupler C2 designed to inject the second signal V2 into the single cavity BDF in a second direction S.sub.2. The cavity BDF comprises: [0101] a first circulator CO1 designed to direct the first signal V1 to a first delay line DL1 comprising the first frequency shifter AOM1, [0102] a second circulator CO2 designed to direct the second signal to the second delay line DL2 comprising the second frequency shifter AOM2, [0103] a bidirectional doped fiber amplifier EDFA (operating in both directions), [0104] an optical bandpass filter BP configured to set said maximum number N of round trips. Thus, the first cavity BDF1 corresponds to the single cavity BDF into which the first signal V1 is injected in the first injection direction S.sub.1 and the second cavity BDF2 corresponds to the single cavity BDF into which the second signal V2 is injected in the second injection direction S.sub.2.

[0105] As an alternative, according to one variant of the fourth embodiment, instead of a single amplifier EDFA and a single optical bandpass filter BP common to the first signal V1 and to the second signal V2, the first delay line DL1 comprises a first amplifier EDFA1 and an optical bandpass BP1 and the second delay line DL1 comprises a second amplifier EDFA2 and an optical bandpass BP2.

[0106] It will be understood that the embodiments of FIGS. 3 to 5 are compatible with the first and second variant of the invention according to the value of the difference f.sub.1.sub.1f.sub.2.sub.2.

[0107] Another subject of the invention is a method for the spectral analysis of a signal of interest Si implemented by the device D of the invention. FIG. 6 schematically illustrates this method, which allows a temporal representation of frequency information (FrFT or power spectrum) of the signal of interest on the basis of the value of the difference f.sub.1.sub.1f.sub.2.sub.2.

[0108] The method of FIG. 6 comprises the following steps: [0109] A. generating the signal of interest Si; [0110] B. spatially splitting the signal of interest into the first signal V1 and the second signal V2; [0111] C. injecting the first signal into the first frequency-shifting optical cavity BDF1 and shifting the optical frequency of the first signal by a first frequency f.sub.1 per round trip in the first cavity [0112] D. injecting the second signal into the second frequency-shifting optical cavity BDF2 and shifting the optical frequency of the second signal by a second frequency f.sub.2 per round trip in the second cavity; [0113] E. coherently detecting the first signal W1 transmitted by the first cavity and the second signal W2 transmitted by the second cavity and generating the photocurrent Tr proportional to the detected luminous intensity, [0114] F. filtering frequencies of the photocurrent that are lower than min [0115] G. computing the square modulus of the filtered photocurrent, based on which a temporal representation of frequency information of said signal of interest is determined.

[0116] The method of FIG. 6 has the notable advantage of being wideband, of having a POI of 100% and good frequency resolution.

[0117] According to a first variant of the method of FIG. 6, the method comprises a step prior to step A), denoted step A0), of adjusting the first or the second cavity in order to cancel out the difference f.sub.1.sub.1f.sub.2.sub.2, so as to compute, in step G, a power spectrum of the signal of interest. This first variant is implemented by the first variant of the device of the invention.

[0118] According to a second variant of the method of FIG. 6, the method comprises a step prior to step A), denoted step A0), of adjusting the first or the second cavity in order to set the difference f.sub.1.sub.1f.sub.2.sub.2 to a desired non-zero value, in order to compute, in step G, a specific order of the real part of a fractional Fourier transform of the signal of interest. This second variant is implemented by the second variant of the device of the invention.

[0119] The equations demonstrating that, depending on the value of the difference f.sub.1.sub.1f.sub.2.sub.2, the trace TT is the real part of a fractional Fourier transform or a power spectrum of the signal of interest, are derived below.

[0120] f.sub.0 denotes the frequency of the injection laser, .sub.1, .sub.2 denote the trip times in the loops, f.sub.1, f.sub.2 denote the shift frequencies per trip and s(t) denotes the signal of interest.

[0121] The electric fields at the output of the first and the second cavity are respectively:

E.sub.1(t)=.sub.n=0.sup.Ns(tn.sub.1).sub.0e.sup.i2f.sup.0.sup.te.sup.i2nf.sup.1.sup.te.sup.if.sup.i.sup..sup.1.sup.n.sup.2(1)

E.sub.2(t)=.sub.n=0.sup.Ns(tn.sub.2).sub.0e.sup.i2f.sup.0.sup.te.sup.i2nf.sup.2.sup.te.sup.if.sup.2.sup..sup.2.sup.n.sup.2(2)

[0122] The intensity detected by PD is:

l(t)=w(t)*E.sub.1(t)E.sub.2*(t)(3)

where w(t) is a time window centered at t (linked to the response of the detection) and * is the convolution product. Omitting the term E.sub.0E.sub.0* gives

l(t)=w(t)*.sub.n,ms(tn.sub.1)s(tm.sub.2)e.sup.i2(nf.sup.1.sup.mf.sup.2.sup.)te.sup.i(f.sup.1.sup..sup.1.sup.n.sup.2.sup.f.sup.2.sup..sup.2.sup.m.sup.2.sup.)(4)

I(t)=w(tt).sub.n,ms(tn.sub.1)s(tm.sub.2)e.sup.i2(nf.sup.1.sup.mf.sup.2.sup.)te.sup.i(f.sup.1.sup..sup.1.sup.n.sup.2.sup.f.sup.2.sup..sup.2.sup.m.sup.2.sup.)dt(5)

[0123] The frequencies f.sub.1 and f.sub.2 are chosen such that f=f.sub.2f.sub.1f.sub.1/2N, f.sub.2/2N. The duration of the window w(t) is chosen to be of the order of 1/Nf. It is therefore greater than 2/f.sub.1 and 2/f.sub.1. Thus, only pairs of integers (n, m=n) have a non-zero contribution to the integral.

[0124] First Variant of the Invention:

[0125] It is assumed here that f.sub.1.sub.1=f.sub.2.sub.2 This then gives:

I(t)=w(tt).sub.ns(tn.sub.1)s(tn.sub.2)e.sup.i2nftdt(6)

[0126] Over the duration of the window, e.sup.i2nft may be assimilated to e.sup.i2nft because the duration of the window is shorter than the period of the functions e.sup.i2nft.

[0127] This thus gives:

I(t)=.sub.ne.sup.i2nftw(tt)s(tn.sub.1)s(tn.sub.2)dt(7)

I(t)=.sub.n<s(tn.sub.1)s(tn.sub.2)>e.sup.i2nft(8)

[0128] where < > represents the average measured over the duration of the window.

[0129] Shifting the time origin gives

I(t)=.sub.n<s(t)s(tn)>e.sup.i2nft(9)

where =.sub.2.sub.1. The correlation function at the time t for the delay T is defined:

C(t,T)=<s(t)s(tT)>.

That is to say I(t)=.sub.nC(t,n)e.sup.i2nft(10)

[0130] It is written that, defining (t) as the Dirac delta function:

[00004]$$\begin{gathered} \begin{array}{cc}I\left(t\right)={.Math.}_{n}C\left(t,t^{}\right)e^{-i2\frac{f}{}{t}^{}t}\left(t^{}-n\right){\mathrm{dt}}^{}= \\ \left(t,t^{}\right)e^{-i2\frac{f}{}{t}^{}t}{.Math.}_n{e}^{i2n\frac{t^{}}{}}{\mathrm{dt}}^{}={.Math.}_n\left(t,t^{}\right)e^{-\mathrm{i2}\left(\frac{f}{}t-\frac{n}{}\right)t^{}}{\mathrm{dt}}^{}& \left(11\right)\end{array} \end{gathered}$$
and finally:

[00005]$\begin{array}{cc}I\left(t\right)={.Math.}_n\tilde{C}\left(t,\frac{f}{}t-\frac{n}{}\right)& \left(12\right)\end{array}$

[0131] where {tilde over (C)}(t,f)=(t,t)e.sup.i2ftdt is the Fourier transform (FT) of .sup.(t, t) with respect to the variable t.

[0132] I(t) is therefore a repetition of the FT of the autocorrelation function at the time t projected over time, that is to say, using the Wiener-Khintchine theorem, the power spectrum of the input signal at the time t. It should be noted here that the signal actually measured by the detector PD is actually 2Re(I(t)). The power spectrum is obtained via a Hilbert transformation of the measured signal.

[0133] The period of the output signal is 1/f. This therefore gives the spectrum of the input signal projected into the time domain (frequency-to-time mapping). The proportionality coefficient is simply /f=(f.sub.2f.sub.1)/(.sub.1.sub.2).

[0134] The temporal resolution is given by the duration of the window w(t), that is to say of the order of 1/Nf. The spectral resolution is therefore:

[00006]$f=\frac{f}{}\frac{1}{Nf}=\frac{1}{N}$

[0135] The spectral width of s(t) able to be measured unambiguously is 1/.

[0136] It should be noted that this result is valid when for all values of f.sub.1.sub.1=f.sub.2 .sub.2 modulo 1 (see equation 5). Modulo 1 is understood here to mean f.sub.1.sub.1=f.sub.2 .sub.2+k, with k.

[0137] Second Variant of the Invention:

[0138] It is now assumed that: f.sub.1.sub.1=f.sub.2 .sub.2+, (||<<1) giving:

[00007]$\begin{array}{cc}I\left(t\right)=\underset{n}{.Math.}C\left(t,n\right)e^{-i2n\mathrm{ft}}{e}^{-i{n}^2}& \left(13\right)\end{array}$

[0139] Following the same process:

[00008]$\begin{array}{cc}\begin{array}{c}I\left(t\right)=\underset{n}{.Math.}C\left(t,t^{}\right)e^{-i2\frac{f}{}{t}^{}t}{e}^{-i\frac{t^{{}^2}}{{}^2}}\left(t^{}-n\right){\mathrm{dt}}^{}\\ =\left(t,t^{}\right)e^{-i2\frac{f}{}{t}^{}t}{e}^{-i\frac{t^{{}^2}}{{}^2}}\underset{n}{.Math.}{e}^{i2n\frac{t^{}}{}}{\mathrm{dt}}^{}\\ =\underset{n}{.Math.}\left(t,t^{}\right)e^{-i2\left(\frac{f}{}t-\frac{n}{}\right)t^{}}{e}^{-i\frac{t^{{}^2}}{{}^2}}{\mathrm{dt}}^{}\end{array}& \left(14\right)\end{array}$

[0140] To give a fractional FT, y=t/ is set. This therefore gives:

[00009]$\begin{array}{cc}I\left(t\right)\underset{n}{.Math.}C\left(t,y\right)e^{-i2\left(\mathrm{ft}-n\right)y}{e}^{-i{y}^2}dy& \left(15\right)\end{array}$

[0141] The expression of the fractional FT (FrFT) of the function C is recognized (apart from a quadratic phase term that would be in the form e.sup.it.sup.2). Here, the order of the FrFT is given by: cot =2. It should be noted that the definition of the parameter y here defines the order of the FrFT.

[0142] Finally, I(t) is therefore the repetition of the FrFT of the autocorrelation function at the time t projected over time, that is to say, using the convolution theorem applied to the FrFT, the square modulus of the FrFT of the input signal.

[0143] As before, the signal actually measured by the detector is actually 2Re(I(t)), which means that the square modulus of the FrFT is obtained via a Hilbert transformation of the measured signal.

[0144] The period of the output signal is 1/f. This therefore gives the FrFT of the input signal projected into the time domain (fractional frequency-to-time mapping).