%?論文仿真原型程序仿真四種（ZFZF-SICMMSEMMSE-SIC）?Vblast接收機(jī)的檢測性能，繪制誤比特率～信噪比曲線。
%?Zhaidewei?2007.4.18?（2007-5-20修正版）
%?Email:zeal0502@yahoo.com.cn
%?qq:402731924

%?發(fā)端初始化===============================================================
%?發(fā)射天線數(shù)tx接收天線數(shù)rx發(fā)射矩陣長度L（幀長）
tx=4;rx=4;L=10000;
Modulation=‘BPSK‘;
EbN0=[0:5:20];
B=30000;Ts=1/24300;
%?建立EbN0與SNR之間的換算關(guān)系
SNR=EbN0-10*log10（Ts*B）;
%?信源A
A=randint（tx*L1）;
%?經(jīng)過BPSK調(diào)制的V-Blast發(fā)射矩陣X
X=zeros（txL）;
for?k=1:tx
X（k:）=（-1）.^（A（k:tx:end）+1）;
end

%?信道傳輸=================================================================

%?快衰落Rayleigh信道H
H=sqrt（1/2）*（randn（rxtxL）+i*randn（rxtxL））;
%?均值為0方差為1的高斯白噪聲n
n=sqrt（1/2）*（randn（rxL）+i*randn（rxL））;
%?未疊加噪聲的接收信號(hào)R
R=zeros（rxL）;
for?k=1:L
????R（:k）=H（::k）*X（:k）;
end

%?檢測

%ZF========================================================================
disp（‘berz‘）;
berz=[];
%?在不同的信噪比下計(jì)算ZF接收機(jī)誤比特率berz
??for?m=SNR
??????m
??????%?每個(gè)子信道的平均信噪比為snr的接受信號(hào)R_noised
??????snr=10^（m/10）;
??????R_noised=R+sqrt（1/snr）*n;
??????x=[];
??????a=zeros（tx*L1）;
??????%?逐時(shí)隙對(duì)接收符號(hào)矢量進(jìn)行檢測，合并得到一幀發(fā)射矩陣X的估計(jì)x
??????for?t=1:L
??????????r=R_noised（:t）;
??????????%?迫零矩陣G
??????????G=pinv（H（::t））;
??????????y=G*r;
??????????xtemp=（y>=0）-（y<0）+0;
??????????x=[xxtemp];
??????end
??????
??????%?從x求A的估計(jì)a
??????for?k=1:tx
??????????a（k:tx:end）=（x（k:tx:end）+1）/2;
??????end
??????%?比較A和a計(jì)算錯(cuò)值率temp_ber
??????[errbittemp_ber]=biterr（Aa）;
??????berz=[berztemp_ber];
??end
??figure
??semilogy（EbN0berz‘o-?b‘）
?
%?ZF-SIC（ordered）==========================================================
disp（‘ZF-SIC‘）;
berzs=[];
??for?m=SNR
??????m
??????snr=10^（m/10）;
??????R_noised=R+sqrt（1/snr）*n;
??????x=[];
??????a=zeros（tx*L1）;
??????for?t=1:L
??????????r=R_noised（:t）;
??????????HH=H（::t）;
??????????G=pinv（HH）;
??????????S=[1:tx];%?S表示一個(gè)時(shí)隙內(nèi)還未檢測的符號(hào)的序號(hào)的集合
??????????xtemp=zeros（tx1）;
??????????%?逐發(fā)射天線進(jìn)行檢測
??????????for?k=1:tx
??????????????%?G的在序號(hào)集合S的非零元素集合內(nèi)范數(shù)最小的行是wki，它是G的第ki行
??????????????[wkiki]=minnorm（GS）;
??????????????%?已經(jīng)檢測過的序號(hào)清零
??????????????S（ki）=0;
??????????????%?判決統(tǒng)計(jì)量y
??????????????y=wki*r;
??????????????%?BPSK判決
??????????????xtemp（ki）=1*（y>=0）-1*（y<0）+0;
??????????????%?SIC串行干擾抵消
??????????????r=r-xtemp（ki）*H（:kit）;
??????????????%?將已經(jīng)檢測的信號(hào)對(duì)應(yīng)的信道矩陣的列清零
??????????????HH（:ki）=zeros（rx1）;
??????????????G=pinv（HH）;
??????????end
??????????x=[xxtemp];
??????end
??????for?k=1:tx
??????????a（k:tx:end）=（x（k:tx:end）+1）/2;
??????end
??????[errbittemp_ber]=biterr（Aa）;
??????berzs=[berzstemp_ber];
??end
??hold?on
??semilogy（EbN0berzs‘o-?r‘）

?%?MMSE====================================================================
?disp（‘MMSE‘）
??berm=[];
??for?m=SNR
??????m
??????snr=10^（m/10）;
??????R_noised=R+sqrt（1/snr）*n;
??????x=[];
??????a=zeros（tx*L1）;
??????for?t=1:L
??????????r=R_noised（:t）;
??????????HH=H（::t）;
??????????xtemp=zeros（tx1）;
??????????w=inv（HH‘*HH+（1/snr）*eye（tx））*HH‘;
?????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2693??2004-10-03?22:04??MYMIMO\demodulation.m

?????文件???????1237??2004-10-03?22:04??MYMIMO\modulation.m

?????文件???????2815??2006-11-21?15:44??MYMIMO\v_blast.m

?????文件???????4573??2007-05-20?12:57??MYMIMO\bd.m

?????文件????????655??2007-05-06?09:18??MYMIMO\minnorm.m

?????文件????????699??2007-05-20?00:36??MYMIMO\說明.txt

?????目錄??????????0??2006-11-21?16:49??MYMIMO

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

????????????????12890????????????????????8