空间多路复用

N = 2;%的发送天线数量M = 2;%的接收天线EbNoVec = 2:3:8;% Eb在dB /不modOrd = 2;%星座大小= 2 ^ modOrd

%创建一个本地随机流使用随机数生成器%可重复性。流= RandStream (“mt19937ar”);%创建相移键控调制和解调系统对象pskModulator = comm.PSKModulator (…“ModulationOrder”2 ^ modOrd,…“PhaseOffset”0,…“BitInput”,真正的);pskDemodulator = comm.PSKDemodulator (…“ModulationOrder”2 ^ modOrd,…“PhaseOffset”0,…“BitOutput”,真正的);从EbNo %计算信噪比为每个独立的传输链路snrIndB = convertSNR (EbNoVec,“ebno”,“BitsPerSymbol”,modOrd);snrLinear = 10 ^ (0.1 * snrIndB);%创建误码率计算系统对象三种不同的接收器zfBERCalc = comm.ErrorRate;mmseBERCalc = comm.ErrorRate;mlBERCalc = comm.ErrorRate;%获得所有钻头和符号组合毫升接收器allBits = int2bit (0:2 ^ (modOrd * N) 1, modOrd * N);allTxSig =重塑(pskModulator (allBits (:)), N, 2 ^ (modOrd * N));%预先分配变量来存储系统结果的速度[BER_ZF, BER_MMSE BER_ML] =交易(0(长度(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 =空间多路复用的;标题(“2 x2未编码的正交相移编码系统”);集(图,“DefaultLegendAutoUpdate”,“关闭”);%循环选择EbNo点为idx = 1:长度(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))%创建随机位向量来调节味精=兰迪(流,[0,1],[N * modOrd, 1]);%调整数据txSig = pskModulator(味精);%平坦瑞利衰落信道与独立的链接rayleighChan = (randn(流,M, N) + 1我* randn(流,M, N)) /√(2);%添加噪声褪色的数据rxSig = awgn (rayleighChan * txSig snrIndB (idx), 0,流);% 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:结束]);%减少订单相应的向量orderVec = orderVec (1, 1: k - 1, k + 1:结束]);%选择下一个符号解码G = (H ' * H)、眼睛(n n + 1);%一样发票(H * H),但速度更快(~ k) = min(诊断接头(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;%的患者开始调零循环为n = 1: n H = H (:, (1: k - 1, k + 1:结束]);orderVec = orderVec (1, 1: k - 1, k + 1:结束]);%的顺序算法(矩阵G计算)是唯一的区别% ZF-SIC接收机G = (H * H + ((n n + 1) / snrLinear (idx)) *眼(n n + 1)) \眼睛(n n + 1);(~ k) = min(诊断接头(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);结束结束%毫升接收机r = rxSig;H = rayleighChan;(~ k) = min(总和(abs (repmat (r, [1、2 ^ (modOrd * N))) - H * allTxSig)。^ 2));estML = allBits (:, k);%更新系统BER_ZF (idx:) = zfBERCalc(味精、estZF);BER_MMSE (idx:) = mmseBERCalc(味精、estMMSE);BER_ML (idx:) = mlBERCalc(味精、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;结束%画线semilogy (EbNoVec BER_ZF (: 1),的r -,…EbNoVec BER_MMSE (: 1),“b -”,…EbNoVec BER_ML (: 1),“g -”);持有从;

openfig (“spatMuxResults.fig”);