信道容量的计算

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

§4.2信道容量的计算

这里,我们介绍一般离散信道的信道容量计算方法,根据信道容量的定义,就是在固定信道的条件下,对所有可能的输入概率分布)(x P 求平均互信息的极大值。前面已知()Y X I ;是输入概率分布的上凸函数,所以极大值一定存在。而);(Y X I 是r 个变量

)}(),(),({21r x p x p x p 的多元函数。并且满足1)(1

=∑=r

i i x p 。所以可用拉格朗日乘子法来

计算这个条件极值。引入一个函数:∑-=i

i

x p Y X I )();(λ

φ解方程组

0)(]

)();([)

(=∑∂-∂∂∂i i

i i x p x p Y X I x p λ

φ

1)(=∑i

i

x p (4.2.1)

可以先解出达到极值的概率分布和拉格朗日乘子λ的值,然后在解出信道容量C 。因为

)

()

(log

)()();(11

i i i i i r

i s

j i y p x y Q x y Q x p Y X I ∑∑===

而)()()(1

i i

r

i i

i x y

Q x p y p ∑==

,所以

e e y p y p i i i i i y p x y Q i x p i x p l o g l o g

))(ln ()(log )

()()()

(==∂∂∂∂。

解(4.2.1)式有

0log )()()()()()(log )(111=--∑∑∑===λe y p x y Q x y Q x p y p x y Q x y Q i

i i i i r i s j i i i i s

j i i (对r i ,,2,1 =都成立) 又因为

)()()(1

j k k

r

k k

y p x y

Q x p =∑=

r

i x y Q s

j i j

,,2,1,1)(1

==∑=

所以(4.2.1)式方程组可以转化为

),,2,1(log )

()(log

)(1r i e y p x y Q x y Q j i j s

j i j =+=∑=λ

1)(1

=∑=r

i i

x p

假设使得平均互信息);(Y X I 达到极值的输入概率分布},,{21r p p p 这样有

e y p x y Q x y Q x p j i j i j r i s

j i log )

()(log

)()(11

+=∑∑==λ

从而上式左边即为信道容量,得 e C log +=λ 现在令

)

()(log

)();(1

j i j s

j i j i y p x y Q x y Q Y x I ∑==

式中,);(Y x I i 是输出端接收到Y 后获得关于i x X =的信息量,即是信源符号i x X =对输出端Y 平均提供的互信息。

一般来讲,);(Y x I i 值与i x 有关。根据(4.2.2)式和(4.2.3)式, C Y x I i =);( ),,2,1(r i = 所以对于一般离散信道有如下定理。

定理 4.2.1 一般离散信道的平均互信息);(Y X I 达到极大值(即等于信道容量)的充要条件是输入概率分布)}(,),({1n x p x p 满足

)(a C Y x I =);(1 对所有的0)(,≠i i x p x )(b C Y x I i ≤);( 对所有的0)(,=i i x p x 这时C 就是所求的信道容量。

对于离散信道来说,其实信道容量还有一个解法:迭代解法。

定理4.2.2 设信道的向前转移概率矩阵为J K i j x y Q Q ⨯=))((,0

P 是任给的输入字母的一个初始概率分布,其所有分量0)(0≠k x P 。按照下式不断地对概率分布进行迭代,更新:

∑=+=K

i r

i

i

r

r k k r

k r P x P P x P x P 1

1)()()

()

()(ββ

其中 r P P k r

k Y x X I P ===)];(exp[)(β

()()()⎪⎪⎭

⎪⎪⎬

⎫⎪⎪⎩⎪⎪⎨⎧=∑∑==J j K i i j r i j k j x y Q P x y Q x y Q 11log exp

由此所得的(

)

Q P I r

,序列收敛于信道容量C 。

我们还可以将上述过程写成算法以便编制程序实现(如图4.2.1) })()(log{

1

∑==K

k k

k

L P x P I β

)}(log{P x ma I k k

U β=

})()(log{1

∑==K

k k k L P x P I β

)}(log{P x ma I k k

U β=

图4.2.1 信道容量的迭代算法

对于一些特殊的离散信道,我们有方便的方法计算其信道容量。

定义4.2.1 设X 和Y 分别表示输入信源与输出信源,则我们称()

Y X H 为损失熵,

()X Y H 为信道噪声熵。

如果信道的损失熵()

0=Y X H ,则次信道容量为

开始 P

P →0

)(P k β

)(P I L )(P I U ε<-L U I I L I C =∑-=1)()()

()

()(P x P P x P x P ββ

相关文档
最新文档