concurrence calculation of qubit system
14 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I need to calcualte the concurrenec for the given density matrix. Pl somebody see my code and solve my issue.
clc
close all
function c=concurrence(rho)
rho=1/2*[0 0 0 0;0 1 -1 0;0 -1 1 0;0 0 0 0];
y=[0 -1i; 1i 0];
R=rho*kron(y,y)*conj(rho)*kron(y,y);
% Real is needed since MATLAB always adds a small imaginary part ...
e=real(sqrt(eig(R)));
e=-sort(-e);
c=max(e(1)-e(2)-e(3)-e(4),0)
end
The concurrenec is not coming. Pl help me to get that.
1 commentaire
Réponses (1)
Ibrahim Issah
le 18 Mar 2021
%% intialize %%
clc
close all
clear
%%%%%%%%%%%%%
%% density matrix%%%%
rho=1/2*[0 0 0 0;0 1 -1 0;0 -1 1 0;0 0 0 0];
%%%%%%% pauli matrix%%%%%%%%%%
y=[0 -1i; 1i 0];
%%% p and ptilda %%%%%%%
R=rho*kron(y,y)*conj(rho)*kron(y,y);
% Real is needed since MATLAB always adds a small imaginary part ...
e=real(sqrt(eig(R)));
e=-sort(-e);
c=max(e(1)-e(2)-e(3)-e(4),0);
disp(c);
%% Apparently you don't need any function for this since your density matrix is already defined. Thank you.
0 commentaires
Voir également
Catégories
En savoir plus sur Creating and Concatenating Matrices dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!