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

%创建一个本地随机流，供随机数生成器使用%可重复性。流= RandStream (“mt19937ar”）;%创建PSK调制器和解调器系统对象pskModulator = comm.PSKModulator (．..“ModulationOrder”2 ^ modOrd,．..“PhaseOffset”0,．..“BitInput”,真正的);pskDemodulator = comm.PSKDemodulator (．..“ModulationOrder”2 ^ modOrd,．..“PhaseOffset”0,．..“BitOutput”,真正的);%为3个不同的接收者创建错误率计算系统对象zfBERCalc = comm.ErrorRate;mmseBERCalc = comm.ErrorRate;mlBERCalc = comm.ErrorRate;获得ML接收器的所有位和符号组合allBits = de2bi (0:2 ^ (modOrd * N) 1,“left-msb”)”;allTxSig =重塑(pskModulator(allBits(:))， N, 2^(modOrd*N));%预分配变量来存储BER结果以提高速度[BER_ZF, BER_MMSE, BER_ML] = deal(0 (length(EbNoVec)， 3));

建立一个可视化的误码率结果的数字无花果=图;网格在；持有在；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:长度(EbNoVec)%重置错误率计算系统对象重置(zfBERCalc);重置(mmseBERCalc);重置(mlBERCalc);%从EbNo计算每个独立传输链路的信噪比snrIndB = EbNoVec(idx) + 10*log10(modOrd);snrLinear = 10 ^ (0.1 * snrIndB);而(BER_ZF (idx, 3) < 1 e5) & & ((BER_MMSE (idx, 2) < 100) | |．..|| (BER_ML(idx, 2) < 100))%创建随机位矢量调制msg = randi(stream， [0 1]， [N*modOrd, 1]);%调整数据txSig = pskModulator(味精);具有独立链路的平坦瑞利衰落信道rayleighChan = (randn(stream, M, N) + 1i*randn(stream, M, N))/根号(2);给褪色的数据添加杂音rxSig = awgn(rayleighChan*txSig, snrIndB, 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 = psk解调器(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:end]);orderVec = orderVec(1， [1:k-1,k+1:end]);%阶算法(矩阵G的计算)是唯一不同的%与ZF-SIC接收机G = (H * H + ((n n + 1) / snrLinear) *眼(n n + 1)) \眼睛(n n + 1);[~， k] = min(diag(G));symNum = orderVec (k);decBits = psk解调器(G(k，:) * H' * r); / /estMMSE(modOrd * (symNum-1) + (1:modOrd)) = decBits;如果n < n r = r - H(:， k) * pskModulator(decBits);结束结束%毫升接收机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);结束%绘制结果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”）;