CN112491546A - Method for extracting random bits shared by atmospheric turbulence optical channels insensitive to transmission power fluctuation - Google Patents

Method for extracting random bits shared by atmospheric turbulence optical channels insensitive to transmission power fluctuation Download PDF

Info

Publication number
CN112491546A
CN112491546A CN202011376004.1A CN202011376004A CN112491546A CN 112491546 A CN112491546 A CN 112491546A CN 202011376004 A CN202011376004 A CN 202011376004A CN 112491546 A CN112491546 A CN 112491546A
Authority
CN
China
Prior art keywords
detector
sampling
sampling value
laser
sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011376004.1A
Other languages
Chinese (zh)
Other versions
CN112491546B (en
Inventor
陈纯毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun University of Science and Technology
Original Assignee
Changchun University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Priority to CN202011376004.1A priority Critical patent/CN112491546B/en
Publication of CN112491546A publication Critical patent/CN112491546A/en
Application granted granted Critical
Publication of CN112491546B publication Critical patent/CN112491546B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
    • H04L9/08—Key distribution or management, e.g. generation, sharing or updating, of cryptographic keys or passwords
    • H04L9/0861—Generation of secret information including derivation or calculation of cryptographic keys or passwords
    • H04L9/0869—Generation of secret information including derivation or calculation of cryptographic keys or passwords involving random numbers or seeds
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
    • H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
    • H04B10/0795—Performance monitoring; Measurement of transmission parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/11—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
    • H04B10/112—Line-of-sight transmission over an extended range
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50—Transmitters
    • H04B10/501—Structural aspects
    • H04B10/503—Laser transmitters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/70—Photonic quantum communication
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
    • H04L9/08—Key distribution or management, e.g. generation, sharing or updating, of cryptographic keys or passwords
    • H04L9/0816—Key establishment, i.e. cryptographic processes or cryptographic protocols whereby a shared secret becomes available to two or more parties, for subsequent use
    • H04L9/085—Secret sharing or secret splitting, e.g. threshold schemes
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
    • H04L9/08—Key distribution or management, e.g. generation, sharing or updating, of cryptographic keys or passwords
    • H04L9/0816—Key establishment, i.e. cryptographic processes or cryptographic protocols whereby a shared secret becomes available to two or more parties, for subsequent use
    • H04L9/0852—Quantum cryptography

Abstract

The invention discloses an atmospheric turbulence optical channel shared random bit extraction method insensitive to transmission power fluctuation. The invention uses the sampling value of the fluctuation of the output power of the laser to multiply the sampling value of the fluctuation of the received optical signal, thereby extracting the random bit. The method of the invention can ensure that the shared random bit can be well extracted from the atmospheric turbulence optical channel when the transmitting power of the laser transceiver A and the laser transceiver B is unstable, namely the extraction of the shared random bit is insensitive to the fluctuation of the transmitting power.

Description

Method for extracting random bits shared by atmospheric turbulence optical channels insensitive to transmission power fluctuation
Technical Field
The invention belongs to the technical field of information security, and relates to an atmospheric turbulence optical channel shared random bit extraction method insensitive to transmission power fluctuation.
Background
A chinese patent application No. 201811370939.1 discloses a method for extracting shared random bits from atmospheric turbulence optical signal fading. A key premise for successfully implementing the method is to ensure that the bidirectional optical transmission channel between the laser transceiver a and the laser transceiver B has good reciprocity, i.e., ensure that the fluctuation of the optical signals received by the laser transceiver a and the laser transceiver B has good correlation. To meet this requirement, in addition to ensuring that the bidirectional transmission channel between the laser transceiver a and the laser transceiver B is reciprocal in the manner described in "Optics Express" 2018, volume 26, page 13, 16422-16441 (i.e., CTSMC transceiving mode), the transmission power of the laser transceiver a and the laser transceiver B is actually required to be stable. In order to obtain stable transmission power, the laser transceiver a and the laser transceiver B need to perform a special control process on the lasers used by them. For example, the master's scientific paper "research and development of control system for stable laser light source output by optical fiber" filed in 2013 at the university of dennan describes a method for controlling the output power of a laser device stably.
For an atmospheric turbulence optical transmission channel, the power transmission time function can be defined as the ratio of the received optical power to the transmitted optical power at time t. Let the power transmission rate of the optical transmission channel from laser transceiver A to laser transceiver B be h at time tAB(t) at time t, let the power transmission rate of the optical transmission channel from laser transceiver B to laser transceiver A be hBA(t) of (d). It is obvious that: h is not less than 0AB(t)≤1，0≤hBA(t) is less than or equal to 1. Due to the effect of turbulence in the atmosphere, hAB(t) and hBA(t) all fluctuate with time. H if the bidirectional optical transmission channels between laser transceiver A and laser transceiver B are completely reciprocalAB(t)＝hBA(t); the fluctuation of the power transmission rate of the bidirectional atmospheric turbulence optical transmission channel can be regarded as a random source shared by the laser transceiver a and the laser transceiver B, and the laser transceiver a and the laser transceiver B can extract shared random bits therefrom to be used as a shared key required for information encryption/decryption operations.
If laser transceiver a and laser transceiver B do not have specific output power stability control over the lasers they use, their transmit power will vary over time. The transmission power of the laser transceiver a at time t is referred to as
Figure BDA0002813558640000011
Wherein
Figure BDA0002813558640000012
Represents the normalized transmit power fluctuation function (maximum value of 1), C, of the laser transceiver AAIs a constant proportionality coefficient (expression)The peak of the transmit power of laser transceiver a). Similarly, the transmission power of the laser transceiver B at time t is denoted as
Figure BDA0002813558640000013
Wherein
Figure BDA0002813558640000014
Represents the normalized transmit power fluctuation function (maximum 1), C, of laser transceiver BBIs a constant scaling factor (representing the peak of the transmitted power of laser transceiver B). Note that the optical power received by laser transceiver a at time t is equal to
Figure BDA0002813558640000021
The optical power received by the laser transceiver B at time t is equal to
Figure BDA0002813558640000022
Even h without special output power stability control of the laserAB(t)＝hBA(t), since the changes of the transmitting power of the laser transceiver a and the laser transceiver B with time are usually irrelevant, the correlation coefficient of the optical power received by the laser transceiver a and the laser transceiver B is obviously smaller than 1; in other words, the variation of the transmitting power of the laser transceiver a and the laser transceiver B with time reduces the correlation between the receiving optical powers of the two. This may further lead to an increased rate of inconsistency of the shared random bits extracted by laser transceiver a and laser transceiver B.
Note that if detected in real time in laser transceiver A, the detection is performed in real time
Figure BDA0002813558640000023
The change curve along with time is detected in real time in the laser transceiver B
Figure BDA0002813558640000024
The variation curve with time is detected, and the variation curve with time of the receiving optical power of the laser transceiver A is detected
Figure BDA0002813558640000025
Detecting the variation curve of the receiving optical power of the laser transceiver B along with the time
Figure BDA0002813558640000026
Then when h isAB(t)＝hBAAt the time of (t), the reaction mixture,
Figure BDA0002813558640000027
and
Figure BDA0002813558640000028
are fully correlated. Therefore, even if the transmission power of the laser transceiver A and the laser transceiver B varies with time, it is possible to transmit the laser beams to the laser transceiver A and the laser transceiver B
Figure BDA0002813558640000029
And
Figure BDA00028135586400000210
seen as a random source shared by laser transceiver a and laser transceiver B to support the extraction of shared random bits. For simplicity of description, the foregoing analysis ignores the detector noise effects. In practice, the presence of detector noise does not have a significant impact on the rationality of the previous analysis, as long as the signal-to-noise ratio of the detection is sufficiently high. Based on the foregoing analysis, the present invention proposes real-time detection in laser transceiver A and laser transceiver B, respectively
Figure BDA00028135586400000211
And
Figure BDA00028135586400000212
and according to
Figure BDA00028135586400000213
And
Figure BDA00028135586400000214
to extract shared random bits. The method of the invention can ensure the laser receiving and transmitting terminal A and the laserWhen the transmitting power of the optical transceiver B is unstable, shared random bits can be well extracted from the atmospheric turbulence optical channel, namely the extraction of the shared random bits is insensitive to the fluctuation of the transmitting power.
Disclosure of Invention
The invention aims to provide a method for extracting shared random bits from an atmospheric turbulence optical channel, which is insensitive to transmission power fluctuation, so that the shared random bits can be well extracted from the atmospheric turbulence optical channel under the condition that the transmission power of a laser transceiver is unstable.
The technical scheme of the method is realized as follows: the atmospheric turbulence optical channel sharing random bit extraction method insensitive to the transmission power fluctuation is characterized in that the required hardware system and the execution operation are as follows:
the laser transceiver A and the laser transceiver B are required to be in mutual sight. The laser transceiver A comprises a laser A, a fiber optic splitter A, a transceiver optical system A, a detector A0, a detector A1 and a computer A. The laser transceiver B comprises a laser B, a fiber splitter B, a transceiving optical system B, a detector B0, a detector B1 and a computer B. As shown in fig. 1, after a laser signal emitted by a laser a is split by a fiber optic splitter a, a part of the laser signal enters a light-emitting and receiving optical system a and is emitted into an atmospheric turbulence channel, and the other part of the laser signal is emitted onto a detector a 0; the laser signal which is sent by the laser A and reaches the light receiving and emitting optical system B after being transmitted through the atmospheric turbulence channel is received by the light receiving and emitting optical system B and then is incident on the detector B1; as shown in fig. 1, after a laser signal emitted by a laser B is split by an optical fiber splitter B, a part of the laser signal enters a light-emitting and receiving optical system B and is emitted into an atmospheric turbulence channel, and the other part of the laser signal is incident on a detector B0; the laser signal which is emitted by the laser B and transmitted through the atmospheric turbulence channel and reaches the light receiving and emitting optical system A is received by the light receiving and emitting optical system A and then is incident on the detector A1; the electrical signals output by the detector A0 and the detector A1 are transmitted to the computer A, and the computer A collects the electrical signal output by the detector A0 and the electrical signal output by the detector A1 in real time. The electrical signals output by the detector B0 and the detector B1 are transmitted to the computer B, and the computer B acquires the electrical signal output by the detector B0 and the electrical signal output by the detector B1 in real time.
The implementation of the method is divided into three parts, a first part, a second part and a third part.
1) The first part of the method enables the laser transceiver A and the laser transceiver B to work normally, and specifically comprises the following operations:
the laser A and the laser B are enabled to work normally, the detector A0, the detector A1, the detector B0 and the detector B1 are enabled to work normally, the computer A and the computer B are enabled to work normally, and the transmitting-receiving optical system A and the transmitting-receiving optical system B are enabled to align with each other and work normally.
2) The second part of the method implements sampling and quantization of the electrical signals output by the detector a0, the detector a1, the detector B0, and the detector B1, and specifically includes operations S1 and S2:
s1 at the slave time tbTo time teIn the time period of (1), the computer A is at intervals of δsThe electrical signals output by detector a0 and detector a1 are sampled and quantized once, respectively. Storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector A0 each time in a memory of the computer A to form a sampling value sequence LISTA0, wherein a first sampling value of the sampling value sequence LISTA0 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector A0, a second sampling value of the sampling value sequence LISTA0 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector A0, and the like; storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector A1 each time in a memory of the computer A to form a sampling value sequence LISTA1, wherein a first sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector A1, a second sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector A1, and the like;
s2 at the slave time tbTo time teIn the time period of (1), the computer B is at intervals of deltasThe electrical signals output by detector B0 and detector B1 are sampled and quantized once, respectively. Storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector B0 each time in a memory of the computer B to form a sampling value sequence LISTB0, wherein a first sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector B0, a second sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector B0, and the like; the sampled values obtained by each sampling and quantization operation performed on the electrical signal output by the detector B1 are stored in the memory of the computer B to form a sampled value sequence LISTB1, the first sampled value of the sampled value sequence LISTB1 is the sampled value obtained by first sampling and quantization of the electrical signal output by the detector B1, the second sampled value of the sampled value sequence LISTB1 is the sampled value obtained by second sampling and quantization of the electrical signal output by the detector B1, and so on.
3) The third part of the method extracts the bit sequence according to the sampling value sequence LITSTA 0, the sampling value sequence LITSTA 1, the sampling value sequence LISTB0 and the sampling value sequence LISTB1, and specifically comprises the following execution steps:
step 301: in the computer a, for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM represents the number of samples in the sample sequence LISTA1, and SA is LA1[ i ═ i]×LA0[i]/Ma0，LA1[i]Representing the ith sample in a sample sequence LISTA1, LA0[ i [ ]]Representing the i-th sample, M, in the sample sequence LISTA0a0Representing the maximum sampling value in the sampling value sequence LITSTA 0, and updating the value of the ith sampling value in the sampling value sequence LITSTA 1 to SA;
step 302: in the computer B, for the j-th sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, and SB is LB1[ j]×LB0[j]/Mb0，LB1[j]Represents the j-th sample value in the sample value sequence LISTB1, LB0[ j [ ]]Representing the jth sample value, M, in the sample value sequence LISTB0b0Represents the maximum sampling value in the sampling value sequence LISTB0, and changes the value of the jth sampling value in the sampling value sequence LISTB1New is SB;
step 303: for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM, if the ith sample in the sample sequence LISTA1 is greater than the decision threshold TaIf not, the ith bit extracted from the sampling value sequence LITSTA 1 by the computer A is 0;
step 304: for the jth sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, if the jth sample value in the sample value sequence LISTB1 is greater than the decision threshold TbIf the bit number j extracted by the computer B from the sample value sequence LISTB1 is 1, otherwise the bit number j extracted by the computer B from the sample value sequence LISTB1 is 0;
step 305: and correcting inconsistent bits in the bit sequence containing NUM bits extracted by the computer A and the bit sequence containing NUM bits extracted by the computer B by utilizing error code estimation, key agreement and error check technologies in the processing after quantum key distribution, so that the computer A and the computer B have the same bit sequence.
The invention has the following positive effects: the present invention multiplies the sampled value of the fluctuation of the output power of the laser by the sampled value of the fluctuation of the received optical signal to obtain the shared random source as described in the background section
Figure BDA0002813558640000041
And
Figure BDA0002813558640000042
the method of the invention can ensure that the shared random bit can be well extracted when the transmitting power of the laser transceiver A and the laser transceiver B is unstable. In other words, the present invention can make the extraction of shared random bits insensitive to the transmit power fluctuations of the laser transceiver.
Drawings
Fig. 1 is a schematic diagram of a system hardware architecture for extracting shared random bits from an atmospheric turbulence optical transmission channel.
Detailed Description
In order that the features and advantages of the method may be more clearly understood, the method is further described below in connection with specific embodiments. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. In this embodiment, the laser transceiver a and the laser transceiver B are respectively located on the roofs of two high-rise buildings, the computer a of the laser transceiver a and the computer B of the laser transceiver B are both connected to the internet, and the computer a and the computer B can communicate with each other through the internet, so as to perform error code estimation, key agreement, and error check operations through an internet channel. The article published in the journal of cryptography 2015, volume 2, and pages 113-121, introduces error code estimation, key agreement, and error checking in detail. The transceiver optical system A and the transceiver optical system B use the method described in Optics Express, 26, paper pages 13, 16422-16441 in 2018 to ensure that the bidirectional transmission channel between the laser transceiver A and the laser transceiver B is reciprocal. Detector a0, detector a1, detector B0, and detector B1 are PIN photodetectors. In fig. 1, a laser a and an optical fiber splitter a are connected by a single mode fiber, the optical fiber splitter a and a transceiver optical system a are connected by a single mode fiber, the optical fiber splitter a and a detector a0 are connected by a single mode fiber, the transceiver optical system a and a detector a1 are connected by a single mode fiber, the laser B and the optical fiber splitter B are connected by a single mode fiber, the optical fiber splitter B and the transceiver optical system B are connected by a single mode fiber, the optical fiber splitter B and a detector B0 are connected by a single mode fiber, and the transceiver optical system B and a detector B1 are connected by a single mode fiber.
The technical scheme of the method is realized as follows: the atmospheric turbulence optical channel sharing random bit extraction method insensitive to the transmission power fluctuation is characterized in that the required hardware system and the execution operation are as follows:
the laser transceiver A and the laser transceiver B are required to be in mutual sight. The laser transceiver A comprises a laser A, a fiber optic splitter A, a transceiver optical system A, a detector A0, a detector A1 and a computer A. The laser transceiver B comprises a laser B, a fiber splitter B, a transceiving optical system B, a detector B0, a detector B1 and a computer B. As shown in fig. 1, after a laser signal emitted by a laser a is split by a fiber optic splitter a, a part of the laser signal enters a light-emitting and receiving optical system a and is emitted into an atmospheric turbulence channel, and the other part of the laser signal is emitted onto a detector a 0; the laser signal which is sent by the laser A and reaches the light receiving and emitting optical system B after being transmitted through the atmospheric turbulence channel is received by the light receiving and emitting optical system B and then is incident on the detector B1; as shown in fig. 1, after a laser signal emitted by a laser B is split by an optical fiber splitter B, a part of the laser signal enters a light-emitting and receiving optical system B and is emitted into an atmospheric turbulence channel, and the other part of the laser signal is incident on a detector B0; the laser signal which is emitted by the laser B and transmitted through the atmospheric turbulence channel and reaches the light receiving and emitting optical system A is received by the light receiving and emitting optical system A and then is incident on the detector A1; the electrical signals output by the detector A0 and the detector A1 are transmitted to the computer A, and the computer A collects the electrical signal output by the detector A0 and the electrical signal output by the detector A1 in real time. The electrical signals output by the detector B0 and the detector B1 are transmitted to the computer B, and the computer B acquires the electrical signal output by the detector B0 and the electrical signal output by the detector B1 in real time.
The implementation of the method is divided into three parts, a first part, a second part and a third part.
1) The first part of the method enables the laser transceiver A and the laser transceiver B to work normally, and specifically comprises the following operations:
the laser A and the laser B are enabled to work normally, the detector A0, the detector A1, the detector B0 and the detector B1 are enabled to work normally, the computer A and the computer B are enabled to work normally, and the transmitting-receiving optical system A and the transmitting-receiving optical system B are enabled to align with each other and work normally.
2) The second part of the method implements sampling and quantization of the electrical signals output by the detector a0, the detector a1, the detector B0, and the detector B1, and specifically includes operations S1 and S2:
s1 at the slave time tbTo time teIn the time period of (1), the computer A is at intervals of δsThe electrical signals output by detector a0 and detector a1 are sampled and quantized once, respectively. Outputs of the detector A0 in each pairThe sampling value obtained by the sampling and quantization operation of the electric signal is stored in a memory of the computer A to form a sampling value sequence LISTA0, the first sampling value of the sampling value sequence LISTA0 is the sampling value obtained by carrying out the first sampling and quantization on the electric signal output by the detector A0, the second sampling value of the sampling value sequence LISTA0 is the sampling value obtained by carrying out the second sampling and quantization on the electric signal output by the detector A0, and the like; storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector A1 each time in a memory of the computer A to form a sampling value sequence LISTA1, wherein a first sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector A1, a second sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector A1, and the like;
s2 at the slave time tbTo time teIn the time period of (1), the computer B is at intervals of deltasThe electrical signals output by detector B0 and detector B1 are sampled and quantized once, respectively. Storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector B0 each time in a memory of the computer B to form a sampling value sequence LISTB0, wherein a first sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector B0, a second sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector B0, and the like; the sampled values obtained by each sampling and quantization operation performed on the electrical signal output by the detector B1 are stored in the memory of the computer B to form a sampled value sequence LISTB1, the first sampled value of the sampled value sequence LISTB1 is the sampled value obtained by first sampling and quantization of the electrical signal output by the detector B1, the second sampled value of the sampled value sequence LISTB1 is the sampled value obtained by second sampling and quantization of the electrical signal output by the detector B1, and so on.
3) The third part of the method extracts the bit sequence according to the sampling value sequence LITSTA 0, the sampling value sequence LITSTA 1, the sampling value sequence LISTB0 and the sampling value sequence LISTB1, and specifically comprises the following execution steps:
step 301: in the computer a, for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM represents the number of samples in the sample sequence LISTA1, and SA is LA1[ i ═ i]×LA0[i]/Ma0，LA1[i]Representing the ith sample in a sample sequence LISTA1, LA0[ i [ ]]Representing the i-th sample, M, in the sample sequence LISTA0a0Representing the maximum sampling value in the sampling value sequence LITSTA 0, and updating the value of the ith sampling value in the sampling value sequence LITSTA 1 to SA;
step 302: in the computer B, for the j-th sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, and SB is LB1[ j]×LB0[j]/Mb0，LB1[j]Represents the j-th sample value in the sample value sequence LISTB1, LB0[ j [ ]]Representing the jth sample value, M, in the sample value sequence LISTB0b0Representing the maximum sampling value in the sampling value sequence LISTB0, and updating the value of the jth sampling value in the sampling value sequence LISTB1 to SB;
step 303: for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM, if the ith sample in the sample sequence LISTA1 is greater than the decision threshold TaIf not, the ith bit extracted from the sampling value sequence LITSTA 1 by the computer A is 0;
step 304: for the jth sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, if the jth sample value in the sample value sequence LISTB1 is greater than the decision threshold TbIf the bit number j extracted by the computer B from the sample value sequence LISTB1 is 1, otherwise the bit number j extracted by the computer B from the sample value sequence LISTB1 is 0;
step 305: and correcting inconsistent bits in the bit sequence containing NUM bits extracted by the computer A and the bit sequence containing NUM bits extracted by the computer B by utilizing error code estimation, key agreement and error check technologies in the processing after quantum key distribution, so that the computer A and the computer B have the same bit sequence.
In bookIn the examples, δs1 ms, tb0 ms, te5000 milliseconds, computer a and computer B are simultaneously from time tbThe sampling is started. Decision threshold TaAnd a decision threshold TbCan be determined as follows: taking each sample value of the sample value sequence LITSTA 1 after Step301 as a random observed value, calculating an empirical cumulative distribution function F of the sample values of the sample value sequence LITSTA 1 in the computer Aa(x) And calculate
Figure BDA0002813558640000071
Taking each sample value of the sample value sequence LISTB1 after Step302 as a random observed value, the empirical cumulative distribution function F of the sample values of the sample value sequence LISTB1 is calculated in the computer Bb(x) And calculate
Figure BDA0002813558640000072
The bit sequence including NUM bits extracted by the computer a in Step305 is extracted in Step303, the ith bit of the bit sequence including NUM bits extracted by the computer a in Step305 is the ith bit extracted by the computer a from the sampled value sequence LISTA1 in Step303, and i is 1,2,3, …, NUM. The bit sequence including NUM bits extracted by the computer B in Step305 is extracted in Step304, the jth bit of the bit sequence including NUM bits extracted by the computer B in Step305 is the jth bit extracted by the computer B from the sampled value sequence LISTB1 in Step304, and j is 1,2,3, …, NUM.
The sample value sequence LISTA1, the sample value sequence LISTB1, the sample value sequence LISTA0, and the sample value sequence LISTB0 all contain NUM sample values. Fiber optic splitter a transmits 3/4 power of the laser signal emitted by laser a to the receiving and emitting optical system a, and transmits the remaining power to detector a 0. The fiber optic splitter B transmits 3/4 power of the laser signal emitted by laser B to the receiving and emitting optical system B and the remaining power to the detector B0.

