clc;
clear;
tic;
data_length?=?10000;%發送數據長度
Nt?=?4;
Nr?=?4;
R_c=1;????%其為糾錯編碼的碼率，無編碼系統中為1
E_s=1;????%發送符號能量。采用歸一化星座，E_s=1
M_c=2;???%每調制符號中的比特數
SNR?=?0:2:14;
l?=?data_length/Nt;%每根天線發送數據長度

for?j=?1:length（SNR）
????SNR（j）
????sigma2=（Nr*E_s）/（10^（SNR（j）/10）*M_c*R_c）;%BSNR=（Nr*E_s）/（M_c*N_0*R_c）;
????Data?=?randint（1data_length）;%要發送的信號，其為01
????Data_mod?=?2*Data-1;%BPSK調制為0：-11:1。
????s?=?reshape（Data_modNtl）;
????MMSE_detect=[];
????ZF_detect=[];
????ML_detect=[];
????Noise_detect=[];
????for?m?=?1:?l
????????
????????H?=?（1/sqrt（2））*（randn（NrNt）+i*randn（NrNt））;%信道滿足循環復高斯分布C~CN（0，1），要求實數和虛數都獨立服從C~N（01/2）
????????N?=?sqrt（sigma2/2）*（randn（Nr1）+i*randn（Nr1））;%噪聲滿足循環復高斯分布C~CN（0，sigma2），要求實數和虛數都獨立服從C~N（0sigma2/2）
????	S?=?s（:?m?）;
????????R?=?H*S+N;
????????
????????G1?=?（（H‘*H+sigma2*eye（Nt））\H‘）*R;%MMSE
????????R1?=?（（sign（real（G1‘）））+1）/2;%解調和檢測
????????MMSE_detect?=?[MMSE_detect?R1];%MMSE檢測數據
????????
????????G2=（（H‘*H）\H‘）;%pinv（H）廣義逆跑程序時間久
????????G2=G2*R;%ZF
????????R2=（（sign（real（G2‘）））+1）/2;%解調和檢測
????????ZF_detect?=?[ZF_detect?R2];%ZF檢測數據
????????
????????G3=[0000000100100011010001010110011110001001101010111100110111101111];
????????G3=reshape（G3416）;%ML
????????G4=G3*2-1;
????????for?n=1:2^Nt
????????????ML_v（n）=（R-H*G4（:n））‘*（R-H*G4（:n））;?%求矩陣2范數
????????end
????????????[XY]=min（ML_v）;%誤差最小
????????????ML_detect=[ML_detect?G3（:Y）‘];
????????????
%??????????G4=（inv（sigma2*eye（Nt）））*R;
%???????????R3?=?（（sign（real（G4‘）））+1）/2;%解調和檢測
%?????????Noise_detect?=?[Noise_detect?R3];%MMSE檢測數據
????end
????MMSE_ber（j）?=?（length（find（Data~=MMSE_detect）））/data_length;%MMSE誤比特率
????ZF_ber（j）?=?（length（find（Data~=ZF_detect）））/data_length;%ZF誤比特率
????ML_ber（j）?=?（length（find（Data~=ML_detect）））/data_length;%ML誤比特率
%?????noise_ber（j）?=?（length（find（Data~=Noise_detect）））/data_length;%ML誤比特率
end
semilogy（SNRMMSE_ber‘b-o‘SNRZF_ber‘r-*‘SNRML_ber‘k-d‘）;%SNRnoise_ber‘g-^‘
xlabel（‘E_b/N_0（dB）‘）;
ylabel（‘BER‘）;
legend（‘MMSE‘‘ZF‘‘ML‘）;
toc;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1528??2013-02-26?21:39??OSIC_QPSK.m

?????文件???????2300??2013-08-23?09:31??BPSK_ML_ZF_MMSE.m

-----------?---------??----------?-----??----

?????????????????3828????????????????????2