本发明涉及一种根据通信数据确定某个偏振态演化历史的方法，尤其是在BB84等协议中演化重构不完全时从两个可能的演化历史中确定出真实的那一个演化曲线的方法。
背景技术：
目前，公知的量子密钥分发协议中，BB84协议、E91协议及在这两者之上的诱骗态等协议均使用两组非正交的基进行编码和解码。对于偏振态这样一个二维的量子态，需要使用三组非正交的基测量才能完全层析。当使用两组非正交基测量时，层析就会出现二义性，所以在BB84等协议的量子密钥分发系统中基于层析的结果进行偏振补偿时就需要采取措施克服这种二义性[专利201510228202.6]。不过，从量子密钥分发的误码率和不确定率中重构的量子态是一种将一段时间内的偏振态近似为相等后获得的结果，从而在还原历史的偏振态和预测当前的偏振态时会出现较大的误差。为准确地再现偏振态随时间变化的曲线，可以采用演化重构的方法。但是，在按照BB84等使用两组非正交基的协议实现的量子密钥分发系统中进行演化重构时，同样会出现二义性的结果，因而无法基于演化重构进行偏振态预测和补偿。
技术实现要素：
本发明的目的是消除在按照BB84等使用两组非正交基的协议实现的量子密钥分发系统中进行偏振演化重构的二义性，确定出实际的演化历史，以便顺利实现基于演化历史进行偏振态预测和补偿。
本发明的技术方案是：
本发明提供一种演化重构不完全时的偏振态二义性消除方法，它包括以下步骤：
(a)发送端和接收端使用两组非正交的基进行量子通信，将待补偿的偏振态记为|H>；
(b)在当前补偿周期获取系统两组基测量|H>的算子平均值函数，分别记为ρ|H>1(t)、ρ|H>2(t)；
(c)从ρ|H>1(t)和ρ|H>2(t)获取出相同时段内与系统使用的两组基非正交的第三组基测量|H>的算子平均值函数，得到两个可能的结果和
(d)分别根据两组算子平均值函数和获取出相同时段内两个可能的时变偏振态|XA(t)>和|XB(t)>；
(e)记录当前补偿周期窃听检测获得的量子误码率QBERp，任选|XA(t)>或|XB(t)>为前述时段内真实的演化历史，据此进行偏振态预测和补偿，将选择的真实演化历史记为|X(t)>，未选择的时变偏振态记为|X′(t)>；
(f)在下一个补偿周期，记前一次补偿后窃听检测获得的量子误码率为QBERn，并按步骤(b)～(d)得到新的两个可能的时变偏振态|YA(t)>和|YB(t)>；
(g)如果QBERn大于QBERp，则将|X′(t)>确定为真实的演化历史，并改记为|X(t)>；否则，|X(t)>不做改动；
(h)比较|YA(t)>和|YB(t)>与更新后的|X(t)>在重合时间段的函数值的重合度，如果|YA(t)>与|X(t)>重合度高，则将|YA(t)>确定为新的演化历史，并改记为|X(t)>，将|YB(t)>改记为|X′(t)>；如果|YB(t)>与|X(t)>重合度高，则将|YB(t)>确定为新的演化历史，并改记为|X(t)>，将|YA(t)>改记为|X′(t)>；
(i)将QBERn改记为QBERp，并根据|X(t)>进行当前补偿周期内的偏振态预测和补偿；
(j)转到步骤(f)继续执行。
进一步地，步骤(b)中，如果ρ|H>1(t)或ρ|H>2(t)获取失败，则转到步骤(a)继续执行。
进一步地，步骤(h)中，比较|YA(t)>和|YB(t)>与|X(t)>在重合时间段[t1，t2]的函数值的重合度方法为：计算积分得到数值σA，计算积分得到数值σB，如果σA大于σB，则判定|YA(t)>与|X(t)>重合度高，否则判定|YB(t)>与|X(t)>重合度高。
进一步地，步骤(h)中，比较|YA(t)>和|YB(t)>与|X(t)>在重合时间段[t1，t2]的函数值的重合度方法为：在[t1，t2]内任选N个时刻τ1…τN，计算求和得到数值计算求和得到数值如果大于判定|YA(t)>与|X(t)>重合度高，否则判定|YB(t)>与|X(t)>重合度高。
进一步地，N≥100。
本发明的有益效果：
本发明中，当量子密钥分发系统按照BB84等使用两组非正交基的协议实现时，本发明的方法可以结合偏振补偿的步骤消除从量子通信的测量结果中重构出的偏振态演化二义性，为实时的偏振态预测和补偿提供确切的历史数据。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
具体实施方式
下面将描述本发明的优选实施方式。虽然实施例中显示了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
一个按照BB84协议实现的量子密钥分发系统，该系统使用两组非正交基{|H>，|V>}、{|+>，|->}进行编解码，在依据演化重构结果进行偏振态预测和补偿的过程中采用演化重构不完全时的偏振态二义性消除方法，它包括以下步骤：
(a)发送端和接收端进行量子通信，将重构的偏振态记为|H>；
(b)在当前补偿周期获取|H><H|和|+><+|两个算子测量|H>的平均值函数，记为ρ|H>1(t)、ρ|H>2(t)；如果ρ|H>1(t)或ρ|H>2(t)获取失败，则转到步骤(a)继续执行；
(c)根据ρ|H>1(t)和ρ|H>2(t)获取出圆偏振基{|R>，|L>}所含的算子|R><R|测量|H>的平均值函数，得到两个可能的结果和
(d)分别根据两组算子平均值函数和获取出相同时段内两个可能的时变偏振态|XA(t)>和|XB(t)>；
(e)记录当前补偿周期窃听检测获得的量子误码率QBERp，任选|XA(t)>或|XB(t)>为前述时段内真实的演化历史，据此进行偏振态预测和补偿，将选择的真实演化历史记为|X(t)>，未选择的时变偏振态记为|X′(t)>；
(f)在下一个补偿周期，记前一次补偿后窃听检测获得的量子误码率为QBERn，并按步骤(b)～(d)得到新的两个可能的时变偏振态|YA(t)>和|YB(t)>；
(g)如果QBERn大于QBERp，则将|X′(t)>确定为真实的演化历史，并改记为|X(t)>；否则，|X(t)>不做改动；
(h)记|YA(t)>和|YB(t)>与|X(t)>重合时间段为[-T，0]，在这个时间区间上每隔T/100取1个点，共100个点{-T+T/100，…，ti，…，0}，计算两个求和和分别得到数值和数值如果大于则将|YA(t)>确定为新的演化历史，并改记为|X(t)>，将|YB(t)>改记为|X′(t)>，否则将|YB(t)>确定为新的演化历史，并改记为|X(t)>，将|YA(t)>改记为|X′(t)>；
(i)将QBERn改记为QBERp，并根据|X(t)>进行当前补偿周期内的偏振态预测和补偿。
(j)转到步骤(f)继续执行。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。