cycle_shift_delay18.m

是自己自创的好程序

clear all;
close all;
tic;
for trial=1:10000
XX=randint(1,128,[0,3]);
signal=qammod(XX,4)/sqrt(2);
x=signal;
C=[x;x(128),x(1:127);x(127:128),x(1:126);x(126:128),x(1:125)];
XX1=randint(1,128,[0,3]);
signal=qammod(XX1,4)/sqrt(2);
x=signal;
deta_C=[x;x(128),x(1:127);x(127:128),x(1:126);x(126:128),x(1:125)];
if(sum(XX~=XX1)>1)
A=(C-deta_C)*(C-deta_C)';
trial_rank(trial)=rank(A);
Gc=(det(A)).^(1/4);
trial_det(trial)=Gc;
else
    trial_rank(trial)=4;
end;

end;
toc;