clear?all;
datasize=100000;
EbN0=0:2:20;
M=4;
x=randsrc（2datasize/2[0:3]）;
x1=pskmod（xMpi/4）;
h=randn（4datasize/2）+j*randn（4datasize/2）;
h=h./sqrt（2）;
for?ii=1:length（EbN0）
????sigma2=sqrt（1/（2*10.^（EbN0（ii）/10）））;
????n=sigma2*（randn（4datasize/2）+j*randn（4datasize/2））;
????F=conj（h）*inv（h*conj（h））;
????beta=sqrt（1/trace（F*conj（F）））;
????W=1/beta*conj（h）*inv（（h*conj（h）））;
????y2=1/beta（h*w*x+n）