stbc_cond.m

是自己自创的好程序

clear all;
close all;
N=10000;
M=4;
tic;
sss=qammod(randint(N,1,4),4);
ss=kron(ones(1,4),sss);
s=ss*diag([1,1,exp(sqrt(-1)*34.4/180*pi),exp(sqrt(-1)*34.4/180*pi)]);
for i=1:N
    X{i}=[s(i,:);
          -conj(s(i,2)),conj(s(i,1)),-conj(s(i,4)),conj(s(i,3));
          s(i,[3,4]),s(i,[1,2]);
          -conj(s(i,4)),conj(s(i,3)),-conj(s(i,2)),conj(s(i,1))];
          con(i)=cond(X{i});
end
cnt=1;
for i=1:N
    cc(i)=cond(X{i}'*X{i});
end
toc;