单径瑞利信道的BPSK相干解调的理论误码率推导与MATLAB分析(1)

一、问题

       最近我继续攻克老师布置的作业，也是通过博客来记录一下自己的解决问题的过程，问题如下：

(1) 请推导出单径瑞利信道中的BPSK相干解调的理论误码率性能，并画出比特信噪比(EbN0)与误码率(SER)的关系曲线。

通过阅读与理解CSDN、相关书籍的资料，我基本理解了题目的理论误码率推导过程，下面是对于单径瑞利信道，采用BPSK调制以及相干解调的理论误比特率推导过程。

二、理论误码率推导过程

1）表示发送信号
       我们设二进制相移键控BPSK中发射的两个信号波形分别表示为$s_1\left(t\right)=g\left(t\right),s_2\left(t\right)=-g\left(t\right)$。其中$g\left(t\right)$是在一个符号间隔$T$内非零的任意实信号脉冲，脉冲能量为${\xi }_g$。这两个发射信号是双极性信号，可由能量表示为$s_1\left(t\right)=\sqrt{{\xi }_g},s_2\left(t\right)=-\sqrt{{\xi }_g}$。两个信号点如图2.1所示。

图2.1: BPSK的两个信号点

2）表示接收信号
       单径瑞利信道可视为平坦慢衰落信道，平坦衰落的信道会使发送信号$s\left(t\right)$发生乘性失真，而慢衰落信道则可将乘性失真过程在至少一个符号间隔 $T$内看作为一个常数的过程。因此我们假设两个信号是等概发送，若发送信号为$s_1\left(t\right)=g\left(t\right)$，则在一个符号间隔$T$内的等效低通接收信号为
$$r_1\left(t\right)=\alpha e^{-j\varphi }{s}_1\left(t\right)+z\left(t\right)\left(0\le t\le T\right)\text{\hspace{1em}(2.1)}$$

其中，$\alpha$为信号衰减，$\varphi$是信号相移，$z\left(t\right)$为复高斯白噪声过程。
       由于单径瑞利信道为慢衰落信道，其信道衰落足够慢，这样会使得相移$\varphi$能够从接收信号中无误差地估计出来。因此可以实现接收信号的理想相干检测，即相干解调。于是对于BPSK调制，可用一个匹配滤波器来处理接收信号。

3）计算固定衰减条件的条件误码率
       对于慢衰落信道，即信号衰减$\alpha$固定，不随时间发生改变。那么由匹配滤波器的解调器得到的接收信号为
$$r=\alpha s_1+n=\alpha \sqrt{{\xi }_g}+n\text{\hspace{1em}(2.2)}$$

其中$n$表示均值为0，噪声方差为${\sigma }_n^2=\frac{N_0}{2}$的加性高斯噪声分量。我们可以根据判决变量来确定误码率。将接收信号$r$与阈值0比较。若$r>0$则判决接收信号为$s_1\left(t\right)$。若$r<0$则判决为$s_2\left(t\right)$。因此接收信号$r$判决为$s_1\left(t\right),s_2\left(t\right)$的概率密度函数分别为
$$p\left(r|s_1\right)=-\frac{1}{\sqrt{\pi N_0}}{e}^{-{\left(r-\alpha \sqrt{{\xi }_g}\right)}^2/N_0}\text{\hspace{1em}(2.3)}$$

$$p\left(r|s_2\right)=-\frac{1}{\sqrt{\pi N_0}}{e}^{-{\left(r+\alpha \sqrt{{\xi }_g}\right)}^2/N_0}\text{\hspace{1em}(2.4)}$$

两个信号的判决的概率密度函数如图2.2所示。

图2.2: 判决的概率密度函数
       在发送信号为 $s_1\left(t\right)=g\left(t\right)$的条件下，错误概率为 $r<0$的概率，即
$$p\left(e|s_1\right)={\int }_{-\infty }^{0}p\left(r|s_1\right)dr=-\frac{1}{\sqrt{\pi N_0}}{\int }_{-\infty }^0\exp \left[-{\left(r-\alpha \sqrt{{\xi }_g}\right)}^2/N_0\right]dr=Q\left(\sqrt{\frac{2{\alpha }^2{\xi }_g}{N_0}}\right)$$

其中$Q\left(x\right)$为Q函数，是用来表示高斯分布的概率密度函数尾部曲线下的面积，定义为$Q\left(x\right)=\frac{1}{\sqrt{2\pi }}{\int }_x^{\infty }{e}^{-t^2/2}dt,x\ge 0$。
       同样假设发送信号为$s_2\left(t\right)=-g\left(t\right),r=-\alpha \sqrt{{\xi }_g}+n$，其$r>0$的错误概率也为$p\left(e|s_2\right)=Q\left(\sqrt{\frac{2{\alpha }^2{\xi }_g}{N_0}}\right)$。因为$s_1\left(t\right),s_2\left(t\right)$是等概发送的，由全概率公式可得平均错误概率为