Claims (1)

1. An atmospheric turbulence optical channel shared random bit extraction method insensitive to transmit power fluctuations, characterized by the hardware system required and the operations performed as follows:
a laser transceiver A and a laser transceiver B are required, and the laser transceiver A and the laser transceiver B are in mutual sight; the laser transceiver A comprises a laser A, an optical fiber splitter A, a transceiving optical system A, a detector A0, a detector A1 and a computer A; the laser transceiver B comprises a laser B, a fiber splitter B, a transceiving optical system B, a detector B0, a detector B1 and a computer B; after a laser signal emitted by the laser A is split by the optical fiber splitter A, one part of the laser signal enters the light-emitting and receiving optical system A and is emitted to an atmospheric turbulence channel, and the other part of the laser signal is emitted to the detector A0; the laser signal which is sent by the laser A and reaches the light receiving and emitting optical system B after being transmitted through the atmospheric turbulence channel is received by the light receiving and emitting optical system B and then is incident on the detector B1; after the laser signal emitted by the laser B is split by the optical fiber splitter B, one part of the laser signal enters the light-emitting and light-emitting optical system B and is emitted into an atmospheric turbulence channel, and the other part of the laser signal is emitted onto a detector B0; the laser signal which is emitted by the laser B and transmitted through the atmospheric turbulence channel and reaches the light receiving and emitting optical system A is received by the light receiving and emitting optical system A and then is incident on the detector A1; the electric signals output by the detector A0 and the detector A1 are transmitted to the computer A, and the computer A collects the electric signal output by the detector A0 and the electric signal output by the detector A1 in real time; the electric signals output by the detector B0 and the detector B1 are transmitted to the computer B, and the computer B acquires the electric signal output by the detector B0 and the electric signal output by the detector B1 in real time;
the implementation of the method is divided into three parts, namely a first part, a second part and a third part;
1) the first part of the method enables the laser transceiver A and the laser transceiver B to work normally, and specifically comprises the following operations:
enabling the laser A and the laser B to normally work, enabling the detector A0, the detector A1, the detector B0 and the detector B1 to normally work, enabling the computer A and the computer B to normally work, and enabling the light receiving and emitting optical system A and the light receiving and emitting optical system B to be aligned with each other and to normally work;
2) the second part of the method implements sampling and quantization of the electrical signals output by the detector a0, the detector a1, the detector B0, and the detector B1, and specifically includes operations S1 and S2:
s1 at the slave time tbTo time teIn the time period of (1), the computer A is at intervals of δsThe electrical signals output by the detector a0 and the detector a1 are sampled and quantized once; storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector A0 each time in a memory of the computer A to form a sampling value sequence LISTA0, wherein a first sampling value of the sampling value sequence LISTA0 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector A0, a second sampling value of the sampling value sequence LISTA0 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector A0, and the like; storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector A1 each time in a memory of the computer A to form a sampling value sequence LISTA1, wherein a first sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector A1, a second sampling value of the sampling value sequence LISTA1 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector A1, and the like;
s2 at the slave time tbTo time teIn the time period of (1), the computer B is at intervals of deltasThe electrical signals output by the detector B0 and the detector B1 are sampled and quantized once; storing sampling values obtained by performing sampling and quantization operations on the electric signal output by the detector B0 each time in a memory of the computer B to form a sampling value sequence LISTB0, wherein a first sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing first sampling and quantization on the electric signal output by the detector B0, a second sampling value of the sampling value sequence LISTB0 is a sampling value obtained by performing second sampling and quantization on the electric signal output by the detector B0, and the like; the sampled values obtained by sampling and quantizing the electrical signal output by the detector B1 at a time are stored in a computerB, a sampling value sequence LISTB1 is formed in the memory, the first sampling value of the sampling value sequence LISTB1 is a sampling value obtained by first sampling and quantizing the electrical signal output by the detector B1, the second sampling value of the sampling value sequence LISTB1 is a sampling value obtained by second sampling and quantizing the electrical signal output by the detector B1, and so on;
3) the third part of the method extracts the bit sequence according to the sampling value sequence LITSTA 0, the sampling value sequence LITSTA 1, the sampling value sequence LISTB0 and the sampling value sequence LISTB1, and specifically comprises the following execution steps:
step 301: in the computer a, for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM represents the number of samples in the sample sequence LISTA1, and SA is LA1[ i ═ i]×LA0[i]/Ma0，LA1[i]Representing the ith sample in a sample sequence LISTA1, LA0[ i [ ]]Representing the i-th sample, M, in the sample sequence LISTA0a0Representing the maximum sampling value in the sampling value sequence LITSTA 0, and updating the value of the ith sampling value in the sampling value sequence LITSTA 1 to SA;
step 302: in the computer B, for the j-th sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, and SB is LB1[ j]×LB0[j]/Mb0，LB1[j]Represents the j-th sample value in the sample value sequence LISTB1, LB0[ j [ ]]Representing the jth sample value, M, in the sample value sequence LISTB0b0Representing the maximum sampling value in the sampling value sequence LISTB0, and updating the value of the jth sampling value in the sampling value sequence LISTB1 to SB;
step 303: for the ith sample in the sample sequence LISTA1, i is 1,2,3, …, NUM, if the ith sample in the sample sequence LISTA1 is greater than the decision threshold TaIf not, the ith bit extracted from the sampling value sequence LITSTA 1 by the computer A is 0;
step 304: for the jth sample value in the sample value sequence LISTB1, j is 1,2,3, …, NUM, if the jth sample value in the sample value sequence LISTB1 is greater than the decision threshold TbThen the jth sample value extracted from the sample value sequence LISTB1 by the computer BThe bit is 1, otherwise the jth bit extracted from the sample value sequence LISTB1 by the computer B is 0;
step 305: and correcting inconsistent bits in the bit sequence containing NUM bits extracted by the computer A and the bit sequence containing NUM bits extracted by the computer B by utilizing error code estimation, key agreement and error check technologies in the processing after quantum key distribution, so that the computer A and the computer B have the same bit sequence.
CN202011376004.1A 2020-12-03 2020-12-03 Method for extracting shared random bits of atmospheric turbulence optical channel Active CN112491546B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011376004.1A CN112491546B (en) 2020-12-03 2020-12-03 Method for extracting shared random bits of atmospheric turbulence optical channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011376004.1A CN112491546B (en) 2020-12-03 2020-12-03 Method for extracting shared random bits of atmospheric turbulence optical channel

