基于matlab模拟量子密钥分发密钥率仿真.zip

clc; %%clear command clear; %%clear variables close all; %%close files %%%%%%%%%%%%%GMCS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha=0.2; %db/km VA=4; %modulation variance V=VA+1; beta=0.95; %reconciliation efficiency eta=0.6; %quantum efficiency vel=0.1; %electronic noise epsilon=0.0001; %excess noise L=0:1:500; KeyRate_hom = zeros(length(L), 1);%Holevo bound using homodyne KeyRate_het = zeros(length(L), 1);%Holevo bound using heterdyne % epsilon=[0.01,0.03,0.05,0.07]; % for i = 1: length(epsilon) % figure(2); % [KeyRate_hom,KeyRate_het] = GMCScvqkd(L, alpha, VA, V, eta, vel, epsilon(i), beta); % semilogy(L, KeyRate_hom, '-b'); % hold on; % % semilogy(L, KeyRate_het, '-b'); % % hold on; % end [KeyRate_hom,KeyRate_het] = GMCScvqkd(L, alpha, VA, V, eta, vel, epsilon, beta); figure(2); semilogy(L, KeyRate_hom, '-b'); grid on; box on; hold on; figure(2); semilogy(L, KeyRate_het, '-r'); hold on; xlabel('Distance[km]'); ylabel('Security Key Rate [bits/pulse]'); title('Comparison of protocols'); %%%%%%%%%%%%%%%%%%DMCS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha=0.2; %db/km VA=0.3; %modulation variance V=VA+1; beta=0.95; %reconciliation efficiency eta=0.6; %quantum efficiency vel=0.1; %electronic noise epsilon=0.00001; %excess noise L=0:1:500; KeyRate_hom_dct = zeros(length(L), 1);%Holevo bound using homodyne KeyRate_het_dct = zeros(length(L), 1);%Holevo bound using heterdyne [KeyRate_hom_dct,KeyRate_het_dct] = DMCScvqkd(L, alpha, VA, V, eta, vel, epsilon, beta); figure(2); semilogy(L, KeyRate_hom_dct, ':b'); hold on; figure(2); semilogy(L, KeyRate_het_dct, ':r'); hold on; %%%%%%%%%%%%%%%%%%decoy%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L = 0:1:500; mu = 0.1; nu = 0.05; etab = 0.4; etad = 0.15; pd = 2*10^(-12); edet = 0.01; fEmu = 1; % edet = 0.01; % edet = [0.01,0.03,0.05,0.07]; % K = decoydvqkd(L0, 0.1, 0.05, 0.4, 0.15, 2*10^(-6), edet, 1); K = decoydvqkd(L, mu, nu, etab, etad, pd, edet, fEmu); figure(2); semilogy(L, K,'--r'); hold on; % for i = 1: length(edet) % K = decoydvqkd(L, mu, nu, etab, etad, pd, edet(i), fEmu); % figure(2); % semilogy(L, K,'--r'); % hold on; % end % for i = 1: length(edet) % K = decoydvqkd(L0, 0.1, 0.05, 0.4, 0.17, 2*10^(-6), edet(i), 1); % figure(2); % h_dvt = semilogy(L0, K); % hold on; % end %%%%%%%%%%%%%%%%%%MDICVQKD%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha=0.2; %db/km VM=80; VA=VM+1; VB=VM+1; V=VB; beta=1; %reconciliation efficiency eta=1; %quantum efficiency vel=0; %electronic noise eA=0.0001; %excess noise eB=0.0001; Lac=0:1:500; %distance between Alice and Bob Lbc=0; K=MDIcvqkd(Lac, Lbc, alpha, VA, VB, V, eta, vel, eA, eB, beta); semilogy(Lac,K,'--r'); hold on; Lac=0:1:100; %distance between Alice and Bob Lbc=0:1:100; K2=MDIcvqkd(Lac, Lbc, alpha, VA, VB, V, eta, vel, eA, eB, beta); semilogy(Lac+Lbc,K2,'--b'); hold on; %%%%%%%%%%%%%%%%%%PLOB%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L=0:1:500; t=power(10,-alpha.*L/10); %channel losses in channel AC % t=0:0.0000001:1; L=50*(-log10(t)); k=-log2(1-t); semilogy(L,k); kZero = 1*1e-15*ones(1,length(L));%Holevo bound using homodyne legend('GMCS Homodyne','GMCS Heterdyne','DMCS Homedyne','DMCS Heterdyne','DecoyDVQKD','MDI-QKD asymmetric','MDI-QKD symmetric','PLOB'); y1_y2=[KeyRate_hom';kZero]; % 是一个2×ls的矩阵，第一行为y1，第二行为y2 maxY1vsY2=max(y1_y2); % 1×ls的是一个行向量，表示y1_y2每一列的最大值，即x相同时、y1与y2的的最大值 minY1vsY2=min(y1_y2); % 1×ls的是一个行向量，表示y1_y2每一列的最小值，即x相同时、y1与y2的的最小值 yForFill=[maxY1vsY2,fliplr(minY1vsY2)]; xForFill=[L,fliplr(L)]; % fill(xForFill,yForFill,'r','FaceColor','r','FaceAlpha',0.5,'EdgeAlpha',1,'EdgeColor','r'); % 填充并设置图形格式 % fill(xForFill,yForFill,'r'); % 填充并设置图形格式