N = 2;%发射天线个数M = 2;%接收天线个数EbNoVec = 2:3:8;% Eb/No in dBmodOrd = 2;%星座大小= 2^modOrd

创建一个本地随机流，供随机数生成器使用%可重复性。流= RandStream(“mt19937ar”）;创建PSK调制器和解调器系统对象pskModulator = com . pskModulator (.．.“ModulationOrder”2 ^ modOrd,.．.“PhaseOffset”0,.．.“BitInput”,真正的);pskDemodulator = com . pskDemodulator (.．.“ModulationOrder”2 ^ modOrd,.．.“PhaseOffset”0,.．.“BitOutput”,真正的);计算每个独立传输链路的EbNo信噪比snrIndB = convertSNR(EbNoVec，“ebno”，“BitsPerSymbol”, modOrd);snrLinear = 10.^(0.1*snrIndB);为3个不同的接收器创建错误率计算系统对象zfBERCalc = com . errorrate;mmseBERCalc = com . errorrate;mlBERCalc = com . errorrate;获取ML接收器的所有位和符号组合allBits = int2bit(0:2^(modOrd*N)-1, modOrd*N);allTxSig =重塑(pskModulator(allBits(:))， N, 2^(modOrd*N));预分配变量存储误码率结果以提高速度[BER_ZF, BER_MMSE, BER_ML] = deal(0 (length(EbNoVec)， 3));

设置一个数字以可视化BER结果图;网格在；持有在；ax = fig.CurrentAxes;斧子。YScale =“日志”；xlim ([EbNoVec (1) -0.01 EbNoVec(结束)]);ylim ([1 e - 3]);包含(“Eb /不(dB)”）;ylabel (“方方面面”）;fig.NumberTitle =“关闭”；fig.Renderer =“zbuffer”；fig.Name =空间多路复用的；标题(“2x2非编码QPSK系统”）;集(图,“DefaultLegendAutoUpdate”，“关闭”）;%循环选定的EbNo点为idx = 1:length(EbNoVec)%重置错误率计算系统对象重置(zfBERCalc);重置(mmseBERCalc);重置(mlBERCalc);而(BER_ZF (idx, 3) < 1 e5) & & ((BER_MMSE (idx, 2) < 100) | |.．.(BER_ZF(idx, 2) < 100) || (BER_ML(idx, 2) < 100))创建随机位向量调制msg = randi(stream， [0 1]， [N*modOrd, 1]);调制数据txSig = pskModulator(msg);%具有独立链路的平坦瑞利衰落信道rayleighChan = (randn(stream, M, N) + 1i*randn(stream, M, N))/√(2);为褪色数据添加噪声rxSig = awgn(rayleighChan*txSig, snrIndB(idx)， 0, stream);% ZF-SIC接收器r = rxSig;H = rayleighChan;假设完美的信道估计%初始化estZF = 0 (N*modOrd, 1);orderVec = 1:N;k = N+1;开始ZF清空循环为n = 1: n缩小H以删除最后一个解码符号的效果H = H(:， [1:k-1,k+1:end]);相应地收缩顺序向量orderVec = orderVec(1， [1:k-1,k+1:end]);选择下一个要解码的符号G = (H'*H) \ eye(N-n+1);%与inv(H'*H)相同，但更快[~， k] = min(diag(G));symNum = orderVec(k);硬解码所选符号decBits = pskDemodulator(G(k，:) * H' * r);estZF(modOrd * (symNum-1) + (1:modOrd)) = decBits;从r中减去最后一个解码符号的效果如果n < n r = r - H(:， k) * pskModulator(decBits);结束结束% MMSE-SIC接收器r = rxSig;H = rayleighChan;%初始化estMMSE = 0 (N*modOrd, 1);orderVec = 1:N;k = N+1;开始MMSE空化循环为n = 1: n H = H(:， [1:k-1,k+1:结束]);orderVec = orderVec(1， [1:k-1,k+1:end]);%阶算法(矩阵G计算)唯一的区别% ZF-SIC接收器G = (H * H + ((n n + 1) / snrLinear (idx)) *眼(n n + 1)) \眼睛(n n + 1);[~， k] = min(diag(G));symNum = orderVec(k);decBits = pskDemodulator(G(k，:) * H' * r);estMMSE(modOrd * (symNum-1) + (1:modOrd)) = decBits;如果n < n r = r - H(:， k) * pskModulator(decBits);结束结束% ML接收器r = rxSig;H = rayleighChan;(~ k) = min(总和(abs (repmat (r,[1、2 ^ (modOrd * N))) - H * allTxSig)。^ 2));estML = allBits(:，k);%更新误码率BER_ZF(idx，:) = zfBERCalc(msg, estZF);BER_MMSE(idx，:) = mmseBERCalc(msg, estMMSE);BER_ML(idx，:) = mlBERCalc(msg, estML);结束%绘图结果semilogy(EbNoVec(1:idx)， BER_ZF(1:idx, 1)，的r *，.．.EbNoVec(1:idx)， BER_MMSE(1:idx, 1)，“波”，.．.EbNoVec(1:idx)， BER_ML(1:idx, 1)，“gs”）;传奇(“ZF-SIC”，“MMSE-SIC”，毫升的）;drawnow;结束%画线符号学(EbNoVec, BER_ZF(:， 1)，的r -，.．.EbNoVec, BER_MMSE(:， 1)，“b -”，.．.EbNoVec, BER_ML(:， 1)，“g -”）;持有从；

openfig (“spatMuxResults.fig”）;