Publications (2)

Publication Number Publication Date
CN112491546A true CN112491546A (en) 2021-03-12
CN112491546B CN112491546B (en) 2022-08-23

Family

ID=74937548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011376004.1A Active CN112491546B (en) 2020-12-03 2020-12-03 Method for extracting shared random bits of atmospheric turbulence optical channel

Country Status (1)

Country Link
CN (1) CN112491546B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114422119A (en) * 2021-12-21 2022-04-29 长春理工大学 Atmospheric optical channel key extraction method for realizing decorrelation of measurement samples based on random modulation
CN115208469A (en) * 2022-09-09 2022-10-18 长春理工大学 Optical power stable receiving method and device based on reciprocity of atmospheric channels

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100080386A1 (en) * 2008-09-30 2010-04-01 The Mitre Corporation Generating Identical Numerical Sequences Utilizing a Physical Property and Secure Communication Using Such Sequences
US20120045053A1 (en) * 2010-08-20 2012-02-23 Bing Qi Random signal generator using quantum noise
WO2016119867A1 (en) * 2015-01-29 2016-08-04 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Cryptographic key distribution based on turbulence-induced beam fluctuations
CN107992285A (en) * 2017-11-22 2018-05-04 长春理工大学 From the method for laser atmospheric turbulence transmission hot spot image zooming-out true random number
CN108880780A (en) * 2018-05-23 2018-11-23 太原理工大学 Key secure distribution system and method based on Chaotic Synchronous common signal channel feature
US20190044712A1 (en) * 2016-11-08 2019-02-07 Microsoft Technology Licensing, Llc Cryptographic key creation using optical parameters
CN109474421A (en) * 2018-11-20 2019-03-15 长春理工大学 From the atmospheric turbulance optical signal method on the downside for extracting shared random bit
CN110890917A (en) * 2019-12-10 2020-03-17 长春理工大学 Atmospheric turbulence optical channel shared random bit extraction method based on frequency spectrum domain
CN110995350A (en) * 2019-12-10 2020-04-10 长春理工大学 Transmit-receive mode regulation and control method of transmit-receive common mode coupling atmosphere optical channel based on spatial filtering

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100080386A1 (en) * 2008-09-30 2010-04-01 The Mitre Corporation Generating Identical Numerical Sequences Utilizing a Physical Property and Secure Communication Using Such Sequences
US20120045053A1 (en) * 2010-08-20 2012-02-23 Bing Qi Random signal generator using quantum noise
WO2016119867A1 (en) * 2015-01-29 2016-08-04 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Cryptographic key distribution based on turbulence-induced beam fluctuations
US20190044712A1 (en) * 2016-11-08 2019-02-07 Microsoft Technology Licensing, Llc Cryptographic key creation using optical parameters
CN107992285A (en) * 2017-11-22 2018-05-04 长春理工大学 From the method for laser atmospheric turbulence transmission hot spot image zooming-out true random number
CN108880780A (en) * 2018-05-23 2018-11-23 太原理工大学 Key secure distribution system and method based on Chaotic Synchronous common signal channel feature
CN109474421A (en) * 2018-11-20 2019-03-15 长春理工大学 From the atmospheric turbulance optical signal method on the downside for extracting shared random bit
CN110890917A (en) * 2019-12-10 2020-03-17 长春理工大学 Atmospheric turbulence optical channel shared random bit extraction method based on frequency spectrum domain
CN110995350A (en) * 2019-12-10 2020-04-10 长春理工大学 Transmit-receive mode regulation and control method of transmit-receive common mode coupling atmosphere optical channel based on spatial filtering

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
KONSTANTIN KRAVTSOV ET AL.: ""Physical layer secret key generation for fiber-optical networks"", 《OPTICS EXPRESS》 *
MARVIN D. DRAKE ET AL.: ""Optical key distribution system using atmospheric turbulence as the randomness generating function: classical optical protocol for information assurance"", 《OPTICAL ENGINEERING》 *
刘中辉 等: ""基于大气湍流传输激光散斑的真随机数提取研究"", 《红外与激光工程》 *
陈纯毅 等: ""双Johnson SB分布湍流扰动OAM光通信接收信号生成研究"", 《通信学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114422119A (en) * 2021-12-21 2022-04-29 长春理工大学 Atmospheric optical channel key extraction method for realizing decorrelation of measurement samples based on random modulation
CN114422119B (en) * 2021-12-21 2023-06-16 长春理工大学 Atmospheric optical channel key extraction method for realizing decorrelation of measurement samples based on random modulation
CN115208469A (en) * 2022-09-09 2022-10-18 长春理工大学 Optical power stable receiving method and device based on reciprocity of atmospheric channels
CN115208469B (en) * 2022-09-09 2022-11-29 长春理工大学 Optical power stable receiving method and device based on reciprocity of atmospheric channels

