1565 words
8 minutes
QAM在AWGN信道中的理论BER计算

引言#

正交幅度调制(QAM)作为现代数字通信的核心技术,通过幅度和相位联合调制实现了高频谱效率传输。从Wi-Fi到5G,从光纤通信到卫星链路,QAM技术始终扮演着关键角色。误码率(BER)作为衡量通信系统可靠性的核心指标,其理论计算方法对系统设计具有重要指导意义。本文将系统性推导4QAM、16QAM的BER表达式。

4-QAM#

1. 信号能量#

在一个4-QAM的星座图中,以A点举例。

SA=E+jES_A = \sqrt{E} + j\sqrt{E}Es=(E)2+(E)2=2EE_s = (\sqrt{E})^2 + (\sqrt{E})^2 = 2E

其中,E\sqrt{E}IIQQ 两个维度的各自的幅度,EE 是每个维度的能量,每个符号在两个幅度上的能量之和是 2E2E

因为是4-QAM,每个符号携带的比特数为 log2(M)=2\log_2(M) = 2,因此:

Eb=Eslog2(M)=Es2=EE_b = \frac{E_s}{\log_2(M)} = \frac{E_s}{2} = E

image.png

2. AWGN信道#

在AWGN信道中,接收到的信号 rr 可以表示为:

r=s+nr = s + n

其中 ss 是发送的信号,nn 是噪声,为复高斯随机变量,其实部 nIn_I 与虚部 nQn_Q 相互独立,每个正交的分量都服从均值为0,方差 σ2=N02\sigma^2 = \frac{N_0}{2} 的分布,其中整个复信号的方差为 N0N_0

3. 判决#

在接收端,采用最小距离判决规则,即选择与接收信号 rr(也就是A点)欧氏距离最小的符号点。由于4-QAM符号点是对称的,判决区域可以用坐标轴划分:

  • II 分量:如果 rI>0r_I > 0,判定 sI=Es_I = \sqrt{E},否则 sI=Es_I = -\sqrt{E}
  • QQ 分量:如果 rQ>0r_Q > 0,判定 sQ=Es_Q = \sqrt{E},否则 sQ=Es_Q = -\sqrt{E}

IIQQ 的判决是相互独立的。

4. 计算BER#

噪声 nIn_InQn_Q 服从高斯分布:

nIN(0,N02),nQN(0,N02)n_I \sim \mathcal{N} \left(0, \frac{N_0}{2}\right), \quad n_Q \sim \mathcal{N} \left(0, \frac{N_0}{2}\right)

Q函数的定义为:

Q(x)=12πxet22dtQ(x) = \frac{1}{\sqrt{2 \pi}} \int_{x}^{\infty} e^{- \frac{t^2}{2}}dt

对于随机变量 XN(μ,σ2)X \sim \mathcal{N} (\mu , \sigma^2),有:

P(X<a)=Q(μaσ)P(X < a) = Q\left(\frac{\mu - a}{\sigma}\right)P(X>a)=Q(aμσ)P(X > a) = Q\left(\frac{a - \mu}{\sigma}\right)

以A点举例,s=(E,E)s = (\sqrt{E}, \sqrt{E})。在接收端,从I分量的角度,计算A落在错误的区域 rI<0r_I < 0 的概率:

r=E+nIr = \sqrt{E} + n_IP(rI<0)=P(E+nI<0)=P(nI<E)=Q(Eσ)P(r_I < 0) = P(\sqrt{E} + n_I < 0) = P(n_I < - \sqrt{E}) = Q\left(\frac{\sqrt{E}}{\sigma}\right)

则落在I分量上正确的概率为 1Q(Eσ)1 - Q(\frac{\sqrt{E}}{\sigma}),同时在两个分量上正确的概率为:

