%?仿真五種（MRCZFMMSE，ZF-SICMMSE-SIC）?Vblast接收的檢測性能，繪制誤比特率～接收天線數曲線。?
%?發端初始化===============================================================?
%?發射天線數tx接收天線數rx發射矩陣長度L（幀長）?
clear?all;
clc;

tx=10;
rx_number=[20:20:200];
L=10000;?
Modulation=‘BPSK‘;
EbN0=5;
B=30000;Ts=1/24300;?
SNR=EbN0-10*log10（Ts*B）;
%?SNR=5;
%?信源A?
A=randint（tx*L1）;?
%?經過BPSK調制的V-Blast發射矩陣X?
X=zeros（txL）;?
for?k=1:tx?
????X（k:）=（-1）.^（A（k:tx:end）+1）;?
end?

%?檢測

%?MRC======================================================================
disp（‘MRC‘）;
berm=[];

for?rx=rx_number
????rx
????%?快衰落Rayleigh信道H?
????H=sqrt（1/2）*（randn（rxtxL）+i*randn（rxtxL））;?
????%?均值為0方差為1的高斯白噪聲n?
????n=sqrt（1/2）*（randn（rxL）+i*randn（rxL））;
????%?未疊加噪聲的接收信號R?
????R=zeros（rxL）;?
????for?k=1:L
???????R（:k）=sqrt（1/tx）*H（::k）*X（:k）;?
????end
????
????snr=10^（SNR/10）;
????R_noised=R+sqrt（1/snr）*n;
????x=[];
????a=zeros（tx*L1）;
????%?逐時隙對接收符號矢量進行檢測，合并得到一幀發射矩陣X的估計x?
????for?t=1:L
????????r=R_noised（:t）;
????????HH=H（::t）;
????????G=HH‘;
????????y=G*r;
????????xtemp=（y>0）-（y<0）;
????????x=[xxtemp];
????end
????%?從x求A的估計a?
????for?k=1:tx
????????a（k:tx:end）=（x（k:tx:end）+1）/2;
????end
????%?比較A和a計算錯值率temp_ber
????[errbittemp_ber]=biterr（Aa）;
????berm=[bermtemp_ber];
end
figure
semilogy（rx_numberberm‘*-?g‘）;
hold?on

%?ZF=======================================================================
disp（‘ZF‘）;
berz=[];

for?rx=rx_number
????rx
????%?快衰落Rayleigh信道H?
????H=sqrt（1/2）*（randn（rxtxL）+i*randn（rxtxL））;?
????%?均值為0方差為1的高斯白噪聲n?
????n=sqrt（1/2）*（randn（rxL）+i*randn（rxL））;
????%?未疊加噪聲的接收信號R?
????R=zeros（rxL）;?
????for?k=1:L
???????R（:k）=sqrt（1/tx）*H（::k）*X（:k）;?
????end
????
????snr=10^（SNR/10）;
????R_noised=R+sqrt（1/snr）*n;
????x=[];
????a=zeros（tx*L1）;
????%?逐時隙對接收符號矢量進行檢測，合并得到一幀發射矩陣X的估計x?
????for?t=1:L
????????r=R_noised（:t）;
????????HH=H（::t）;
????????G=pinv（HH）;
????????y=G*r;
????????xtemp=（y>0）-（y<0）;
????????x=[xxtemp];
????end
????%?從x求A的估計a?
????for?k=1:tx
????????a（k:tx:end）=（x（k:tx:end）+1）/2;
????end
????%?比較A和a計算錯值率temp_ber
????[errbittemp_ber]=biterr（Aa）;
????berz=[berztemp_ber];
end
semilogy（rx_numberberz‘o-?b‘）;

%?MMSE=====================================================================
disp（‘MMSE‘）;
bermm=[];

for?rx=rx_number
????rx
????%?快衰落Rayleigh信道H?
????H=sqrt（1/2）*（randn（rxtxL）+i*randn（rxtxL））;?
????%?均值為0方差為1的高斯白噪聲n?
????n=sqrt（1/2）*（randn（rxL）+i*randn（rxL））;
????%?未疊加噪聲的接收信號R?
????R=zeros（rxL）;?
????for?k=1:L
???????R（:k）=sqrt（1/tx）*H（::k）*X（:k）;?
????end
????
????snr=10^（SNR/10）;
????R_noised=R+sqrt（1/snr）*n;
????x=[];
????a=zeros（tx*L1）;
????%?逐時隙對接收符號矢量進行檢測，合并得到一幀發射矩陣X的估計x?
????for?t=1:L
????????r=R_noised（:t）;
????????HH=H（::t）;
????????G=inv（HH‘*HH+（1/snr）*eye（tx））*HH‘;
????????y=G*r;
????????xtemp=（y>0）-（y<0）;
????????x=[xxtemp];
????end
????%?從x求A的估計a?
????for?k=1:tx
????????a（k:tx:end）=（x（k:

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

?????文件???????5951??2017-12-14?21:49??MIMO和Massive?MIMO信號檢測算法仿真程序\bpsk_large_linear_nolinear.m

?????文件???????5263??2017-12-11?19:08??MIMO和Massive?MIMO信號檢測算法仿真程序\bpsk_test_linear_nolinear.m

?????文件????????205??2017-12-10?19:06??MIMO和Massive?MIMO信號檢測算法仿真程序\minnorm.m

?????文件???????7908??2017-12-14?23:15??MIMO和Massive?MIMO信號檢測算法仿真程序\qam_large_linear_nolinear1.m

?????文件???????7189??2017-12-14?19:45??MIMO和Massive?MIMO信號檢測算法仿真程序\qam_test_linear_nolinear1.m

?????文件???????8193??2017-12-14?22:58??MIMO和Massive?MIMO信號檢測算法仿真程序\qpsk_larage_linear_nolinear.m

?????文件???????7468??2017-12-14?20:01??MIMO和Massive?MIMO信號檢測算法仿真程序\qpsk_test_linear_nolinear.m

?????文件?????600700??2019-04-19?02:21??MIMO和Massive?MIMO信號檢測算法仿真程序\畢設仿真.docx

?????目錄??????????0??2019-04-24?10:53??MIMO和Massive?MIMO信號檢測算法仿真程序

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

???????????????642877????????????????????9