%%?
%?ML-ZF-MMSE?檢測算法比較
%?bpsk調制-瑞利衰落信道-2個發送端，2個接收端

%%
clc
clear
N?=?10^4;?%?發送的符號數目
Eb_N0_dB?=?0:20;?%?信噪比范圍
nTx?=?2;
nRx?=?2;
%%
%ML
for?ii?=?1:length（Eb_N0_dB）

????%?發送端
????ip?=?rand（1N）>0.5;?%?等概率產生0和1
????s?=?2*ip-1;?%?BPSK?調制?0?->?-1;?1?->?0
????sMod?=?kron（sones（nRx1））;?%?
????sMod?=?reshape（sMod[nRxnTxN/nTx]）;?%?將矩陣轉換為[nRxnTxN/nTx?]形式?

????h?=?1/sqrt（2）*[randn（nRxnTxN/nTx）?+?j*randn（nRxnTxN/nTx）];?%?瑞利衰落信道
????n?=?1/sqrt（2）*[randn（nRxN/nTx）?+?j*randn（nRxN/nTx）];?%?0均值高斯白噪聲

????%?加入噪聲后在信道中傳輸
%?????a?=?h.*sMod;?%將兩個矩陣的第三維分別相乘
%?????b?=?sum（h.*sMod2）;%將第三維加為一列
%?????c?=?squeeze（sum（h.*sMod2））;%變換為2x10矩陣
%?????d?=?10^（-Eb_N0_dB（ii）/20）*n;%噪聲為2x10矩陣
????
????y?=?squeeze（sum（h.*sMod2））?+?10^（-Eb_N0_dB（ii）/20）*n;

????%?最大似然接收
????%?-

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????7213??2011-09-22?16:46??mimo_ml_zf_mmse_bpsk_rayleigh_channel.m