欢迎光临散文网 会员登陆 & 注册

AWGN信道下的误码率分析-BPSK

2023-04-20 14:24 作者:乐吧的数学  | 我要投稿

录制的视频在:https://www.bilibili.com/video/BV1bs4y1d7zY/


这篇文章主要是翻译了一篇英文博客的,自己略加修改。

原文链接: https://www.gaussianwaves.com/2012/07/intuitive-derivation-of-performance-of-an-optimum-bpsk-receiver-in-awgn-channel/



信道是 AWGN 信道,没有衰落,则可以表示为:   r = x + w

其中 x 是 BPSK 调制, w 是高斯白噪声,符合均值为 0 ,方差为  $$\frac{N_0}{2}$$ 的 高斯分布。

因为 x 是 BPSK  调制,我们假定其信号能量为 Es (因为是 BPSK,所以 Eb = Es).


当传输的比特 是 0 时,则 S_0%20%3D%20-%5Csqrt%7BE_s%7D, 因为 w 是均值为 0 ,方差为  %5Cfrac%7BN_0%7D%7B2%7D 的 高斯分布,所以, r  是均值为 %20-%5Csqrt%7BE_s%7D ,方差为  %5Cfrac%7BN_0%7D%7B2%7D 的 高斯分布.

当传输的比特 是 1 时,则 S_1%20%3D%20%2B%5Csqrt%7BE_s%7D, 因为 w 是均值为 0 ,方差为  %5Cfrac%7BN_0%7D%7B2%7D 的 高斯分布,所以, r  是均值为 %2B%5Csqrt%7BE_s%7D ,方差为  %5Cfrac%7BN_0%7D%7B2%7D 的 高斯分布.

表示成数学公式为:

p(r%7C0_T)%20%3D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20e%5E%7B-%5Cfrac%7B(r-(-%5Csqrt%7BE_s%7D))%5E2%7D%7BN_0%7D%7D%20%20%20%5Ctag%20%7B1A%7D
p(r%7C1_T)%20%3D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20e%5E%7B-%5Cfrac%7B(r-(%2B%5Csqrt%7BE_s%7D))%5E2%7D%7BN_0%7D%7D%20%20%20%5Ctag%20%7B1B%7D
那么做最优判决的时候,如下图所示:


那么发生比特错误的情况为:

* 发送的是 0,但是被判决为 1
* 发送的是 1,但是被判决为 0

用数学公式表示为:

P(error%20)%20%3D%20P(%E5%88%A4%E5%86%B3%E4%B8%BA1%2C%20%E5%8F%91%E9%80%810)%20%2B%20%20P(%E5%88%A4%E5%86%B3%E4%B8%BA0%2C%20%E5%8F%91%E9%80%811)%20%20%5Ctag%202


即:

P(e)%20%3D%20P(1_D%2C%200_T)%20%2B%20P(0_D%2C1_T)%20%20%5Ctag%203


其中 D 表示 Decided,即被判决出来的; T 表示 transmit,即发送的。

用贝叶斯定理,公式 (3) 可以表示为:

P(e)%20%3D%20P(1_D%7C%200_T)P(0_T)%20%2B%20P(0_D%7C1_T)P(1_T)%20%20%5Ctag%204


从下图的示意,我们可以知道 P(1_D%7C%200_T) 和 P(0_D%7C1_T) 的数学表达式,实际上表示的就是下图中红色的区域的概率。




则:

P(1_D%7C%200_T)%20%3D%20%5Cint_0%5E%7B%2B%5Cinfty%7D%20p(r%7C0_T)%20dr%20%20%20%5Ctag%7B5A%7D


以及:

P(1_D%7C%200_T)%20%3D%20%5Cint_%7B%2B%5Cinfty%7D%5E0%20p(r%7C1_T)%20dr%20%20%20%5Ctag%7B5B%7D


把公式 (1A) 代入 (5A),把公式 (1B) 代入 (5B) 有:

