%%?mimo?stbc?qpsk?
%該程序?yàn)檫M(jìn)行2發(fā)1收和2發(fā)2收的的性能仿真輸出誤碼率和信噪比曲線
%用的是復(fù)正交編碼，編碼速率為1.


clc;
clear?all;
Nt=2;%發(fā)射天線數(shù)
Num_sym=500;?%要發(fā)射的符號(hào)數(shù)
EN=20;??%最大信噪比
sym=randint（1Num_sym4）;????????%產(chǎn)生要發(fā)射的符號(hào)0123
s_qpsk=pskmod（sym4）;????%qpsk調(diào)制
sigPower=sum（abs（s_qpsk（:））.^2）/Num_sym;??%信號(hào)功率
for?Nr=1:2???%接收天線數(shù)
????for?snr=1:EN???%信噪比1到20dB
????????for?cc=1:100???%每個(gè)信噪比下進(jìn)行100次仿真取誤碼率的平均
????????????sigma=10^（snr/10）;
????????????noisePower=sigPower/sigma;?%噪聲功率
????????????channel=?sqrt（1/2）*（randn（2Nr）+j*randn（2Nr））;???%產(chǎn)生個(gè)瑞利信道
????????????for?i=1:2:Num_sym????????????%stbc編碼，復(fù)正交編碼
????????????????count=0;??%誤碼數(shù)
????????????????Stbc=[s_qpsk（i）?s_qpsk（i+1）;-conj（s_qpsk（i+1））?conj（s_qpsk（i））];
????????????????noise=sqrt（noisePower）*（randn（2Nr）+j*randn（2Nr））;??%高斯白噪聲
????????????????r=Stbc*channel+noise;???%接收信號(hào)
????????????????xc1=[00];xc2=[00];
????????????????if?Nr==1
????????????????????xc1（Nr）=conj（channel（1））*r（1）+channel（2）*conj（r（2））;%1個(gè)接收天線用于檢驗(yàn)的組合信號(hào)
????????????????????xc2（Nr）=conj（channel（2））*r（1）-channel（1）*conj（r（2））;
????????????????elseif?Nr==2
????????????????????for?k=1:Nr
?????????????????????????
?????????????????????????xc1（Nr）=xc1