利用Matlab计算多个BSC信道串联后的信道容量【P22214105宣烨 P22214107程瀚 P22214103蒋强理 P22214136王允闲 P22214117韩文静】

本文为安徽大学电子信息工程学院P22214105宣烨 、P22214107程瀚、 P22214103蒋强理、 P22214136王允闲、 P22214117韩文静的信息论与编码作业，能力不佳，望请指正。

---

前言

本文主要利用Matlab的矩阵计算的功能大大减少了计算量，并画出信道容量随
的变化情况

---

一、什么是BSC串联信道

1.BSC信道

二进制对称信道（Binary Symmetric Channel）是离散无记忆信道（discrete memoryless channel）在J=K=2时的特例。

它的输入和输出都只有0和1两种符号，并且发送0而接收到1，以及发送1而收到0（即误码）的概率相同，所以称信道是对称的。此时条件差错概率（conditional probability）由p表示。二进制对称信道的转移概率如下图所示

2.BSC信道的串联

在实际通信系统中，信号往往要通过几个环节的传输，或
多步处理，这些传输和处理都可以看成信道，他们串联而成一个串联信道，如下图所示，

二、Matlab对多个串联BSC信道容量的求解

1.从一个BSC信道开始

求解一个BSC信道的信道容量

代码如下（示例）：

clc
x=0:0.0001:1;
C=1+x.*log2(x)+(1-x).*log2(1-x);
plot(x,C)

得到图形

2.计算n个BSC信道串联的信道容量

代码如下（示例）：

syms n s;
% x=0:0.0001:1;
x=s;
X0=[1-x x
    x 1-x]; 
X1=X0^n;        %n次串联,即矩阵相乘n次
X1
C=1+X1(1,1).*log2(X1(1,1))+X1(1,2).*log2(X1(1,2)); %带入公式求解
simplify(C)

得到结果为
ans =
(log((1 - 2s)^n/2 + 1/2)((1 - 2s)^n/2 + 1/2))/log(2) - (log(1/2 - (1 - 2s)^n/2)((1 - 2s)^n/2 - 1/2))/log(2) + 1
化简后可以得到

三、串级BSC信道容量与单级BSC信道容量的联系

下面给出求解n等于1，2，3的BSC信道容量对比图，代码如下
我们可以利用上述所求的代码求解
n=2时的信道容量为

1+((1-x).^2+x.^2).*log2((1-x).^2+x.^2)+(1-((1-x).^2+x.^2)).*log2(1-((1-x).^2+x.^2))

n=3时的信道容量为

1+((1-(1-2*x).^3)/2).*log2((1-(1-2*x).^3)/2)+(1-((1-(1-2*x).^3)/2)).*log2(1-((1-(1-2*x).^3)/2))

因此我们可以求解其信道容量随变化的图像

1.代码

clc
x=0:0.0001:1;
C=1+x.*log2(x)+(1-x).*log2(1-x);
plot(x,C)
% simplify(C)
C2=1+((1-x).^2+x.^2).*log2((1-x).^2+x.^2)+(1-((1-x).^2+x.^2)).*log2(1-((1-x).^2+x.^2))
C3=1+((1-(1-2*x).^3)/2).*log2((1-(1-2*x).^3)/2)+(1-((1-(1-2*x).^3)/2)).*log2(1-((1-(1-2*x).^3)/2))
plot(x,C,x,C2,x,C3)

2.图像

四、总结

由图可得
串接的信道越多，其信道容量可能会越小，当串接信道数无限大时，信道容量就有可能趋于零。