P(1_D%7C%200_T)%20%3D%20%5Cint_0%5E%7B%2B%5Cinfty%7D%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20%20e%5E%7B-%5Cfrac%7B(r-(-%5Csqrt%7BE_s%7D))%5E2%7D%7BN_0%7D%7D%20dr%0A%3D%20%5Cint_0%5E%7B%2B%5Cinfty%7D%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20%20e%5E%7B-%5Cfrac%7B(r%2B%5Csqrt%7BE_s%7D)%5E2%7D%7BN_0%7D%7D%20dr%0A%5Ctag%7B6A%7D


以及
P(1_D%7C%200_T)%20%3D%20%5Cint_%7B%2B%5Cinfty%7D%5E0%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20e%5E%7B-%5Cfrac%7B(r-(%2B%5Csqrt%7BE_s%7D))%5E2%7D%7BN_0%7D%7D%20%20dr%20%20%0A%3D%20%5Cint_%7B%2B%5Cinfty%7D%5E0%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20e%5E%7B-%5Cfrac%7B(r-%5Csqrt%7BE_s%7D)%5E2%7D%7BN_0%7D%7D%20%20dr%0A%5Ctag%7B6B%7D

因为这两个分布具有对称性,所以,公式 (4) 可以推导为(把公式 6 A  和 B 代入公式 (4) ,并利用对称性 , 且发送的符号是等概率分布的):

%5Cbegin%7Baligned%7D%0AP(e)%20%26%3D%20%5Cfrac%7B1%7D%7B2%7DP(1_D%7C%200_T)%20%2B%20%5Cfrac%7B1%7D%7B2%7DP(0_D%7C1_T)%5C%5C%20%20%5C%5C%0A%26%3D%20P(1_D%7C%200_T)%5C%5C%20%20%5C%5C%0A%26%3D%20%20%5Cint_0%5E%7B%2B%5Cinfty%7D%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20%20e%5E%7B-%5Cfrac%7B(r%2B%5Csqrt%7BE_s%7D)%5E2%7D%7BN_0%7D%7D%20dr%0A%5Cend%7Baligned%7D%0A%5Ctag%207


我们继续来分析公式 (7) 中的积分公式,做积分变量的变量代换,令:

t%20%3D%20%5Cfrac%7Br%2B%5Csqrt%7BE_s%7D%7D%7B%5Csqrt%7BN_0%2F2%7D%7D


 则积分上下限就变成:

r%3D0%2C%20t%20%3D%5Cfrac%7B%5Csqrt%7BE_s%7D%7D%7B%5Csqrt%7BN_0%2F2%7D%7D%20%3D%20%5Csqrt%7B%5Cfrac%7BE_s%7D%7BN_0%2F2%7D%7D%20%20%5C%5C%0Ar%20%3D%2B%5Cinfty%2C%20t%3D%2B%5Cinfty


则公式(7) 的积分就变成:

%5Cint_0%5E%7B%2B%5Cinfty%7D%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cpi%20N_0%7D%7D%20%20e%5E%7B-%5Cfrac%7B(r%2B%5Csqrt%7BE_s%7D)%5E2%7D%7BN_0%7D%7D%20dr%0A%3D%5Cint_%7B%5Csqrt%7B%5Cfrac%7BE_s%7D%7BN_0%2F2%7D%7D%20%7D%5E%7B%2B%5Cinfty%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%7D%20%20%20e%5E%7B-%5Cfrac%7Bt%5E2%7D%7B2%7D%7D%20dt%20%20%5Ctag%208


这里就刚好是推导出来了 Q 函数的模样, 均值为 0 ,方差为 1 的标准正态分布,从 x 开始一直到无穷,计算其概率的大小,这是一个关于 x 的函数,即 Q(x),定义如下:

Q(x)%20%3D%5Cint_x%5E%7B%2B%5Cinfty%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%7D%20%20%20e%5E%7B-%5Cfrac%7Bt%5E2%7D%7B2%7D%7D%20dt%20%20%5Ctag%209

把公式 (8) 和 (9) 逐级带回到 (7) 有:

P(e)%20%3D%20Q%5Cleft%20(%5Csqrt%7B%5Cfrac%7BE_s%7D%7BN_0%2F2%7D%7D%20%5Cright%20)


画上图的 Python 代码


AWGN信道下的误码率分析-BPSK的评论 (共 条)

分享到微博请遵守国家法律