MIMO的ML检测matlab程序

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