[1Q(Eσ)]2\left[1 - Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2

在两个分量上出错的概率(符号错误率)为:

Pser=1[1Q(Eσ)]2=2Q(Eσ)Q(Eσ)2P_{ser} = 1 - \left[1 - Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2 = 2Q\left(\frac{\sqrt{E}}{\sigma}\right) - Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

代入 σ=N0/2\sigma = \sqrt{N_0/2}E=EbE = E_b

Pser=2Q(EsN0)Q(EsN0)2P_{ser} = 2Q\left(\sqrt{\frac{E_s}{N_0}}\right) - Q\left(\sqrt{\frac{E_s}{N_0}}\right)^2

EbN01\frac{E_b}{N_0} \gg 1 时,可以忽略掉高阶项 Q()2Q(\cdot)^2

Pser2Q(EsN0)P_{ser} \approx 2Q\left(\sqrt{\frac{E_s}{N_0}}\right)BERPserBits per Symbol=Q(EsN0)=Q(2EbN0)BER \approx \frac{P_{ser}}{\text{Bits per Symbol}} = Q\left(\sqrt{\frac{E_s}{N_0}}\right) = Q\left(\sqrt{\frac{2E_b}{N_0}}\right)

16-QAM#

image.png

1. 信号能量#

在一个16-QAM的星座图中,有三种类型的点:角点,边缘点,内点。

其中:

  • 角点:J、G、P、M
  • 边缘点:I、H、E、F、O、N、K、L
  • 内点:A、B、C、D

本文的角点、边缘点、内点分别使用 P:(3E,3E)P:(3 \sqrt{E}, -3 \sqrt{E})H:(E,3E)H:(\sqrt{E}, 3 \sqrt{E})A:(E,E)A:(\sqrt{E}, \sqrt{E}) 来举例。

SA=E+jES_A = \sqrt{E} + j\sqrt{E}SH=E+j3ES_H = \sqrt{E} + j3\sqrt{E}SP=3Ej3ES_P = 3\sqrt{E} - j3\sqrt{E}

其中,E\sqrt{E}IIQQ 两个维度的各自的幅度,EE 是每个维度的能量。

因为是16-QAM,每个符号携带的比特数为 log2(M)=4\log_2(M) = 4,因此:

Eb=Eslog2(M)=Es4E_b = \frac{E_s}{\log_2(M)} = \frac{E_s}{4}

其中16个点的平均符号能量 EsE_s 为:

Es=18E4+10E8+2E416=10EE_s = \frac{18E \cdot 4 + 10E \cdot 8 + 2E \cdot 4}{16} = 10EE=Es10=4Eb10=25EbE = \frac{E_s}{10} = \frac{4E_b}{10} = \frac{2}{5} E_b

2. AWGN信道#

在AWGN信道中,接收到的信号 rr 可以表示为:

r=s+nr = s + n

其中 ss 是发送的信号,nn 是噪声,为复高斯随机变量,其实部 nIn_I 与虚部 nQn_Q 相互独立,每个正交的分量都服从均值为0,方差 σ2=N02\sigma^2 = \frac{N_0}{2} 的分布,其中整个复信号的方差为 N0N_0

3. 判决#

3.1 内点#

以A点为例,需要满足在 I:(0,2E)I:(0, 2\sqrt{E})Q:(0,2E)Q:(0, 2\sqrt{E}) 的范围内,即判定为正确的信号。

在接收端,从I分量的角度,计算A落在错误的区域 rI<0r_I < 0rI>2Er_I > 2\sqrt{E} 的概率:

r=E+nIr = \sqrt{E} + n_IP(rI<0)=P(E+nI<0)=P(nI<E)=Q(Eσ)P(r_I < 0) = P(\sqrt{E} + n_I < 0) = P(n_I < - \sqrt{E}) = Q\left(\frac{\sqrt{E}}{\sigma}\right)P(rI>2E)=P(E+nI>2E)=P(nI>E)=Q(Eσ)P(r_I > 2\sqrt{E}) = P(\sqrt{E} + n_I > 2\sqrt{E}) = P(n_I > \sqrt{E}) = Q\left(\frac{\sqrt{E}}{\sigma}\right)

在I分量上出错的概率为:

Q(Eσ)+Q(Eσ)=2Q(Eσ)Q\left(\frac{\sqrt{E}}{\sigma}\right) + Q\left(\frac{\sqrt{E}}{\sigma}\right) = 2Q\left(\frac{\sqrt{E}}{\sigma}\right)

则落在I分量上正确的概率为 12Q(Eσ)1 - 2Q(\frac{\sqrt{E}}{\sigma})。在这个例子中,I分量和Q分量相互独立,概率相同,同时在两个分量上正确的概率为:

[12Q(Eσ)]2\left[1 - 2Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2

在两个分量上出错的概率为:

Pser1=1[12Q(Eσ)]2=4Q(Eσ)4Q(Eσ)2P_{ser1} = 1 - \left[1 - 2Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2 = 4Q\left(\frac{\sqrt{E}}{\sigma}\right) - 4Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

3.2 边缘点#

以H点为例,需要满足在 I:(0,2E)I:(0, 2\sqrt{E})Q:(2E,+)Q:(2\sqrt{E}, +\infty) 的范围内,即判定为正确的信号。

在接收端,从I分量的角度,计算H落在错误的区域 rI<0r_I < 0rI>2Er_I > 2\sqrt{E} 的概率:

同内点的情况,在I分量上出错的概率为 2Q(Eσ)2Q(\frac{\sqrt{E}}{\sigma})。 则落在I分量上正确的概率为 12Q(Eσ)1 - 2Q(\frac{\sqrt{E}}{\sigma})

从Q分量的角度看,计算H落在错误的区域 rQ<2Er_Q < 2\sqrt{E} 的概率:

r=3E+nQr = 3\sqrt{E} + n_QP(rQ<2E)=P(3E+nQ<2E)=P(nQ<E)=Q(Eσ)P(r_Q < 2\sqrt{E}) = P(3\sqrt{E} + n_Q < 2\sqrt{E}) = P(n_Q < - \sqrt{E}) = Q\left(\frac{\sqrt{E}}{\sigma}\right)

则落在Q分量上正确的概率为 1Q(Eσ)1 - Q(\frac{\sqrt{E}}{\sigma})

同时在两个分量上正确的概率为:

Pcorrect=[12Q(Eσ)][1Q(Eσ)]P_{correct} = \left[1 - 2Q\left(\frac{\sqrt{E}}{\sigma}\right)\right] \cdot \left[1 - Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]

在两个分量上出错的概率为:

Pser2=1Pcorrect=3Q(Eσ)2Q(Eσ)2P_{ser2} = 1 - P_{correct} = 3Q\left(\frac{\sqrt{E}}{\sigma}\right) - 2Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

3.3 角点#

以P点为例,需要满足在 I:(2E,+)I:(2\sqrt{E}, +\infty)Q:(2E,)Q:(- 2\sqrt{E}, -\infty) 的范围内,即判定为正确的信号。

在接收端,从I分量的角度,计算P落在错误的区域 rI<2Er_I < 2\sqrt{E} 的概率:

r=3E+nIr = 3\sqrt{E} + n_IP(rI<2E)=P(3E+nI<2E)=P(nI<E)=Q(Eσ)P(r_I < 2\sqrt{E}) = P(3\sqrt{E} + n_I < 2\sqrt{E}) = P(n_I < - \sqrt{E}) = Q\left(\frac{\sqrt{E}}{\sigma}\right)

则落在I分量上正确的概率为 1Q(Eσ)1 - Q(\frac{\sqrt{E}}{\sigma})。在这个例子中,I分量和Q分量相互独立,概率相同,同时在两个分量上正确的概率为:

[1Q(Eσ)]2\left[1 - Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2

在两个分量上出错的概率为:

Pser3=1[1Q(Eσ)]2=2Q(Eσ)Q(Eσ)2P_{ser3} = 1 - \left[1 - Q\left(\frac{\sqrt{E}}{\sigma}\right)\right]^2 = 2Q\left(\frac{\sqrt{E}}{\sigma}\right) - Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

4. 计算BER#

三类符号点出错的概率分别为:

  • 内点:Pser1=4Q(Eσ)4Q(Eσ)2P_{ser1} = 4Q\left(\frac{\sqrt{E}}{\sigma}\right) - 4Q\left(\frac{\sqrt{E}}{\sigma}\right)^2
  • 边缘点:Pser2=3Q(Eσ)2Q(Eσ)2P_{ser2} = 3Q\left(\frac{\sqrt{E}}{\sigma}\right) - 2Q\left(\frac{\sqrt{E}}{\sigma}\right)^2
  • 角点:Pser3=2Q(Eσ)Q(Eσ)2P_{ser3} = 2Q\left(\frac{\sqrt{E}}{\sigma}\right) - Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

总的符号错误概率为:

Pser=4Pser1+8Pser2+4Pser316=3Q(Eσ)94Q(Eσ)2P_{ser} = \frac{4 \cdot P_{ser1} + 8 \cdot P_{ser2} + 4 \cdot P_{ser3}}{16} = 3 Q\left(\frac{\sqrt{E}}{\sigma}\right) - \frac{9}{4} Q\left(\frac{\sqrt{E}}{\sigma}\right)^2

根据 σ=N02\sigma = \sqrt{\frac{N_0}{2}}E=25EbE = \frac{2}{5}E_b

Eσ=EN0/2=2EN0=225EbN0=45EbN0\frac{\sqrt{E}}{\sigma} = \sqrt{\frac{E}{N_0/2}} = \sqrt{\frac{2E}{N_0}} = \sqrt{\frac{2 \cdot \frac{2}{5}E_b}{N_0}} = \sqrt{\frac{4}{5} \frac{E_b}{N_0}}

代入得:

Pser=3Q(45EbN0)94Q(45EbN0)2P_{ser} = 3 Q\left(\sqrt{\frac{4}{5} \cdot \frac{E_b}{N_0}}\right) - \frac{9}{4} Q\left(\sqrt{\frac{4}{5} \cdot \frac{E_b}{N_0}}\right)^2

EbN01\frac{E_b}{N_0} \gg 1 时,可以忽略掉高阶项:

Pser3Q(45EbN0)P_{ser} \approx 3 Q\left(\sqrt{\frac{4}{5} \cdot \frac{E_b}{N_0}}\right)BERPserBits per Symbol=34Q(45EbN0)BER \approx \frac{P_{ser}}{\text{Bits per Symbol}} = \frac{3}{4} Q\left(\sqrt{\frac{4}{5} \cdot \frac{E_b}{N_0}}\right)
QAM在AWGN信道中的理论BER计算
https://fuwari.vercel.app/posts/mqam/
Author
Akatsuki Sky
Published at
2025-03-02
License
CC BY-NC-SA 4.0