我目前正在处理一个棘手的通信链，我想在 Matlab 中进行模拟。

该链类似于经典的 BPSK 调制器/解调器链，但 BPSK 调制信号已使用 PM（相位调制）调制方法与其他载波合并。

因此，在 PM 调制步骤之后，信号可以表示为（$\mathrm{h}$是调制指数和$d(t)$是 NRZ 数据信号）：

$$\cos\left( \omega_c t + \mathrm{h} d(t) \cos(\omega_{sc} t + \varphi_{sc}) + \varphi_c \right)$$

我的目标是找到曲线$BER = f\left(\frac{E_b}{N_0}\right)$（这应该符合理论上的 BPSK BER 曲线）。

好吧，我曾经在没有 PM 步骤的情况下使用这种链。为了模拟它，我通常只使用理论噪声功率$P_n = \frac{\sigma_c^2}{\frac{E_b}{N_0} * \mathrm{BPS}}$创建标准偏差的高斯噪声$\sqrt{\frac{Pn}{2}}$. 但我目前的情况有一些问题：我在 Matlab 中制作了这个链，但我无法模拟它，因为当我的信号介于两者之间时，我不知道如何找到我必须应用于高斯噪声的标准偏差PM 步骤。

实际上，我不知道如何找到它...由于 PM 调制/解调步骤的非线性，我只是精神错乱，我不知道该怎么做$P_n$...

我已经在 PM 调制步骤之后计算了理论信号自相关/PSD，它使用第一类 Bessel 函数（类似于本文的表达式（9） ）。但老实说，我不知道该怎么办...