%%?mimo?stbc?qpsk?
%該程序為進行2發1收和2發2收的的性能仿真輸出誤碼率和信噪比曲線
%用的是復正交編碼，編碼速率為1.


clc;
clear?all;
Nt=2;%發射天線數
Num_sym=500;?%要發射的符號數
EN=20;??%最大信噪比
sym=randint（1Num_sym4）;????????%產生要發射的符號0123
s_qpsk=pskmod（sym4）;????%qpsk調制
sigPower=sum（abs（s_qpsk（:））.^2）/Num_sym;??%信號功率
for?Nr=1:2???%接收天線數
????for?snr=1:EN???%信噪比1到20dB
????????for?cc=1:100???%每個信噪比下進行100次仿真取誤碼率的平均
????????????sigma=10^（snr/10）;
????????????noisePower=sigPower/sigma;?%噪聲功率
????????????channel=?sqrt（1/2）*（randn（2Nr）+j*randn（2Nr））;???%產生個瑞利信道
????????????for?i=1:2:Num_sym????????????%stbc編碼，復正交編碼
????????????????count=0;??%誤碼數
????????????????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;???%接收信號
????????????????xc1=[00];xc2=[00];
????????????????if?Nr==1
????????????????????xc1（Nr）=conj（channel（1））*r（1）+channel（2）*conj（r（2））;%1個接收天線用于檢驗的組合信號
????????????????????xc2（Nr）=conj（channel（2））*r（1）-channel（1）*conj（r（2））;
????????????????elseif?Nr==2
????????????????????for?k=1:Nr
?????????????????????????
?????????????????????????xc1（Nr）=xc1