Also Published As

Publication number Publication date
CN112491546B (en) 2022-08-23

Similar Documents

Publication Publication Date Title
CN112491546B (en) Method for extracting shared random bits of atmospheric turbulence optical channel
CN105337730A (en) Single-photon polarization control method based on phase encoding quantum key distribution (QKD) system and single-photon polarization control device based on phase encoding quantum key distribution (QKD) system
CN109474421B (en) Method for extracting shared random bit from atmospheric turbulence optical signal fading
CN110890917B (en) Atmospheric turbulence optical channel shared random bit extraction method based on frequency spectrum domain
US11057202B2 (en) Pulsed light high-speed polarization locking method of continuous-variable quantum key distribution system
CN112039919B (en) Atmospheric turbulence optical channel shared random bit extraction method based on frequency equalization
CN112929160B (en) Plug-and-play reference system and measuring equipment independent quantum key distribution system and method
CN114422119A (en) Atmospheric optical channel key extraction method for realizing decorrelation of measurement samples based on random modulation
TW349301B (en) Methods for generating side information in the presence of time-selective fading
CN112953973B (en) Hybrid attack detection method for continuous variable quantum key distribution system
CN101403824A (en) Real-time polarization control method suitable for quantum secret communication
CN113541819B (en) Time synchronization system for quantum key distribution
CN108011672A (en) The detection device and method that compatible polarization demodulation number of photons is differentiated in laser communication
CN112929163B (en) Measuring device-independent continuous variable quantum key distribution method and system
CN110620619B (en) Quantum communication system, transmitting end thereof and quantum communication method
CN111769881B (en) Method and system for improving phase compensation precision and communication efficiency of CVQKD system
Zhu et al. Optical performance monitoring via domain adversarial adaptation in few-mode fiber
CN112491547B (en) Atmospheric turbulence optical channel shared random bit extraction system
CN111756526B (en) Quantum key distribution system, transmitting end, receiving end and communication method thereof
CN102104430B (en) 1,550nm-waveband light beam tracking and communication integrated optical detection device
CN109361456B (en) Method for measuring signal fluctuation correlation coefficient of bidirectional atmospheric turbulence optical channel
Timofeev The influence of the time of single photon transmission of information on the reliability of its reception in a quantum cryptographic communication channel
CN109429508A (en) A kind of photon detection system
CN207460200U (en) A kind of quantum key dispatching system based on Waveguide interference ring
CN112787597B (en) Noiseless linear amplification method for three-degree-of-freedom coded single photon bit

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant