MIMO信道的信道容量
 
//受工作时间所限，本文中的仿真结果均未在博客中做出整理，如对仿真结果感兴趣，本文附带了本文中所有结果的仿真程序，请自行修改参数后使用Matlab仿真查看
假设信道容量分析模型的信道的衰落类型是频率非选择性衰落，信道系数为0均值单位方差的复高斯随机变量。假设系统的发射端配有NT根天线，接收端配有NR根天线，总的发射功率为P。假设发射端未知信道的瞬时状态信息，因此只能将所有的功率平分给所有的天线，所以每根发射天线的功率为P/NT，且每根接收天线上的噪声功率为sigma^2，于是每根接受天线上的信噪比（SNR）为

    这样就可以用NR*NT的复数矩阵表示信道矩阵，H的第ji个元素H(j,i)表示第i根发射天线到第j根接收天线的信道衰落系数。
    下面具体分析单输入单输出（SISO）、单输入多输出（SIMO）、多输入单输出（MISO）和多输入多输出（MIMO）四种情况下的信道容量。
1. SISO信道的信道容量
 
    采用单根天线发射和单根天线接收的通信系统称为SISO系统。对于加性高斯白噪声SISO信道，NT=NR=1，信道矩阵H=H(1,1)=1，假设接收端的信噪比为$\xi$（希腊字母），则根据香农公式，该信道的归一化容量可以表示为：

    该容量的取值不受编码或者信号设计复杂性的限制，即只要信噪比每增加3dB（2倍），信道容量增加1bps/Hz。
    实际的无线信道是时变的，会受到衰落的影响，如果用H表示信道系数的瞬时值，则信道容量可以表示为

    上式的结果是一个随机变量，可以计算其概率分布，瑞利平衰落的信道容量的累计概率分布的仿真结果如图XX所示。从图中可以看出，由于受到衰落的影响，其与加性高斯白噪声信道的信道容量有比较大的差距。
    从信道容量的累计概率分布图中可以提取出两个与基带系统设计相关的统计参数，一个是平均信道容量Cav，即C的所有样本的平均值，它表示该无线链路能够提供的平均数据传输速率；另一个参数是终端容量（Outage Capacity）Cout，它定义了确保传输的高可靠性的数据传输速率，即prob{C>Cout}。
2. MISO信道的信道容量
 
    对于多输入但出书信道，发射端配有NT根天线，接收端只有一根天线，即NR=1，信道矩阵H也随之变成一个矢量H=[H(1,1), H(1,2), … , H(1,NT)]，其中H(1,i)表示从发射端的第i根天线到达接收端的信道的信道系数。
    如果信道系数固定，则该信道的信道容量可以表示为：

其中。
    由于假定信道的系数固定，则受到功率归一化的影响，该信道容量不会随着发射天线数目的增长而增长。
    如果信道系数随机变化，则该信道容量的表达式同上，但是变成一个随H变化的随机变量，norm(H, ‘fro’)^2服从自由度为2NT的$\chi^{2}$（希腊字母）分布。
    图XX给出了信道容量随发射天线数变化的累计概率分布函数，仿真假定信道系数是复高斯随机变量，发射天线数分别取1，2，4，8，16，信噪比取为10dB和100dB，从图中可以看出，信道容量会随着发射天线数的增加而增加，增长的极限是加性高斯白噪声信道的信道容量，这就从信道容量的角度解释了分集能够带来的系统可靠性的提升。
3. SIMO信道的信道容量
 
    对于单输入多输出信道，发射端只有一根天线，接收端配有NR，其信道可以看成是一个NR维列矢量：H=[H(1,1), H(2,1), … , H(NT,1)]，其中H(j,1)是发射端的天线到接收端的第j根天线的信道系数。
    如果信道系数是固定的，则该信道的信道容量可以表示为：

其中。
    由于信道系数被归一化，从信道容量的计算公式可以看出，SIMO信道与SISO信道相比获得了大小为NR的阵列增益。
    类似地，信道系数是复高斯随机变量情况下的信道容量公式同上，C成为一个随机变量。
    图XX给出了信道容量随接收天线数变化的累计概率分布函数，仿真假定信道系数是复高斯随机变量，发射天线数分别取1，2，4，8，16，信噪比取为10dB和100dB。
4. MIMO信道的信道容量
 
    对于配备有NT根发射天线和NR根接收天线的MIMO信道，发射端在不知道传输信道的状态信息的条件下，其信道容量的公式可以表示为：

5. SISO、MISO、SIMO和MIMO信道的信道容量比较
    通过上述比较可以看出，MIMO方式能够通过多种手段改善平均信道容量和中断概率，其改善平均信道容量主要是通过阵列增益和复用增益实现的，改善中断概率主要是利用了分集增益，这与前文中提到的MIMO系统的三种主要增益是吻合的。
    图XX给出了四种信道情况下的比较累计概率分布曲线，可以从中看出分集增益、阵列增益和复用增益对信道容量的改善作用。
6. 天线相关对信道容量的影响
   MIMO信道的信道容量在很大程度上取决于天线之间的相关性和MIMO信道矩阵的满秩情况。如前所述，对于收发端都具有相关性的MIMO信道，其信道矩阵可以表示为：

    为了分析和比较天线的相关性对信道容量的具体影响，这里以2*2的信道为例进行仿真，相关系数分别取为0，0.2，0.5，0.9，结果如图XX所示。从图中可以看出，天线的相关系数越大，对信道容量的影响越大，而且信噪比越大，造成的损失越大。
7. 莱斯因子对信道容量的影响
    在同时存在直射路径和散射分量的传播环境中，MIMO系统一般采用莱斯信道模型来近似，其中莱斯银子K的选择决定了该模型中直射分量的大小。实际情况中，MIMO莱斯信道矩阵可以表示为：

    式中HLOS表示无散射，仅有直射路径的信道矩阵；HRay表示仅有散射、没有直射路径的信道矩阵。D^2和2sigma^2分别为直射分量和散射分量的功率，两者之比就是莱斯因子K。
    为了分析和比较莱斯因子对MIMO信道容量的具体影响，本文以2*2的MIMO信道为例进行计算机仿真，选择HRay为独立同分布的瑞利平衰落，HLOS=[1, 1; 1, 1]，而莱斯因子分别取0、1、10和100四种情况，得到的仿真结果如图XX所示，从图中可以看出当莱斯因子较小，如K=1时，对信道偶然那个亮的影响不大，而当K因子较大，如K=100时，对信道容量的影响较大，而且随着信噪比的增大，损害的程度不断加大。
    从这里我们可以看出，直达径对MIMO系统的信道容量起到了破坏作用，这与前文中指出的MIMO系统适合工作在富散射环境中的结论是吻合的。
8. 自适应功率分配的信道容量
    如果发射端获得了信道的状态信息，那么它可以利用信道状态信息对发射信息进行预处理，以获得最大的信道容量。
    对信道矩阵H执行SVD分解，得到

    在发射端预乘右矩阵V，在接收端乘上左矩阵U^H，则MIMO信道等效为对角矩阵Sigma，也就是说多条相互独立的链路可以分别传输数据。在这样的条件下，信道容量的公式可以改写为：

其中。
    在总功率受限的情况下找到C的最大值在数学上是一个拉格朗日乘数问题，在信息论上可以用最佳能量分配的注水算法[Thomas Cover, 1991]解决，其原理可以用图XX来描述.
 

    注水的原理是，由1/λ和子信道的状态来决定在各子信道中如何分配功率。（图里面P*就表示对那个子信道分配的功率多少，也就是1/λ和阶梯层的差值）。信道状态越好，分配的功率越多，越差就越少，当状态差的程度超过阈值1/λ就不在该子信道上分配功率了。
 
    注水算法的求解是一个迭代的过程，其复杂度相对较高，可以考虑寻找注水原理的目标函数的多项式逼近函数。
 
附：MIMO信道容量的Matlab仿真程序
 
%% Simulation of Capacity of MIMO channel
% Author: alpswy
% Email: alpswy@gmail.com
% Revise History:
%  2009-08-22 Created
 
%% Parameters for normal simulation
ChanType = 'Rayleigh'; % AWGN, Rayleigh,CorrRayleigh, Rice
LOOP = 10000; % number of channel tapsamples
const = 1/sqrt(2); % constant number
NT = 4; % number of receive antennas
NR = 4; % number of transmit antennas
SNRdB = 10; % SNR in dB
SNRdec = 10^(0.1*SNRdB); % SNR in decibel
 
%% Parameters for displaying
numGroup = 100; % variable used to countsamples
Range_MIN = 0; % minimum value of figure
Range_MAX = 6; % maximum value of figure
 
%% Parameters for simulation of correlationof channel
rho_t = 0;
rho_r = 0;
 
%% Parameters for simulation of Ricechannel
K = 1; % Rice factor
D = sqrt(K);
 
switch ChanType
   case 'AWGN'
       H = ones(NR, NT);
       C = log2(det(eye(NR) + SNRdec/NT*H*H'));
       C_1 = [C, C];
       plot(C_1, 0:1);
       axis([Range_MIN, Range_MAX, 0, 1]);
   case 'Rayleigh'
        C = zeros(1, LOOP);
       C_2 = zeros(1, numGroup);
       for idx_LOOP = 1:LOOP
           H = const*(randn(NR, NT) + 1j*randn(NR, NT));
           C(idx_LOOP) = log2(det(eye(NR) + SNRdec/NT*H*H'));
       end
       C_min = min(C);
       C_max = max(C);
       C_1 = hist(C, numGroup);
       C_2(1) = C_1(1);
       for idx = 2:numGroup
           C_2(idx) = C_2(idx-1) + C_1(idx);
       end
       C_2 = C_2 / LOOP;
       Range_MIN = C_min;
       Range_MAX = C_max;
       X_axe = Range_MIN + [0:(numGroup-1)]/numGroup*(Range_MAX-Range_MIN);
       plot(X_axe, C_2);
   case 'CorrRayleigh'
       C = zeros(1, LOOP);
       C_2 = zeros(1, numGroup);
       R_tx = [1, rho_t; rho_t, 1];% for 2*2 only
       R_rx = [1, rho_r; rho_r, 1];
       for idx_LOOP = 1:LOOP
           H = const*(randn(NR, NT) + 1j*randn(NR, NT));
           H_Corr = R_rx * H * R_tx;
           C(idx_LOOP) = log2(det(eye(NR) + SNRdec/NT*H_Corr*H_Corr'));
       end
       C_min = min(C);
       C_max = max(C);
       C_1 = hist(C, numGroup);
       C_2(1) = C_1(1);
       for idx = 2:numGroup
           C_2(idx) = C_2(idx-1) + C_1(idx);
       end
       C_2 = C_2 / LOOP;
       Range_MIN = C_min;
       Range_MAX = C_max;
       X_axe = Range_MIN + [0:(numGroup-1)]/numGroup*(Range_MAX-Range_MIN);
       plot(X_axe, C_2);
   case 'Rice'
       C = zeros(1, LOOP);
       C_2 = zeros(1, numGroup);
       for idx_LOOP = 1:LOOP
           H = const*(randn(NR, NT) + 1j*randn(NR, NT)) + D * ones(NR, NT);
           C(idx_LOOP) = log2(det(eye(NR) + SNRdec/NT*H*H'));
       end
       C_min = min(C);
       C_max = max(C);
       C_1 = hist(C, numGroup);
       C_2(1) = C_1(1);
       for idx = 2:numGroup
           C_2(idx) = C_2(idx-1) + C_1(idx);
       end
       C_2 = C_2 / LOOP;
       Range_MIN = C_min;
       Range_MAX = C_max;
       X_axe = Range_MIN + [0:(numGroup-1)]/numGroup*(Range_MAX-Range_MIN);
       plot(X_axe, C_2);
   otherwise
       disp('Other Channel Type is still under construction. Only AWGN andRayleigh fading is available.')
       keyboard
end

可以严格地证明在被高斯白噪声干扰的信道中传送的最大信息速率C由下述

公式确定

  
C=B*log2(1+S/N) =B*log2(1+S/n0*B)， log2表示以2为底的对数,n0为噪声单边功率谱密度，所以带宽无限大，信道容量并不能无限大！！！

上式表明，当给定S/ n0时，若带宽B趋于无穷大，信道容量不会趋于无限大，而只是S/
n0的1.44倍。这是因为当带宽B增大时，噪声功率也随之增大。

  该式通常称为香农公式。B是信道带宽赫S是信号功率瓦N是噪声

功率瓦。  

  香农公式中的S/N 为无量纲单位。如S/N=1000即信号功率是噪声功率

的1000倍

  但是当讨论信噪比S/N时常以分贝dB为单位。公式如下：

   S/N = 10lgS/N （dB）(分贝与信噪比之间的关系为：dB=10lgS/N) 

  公式表明:信道带宽限制了比特率的增加信道容量还取决于系统信噪比以及编

码技术种类  
连续信道的容量Ct和信道带宽B、信号功率S及噪声功率谱密度n0三个因素有关。 

        在写这篇文章之前我一直对“带宽”这个用语比较迷惑，比如我们通常说“带宽10M的网络”，“网络带宽是10M”等等，很容易跟信道的带宽比如10MHZ、100MHZ相混淆。我查阅了相关的书籍，终于在通信和电路上将这两个概念彻底搞清楚了。

1、网络带宽

        我们通常口头上说的网络带宽，用标准的术语来说应该是信道容量，单位为比特每秒(bit/s)。在通信原理中用香农公式计算出来的信道容量，就是我们口头上说的网络带宽单位为bit/s；同时我们应该注意这个公式中还有个带宽，是信道通频带宽，由电路频率特性决定，单位为HZ，也就是说只要电路设计确定了，那么这个通频带宽就确定了，就是我们通常说的高通、低通、带通滤波器。





2、信道带宽

        信道带宽就是我们经常说的高通、低通、带通，是由电路决定了的。由于波速跟媒质相关，根据c=λf，可知在同一媒质λ（波长）越大，f（频率）越低；而f就是带宽决定了信道的容量，即传输速率。所以在无线通信中，在功率一定的情况下，带宽越大传输的数据速率越高，传输距离就近；带宽越小传输的数据速率越低，传输距离就远。





 

 



3、奈奎斯特定理和香农定理

        由时候会将这两个定理搞混淆，所以这里简单说明一下，奈奎斯特定理是理想情况下的用于计算信道容量，香农定理就是实际工程中计算信道容量。

 

 

 

一、信道容量和吞吐量

信道容量是指互信息量的最大值。具体地说，在一定带宽和信噪比（信干比）下，借助某种编码方式实现无差错传输时可以达到的最大速率，是一个上界，在实际应用中，传输速率小于等于信道容量。考察的对象主要是实际物理信道传输信息的能力，不针对特定的传输技术及编码方案。
	
吞吐量是指某个系统在单位时间内正确传输的信息比特数，并非一个上界。如果引入了某种编码方案，那么在计算吞吐量的时候，一定要把冗余比特（如校验位）除去。另外，吞吐量更侧重于系统层面，例如某个小区的吞吐量。吞吐量考量的主要对象是传输技术（如接入方式，复用方式；双工方式，协作方式等）。在MAC层协议设计时，经常假设信道传输无误码，数据包不能被正确接收的原因在于丢包和碰撞，相应得到的吞吐量比实际吞吐量大，所以一般文章在计算吞吐量之前都会给出定义及相应的假设。
 

二、链路和信道

链路：

    链路指无源的点到点的物理连接。

    有线通信时，链路指两个节点之间的物理线路，如电缆或光纤。无线电通信时，链路指基站和终端之间传播电磁波的路径空间。水声通信时链路指换能器和水听器之间的传播声波的路径空间。

信道：

    信息是抽象的，但传送信息必须通过具体的媒质。例如二人对话，靠声波通过二人间的空气来传送，因而二人间的空气部分就是信道。邮政通信的信道是指运载工具及其经过的设施。无线电话的信道就是电波传播所通过的空间，有线电话的信道是电缆。每条信道都有特定的信源和信宿。在多路通信，例如载波电话中，一个电话机作为发出信息的信源，另一个是接收信息的信宿，它们之间的设施就是一条信道，这时传输用的电缆可以为许多条信道所共用。在理论研究中，一条信道往往被分成信道编码器、信道本身和信道译码器。人们可以变更编码器、译码器以获得最佳的通信效果，因此编码器、译码器往往是指易于变动和便于设计的部分，而信道就指那些比较固定的部分。但这种划分或多或少是随意的，可按具体情况规定。例如调制解调器和纠错编译码设备一般被认为是属于信道编码器、译码器的，但有时把含有调制解调器的信道称为调制信道；含有纠错编码器、译码器的信道称为编码信道。

无线信道：

    无线信道也就是常说的无线的“频段（Channel）”，其是以无线信号作为传输媒体的数据信号传送通道。

 

参考文献：信道容量与吞吐量之间的关系与区别？ - hikinglau的回答 - 知乎 https://www.zhihu.com/question/36364638/answer/104257616

                  计算机网络（第七版），谢希仁，电子工业出版社，2017.

                  链路-百度百科-https://baike.baidu.com/item/%E9%93%BE%E8%B7%AF/9410314?fr=aladdin

                  信道-百度百科-https://baike.baidu.com/item/%E4%BF%A1%E9%81%93/499862?fr=aladdin