Matlab仿真，基于诱骗态量子密钥BB84协议的三强度和无穷强度不同光子数下密钥率对比

function [Y_1,SF_YZL_1,e_1,SF_eXU_1,SF_R,R,Y_1_inf,e_1_inf,R_inf] = process(mu,q,alpha,e_0,e_d,Y_0,eta_bob,nu,Y_1_inf,e_1_inf,R_inf) for cnt = 1:1:250 %% 效率 eta_channel = power(10, -alpha*cnt/10); eta = eta_channel * eta_bob; % 单光子状态的效率 eta_1 = 1 - power(1-eta, 1); %% 诱骗态 % 增益 Q_mu = Y_0 + 1 - exp(-eta*mu); % 量子比特错误率 T_mu = e_0*Y_0 + e_d*(1-exp(-eta*mu)); E_mu = T_mu / Q_mu; %% 单光子率 Y_1_inf(1,cnt) = Y_0 + eta_1 - Y_0*eta_1; %% 单光子错误率 e_1_inf(1,cnt) = (e_0*Y_0 + e_d*eta_1) / Y_1_inf(1,cnt); %% 单光子密钥生成率 % 1-photon gain Q_1 = Y_1_inf(1,cnt)*poisspdf(1,mu); R_inf(1,cnt) = q * (-Q_mu*error_correction_rate(E_mu)*binary_shannon_entropy(E_mu) + Q_1*(1-binary_shannon_entropy(e_1_inf(1,cnt)))); end % %%-------------------------------------------------------------------------------------------- %% 3强度诱骗态BB84 %% 真空状态 % 增益 Q_0 = Y_0; % 量子比特错误率 T_0 = Y_0 * e_0; %% probability of i-photon(i=0,1,2) with strength mu or nu(强度为mu or nu的i光子（i=0，1，2）的概率%%) Pmu_0 = poisspdf(0, mu); %% (P_0)^mu Pmu_1 = poisspdf(1, mu); %% (P_1)^mu Pmu_2 = poisspdf(2, mu); %% (P_2)^mu Pnu_0 = poisspdf(0, nu); %% (P_0)^nu Pnu_1 = poisspdf(1, nu); %% (P_1)^nu Pnu_2 = poisspdf(2, nu); %% (P_2)^nu %% 初始化空向量 Y_1 = zeros(1,250); e_1 = zeros(1,250); R = zeros(1,250); for cnt = 1:1:250 %% 效率 eta_channel = power(10, -alpha*cnt/10); eta = eta_channel * eta_bob; %% 诱骗态 % 增益 Q_nu = overall_qubit_gain(nu, Y_0, eta); % qubit error rate T_nu = overall_qubit_error_rate(nu, Y_0, eta); %% 信号状态 % 增益 Q_mu = overall_qubit_gain(mu, Y_0, eta); % 量子比特错误率 T_mu = overall_qubit_error_rate(mu, Y_0, eta); E_mu = average_qubit_error_rate(Q_mu, T_mu); %% 单光子率 Y_1(1,cnt) = (Pmu_2*(Q_nu - Pnu_0*Y_0) - Pnu_2*(Q_mu - Pmu_0*Y_0)) / (Pnu_1*Pmu_2 - Pmu_1*Pnu_2); % Y_1_alt(1,cnt) = mu/(mu*nu - nu*nu) * (Q_nu*exp(nu) - Q_mu*exp(mu)*(nu*nu/(mu*mu)) - ((mu*mu - nu*nu)/(mu*mu))*Y_0) %% 单光子错误率 e_1(1,cnt) = (T_nu - Pnu_0*Y_0*e_0) / (Pnu_1 * Y_1(1,cnt)); % e_1_alt(1,cnt) = (average_qubit_error_rate(Q_nu, T_nu)*Q_nu*exp(nu) - e_0*Y_0)/(Y_1_alt(1,cnt)*nu) %%单光子密钥生成率 % 单光子增益 Q_1 = Y_1(1,cnt)*poisspdf(1,mu); R(1,cnt) = q * (Q_1*(1-binary_shannon_entropy(e_1(1,cnt))) - Q_mu*error_correction_rate(E_mu)*binary_shannon_entropy(E_mu)); end % %%-------------------------------------------------------------------------------------------- %%具有统计波动的3强度诱骗态 % 基本脉冲数的参数 N = 10^10; % 单光子中的参数响应率，包括统计波动 gamma = 5.3; % 概率 % 光子强度选择 Pr_mu = 0.6; Pr_nu = 0.3; Pr_o = 0.1; %% 不同诱饵态的总光子数 N_mu = Pr_mu * N; N_nu = Pr_nu * N; N_o = Pr_o * N; % 基础选择 Pr_Z = 0.5; Pr_X = 0.5; %% 特殊基础和诱饵状态的数量 NZ_mu = N_mu * Pr_Z * Pr_Z; NZ_nu = N_nu * Pr_Z * Pr_Z; NX_nu = N_nu * Pr_X * Pr_X; %%Y_0的上下限，背景率 YL_0 = Y_0 * (1 - gamma/sqrt(N_o*Y_0)); YU_0 = Y_0 * (1 + gamma/sqrt(N_o*Y_0)); %% probability of i-photon(i=0,1,2) with strength mu or nu(强度为mu or nu的i光子（i=0，1，2）的概率%%) Pmu_0 = poisspdf(0, mu); %% (P_0)^mu Pmu_1 = poisspdf(1, mu); %% (P_1)^mu Pmu_2 = poisspdf(2, mu); %% (P_2)^mu Pnu_0 = poisspdf(0, nu); %% (P_0)^nu Pnu_1 = poisspdf(1, nu); %% (P_1)^nu Pnu_2 = poisspdf(2, nu); %% (P_2)^nu %% 初始化空向量 SF_YZL_1 = zeros(1,250); SF_eXU_1 = zeros(1,250); SF_R = zeros(1,250); for cnt = 1:1:250 %% 效率 eta_channel = power(10, -alpha*cnt/10); eta = eta_channel * eta_bob; %% 诱骗态 % gain Q_nu = overall_qubit_gain(nu, Y_0, eta); % qubit error rate T_nu = overall_qubit_error_rate(nu, Y_0, eta); %% 信号状态 % 增益 Q_mu = overall_qubit_gain(mu, Y_0, eta); % 量子比特错误率 T_mu = overall_qubit_error_rate(mu, Y_0, eta); E_mu = average_qubit_error_rate(Q_mu, T_mu); DELTA_Z_nu = gamma / sqrt(NZ_nu * Q_nu); DELTA_Z_mu = gamma / sqrt(NZ_mu * Q_mu); DELTA_PRIME_X_nu = gamma / sqrt(NZ_nu * T_nu); %% 单光子率 SF_YZL_1(1,cnt) = (Pmu_2*Q_nu*(1 - DELTA_Z_nu) - Pnu_2*Q_mu*(1 - DELTA_Z_mu) + (Pnu_2*Pmu_0 - Pmu_2*Pnu_0)*YL_0) / (Pnu_1*Pmu_2 - Pmu_1*Pnu_2); %% 单光子错误率 SF_eXU_1(1,cnt) = (T_nu*(1 + DELTA_PRIME_X_nu) - Pnu_0*e_0*YL_0) / (Pnu_1 * SF_YZL_1(1,cnt)); %% 单光子密钥生成率 SF_R(1,cnt) = q * (Pmu_1*SF_YZL_1(1,cnt)*(1-binary_shannon_entropy(SF_eXU_1(1,cnt))) - Q_mu*error_correction_rate(E_mu)*binary_shannon_entropy(E_mu)); end % %%-------------------------------------------------------------------------------------------- end