《MIMO系统检测算法仿真设计.docx》由会员分享，可在线阅读，更多相关《MIMO系统检测算法仿真设计.docx（19页珍藏版）》请在三一文库上搜索。
1、下载可编辑MIMO系统检测算法仿真姓名：学号：目录.专业 .整理 .下载可编辑摘要 - -3一、 MIMO系统简介 - -4二、 MIMO系统容量分析 - -5三、 MIMO多用户系统检测技术 - -83.1 多用户检测技术原理- -83.2 V-BLAST 的迫零检测 - -93.3 ZF-SIC 算法 - -93.4 MMSE算法 - -103.5 MMSE-SIC算法 - -113.6最大似然 (ML)检测算法 -11四、仿真结果-12五、分析与总结-13六、参考文献-15七、附录 -16.专业 .整理 .下载可编辑摘要关于如何在有限的频谱资源上实现高速率和大容量是未来移动通信系统的关键
2、因素之一。 MIMO技术利用多个天线实现多发多收 , 充分开发空间资源 , 可以成倍地提高信道容量 , 同时也可以提高信道的可靠性 , 降低误码率。多用户检测技术与 MIMO结合可以进一步提高传输效率 , 提高系统性能。报告首先介绍了MIMO系统以及多用户检测理论, 然后着重解释了 BLAST系统中一类垂直分层空时码的检测算法 , 依据信号模型 , 分析推导了基于迫零准则和最小均方误差准则的估计算法 , 并在此基础上采用了以上算法与判决反馈及最佳排序思想结合的方法 , 使系统的误码率性能得到了提高。此外，还进行了最大似然算法的仿真。最后通过仿真实验比较了各种算法的性能和特点 , 结果表明分层空
3、时码用于无线通信具有极大优势。.专业 .整理 .下载可编辑一、 MIMO系统简介多输入多输出(MIIMO,Multiple-Input Multiple-Output)技术最早由 Marconi 于 1908 年提出的 , 是利用多天线来抵御信道衰落。 可以简单定义为 : 在一个任意的无线系统中 , 链路的发端和收端都使用多天线 , 也包括单入多出系统和多入单出系统。 MIMO的核心思想是 : 将发送端的信号分开而将接收端天线的信号合并 , 使每个 MIMO用户的传输质量比特误码率 (BER)或数据速率得到改进 , 提高网络服务质量。 MIMO技术可以把多径效应作为有利因素来加以利用 , 在不
4、增加带宽的情况下成倍地提高通信系统的容量和频谱利用率。 在理想情况下 ,MIMO系统的信道容量随发送天线个数和接收天线个数的最小值线性增加 , 从而提供目前其他技术无法达到的容量潜力。多天线发送和接收技术是空间分集与时间分集技术的结合, 有很好的抗干扰能力 ,如果能进一步将多天线发送和接收技术与信道编码技术结合 , 就可以更大程度地提高系统性能。MIMO技术领域另一个研究热点就是空时编码。常见的空时码有分层空时码、空时网格码、 空时分组码。空时码的主要思想是利用空间和时间上的编码实现一定的空间分集和时间分集 , 从而降低信道误码率。 MIMO无线通信技术采用空时处理技术进行信号处理 , 在多径
5、环境下 , 无线 MIMO系统可以极大地提高频谱利用率 , 增加系统的数据传输速率。MIMO技术利用多个天线实现多发多收, 充分开发空间资源 , 可以.专业 .整理 .下载可编辑成倍地提高信道容量 , 同时也可以提高信道的可靠性 , 降低误码率。多用户检测技术与 MIMO结合可以进一步提高传输效率。 但是在 MIMO系统多用户检测的研究上 , 尽管全空间最大似然 (ML)检测能获取最优的检测性能 , 但由于其算法复杂度太高, 一直无法实际应用 , 于是研究复杂度较低 , 性能接近ML 的多用户检测算法具有重要意义,MMSE算法,ZF 算法 , 以及对两者进行改进的 MMSE-SIC算法和 ZF
6、-SIC 算法。二、 MIMO系统容量分析假定一个点对点 MIMO系统有 Mr 个发射天线、 MR 个接收天线。集中于用离散时间描述的复基带线性系统模型 , 系统框图如图 2-1 所示。用 Mr×1 列矩阵 s 表示每个符号周期内的发射信号 , 其中第 i 个元素 si 表示第 i 根天线发射的信号。图 2-1 MIMO 系统框图对于高斯信道 , 按照信息论可知 , 发射信号的最佳分布也应该是.专业 .整理 .下载可编辑高斯分布的。因此 ,s 的元素是零均值独立同分布的高斯变量。发射信号的协方差矩阵为 :R ss E SSH (2-1)E ? 代表均值； S 表示矩阵的复共扼转置矩阵
7、。不管发射天线数是多少 , 总的发射功率限制为 P。如果信道在发射端未知 , 则假定从各个天线发射的信号都有相等的功率P / M T 。这样 , 发射信号的协方差矩阵为 :RssPI M r(2-2)M T由于发射信号的带宽足够窄, 因此可以认为它的频率响应是平坦的。换句话说 , 假定信道是无记忆的。用 M RM T 的复矩阵H描述信道。 hji 表示矩阵 H的第 ji个元素 ,代表从第 i 根发射天线到第j 根接收天线之间的信道衰落系数。 为了规范 , 假定 M R 根接收天线中每一根天线的接收功率都等于总的发射功率。这种假定 , 实际上忽略了信号传播过程中的信号衰减和放大 , 包括阴影、天
8、线增益。于是得到了具有确定参数的信道矩阵 H中的各元素的规范限定 , 如下式所示 :M r2(2-3)hjiM T , j 1,2,., M Ri 1当信道矩阵元素为随机变量时 , 规范就是对上述表达式取期望值。假定己知接收端信道矩阵 , 但发射端不确定 , 那么可以通过向接收端发射训练序列来估计信道矩阵 , 然后再通过可靠的反馈信道将估计的.专业 .整理 .下载可编辑信道状态信息 (CSI) 发送到发射端。 用 M R 1 的列矩阵描述接收端的噪声, 表示为 n。它的元素是统计独立的复零均值高斯变量 , 它具有独立且方差相等的实部与虚部。接收噪声的协方差矩阵为 :Rnn E nnH (2-4
9、)如果 n 的元素之间没有相关性 , 则接收噪声的协方差矩RnnN0 I M R(2-5)M R 个接收分支中每一个都有相同的噪声功率N 0 。接收端基于最大似然准则, 在 M R 根天线上进行联合操作。用M R1的列矩阵描述接收信号, 表示为r, 其中每个复元素代表一根接收天线。 Pr 表示每根接收天线输出端的平均功率。则每根接收天线处的平均信噪比 (SNR)定义为 :Pr(2-6)N 0前面己假定每根天线的总接收功率都等于总发射功率, 这样 ,SNR就等于总的发射功率和每根接收天线的噪声功率的比值, 而且它独立于MT,写为:P(2-7)N 0使用线性模型 , 可将接收矢量表示为rHsn 接
10、收信号的协方差矩阵定义为 E rr H , 利用式 (2-5) 可以得出Rrr HRssH H(2-8)而总接收信号功率可表示为tr (Rss ) 。.专业 .整理 .下载可编辑三、 MIMO多用户系统检测技术3.1 多用户检测技术原理传统的接收机进行的都是单一用户的检测, 它的缺点是在对一个用户解调时没有利用已知的其它用户的信息, 而是将多址干扰当作高斯白噪声处理 , 按照经典的直接扩频理论对每个用户信号分别进行扩频码匹配处理 , 因此抗 MAI 干扰能力较差 , 大大降低了系统容量。MUD技术是在传统的检测技术基础上, 充分利用造成干扰的所有用户的信息 , 对单个用户进行检测, 从而解决抗
11、干扰问题和抗远近效应的问题 , 降低了系统对功率控制精度的要求, 因而能够更加有效地利用上行链路的频谱资源, 从而显著提高系统的容量。MUD的主要思想是 : 把同时占用某个信道的所有用户或某些用户的信号当作有用信号而不是当作干扰信号处理, 利用多个用户的码元、时间、信号幅度以及相位等信息联合检测单个用户的信号 , 即综合利用各种信息及信号处理手段 , 对接收信号进行处理 , 从而达到对多个用户信号的最佳联合检测。图 3-1 多用户检测算法分类.专业 .整理 .下载可编辑3.2 V-BLAST 的迫零检测迫零检测器是一种线性检测器， 它的作用是通过是通过线性滤波将数据流分隔开来 , 再对其进行独
12、立译码。我们假定信道矩阵是可逆的, 则发送信号向量的估计为：x (H H H )-1 Hx H x(3-1)其中 H 表示矩阵 H的伪逆。在迫零 ZF 检测器中，每个数据流的分集增益为nRnT1。因此，迫零检测器是次优的。 迫零检测是 V-BLAST系统最早提出的经典检测算法。迫零算法的目的是找到迫零向量以消除各数据层之间的干扰。3.3 ZF-SIC算法串行干扰抵消算法 (SIC) 的思想是采用过逐步检测相消法来实现子流的区分和检测，考虑MIMO系统的信号检测模型xh1s1h2 s2.hn snn ，之前所讲的ZF和 MMSE算法是将 s 各分量的估计值都得到后再作判决， 而 SIC 算法是每
13、得到一个分量就作判决得:到 si （即根据某种准则找出星座点集合中的一点与之对应），然后得到一个新的接收信号x':'中就消除了 si 的干扰，重复上ix hi s，这时 x述步骤得到 s 所有分量的估计即可。 然而如果某个分量估计出错，则会大大影响剩余分量估计的准确度，即差错传播问题。 引入排序机制的 OSIC算法就在一定程度上削弱了这个影响。ZF 的协方差矩阵 CZF2 (H H H ) 12WZF WZFH ，其中 wi hj0, ij ，1, ij.专业 .整理 .下载可编辑wi (i1,2,., nr ) 作为行向量组成一个矩阵WZF ，这说明s 各分量估计值的方差正比
14、于WZF 相应行向量的范数Wi ，显然应该先选择方差小的分量进行估计，即最先选择 s 的第 k1arg minwk1个分量进行估计，1 k1 nr得 sw x，而后对该结果进行判决得到s 。记 x(1)x H(1)Hs(1)s，k1k1k1，经过第一次判决后， 接收信号矢量更新为 x(2)x(1):，去掉 H 的第hk1 sk1k1 列 hk1得到 H(2)， 去 掉 s 中 已 判 决 出 的 分 量 sk1得 到s(2) =(s1,.,sk1,sk1,.sn )T，计算新的检测子WZF(2)， 选 择 WZF(2)的 第11rk2 arg minwk2行与 x(2) 作内积然后判决即可得到
15、s(2) 中的第 k2 个分量1 k2 nr的估计结果。重复上述步骤直到得到s 所有分量的估计结果。3.4 MMSE 算法MMSE算法在一定程度上改善了迫零算法的性能，它是用一个新的矩阵 H MMSE 来代替迫零算法中的H。利用一下原则得到加权矩阵2H MMSE arg min ( E H MMSE r s )H MMSE其中 E 表示期望值。算法步骤如下：1) 先根据上述原则得到并计算加权矩阵 H MMSE ( H H H2I) 1HH，其中：2 为噪声方差；为发送信号功率。如果每根发射天线信号能量进行归一化，则在公式中的可以省略，此时1 。2)将加权矩阵左乘接收信号，式子变为 H MMSE
16、 r H MMSE H sH MMSE n 。3)直接利用公式 s Q( H MMSE r ) 进行量化，从而对信号进行译码。.专业 .整理 .下载可编辑3.5 MMSE-SIC 算法1) 首先进行初始化 G1 H MMSE ；2)选取信噪比最大的信号 k1arg min (Gi ) j ；j3)第 i 次迭代 Wki(Gi )ki 取第 ki 行；4)计算判决变量 skiWki ri ；5) 进行硬判决 ski Wki ri ；6)抵消这个信号对下面要进行判决的影响 ri 1riski ( H ) ki;7)计算 Gi1H ki ;8)再次选取信噪比最大的信号 ki 1 arg min (G
17、i1 ) j;j9)i i 1。3.6最大似然 (ML) 检测算法最大似然检测是 MIMO检测中的最优检测方式，能获得最高的分集增益，检测出的信号是星座图上最接近发射点的信号。对于MIMO来说，对第 k 个子载波上的信号检测表示如下：2dk (n)ML arg min N zk ( n) H k (n)dk (n)(3-2)dk ( n) ( S)式中， H k (n) 为 MIMO信道估计所估计出的信道频率响应矩阵，S 为BPSK调制中所有发出符号的星座集合，在进行了信道估计后，产生一个对角矩阵， 因此，对各个子载波上的数据 dk ( n) 进行 ML检测时，可以分离开来进行，这时，多载波下
18、的检测就相当于单个载波下的检测。.专业 .整理 .下载可编辑四、仿真结果0BPSK 调 制 的 2x2 型 MIMO 检 测 算 法 仿 真 图10ZFZF-SIC-1MMSE10MMSE-SICML10率特 10比误-2-3-410-510-6100510152025010EbN0(dB)图 4-1不同收发天线数MIMO检测算法仿真图2x2 型4x4 型8x8 型-110率特 10比误-2-310-4100510152025EbN0(dB)图 4-2.专业 .整理 .下载可编辑0不同收发天线数MIMO检测算法仿真图102x2 型2x4 型4x2 型-110-210率特比误-310-410-5
19、100510152025EbN0(dB)图 4-3五、分析与总结由仿真结果可知， 在采用相同的调制方式和相同的信噪比时，非线性检测的ML 算法的误比特率远小于其他算法。线性检测算法中MMSE算法要优于 ZF 算法，这是因为MMSE算法在处理消除符号间干扰的时候考虑了噪声的影响， 从而使得噪声和同频信号相互干扰造成的错误最小。 ZF 算法因为要根除同频信号的相互干扰而忽略了噪声的影响，因此，检测效果差些。在线性检测 ZF 和 MMSE中，当采用了顺序干扰消除， 一次一个地将干扰从接收信号中除去时，性能得到了较好的提升， 从图 4-1 中可.专业 .整理 .下载可编辑以看出， ZF-SIC 算法性
20、能和 MMSE算法性能差不多，但ZF-SIC 算法和 MMSE-SIC算法性能仍不及 ML算法的性能。图 4-2 和图 4-3 是研究在 ZL 算法种发射天线与接收天线不同组合时检测性能的优劣。可以看出，当发收天线数一致时，不影响ZF算法的检测性能。 当发射天线数大于接收天线数时，算法检测性能不受信噪比的影响。 发射天线数小于接收天线数时，算法有更佳的检测性能。.专业 .整理 .下载可编辑六、参考文献1 GB/T 7714 赵莉 . MIMO 检测算法 J. 2010.2 龚兵 , 邓飞其 . LTE 系统中的 MIMO检测算法 J. 重庆理工大学学报自然科学版 , 2009, 23(9):1
21、15-117.3 王利勇 . MIMO 检测算法的研究与实现 D. 西安电子科技大学, 2010.4 肖莉 . LTE 系统中 MIMO检测算法的研究 D. 安徽大学 ,2012.5 刘文龙 , 裴莹莹 , 金明录 . BPSK 通信系统的部分最优 MIMO检测算法 J.信号处理 , 2013, 29(10):1315-1322.6 陈雯柏 , 李卫 , 张小频 . 一种低复杂度的近似最大似然MIMO检测算法 J.哈尔滨工业大学学报 , 2012, 44(5):140-143.7 王灵光 , 刘皓 , 王军 , 等. MIMO 检测算法研究与实现 C/中国西部青年通信学术会议. 2005.专业
22、 .整理 .下载可编辑主要代码% 经过 BPSK调制的 V-Blast 发射矩阵 X七、附录X=zeros(tx,L);for k=1:txX(k,:)=(-1).(A(k:tx:end)+1);End% Rayleigh 信道 H H=sqrt(1/2)*(randn(rx,tx,L)+i*randn(rx,tx,L) n=sqrt(1/2)*(randn(rx,L)+i*randn(rx,L); R=zeros(rx,L);%ZF检测disp('berz');berz=;for m=SNRsnr=10(m/10);R_noised=R+sqrt(1/snr)*n;x=;a=
23、zeros(tx*L,1);for t=1:Lr=R_noised(:,t);G=pinv(H(:,:,t);y=G*r;xtemp=(y>=0)-(y<0)+0;x=x,xtemp;endfor k=1:txa(k:tx:end)=(x(k:tx:end)+1)/2;enderrbit,temp_ber=biterr(A,a);berz=berz,temp_ber;end% ZF-SIC(ordered)算法disp('ZF-SIC');berzs=;for m=SNRmsnr=10(m/10);R_noised=R+sqrt(1/snr)*n;x=;a=zero
24、s(tx*L,1);.专业 .整理 .for t=1:r=RHHG=S=xteforenx=endfor k=1a(enderrbit,berzs=end% MMSE算法disp('MMSE')berm=;for m=SNRmsnr=10R_noisx=;a=zerofor t=1r=Hxtw=inv(HH'*HHy=xtx=end下载可编辑for k=1:txa(k:tx:end)=(x(k:tx:end)+1)/2;enderrbit,temp_ber=biterr(A,a);berm=berm,temp_ber;end% %MMSE算法disp('MMSE
25、')berm=;for m=SNRsnr=10(m/10);R_noised=R+sqrt(1/snr)*n;x=;a=zeros(tx*L,1);for t=1:Lr=R_noised(:,t);HH=H(:,:,t);xtemp=zeros(tx,1);w=inv(HH'*HH+(1/snr)*eye(tx)*HH'y=w*r;xtemp=(y>=0)-(y<0)+0;x=x,xtemp;end%ML算法for ii = 1:length(Eb_N0_dB)ip = rand(1,N)>0.5;s = 2*ip-1;sMod=kron(s,ones
26、(nRx,1);sMod=reshape(sMod,nRx,nTx,N/nTx);h = 1/sqrt(2)*randn(nRx,nTx,N/nTx) + j*randn(nRx,nTx,N/nTx); % Rayleigh channeln = 1/sqrt(2)*randn(nRx,N/nTx) + j*randn(nRx,N/nTx); % white gaussian noise, 0dB variancey = squeeze(sum(h.*sMod,2) + 10(-Eb_N0_dB(ii)/20)*n;sHat1 = 1 1;.专业 .整理 .sHat1 = repmat(sHa
27、t1,1 ,N/2);sHat1Mod = kron(sHat1,ones(nRx,1);sHat2 = 1 -1;sHat2 = repmat(sHat2,1 ,N/2);sHat2Mod = kron(sHat2,ones(nRx,1);sHat2Mod = reshape(sHat2Mod,nRx,nTx,N/nTx);J10 = sum(abs(y - zHat2),1);sHat3 = -1 1;sHat3 = repmat(sHat3,1 ,N/2);sHat3Mod = kron(sHat3,ones(nRx,1);sHat3Mod = reshape(sHat3Mod,nRx,
28、nTx,N/nTx);zHat3 = squeeze(sum(h.*sHat3Mod,2) ;J01 = sum(abs(y - zHat3),1);sHat4 = -1 -1;sHat4 = repmat(sHat4,1 ,N/2);sHat4Mod = kron(sHat4,ones(nRx,1);sHat4Mod = reshape(sHat4Mod,nRx,nTx,N/nTx);zHat4 = squeeze(sum(h.*sHat4Mod,2) ;J00 = sum(abs(y - zHat4),1); function wk,kk = minnorm(G,MM)if nargin=1M=1:size(G,1);elseM=;for i=1:length(MM)if MM(i)=0temp=MM(i);M=M,temp;endendendGG=;for i=1:size(G,1)GG=GG;norm(G(i,:);