诱骗态量子密钥BB84协议的三强度和无穷强度不同光子数下密钥率对比附Matlab代码

%% 初始化
clc
clear all
close all
%  协议相关效率 
q = 0.5;
%  衰减系数（在1550nm光纤中dB/km)



%% figures

L1 = 1:1:250;
p = ['r','.r';'b','.b';'g','.p'];
% Y_1 figure
%figure
%for i = 1 : size(mu,2)
%L12 = plot(L1,Y_1_inf{1,i},p(i,1));hold on
%L11 = plot(L1,Y_1{1,i},p(i,2));hold on
%L13 = plot(L1,SF_YZL_1{1,i},p(i,3));
%xlabel('Distance (km)');
%ylabel('Y1');
% title('Single Photon response rate')单光子响应率
%end

% e_1 figure
%figure

%for i = 1 : size(mu,2)
 %   L22 = plot(L1,e_1_inf{1,i},p(i,1));hold on
%L21 = plot(L1,e_1{1,i},p(i,2));hold on
%L23 = plot(L1,SF_eXU_1{1,i},p(i,3));
%xlabel('Distance (km)');
%ylabel('e1');
%end
% title('Single Photon error rate')单光子错误率

% R figure
%figure

%for i = 1 : size(mu,2)
 %   L32 = semilogy(L1,R_inf{1,i},p(i,1));hold on
%L31 = semilogy(L1,R{1,i},p(i,2));hold on
%L33 = semilogy(L1,SF_R{1,i},p(i,3));
%xlabel('Distance (km)');
%ylabel('Key generation rate');
%end

figure
 for i = 1 : size(mu,2)
    L32 = semilogy(L1,R_inf{1,i},p(i,1));hold on
xlabel('Distance (km)');
ylabel('Key generation rate');
legend('0.2','0.8','0.5');
%title('不同平均光子数无穷多诱骗态');
end


figure
 for i = 1 : size(mu,2)
L31 = semilogy(L1,R{1,i},p(i,2));hold on
xlabel('Distance (km)');
ylabel('Key generation rate');
legend('0.2','0.8','0.5');
%title('不同平均光子数3强度多诱骗态');
end
% %%--------------------------------------------------------------------------------------------