minimal_eulidean0110a.asv

是自己自创的好程序

clear all;
close all;
tic;
t=0;
t1=0;
t2=0;
for trial=1:10000
s=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
%s=qammod(randint(1,4,[0,15]),16)/sqrt(10);
x1=s(1,1); x2=s(1,2); x3=s(1,3); x4=s(1,4);
A=[x1,x2;-x2',x1'];
B=[x3,x4;-x4',x3'];
S=[A,B;B,A];
s1=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
%s1=qammod(randint(1,4,[0,15]),16)/sqrt(10);
x1=s1(1,1); x2=s1(1,2); x3=s1(1,3); x4=s1(1,4);
A=[x1,x2;-x2',x1'];
B=[x3,x4;-x4',x3'];
S1=[A,B;B,A];
if(sum(s~=s1)>0);
    t=t+1;
    B=S-S1;
    A=B*B';
    deta_a=rank(A);
    %deta_a=det(A);
    tr_deta_a=trace(A);
    deta_trail(t)=deta_a;
  tr_deta_trail(t)=tr_deta_a;
end
% deta_trail(trial)=deta_a;
% tr_deta_trail(trial)=tr_deta_a;
%add the 45 degree rotation 20080108.
% s=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
% %s=qammod(randint(1,4,[0,15]),16)/sqrt(10);
% x1=s(1,1); x2=s(1,2); x3=exp(sqrt(-1)*(45)/180*pi)*s(1,3); x4=exp(sqrt(-1)*45/180*pi)*s(1,4);
% A=[x1,x2;-x2',x1'];
% B=[x3,x4;-x4',x3'];
% S=[A,B;B,A];
% s1=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
% %s1=qammod(randint(1,4,[0,15]),16)/sqrt(10);
% x1=s1(1,1); x2=s1(1,2); x3=exp(sqrt(-1)*(45)/180*pi)*s1(1,3); x4=exp(sqrt(-1)*45/180*pi)*s1(1,4);
% A=[x1,x2;-x2',x1'];
% B=[x3,x4;-x4',x3'];
% S1=[A,B;B,A];
% if(sum(s~=s1)>0);
%     t1=t1+1;
%     B=S-S1;
%     A=B*B';
%     deta_a=det(A);
%     tr_deta_a=trace(A);
%     x_deta_trail(t1)=deta_a;
%   x_tr_deta_trail(t1)=tr_deta_a;
% end

%add the 45 degree rotation 20080108.
%add the s3 and s4  rotation degree 20080809.
% s=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
% %s=qammod(randint(1,4,[0,15]),16)/sqrt(10);
% x1=s(1,1); x2=s(1,2); x3=s(1,3); x4=s(1,4);
% A=[x1,x2;-x2',x1'];
% B=[x3,x4;-x4',x3'];
% S=[A,B;B,A];
% S(:,[3,4])=S(:,[3,4])*exp(sqrt(-1)*(45)/180*pi);
% s1=(1-2*randint(1,4)+sqrt(-1)*(1-2*randint(1,4)))/sqrt(2);
% %s1=qammod(randint(1,4,[0,15]),16)/sqrt(10);
% x1=s1(1,1); x2=s1(1,2); x3=s1(1,3); x4=s1(1,4);
% A=[x1,x2;-x2',x1'];
% B=[x3,x4;-x4',x3'];
% S2=[A,B;B,A];
% S2(:,[3,4])=S2(:,[3,4])*exp(sqrt(-1)*(45)/180*pi);
% if(sum(s~=s1)>0);
%     t2=t2+1;
%     B=S-S2;
%     A=B*B';
%     deta_a=det(A);
%     tr_deta_a=trace(A);
%     xx_deta_trail(t2)=deta_a;
%     xx_tr_deta_trail(t2)=tr_deta_a;
% end

%add the s3 and s4 rotation degree 20080809.
end
toc;