MIMO通信系统的信道估计与信号检测

MIMO通信系统的信道估计与信号检测项目意义义一项目意多输入多输出（MIMO）技术由于能够在不增加传输带宽的条件下成倍的提高无线信道的信道容量,因而被认为是下一代移动通信系统4G的关键技术之一。

MIMO技术是未来无线通信系统中实现高数据速率传输、改善传输质量、提高系统容量的重要途径。

MIMO信道模型无论是在MIMO技术的理论研究阶段还是在MIMO系统的应用阶段都是必需的。

因此，MIMO信道的建模是MIMO理论研究中的重要内容。

多输入多输出（MIMO）衰落信道是迄今为止所考虑的单输入单输出（SISO）随机信道的多变量推广。

从SISO入手，逐步增加天线数，通过对MIMO 信道的建模和仿真，深刻理解MIMO的系统的内涵。

二项目内容1.MIMO信道的建模。

搭建1*1，2*2,4*4,8*8，MIMO－任一路的信道符合Rayleigh Fading。

2.在接收端基于导频的信道估计。

3.利用估计的信道分别进行MLD和Zero-forcing信号检测。

4.1×1,2×2,4×4,8×8，（理想信道）模型的传输性能比较。

5.1×1,2×2,4×4,8×8，（估计信道）模型的传输性能比较。

6.估计信道和理想信道（4×4）之间的传输性能比较。

三项目原理(1)MIMO系统模型以2×2MIMO为例：r1=H11*S1+H21*S2+n1n2r2=H12*S1+H22*S2+说明:H信道符合Rayleigh衰落。

n为信道的高斯白噪声。

S为发射信号，r为接收端接收信号。

(2)基于导频的信道估计在2×2MIMO信道模型中，导引信号的数量可以是2当导引信号时p1p2=[10],r1=H11*p1+H21*p2+n1(p1=0),不考虑噪声的影响n2(p1=0),不考虑噪声的影响。

r2=H12*S1+H22*S2+则有：H11=r1/p1；H12=r2/p1；当导引信号时p1p2=[01],r1=H11*p1+H21*p2+n1(p1=0),不考虑噪声的影响r2=H12*S1+H22*S2n2(p1=0),不考虑噪声的影响。

+则有：H21=r1/p2；H22=r2/p2；因此，采用导引结构，可以在接收端估计出信道H为：H=[H11H12H21H21]（3）最大似然检测（MLD）这是理论最优的检测方法，能够完全获得接收分集增益。

由估计理论可知，对于Y=HX+Z中的H的估计，最大似然估计方法构造了一种代价函数p(H|Y,X)，使该代价函数取得最大值的H为最终的估计值：本项目主要针对4×4的MIMO系统进行仿真，星座映射图选择QPSK，在接收端采用最大似然检测算法，主要是考察信道估计算法的性能，所以在检测时采用对发射信号空间进行了穷尽搜索的方法来进行最大似然检测，以充分展示信道估计算法的BER性能。

（4）迫零检测（ZF）迫零检测算法是给接收信号乘以信道矩阵的逆，其他用户对它的的干扰可以消除，但同时噪声也乘以信道矩阵的逆，一般来说，信道矩阵的系数都小于1，那么它的逆就是大于1的，也就是说给噪声乘了一个大于1的因子，必然是放大了噪声。

四项目分析图一图二图三图四结果分析：从图一图三中可以看出，在MLD检测中，当发射和接收的天线数越多，误码率越低，从而证明了MIMO通信系统可以克服多径衰落带来的不良影响，实现了信号传输的可靠性，而且还可以增加系统容量，提高频谱利用率。

从图二图四中可以看出，在ZF检测中，当发射和接收的天线数越多，误码率越高，迫零检测算法是给接收信号乘以信道矩阵的逆，其他用户对它的的干扰可以消除，但同时噪声也乘以信道矩阵的逆，一般来说，信道矩阵的系数都小于1，那么它的逆就是大于1的，也就是说给噪声乘了一个大于1的因子，必然是放大了噪声。

所以误码率会随着天线数的增多而增大。

图五图六结果分析从结果可以看出，使用MLD检测时，理想信道和估计信道误码率总体趋势一样，但理想信道的误码率比估计信道的误码率要低，说明基于导频进行信道估计在MIMO通信系统是可行的。

使用ZF检测时，理想信道和估计信道误码率总体趋势一样，但理想信道的误码率比估计信道的误码率要高，所以一般不建议在基于导频的信道估计中使用ZF检测算法。

