function?[Estimate_Bits]?=?Detector_ML（Rx_data?H?nTx?M）
QAMOrder=log2（M）;????%調制位數?????以下情況都建立在QPSK調制
N_sym=size（Rx_data2）;???%幀數?（每根天線的接受比特）?
Estimate_Bits_Matrix=zeros（nTxN_sym*QAMOrder）;%判決矩陣預分配內存
Estimate_Bits=zeros（1nTx*N_sym*QAMOrder）;%判決結果預分配內存
bits?=?de2bi（0:2^（nTx*QAMOrder）-1?‘left-msb‘）‘;???%十進制轉二進制??一共的bit情況（調制+天線數的結果。一共有8位，調制后變成4個符號，每個符號有2位bit）????為nTx*QAMOrder?x?2^（nTx*QAMOrder）??從上到下排列為高位到低位??注意：有一個轉置的符號
for?k_1=1:N_sym??%判決的符號個數??（判決次數）?列
????rx_data=Rx_data（:k_1）;??%取出每一列的數據??（四根天線的數據都有?4*1）
????for?k_2=1:length（bits）??????%取出對比距離的次數
????????bits_i=bits（:k_2）;????%提取出每一種情況的數值
????????antenna_data_bits=reshape（bits_i[]nTx）‘;??%？？？？？將這種數據分配成發射在天線上的情況。?矩陣為：天線數*調制數
????????tx_m