%ex1.m
%仿真Alamouti?2發(fā)1收空時(shí)編碼性能，調(diào)制方式為QPSK
clear?all
datasize=100000;????????????????????%?仿真的符號(hào)數(shù)
EbNo=0:2:20;?????????????????????????%?信噪比
M=4;????????????????????????????????%?QPSK?modulation
x=randsrc（2datasize/2[0:3]）;??????%?數(shù)據(jù)源符號(hào)
x1=pskmod（xMpi/4）;???????
h=randn（2datasize/2）+j*randn（2datasize/2）;???%Rayleigh衰落信道
h=h./sqrt（2）;?????????????????????????????
for?indx=1:length（EbNo）
????sigma1=sqrt（1/（4*10.^（EbNo（indx）/10）））;??????????????%?SISO信道高斯白噪聲標(biāo)準(zhǔn)差
????n=sigma1*（randn（2datasize/2）+j*randn（2datasize/2））;
????y=x1+n;?????????????????????????????????????????????%?通過AWGN信道
????y1=x1+n./h;?????????????????????????????????????????%?通過SISO瑞利衰落信道后的判決變量
????x2=pskdemod（yMpi/4）;
????x3=pskdem