图七图八结果分析从结果可以看出，使用理想信道时，MLD算法比ZF算法的误码率小，且总体趋势一样，但在高信噪比部分MLD算法误码率比ZF算法的误码率明显低了很多，说明MLD算法全局搜索的性能极佳。

使用估计信道时，MLD算法和ZF算法的误码率相差不大，其中MLD算法的信噪比稍低一些。

从运算量和运行时间来看，当num=256，loop=1000时，使用MLD算法一共大约1200s，而使用ZF算法一共需要时间大约180s。

MLD算法的运算量远大于ZF算法，所用时间也长于ZF算法。

五项目总结通过本课程的学习，我们知道MIMO将多径无线信道与发射、接收视为一个整体进行优化，从而实现高的通信容量和频谱利用率。

这是一种近于最优的空域时域联合的分集和干扰对消处理。

而本项目通过建立MIMO通信系统，并基于导频估计出信道H，再利用估计的信道分别进行MLD和Zero-forcing信号检测，最后应用MATLAB对本系统进行仿真，画图，观察图形，通过仿真结果进行误码率分析。

从上面进行到的项目过程，我们进一步了解到MIMO系统突出的性能优势，为未来学习4G通信打下基础。

六参考文献[1].樊昌信通信原理北京国防工业出版社2001[2].魏巍MATLAB信息工程工具箱技术手册北京国防工业出版社2004[3].陈永春MATLAB M语言高级编程北京清华大学出版社2004[4].肖扬MIMO多天线无线通信系统北京人民邮电出版社2009[5].章坚武移动通信陕西西安电子科技大学出版社2003七、附录%4x4天线MLD检测，理想信道与估计信道的MIMO系统误码率分析%clear all;close;Eb_N0=-10:1:10;tx=4;%发送天线数目%rx=4;%接收天线数目%num=256;loop=1000;%循环次数（100/BER=1000）%SNR_num=length(Eb_N0);%number of SNRerror=zeros(2,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列modSig=modBaseBandQPSK(Sig);%QPSK调制%估计信道发送矩阵初始化%transmit_array_estimated=zeros(tx,num/tx/2+tx);transmit_array_estimated(:,1:tx)=eye(tx,tx);transmit_array_estimated(:,tx+1:num/tx/2+tx)=series2array(modSig,tx)./sqrt(tx);%将序列转换为数组，用于送给不同天线，有功率分配%%理想信道发送矩阵初始化%transmit_array_ideal=series2array(modSig,tx)./sqrt(tx);%4x4MIMO信道%Channel=zeros(tx,rx,num);for n=1:txfor m=1:rxChannel(n,m,:)=wgn(1,1,0,'complex');endendH=Channel;%通过MIMO理想信道%for n=1:num/tx/2recSig_ideal(:,n)=H(:,:,n)*transmit_array_ideal(:,n);endfor n=1:txrecSig_ideal(n,:)=awgn(recSig_ideal(n,:),Eb_N0(p));end%通过MIMO估计信道%for n=1:num/tx/2+txrecSig_estimated(:,n)=H(:,:,n)*transmit_array_estimated(:,n);%信号通过4x4MIMO信道end%加awgn噪声%for n=1:txrecSig_estimated(n,:)=awgn(recSig_estimated(n,:),Eb_N0(p));end%估计信道H_estimated=recSig_estimated(:,1:tx);%信道估计%receive=recSig_estimated(:,tx+1:num/tx/2+tx);%有用信息（去除信道估计用的信息后）%%最大似然检测%S_estimated=MLD1(receive,H_estimated);%估计信道%S_ideal=MLD1(recSig_ideal,H(:,:,num));%理想信道%%将数组转换为序列%x_estimated=conj(S_estimated(:)');%估计信道%x_ideal=conj(S_ideal(:)');%理想信道%%QPSK解调%demodSig_estimated=demodQPSK(x_estimated);%估计信道%demodSig_ideal=demodQPSK(x_ideal);%理想信道%%错误点数%error(1,p)=error(1,p)+length(find((demodSig_estimated-Sig)~=0));%估计信道%error(2,p)=error(2,p)+length(find((demodSig_ideal-Sig)~=0));%理想信道% enderror(:,p)=error(:,p)/loop/num%误码率%tocend%画误码率曲线semilogy(Eb_N0,error(1,:),'r*-',Eb_N0,error(2,:),'b*-');grid on;title('4×4天线，MLD检测，理想信道和估计信道的性能比较(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('估计信道误码率曲线','理想信道误码率曲线');====================================================================== ===========================%4x4天线,理想信道vs估计信道,ZF检测%clear all;close;Eb_N0=-10:1:10;tx=4;%发送天线数目%rx=4;%接收天线数目%num=256;loop=1000;%循环次数（100/BER=1000）%SNR_num=length(Eb_N0);%number of SNRBER=zeros(2,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列modSig=modBaseBandQPSK(Sig);%QPSK调制%估计信道发送矩阵初始化%transmit_array_estimated=zeros(tx,num/tx/2+tx);transmit_array_estimated(:,1:tx)=eye(tx,tx);transmit_array_estimated(:,tx+1:num/tx/2+tx)=series2array(modSig,tx)./sqrt(tx);%将序列转换为数组，用于送给不同天线，有功率分配%%理想信道发送矩阵初始化%transmit_array_ideal=series2array(modSig,tx)./sqrt(tx);%4x4MIMO信道%Channel=zeros(tx,rx,num);for n=1:txfor m=1:rxChannel(n,m,:)=wgn(1,1,0,'complex');endendH=Channel;h=(1/sqrt(2))*H(:,:,num)/(sqrt(tx));%对信道进行归一化处理%%通过MIMO理想信道%for n=1:num/tx/2recSig_ideal(:,n)=H(:,:,n)*transmit_array_ideal(:,n);endfor n=1:txrecSig_ideal(n,:)=awgn(recSig_ideal(n,:),Eb_N0(p));end%通过MIMO估计信道%for n=1:num/tx/2+txrecSig_estimated(:,n)=H(:,:,n)*transmit_array_estimated(:,n);%信号通过4x4MIMO信道end%加awgn噪声%for n=1:txrecSig_estimated(n,:)=awgn(recSig_estimated(n,:),Eb_N0(p));end%估计信道H_estimated=recSig_estimated(:,1:tx);%信道估计%receive=recSig_estimated(:,tx+1:num/tx/2+tx);%有用信息（去除信道估计用的信息后）%%迫零检测%%估计信道%y_estimated=pinv(h)*receive;x_estimated=conj(y_estimated(:)');demodSig_estimated_ZF=demodQPSK(x_estimated);%QPSK解调%%理想信道%y_ideal=pinv(h)*recSig_ideal;x_ideal=conj(y_ideal(:)');demodSig_ideal_ZF=demodQPSK(x_ideal);%QPSK解调%%错误点数%BER(1,p)=BER(1,p)+length(find((demodSig_estimated_ZF-Sig)~=0));%估计信道%BER(2,p)=BER(2,p)+length(find((demodSig_ideal_ZF-Sig)~=0));%理想信道%endBER(:,p)=BER(:,p)/loop/num%误码率%tocend%画误码率曲线semilogy(Eb_N0,BER(1,:),'r*-',Eb_N0,BER(2,:),'b*-');grid on;title('4×4天线，ZF检测，通过理想信道和估计信道的性能比较(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('估计信道误码率曲线','理想信道误码率曲线');====================================================================== =========================%估计信道，MLD检测的1x1,2x2,4x4天线的MIMO系统误码率曲线比较分析%clear all;close;tx1=1;%发送天线数目为1%rx1=1;%接收天线数目为1%tx2=2;%发送天线数目为2%rx2=2;%接收天线数目为2%tx3=4;%发送天线数目为4%rx3=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数%100/BER=1000SNR_num=length(Eb_N0);BER=zeros(3,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%%发送信号矩阵初始化%transmit_array1=zeros(tx1,num/tx1/2+tx1);transmit_array1(:,1:tx1)=eye(tx1,tx1);transmit_array2=zeros(tx2,num/tx2/2+tx2);transmit_array2(:,1:tx2)=eye(tx2,tx2);transmit_array3=zeros(tx3,num/tx3/2+tx3);transmit_array3(:,1:tx3)=eye(tx3,tx3);%发送序列转换为数组，用于送给不同天线，有功率分配%transmit_array1(:,tx1+1:num/tx1/2+tx1)=series2array(modSig,tx1)./sqrt(tx1);transmit_array2(:,tx2+1:num/tx2/2+tx2)=series2array(modSig,tx2)./sqrt(tx2);transmit_array3(:,tx3+1:num/tx3/2+tx3)=series2array(modSig,tx3)./sqrt(tx3);%1x1MIMO信道%Channel=zeros(tx1,rx1,num);for n1=1:tx1for m1=1:rx1Channel(n1,m1,:)=wgn(1,1,0,'complex');endendH1=Channel;%2x2MIMO信道%Channel2=zeros(tx2,rx2,num);for n2=1:tx2for m2=1:rx2Channel2(n2,m2,:)=wgn(1,1,0,'complex');endendH2=Channel2;%4x4MIMO信道%Channel3=zeros(tx3,rx3,num);for n3=1:tx3for m3=1:rx3Channel3(n3,m3,:)=wgn(1,1,0,'complex');endendH3=Channel3;%信号通过MIMO信道%for n=1:num/tx1/2+tx1recSig1(:,n)=H1(:,:,n)*transmit_array1(:,n);%信号通过1x1MIMO信道endfor n=1:num/tx2/2+tx2recSig2(:,n)=H2(:,:,n)*transmit_array2(:,n);%信号通过2x2MIMO信道endfor n=1:num/tx3/2+tx3recSig3(:,n)=H3(:,:,n)*transmit_array3(:,n);end%加awgn噪声%for n=1:tx1%1x1recSig1(n,:)=awgn(recSig1(n,:),Eb_N0(p));endfor n=1:tx2%2x2recSig2(n,:)=awgn(recSig2(n,:),Eb_N0(p));endfor n=1:tx3recSig3(n,:)=awgn(recSig3(n,:),Eb_N0(p));end%信道估计%H_estimated1=recSig1(:,1:tx1);%1x1H_estimated2=recSig2(:,1:tx2);%2x2H_estimated3=recSig3(:,1:tx3);%4x4%有用信息（去除信道估计用的信息后）%receive1=recSig1(:,tx1+1:num/tx1/2+tx1);receive2=recSig2(:,tx2+1:num/tx2/2+tx2);receive3=recSig3(:,tx3+1:num/tx3/2+tx3);%最大似然检测%S1=MLD1(receive1,H_estimated1);%1x1S2=MLD1(receive2,H_estimated2);%2x2S3=MLD1(receive3,H_estimated3);%4x4%将数组转换为序列%x1=conj(S1(:)');x2=conj(S2(:)');x3=conj(S3(:)');%QPSK解调%demodSig1=demodQPSK(x1);demodSig2=demodQPSK(x2);demodSig3=demodQPSK(x3);BER(1,p)=BER(1,p)+length(find((demodSig1-Sig)~=0));%1x1错误点数BER(2,p)=BER(2,p)+length(find((demodSig2-Sig)~=0));%2x2错误点数BER(3,p)=BER(3,p)+length(find((demodSig3-Sig)~=0));%4x4错误点数endBER(:,p)=BER(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,BER(1,:),'r*-',Eb_N0,BER(2,:),'g*-',Eb_N0,BER(3,:),'b*-');grid on;title('导频估计信道,MLD检测，MIMO系统1x1,2x2,4x4天线性能比较(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('1x1天线误码率曲线','2x2天线误码率曲线','4x4天线误码率曲线');====================================================================== ============================%估计信道,MLD和ZF检测下的4x4天线MIMO系统误码率曲线比较分析%clear all;close;tx=4;%发送天线数目为4%rx=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数%100/BER=1000SNR_num=length(Eb_N0);BER=zeros(2,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%transmit_array=zeros(tx,num/tx/2+tx);transmit_array(:,1:tx)=eye(tx,tx);%4x4天线导频信号%%将序列转换为数组，用于送给不同天线，有功率分配%transmit_array(:,tx+1:num/tx/2+tx)=series2array(modSig,tx)./sqrt(tx);%4x4MIMO信道%Channel=zeros(tx,rx,num);for n=1:txfor m=1:rxChannel(n,m,:)=wgn(1,1,0,'complex');endendH=Channel;%信号通过MIMO信道%for n=1:num/tx/2+txrecSig(:,n)=H(:,:,n)*transmit_array(:,n);end%加awgn噪声%for n=1:txrecSig(n,:)=awgn(recSig(n,:),Eb_N0(p));end%信道估计%H_estimated=recSig(:,1:tx);%4x4receive=recSig(:,tx+1:num/tx/2+tx);%有用信息（去除信道估计用的信息后）%%最大似然检测%S=MLD1(receive,H_estimated);%4x4x_MLD=conj(S(:)');%将数组转换为序列%demodSig_MLD=demodQPSK(x_MLD);%QPSK解调%BER(1,p)=BER(1,p)+length(find((demodSig_MLD-Sig)~=0));%错误点数%迫零检测%h=(1/sqrt(2))*H(:,:,num)/(sqrt(tx));%对信道进行归一化处理%y=pinv(h)*receive;x_ZF=conj(y(:)');demodSig_ZF=demodQPSK(x_ZF);BER(2,p)=BER(2,p)+length(find((demodSig_ZF-Sig)~=0));%错误点数endBER(:,p)=BER(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,BER(1,:),'r*-',Eb_N0,BER(2,:),'g*-');grid on;title('估计信道,MLD和ZF检测下的4x4天线MIMO系统误码率曲线比较分析(11组)'); xlabel('E_b/N_0(dB)');ylabel('误码率');legend('MLD检测天线误码率曲线','ZF检测天线误码率曲线');====================================================================== ================%导频序列估计信道,ZF检测下,MIMO系统1x1,2x2,4x4天线性能比较%clear all;close;tx1=1;%发送天线数目为1%rx1=1;%接收天线数目为1%tx2=2;%发送天线数目为2%rx2=2;%接收天线数目为2%tx3=4;%发送天线数目为4%rx3=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数%100/BER=1000SNR_num=length(Eb_N0);BER=zeros(3,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%transmit_array1=zeros(tx1,num/tx1/2+tx1);transmit_array1(:,1:tx1)=eye(tx1,tx1);%1x1天线导频信号%transmit_array2=zeros(tx2,num/tx2/2+tx2);transmit_array2(:,1:tx2)=eye(tx2,tx2);%2x2天线导频信号%transmit_array3=zeros(tx3,num/tx3/2+tx3);transmit_array3(:,1:tx3)=eye(tx3,tx3);%4x4天线导频信号%%将序列转换为数组，用于送给不同天线，有功率分配%transmit_array1(:,tx1+1:num/tx1/2+tx1)=series2array(modSig,tx1)./sqrt(tx1);transmit_array2(:,tx2+1:num/tx2/2+tx2)=series2array(modSig,tx2)./sqrt(tx2);transmit_array3(:,tx3+1:num/tx3/2+tx3)=series2array(modSig,tx3)./sqrt(tx3);%1x1MIMO信道%Channel=zeros(tx1,rx1,num);for n1=1:tx1for m1=1:rx1Channel(n1,m1,:)=wgn(1,1,0,'complex');endendH1=Channel;%2x2MIMO信道%Channel2=zeros(tx2,rx2,num);for n2=1:tx2for m2=1:rx2Channel2(n2,m2,:)=wgn(1,1,0,'complex');endendH2=Channel2;%4x4MIMO信道%Channel3=zeros(tx3,rx3,num);for n3=1:tx3for m3=1:rx3Channel3(n3,m3,:)=wgn(1,1,0,'complex');endendH3=Channel3;%信号通过MIMO信道%for n=1:num/tx1/2+tx1recSig1(:,n)=H1(:,:,n)*transmit_array1(:,n);%信号通过1x1MIMO信道endfor n=1:num/tx2/2+tx2recSig2(:,n)=H2(:,:,n)*transmit_array2(:,n);%信号通过2x2MIMO信道endfor n=1:num/tx3/2+tx3recSig3(:,n)=H3(:,:,n)*transmit_array3(:,n);%信号通过4x4MIMO信道end%加awgn噪声%for n=1:tx1%1x1recSig1(n,:)=awgn(recSig1(n,:),Eb_N0(p));endfor n=1:tx2%2x2recSig2(n,:)=awgn(recSig2(n,:),Eb_N0(p));endfor n=1:tx3%4x4recSig3(n,:)=awgn(recSig3(n,:),Eb_N0(p));end%信道估计%H_estimated1=recSig1(:,1:tx1);%1x1H_estimated2=recSig2(:,1:tx2);%2x2H_estimated3=recSig3(:,1:tx3);%4x4%有用信息（去除信道估计用的信息后）%receive1=recSig1(:,tx1+1:num/tx1/2+tx1);receive2=recSig2(:,tx2+1:num/tx2/2+tx2);receive3=recSig3(:,tx3+1:num/tx3/2+tx3);%迫零检测%h1=(1/sqrt(2))*H1(:,:,num)/(sqrt(tx1));h2=(1/sqrt(2))*H2(:,:,num)/(sqrt(tx2));h3=(1/sqrt(2))*H3(:,:,num)/(sqrt(tx3));y1=pinv(h1)*receive1;y2=pinv(h2)*receive2;y3=pinv(h3)*receive3;%将数组转换为序列%x1=conj(y1(:)');x2=conj(y2(:)');x3=conj(y3(:)');%QPSK解调%demodSig1=demodQPSK(x1);demodSig2=demodQPSK(x2);demodSig3=demodQPSK(x3);BER(1,p)=BER(1,p)+length(find((demodSig1-Sig)~=0));%1x1错误点数BER(2,p)=BER(2,p)+length(find((demodSig2-Sig)~=0));%2x2错误点数BER(3,p)=BER(3,p)+length(find((demodSig3-Sig)~=0));%4x4错误点数endBER(:,p)=BER(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,BER(1,:),'r*-',Eb_N0,BER(2,:),'g*-',Eb_N0,BER(3,:),'b*-');grid on;title('导频序列估计信道,ZF检测下,MIMO系统1x1,2x2,4x4天线性能比较(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('1x1天线误码率曲线','2x2天线误码率曲线','4x4天线误码率曲线');====================================================================== ===============%理想信道,ZF检测的1x1,2x2,4x4天线的MIMO系统误码率比较%clear all;close;tx1=1;%发送天线数目为1%rx1=1;%接收天线数目为1%tx2=2;%发送天线数目为2%rx2=2;%接收天线数目为2%tx3=4;%发送天线数目为4%rx3=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数（100/BER=1000）%SNR_num=length(Eb_N0);error=zeros(3,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%%将序列转换为数组，用于送给不同天线，有功率分配%transmit_array1=series2array(modSig,tx1)./sqrt(tx1);transmit_array2=series2array(modSig,tx2)./sqrt(tx2);transmit_array3=series2array(modSig,tx3)./sqrt(tx3);%产生MIMO1x1信道%Channel1=zeros(tx1,rx1,num);for n1=1:tx1for m1=1:rx1Channel1(n1,m1,:)=wgn(1,1,0,'complex');endendH1=Channel1;h1=(1/sqrt(2))*H1(:,:,num)/(sqrt(tx1));%产生MIMO2x2信道%Channel2=zeros(tx2,rx2,num);for n2=1:tx2for m2=1:rx2Channel2(n2,m2,:)=wgn(1,1,0,'complex');endendH2=Channel2;h2=(1/sqrt(2))*H2(:,:,num)/(sqrt(tx2));%产生MIMO4x4信道%Channel3=zeros(tx3,rx3,num);for n3=1:tx3for m3=1:rx3Channel3(n3,m3,:)=wgn(1,1,0,'complex');endendH3=Channel3;h3=(1/sqrt(2))*H3(:,:,num)/(sqrt(tx3));%信号通过MIMO信道%for n=1:num/tx1/2recSig1(:,n)=H1(:,:,n)*transmit_array1(:,n);endfor n=1:num/tx2/2recSig2(:,n)=H2(:,:,n)*transmit_array2(:,n);endfor n=1:num/tx3/2recSig3(:,n)=H3(:,:,n)*transmit_array3(:,n);end%加噪声%for n=1:tx1recSig1(n,:)=awgn(recSig1(n,:),Eb_N0(p));endfor n=1:tx2recSig2(n,:)=awgn(recSig2(n,:),Eb_N0(p));endfor n=1:tx3recSig3(n,:)=awgn(recSig3(n,:),Eb_N0(p));end%迫零检测%y1=pinv(h1)*recSig1;y2=pinv(h2)*recSig2;y3=pinv(h3)*recSig3;%将数组转换为序列%x_ZF1=conj(y1(:)');x_ZF2=conj(y2(:)');x_ZF3=conj(y3(:)');%QPSK解调%demodSig_ZF1=demodQPSK(x_ZF1);demodSig_ZF2=demodQPSK(x_ZF2);demodSig_ZF3=demodQPSK(x_ZF3);%错误点数%error(1,p)=error(1,p)+length(find((demodSig_ZF1-Sig)~=0));error(2,p)=error(2,p)+length(find((demodSig_ZF2-Sig)~=0));error(3,p)=error(3,p)+length(find((demodSig_ZF3-Sig)~=0));enderror(:,p)=error(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,error(1,:),'r*-',Eb_N0,error(2,:),'g*-',Eb_N0,error(3,:),'b*-');grid on;title('理想信道下ZF检测,MIMO系统1x1,2x2,4x4天线性能(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('1x1天线误码率曲线','2x2天线误码率曲线','4x4天线误码率曲线');====================================================================== ===============%理想信道MLD的1x1,2x2,4x4天线的MIMO系统误码率比较%clear all;close;tx1=1;%发送天线数目为1%rx1=1;%接收天线数目为1%tx2=2;%发送天线数目为2%rx2=2;%接收天线数目为2%tx3=4;%发送天线数目为4%rx3=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数（100/BER=1000）%SNR_num=length(Eb_N0);error=zeros(3,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%%将序列转换为数组，用于送给不同天线，有功率分配%transmit_array1=series2array(modSig,tx1)./sqrt(tx1);transmit_array2=series2array(modSig,tx2)./sqrt(tx2);transmit_array3=series2array(modSig,tx3)./sqrt(tx3);%产生MIMO1x1信道%Channel1=zeros(tx1,rx1,num);for n1=1:tx1for m1=1:rx1Channel1(n1,m1,:)=wgn(1,1,0,'complex');endendH1=Channel1;%产生MIMO2x2信道%Channel2=zeros(tx2,rx2,num);for n2=1:tx2for m2=1:rx2Channel2(n2,m2,:)=wgn(1,1,0,'complex');endendH2=Channel2;%产生MIMO4x4信道%Channel3=zeros(tx3,rx3,num);for n3=1:tx3for m3=1:rx3Channel3(n3,m3,:)=wgn(1,1,0,'complex');endendH3=Channel3;%信号通过MIMO信道%for n=1:num/tx1/2recSig1(:,n)=H1(:,:,n)*transmit_array1(:,n);endfor n=1:num/tx2/2recSig2(:,n)=H2(:,:,n)*transmit_array2(:,n);endfor n=1:num/tx3/2recSig3(:,n)=H3(:,:,n)*transmit_array3(:,n);end%加噪声%for n=1:tx1recSig1(n,:)=awgn(recSig1(n,:),Eb_N0(p));endfor n=1:tx2recSig2(n,:)=awgn(recSig2(n,:),Eb_N0(p));endfor n=1:tx3recSig3(n,:)=awgn(recSig3(n,:),Eb_N0(p));end%最大似然检测%S1=MLD1(recSig1,H1);S2=MLD1(recSig2,H2);S3=MLD1(recSig3,H3);%将数组转换为序列%x1=conj(S1(:)');x2=conj(S2(:)');x3=conj(S3(:)');%QPSK解调%demodSig1=demodQPSK(x1);demodSig2=demodQPSK(x2);demodSig3=demodQPSK(x3);%错误点数%error(1,p)=error(1,p)+length(find((demodSig1-Sig)~=0));error(2,p)=error(2,p)+length(find((demodSig2-Sig)~=0));error(3,p)=error(3,p)+length(find((demodSig3-Sig)~=0));enderror(:,p)=error(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,error(1,:),'r*-',Eb_N0,error(2,:),'g*-',Eb_N0,error(3,:),'b*-'); grid on;title('理想信道，MLD检测,MIMO系统1x1,2x2,4x4天线性能比较(11组)'); xlabel('E_b/N_0(dB)');ylabel('误码率');legend('1x1天线误码率曲线','2x2天线误码率曲线','4x4天线误码率曲线');====================================================================== ======%理想信道下,MIMO系统4x4天线MLD检测与ZF检测的性能比较%clear all;close;tx=4;%发送天线数目为4%rx=4;%接收天线数目为4%Eb_N0=-10:1:10;num=256;loop=1000;%循环次数（100/BER=1000）%SNR_num=length(Eb_N0);error=zeros(2,SNR_num);for p=1:SNR_numticfor w=1:loopSig=randint(1,num);%随机产生二进制序列%modSig=modBaseBandQPSK(Sig);%QPSK调制%%将序列转换为数组，用于送给不同天线，有功率分配%transmit_array=series2array(modSig,tx)./sqrt(tx);%产生MIMO4x4信道%Channel=zeros(tx,rx,num);for n=1:txfor m=1:rxChannel(n,m,:)=wgn(1,1,0,'complex');endendH=Channel;h=(1/sqrt(2))*H(:,:,num)/(sqrt(tx));%对信道进行归一化处理%%信号通过MIMO信道%for n=1:num/tx/2recSig(:,n)=H(:,:,n)*transmit_array(:,n);end%加噪声%for n=1:txrecSig(n,:)=awgn(recSig(n,:),Eb_N0(p));end%最大似然检测%S=MLD1(recSig,H);%将数组转换为序列%x_MLD=conj(S(:)');demodSig_MLD=demodQPSK(x_MLD);%QPSK解调%error(1,p)=error(1,p)+length(find((demodSig_MLD-Sig)~=0));%错误点数%%迫零检测%y=pinv(h)*recSig;%将数组转换为序列%x_ZF=conj(y(:)');%QPSK解调%demodSig_ZF=demodQPSK(x_ZF);error(2,p)=error(2,p)+length(find((demodSig_ZF-Sig)~=0));enderror(:,p)=error(:,p)/loop/num%误码率%tocend%画误码率曲线%semilogy(Eb_N0,error(1,:),'r*-',Eb_N0,error(2,:),'g*-');grid on;title('理想信道下,MIMO系统4x4天线MLD检测与ZF检测的性能比较(11组)');xlabel('E_b/N_0(dB)');ylabel('误码率');legend('MLD检测误码率曲线','ZF检测误码率曲线');====================================================================== =====function[array]=series2array(series,M)%This function is used to change series in to arrayL=length(series);array=zeros(M,L/M);for m=1:Mfor n=1:L/M;array(m,n)=series((n-1)*M+m);endend====================================================================== =======function[S]=MLD1(r,H)%This function is maximum likelihood detection function.%%r：Received signal%%H：Estimated channel or ideal channel%[N,len]=size(r);%len:length of Signal N:the num of antenna%S=zeros(N,len);%Each antenna sent p signal(Power allocation)%s=[exp(1j*pi*5/4)exp(1j*pi*3/4)exp(1j*pi*7/4)exp(1j*pi/4)]./sqrt(N); X=zeros(N,4^N);k=1;if N==1%antenna1x1X=s;else if N==2%antenna2x2for t=1:4for p=1:4X(1,k)=s(1,t);X(2,k)=s(1,p);k=k+1;endendelse if N==4%antenna4x4for t=1:4for p=1:4for u=1:4for q=1:4X(1,k)=s(1,t);X(2,k)=s(1,p);X(3,k)=s(1,u);X(4,k)=s(1,q);k=k+1;endendendendelse if N~=1&&N~=2&&N~=4disp('error');endendendend %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% matrix=zeros(N,4^N);distance=zeros(1,length(X));for n=1:lenif(ndims(H)~=2)%H is not the estimated channel%H(:,:,1)=H(:,:,n);endmatrix=H(:,:,1)*X;%HS matrix%for m=1:length(X)distance(1,m)=norm(matrix(:,m)-r(:,n));end[~,index]=min(distance);%Returns the index of the minimum range%S(:,n)=X(:,index);end====================================================================== ========================function y=modBaseBandQPSK(sig)%mod BaseBandQPSK modulate the signal of basebank in QPSK method(carrrier%frequency is0Hz.%y=modBaseBandQPSK(sig)modulated the signals specified by sig.y is%modulated signals.%---change the signals from binary to quaternariant and modulate[signal_length_row,signal_length_col]=size(sig);%Qsig=zeros(signal_length_row,signal_length_col);y=zeros(signal_length_row,signal_length_col/2);%----initialfor n=1:signal_length_rowfor m=1:signal_length_col/2Qsig=2*sig(n,2*m-1)+sig(n,2*m);switch(Qsig)case0%------00-->pi/4y(n,m)=exp(1i*pi/4);case1%------01-->7pi/4y(n,m)=exp(1i*7*pi/4);case2%------10-->3pi/4y(n,m)=exp(1i*3*pi/4);case3%------11-->5pi/4y(n,m)=exp(1i*5*pi/4);end;endendend====================================================================== =========================function y=demodQPSK(rec_Sig)%This function is used for QPSK signal demodulation.%[signal_length_row,signal_length_col]=size(rec_Sig);y=zeros(signal_length_row,signal_length_col*2);for n=1:signal_length_rowfor m=1:signal_length_colif real(rec_Sig(n,m))>0%'00''01'-->R[x]>0%y(n,2*m-1)=0;else%'10''11'-->R[x]<0%y(n,2*m-1)=1;endif imag(rec_Sig(n,m))>0%'00''01'-->I[x]>0%y(n,2*m)=0;else%'10''11'-->I[x]<0%y(n,2*m)=1;endendendend31/31。