%ex1.m
%仿真Alamouti?2發1收空時編碼性能，調制方式為QPSK
clear?all
datasize=100000;????????????????????%?仿真的符號數
EbNo=0:2:20;?????????????????????????%?信噪比
M=4;????????????????????????????????%?QPSK?modulation
x=randsrc（2datasize/2[0:3]）;??????%?數據源符號
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信道高斯白噪聲標準差
????n=sigma1*（randn（2datasize/2）+j*randn（2datasize/2））;
????y=x1+n;?????????????????????????????????????????????%?通過AWGN信道
????y1=x1+n./h;?????????????????????????????????????????%?通過SISO瑞利衰落信道后的判決變量
????x2=pskdemod（yMpi/4）;
????x3=pskdemod（y1Mpi/4）;
????sigma2=sqrt（1/（2*10.^（EbNo（indx）/10）））;??????????????????????%?Alamouti方案每個子信道高斯白噪聲標準差
????n=sigma2*（randn（2datasize/2）+j*randn（2datasize/2））;
????n1（1:）=（conj（h（1:））.*n（1:）+h（2:）.*conj（n（2:）））./（sum（abs（h）.^2））;????%?Alamouti方案的判決變量
????n1（2:）=（conj（h（2:））.*n（1:）-h（1:）.*conj（n（2:）））./（sum（abs（h）.^2））;
????y3=x1+n1;
????x4=pskdemod（y3Mpi/4）;
????[tempber1（indx）]=biterr（xx2log2（M））;
????[tempber2（indx）]=biterr（xx3log2（M））;
????[tempber3（indx）]=biterr（xx4log2（M））;
????
end
semilogy（EbNober1‘-k*‘EbNober2‘-ko‘EbNober3‘-kd‘）
grid?on
legend（‘AWGN信道‘‘SISO瑞利衰落信道‘‘2發1收Alamouti方案‘）
xlabel（‘信噪比EbNo（dB）‘）
ylabel（‘誤比特率（BER）‘）
title（‘2發1收Alamouti方案在瑞利衰落信道下的性能‘）