砸饭碗系列一、BPSK高斯、瑞利、莱斯信道误码率对比（手把手推导瑞利信道理论误码率公式）

科研Beatles

已于 2024-05-13 08:13:27 修改

阅读量3.2k

点赞数 12

文章标签：信息与通信

于 2024-05-11 20:40:56 首次发布

本文链接：https://blog.csdn.net/weixin_43702241/article/details/138728909

版权

本篇文章整合了所有可以用的代码，需要代码请添加微信公众号，有问题请与作者联系qq:286820790

一、砸csdn上一些过于贵的代码人的饭碗

本文是做了BPSK高斯信道理论以及仿真的误码率曲线、瑞利信道理论以及仿真的误码率曲线、莱斯信道仿真的误码率曲线。（注意没有莱斯信道理论的误码率曲线公式，为什么呢，因为我没做出来！）

做这个的原因是，想学习瑞利跟莱斯信道，结果发现csdn 上拥有这个曲线的博主，居然要收80多块钱！于是果断决定砸他饭碗！嗯哼，搞了五天左右几乎下载csdn上所有代码整合出了，这三个信道的误码率仿真图。但是美中不足的是缺少莱斯信道的理论误码率曲线。这个日后如果有机会搞清楚后，会补上。

二、瑞利信道

对于 $y=h x+n$ ,h为信道衰落。y为接收数据，x为发射数据，n服从均值0方差 $\sigma^2$ 的高斯分布。

高斯信道下BPSK理论误码率公式,其中 $E_bN_o$ 是信噪比dB形式，也就是我们设置的0:1:20;

$P_b=\frac{1}{2} \operatorname{erfc}\left(\sqrt{{snr}}\right)=Q(\sqrt{2snr})$

$snr=10^{E_bN_o/10}$

当信号经过瑞利衰落信道后，接收信噪比为 $\gamma=|h|^{2}snr$ ，因此，对于给定ℎ，理论误比特率为

$P_b=\frac{1}{2} \operatorname{erfc}\left(\sqrt{{\gamma}}\right)=Q(\sqrt{2\gamma})$

由于 $|h|^{2}$ 服从自由度为2的卡方分布，因此 $\gamma$ 也服从卡方分布只不过多个倍数snr而已，卡方分布这个大家可以自己去查查，我们直接给出公式

$f(x)=\left\{\begin{array}{c}\frac{1}{2 \sigma^2} e^{-\frac{x}{2 \sigma^2}}, x>0 \\ \text { b, otherwise }\end{array}\right.$

$p(\gamma)=\frac{1}{2 snr} e^{\frac{-\gamma}{2snr}}, \gamma \geq 0$

但是对于h来说他是个随机变量，因此必须以积分的形式求出误码率，积分公式如下

公式一、 $P_b=\int_0^{\infty} \frac{1}{2} \operatorname{erfc}(\sqrt{\gamma}) p(\gamma) d \gamma$

通过积分结果可以得到

公式二、 $P_b=\frac{1}{2}\left(1-\sqrt{\frac{snr}{snr+1}}\right)$

但是很多人疑问如何从公式一的积分形式得到公式二、最终的结果的，这里我给大家推导一下，其实就是用了高数中分部积分的原理而已。积分过程如下：

$P_b=\int_0^{\infty} \frac{1}{2} \operatorname{erfc}(\sqrt{\gamma}) p(\gamma) d \gamma \\=\int_0^{\infty}[ \frac{1}{\sqrt{\pi}}\int_\gamma^{\infty}\exp(-t^2)d\gamma ]*\frac{1}{2 snr} e^{\frac{-\gamma}{2snr}}d \gamma\\=\frac{-1}{\sqrt{\pi}}\int_0^{\infty}[ \int_\gamma^{\infty}\exp(-t^2)d\gamma ]*de^{\frac{-\gamma}{2snr}} \\=\frac{-1}{\sqrt{\pi}}*([ \int_{\sqrt{\gamma}}^{\infty}\exp(-t^2)d\gamma *e^{\frac{-\gamma}{2snr}}]_0^{\infty} -\int_0^{\infty}e^{-\gamma}*e^{\frac{-\gamma}{2snr}})d\gamma$