$$P_2\left({\gamma }_b\right)=p\left(s_1\right)p\left(e|s_1\right)+p\left(s_2\right)p\left(e|s_2\right)=\frac{1}{2}p\left(e|s_1\right)+\frac{1}{2}p\left(e|s_2\right)=Q\left(\sqrt{\frac{2{\alpha }^2{\xi }_g}{N_0}}\right)=Q\left(\sqrt{2{\gamma }_b}\right)$$

其中${\gamma }_b=\frac{{\alpha }^2{\xi }_g}{N_0}$。上式$P_2\left({\gamma }_b\right)$即为条件误比特率。

4）计算衰减变化的误比特率
       上式得到的条件误比特率$P_2\left({\gamma }_b\right)$的条件是信号衰减$\alpha$固定不变。为了得到信号衰减$\alpha$随机变化的误码率，需要将条件误比特率$P_2\left({\gamma }_b\right)$对${\gamma }_b$的概率密度函数求平均，即衰减变化的误比特率为
$$P_2={\int }_0^{\infty }{P}_2\left({\gamma }_b\right)p\left({\gamma }_b\right)d{\gamma }_b\text{\hspace{1em}(2.5)}$$

其中$p\left({\gamma }_b\right)$是$\alpha$为随机变量${\gamma }_b$时的概率密度函数。因为作为随机变量的信号衰减$\alpha$服从瑞利分布，所以${\alpha }^2$是自由度为2的${\chi }^2$分布，即${\alpha }^2\backsim {\chi }^2\left(2\right)$。因此，${\gamma }_b$也是${\chi }^2$分布，则其概率密度函数为
$$p\left({\gamma }_b\right)=\frac{1}{{\bar{\gamma }}_b}{e}^{-{\gamma }_b/{\bar{\gamma }}_b}\text{\hspace{1em}(2.6)}$$

其中${\bar{\gamma }}_b$为平均信噪比，定义为${\bar{\gamma }}_b=\frac{{\xi }_g}{N_0}E\left({\alpha }^2\right)=2\frac{{\xi }_g}{N_0}$。其中$E\left({\alpha }^2\right)$为${\alpha }^2$的均值，为2。
       因此，将式(2.6)的${\gamma }_b$的概率密度函数$p\left({\gamma }_b\right)$代入衰减变化的误比特率公式(2.5)中，计算积分，即可得到误比特率为
$$P_2=\frac{1}{2}\left(1-\sqrt{\frac{{\bar{\gamma }}_b}{1+{\bar{\gamma }}_b}}\right)=\frac{1}{2}\left(1-\sqrt{\frac{2\frac{{\xi }_g}{N_0}}{1+2\frac{{\xi }_g}{N_0}}}\right)\text{\hspace{1em}(2.7)}$$

定义比特信噪比为$BSNR=2\frac{{\xi }_g}{N_0}$。因此式(2.7)的误比特率可进一步用比特信噪比表示为
$$P_2=\frac{1}{2}\left(1-\sqrt{\frac{BSNR}{1+BSNR}}\right)\text{\hspace{1em}(2.8)}$$

式(2.8)即为在单径瑞利信道中使用BPSK相干解调的理论误比特率。需要注意的是，在推导过程中，我们是假定在慢衰落信道中得到的相移估计值是无噪的，而这种理想条件在实际中并不成立。

三、使用MATLAB绘制理论误码率曲线

       首先MATLAB软件自带有绘制误比特率曲线的工具，即bertool工具。我们先使用bertool工具来绘制理论误比特率曲线，设置比特信噪比(Eb/N0)范围为-30到30dB，信道选用单径瑞利信道，调制类型选用BPSK，并使用相干解调，设置如图3.1所示：

图3.1: bertool工具的参数设置

       使用bertool工具绘制出的理论误比特率曲线如图3.2所示，可以得出随着比特信噪比的增加，使用BPSK相干解调的理论误比特率逐渐减小。在比特信噪比为10dB时误比特率约为0.02327。

图3.2: 使用bertool工具绘制出的理论误码率曲线

       接着我根据推导出的单径瑞利信道中的使用BPSK相干解调的理论误比特率公式，即式(2.8)，使用MATLAB编写代码来实现理论误比特率曲线的绘制。使用推导出的公式得出的理论误比特率曲线如图3.3所示。

图3.3: 使用推导出的公式得出的理论误码率曲线

根据图形也可以得出随着比特信噪比的增加，理论误码率在逐渐减小，当比特信噪比为10dB时误码率约为0.02327，该数值基本与使用bertool工具绘制出的理论误码率相同。

四、总结

       谢谢大家看到这里！这次的作业时隔将近一个月才写好，期间主要在搞其他任务，哈哈不过总算是写完一半了，这次主要是通过看经典书籍来消化知识的，对于第二问的信道估计也快解决了，应该很快就能做出来，还请读者们多多指教，谢谢！世上无难事，只要肯登攀。加油！

这里是本文的参考文献：
[1]:数字通信（第四版）
[2]:单径瑞利信道中的BPSK相干解调的（理论）误码率性能
[3]:MATLAB的bertool绘制误码率理论值与仿真值对比曲线
[4]:dB的换算