1.一种连续变量量子密钥分发安全防御方法,其特征在于,包括:
步骤A:散粒噪声监控,通过实时监控接收端本振光强来评估散粒噪声方差;
步骤B:密钥率计算,根据步骤A所估计的实时散粒噪声方差,计算运行密钥率,计算方法如下:
通信时,Alice制备N个分别服从均匀及瑞利分布的相角和幅值数据序列Θ={θ1,θ2,LθN},Σ={α1,α2,LαN},根据变换关系其对应的高斯随机数序列为将信号Θ,Σ分别经过线性变换,得到适合加载在幅度以及相位调制器上的电压数据信号Θ'={θ1',θ2',LθN'},Σ'={α1',α2',LαN'},其中θi'=k1θi,αi'=k2αi,当对相干态进行幅度和相位调制后,编码后的相干态的光场相角和幅度分别为Φ={φ1,φ2,LφN},Ω={β1,β2,LβN},满足关系φi=f1(k1θi),βi=f2(k2αi),根据变换关系其对应的正则分量的数值也服从高斯分布,序列为以及此时存在对应关系
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>A</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>P</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
编码后的相干态经过量子信道发送给Bob,Bob通过零差检测器随机选择测量收到相干态的一个正则分量,假设Bob检测器对不同正则分量的响应值分别为 根据检测器的线性关系,可以根据检测器输出信号推断得到Bob端收到的光场信号的对应正则分量值为以及满足关系
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>y</mi>
<mrow>
<mi>A</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>Y</mi>
<mrow>
<mi>A</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>y</mi>
<mrow>
<mi>P</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<msubsup>
<mi>Y</mi>
<mrow>
<mi>P</mi>
<mi>M</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
通过测量基比对,Alice和Bob获得一串关于光场正则分量的相关联的数据,假设Alice对比基后的数据为X={x1,x2,L xN},Bob端的数据为Y={y1,y2,L yN},考虑高斯模型,Alice和Bob交换m组数据{xi,yi}i=1,K,m用于评估通信参数,则y=tx+z,其中z服从一个均值为零,方差为σ2=(1+ηTεc+νel)N0,其中vel是归一化了的检测器电噪声,得到正态线性模型下的最大似然估计及分别为:
<mrow>
<mover>
<mi>t</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
<msup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>t</mi>
<mo>^</mo>
</mover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>,</mo>
<msub>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
<mi>A</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
它们为独立分量,对应地服从如下分布
<mrow>
<mover>
<mi>t</mi>
<mo>^</mo>
</mover>
<mo>:</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mfrac>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mfrac>
<mrow>
<mi>m</mi>
<msup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>,</mo>
<mfrac>
<mrow>
<mi>m</mi>
<msub>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
<mi>A</mi>
</msub>
</mrow>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
</mfrac>
<mo>:</mo>
<msup>
<mi>&chi;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中t,σ2,VA为参数的真实值,通过上面的参数评估方法,可以得到信道的传输透过率总的归一化噪声参数为信道过噪声参数为调制方差参数为
检测器相关参数η,vel都是预先获取的稳定值,并根据密钥分发过程中散粒噪声方差N0值,计算得到系统信道过噪声运行参数,并通过公式
R=βIAB-χBE (5)
来计算系统运行密钥率,其中VA为调制方差,β为数据协商效率,χtot=χline+χh/T,χline=1/T-1+εc,χh=[(1-η)+υel]/η,且
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>A</mi>
<mo>&PlusMinus;</mo>
<msqrt>
<mrow>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<mi>B</mi>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>T</mi>
<mo>+</mo>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>+</mo>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>&PlusMinus;</mo>
<msqrt>
<mrow>
<msup>
<mi>C</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<mi>D</mi>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A&chi;</mi>
<mi>h</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msqrt>
<mi>B</mi>
</msqrt>
<mo>+</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>t</mi>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mo>=</mo>
<msqrt>
<mi>B</mi>
</msqrt>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mi>B</mi>
</msqrt>
<msub>
<mi>&chi;</mi>
<mi>h</mi>
</msub>
</mrow>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>+</mo>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>t</mi>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
根据式(5)-(7)计算,即可得到密钥率的值;
上述两个步骤同时并行执行。