%%仿真Alamouti?2發多收空時編碼性能，調制方式為BPSK
%%
%%
clear?all
datasize=100000;????????????????????%?仿真的符號數
EbNo=0:2:20;?????????????????????????%?信噪比
M=2;????????????????????????????????%?BPSK?modulation
x=randsrc（2datasize/2[0:1]）;??????%?數據源符號//行列不變
x1=pskmod（xMpi）;???????
h=randn（8datasize/2）+j*randn（8datasize/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（1:2:）;?????????????????????????????????????????%?通過SISO瑞利衰落信道后的判決變量
????x2=pskdemod（yMpi）;
????x3=pskdemod（y1Mpi）;
????sigma2=sqrt（1/（2*10.^（EbNo（indx）/10）））;??????????????????????%?Alamouti方案每個子信道高斯白噪聲標準差
????n=sigma2*（randn（8datasize/2）+j*randn（8datasize/2））;
????n1（1:）=（conj（h（1:））.*n（1:）+h（2:）.*conj（n（2:）））./（sum（abs（h（1:2:））.^2））;????%?2發1收Alamouti方案判決變量中的噪聲項
????n1（2:）=（conj（h（2:））.*n（1:）-h（1:）.*conj（n（2:）））./（sum（abs（h（1:2:））.^2））;
????y=x1+n1;
????x4=pskdemod（yMpi）;
????n2（1:）=（conj（h（1:））.*n（1:）+h（2:）.*conj（n（2:））+conj（h（3:））.*n（3:）+h（4:）.*conj（n（4:）））./（sum（abs（h（1:4:））.^2））;????%?2發2收Alamouti方案判決變量中的噪聲項
????n2（2:）=（conj（h（2:））.*n（1:）-h（1:）.*conj（n（2:））+conj（h（4:））.*n（3:）-h（3:）.*conj（n（4:）））./（sum（abs（h（1:4:））.^2））;????%第二時隙，
????y1=x1+n2;
????x5=pskdemod（y1Mpi）;
????n3（1:）=（conj（h（1:））.*n（1:）+h（2:）.*conj（n（2:））+conj（h（3:））.*n（3:）+h（4:）.*conj（n（4:））+?conj（h（5:））.*n（5:）+h（6:）.*conj（n（6:）））./（sum（abs（h（1:6:））.^2））;????%?2發2收Alamouti方案判決變量中的噪聲項
????n3（2:）=（conj（h（2:））.*n（1:）-h（1:）.*conj（n（2:））+conj（h（4:））.*n（3:）-h（3:）.*conj（n（4:））+?conj（h（6:））.*n（5:）-h（5:）.*conj（n（6:）））./（sum（abs（h（1:6:））.^2））;????%第二時隙，
????y1=x1+n3;
????x6=pskdemod（y1Mpi）;
????n4（1:）=（conj（h（1:））.*n（1:）+h（2:）.*conj（n（2:））+conj（h（3:））.*n（3:）+h（4:）.*conj（n（4:））+?conj（h（5:））.*n（5:）+h（6:）.*conj（n（6:））+conj（h（7:））.*n（7:）+h（8:）.*conj（n（8:）））./（sum（abs（h）.^2））;????%?2發2收Alamouti方案判決變量中的噪聲項
????n4（2:）=（conj（h（2:））.*n（1:）-h（1:）.*conj（n（2:））+conj（h（4:））.*n（3:）-h（3:）.*conj（n（4:））+?conj（h（6:））.*n（5:）-h（5:）.*conj（n（6:））+conj（h（8:））.*n（7:）-h（7:）.*conj（n（8:）））./（sum（abs（h）.^2））;????%第二時隙，
????y1=x1+n4;
????x7=pskdemod（y1Mpi）;
????[tempber1（indx）]=biterr（xx2log2（M））;
????[tempber2（indx）]=biterr（xx3log2（M））;
????[tempber3（indx）]=biterr（xx4log2（M））;
????[tempber4（indx）]=biterr（xx5log2（M））;
????[tempber5（indx）]=biterr（xx6log2（M））;
????[tempber6（indx）]=biterr（xx7log2（M））;????
????
end

figure;
semilogy（EbNober2‘-o‘EbNober3‘-<‘EbNober4‘->‘EbNober5‘-d‘EbNober6‘-^‘）
ber4
grid?on
legend（‘未采用編碼‘‘2發1收Alamouti方案‘‘2發2收Alamouti方案‘‘2發3收Alamouti方案‘‘2發4收Alamouti方案‘）
xlabel（‘信噪比EbNo（dB）‘）
ylabel（‘誤比特率（BER）‘）
title（‘二進制調制下的性能曲線（BPSK）‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-04-16?18:37??STBCv2\
?????文件????????3109??2016-04-14?17:56??STBCv2\BPSK.m
?????文件????????3134??2016-04-14?17:28??STBCv2\PSK16.m
?????文件????????3136??2016-04-14?17:32??STBCv2\PSK32.m
?????文件????????3126??2016-04-14?17:25??STBCv2\PSK8.m
?????文件????????3132??2016-04-14?17:57??STBCv2\QPSK.m
?????文件????????1992??2016-04-14?17:51??STBCv2\xPSK.m
?????文件?????????854??2016-04-14?18:02??STBCv2\xPSK2.m