MSK接收

---

前言

MSK解调 M序列相关 峰值检测 抽样判决 RS解码 等流程记录

一、MSK解调

MSK信号调制好后表达式如下所示：
$$S_{MSK}(t)=\cos \left({\omega }_{c}t+\frac{\pi a_k}{2T_s}t+{\phi }_k\right),(k-1)T_s\le t\le k{T}_s,k=1,2,3,\cdots \cdots$$

在移动或是室内的无线通信中,由多径引起的衰落将使相干解调的性能严重下降,误码率随之上升。此时,非相干解调更为合适。同时非相干解调结构简单,不需要载波同步，也不需要频偏、相偏的估计补偿，易于实现，后续采用非相干解调。

与I路本地载波相乘得（存在相偏）：
$$\begin{gathered} S_{MSK}(t)\cos \left({\omega }_{c}t+\psi \right)=\cos \left({\omega }_{c}t+\frac{\pi t}{2T_s}{a}_k+{\phi }_k\right)\cos \left({\omega }_{c}t+\psi \right) \\ =\frac{1}{2}\cos \left(2{\omega }_{c}t+\frac{\pi t}{2T_s}{a}_k+{\phi }_k+\psi \right)+\frac{1}{2}\cos \left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right) \end{gathered}$$

与Q路本地载波相乘得（存在相偏）：
$$\begin{gathered} S_{MSK}(t)\left(-\sin \left({\omega }_{c}t+\psi \right)\right)=\cos \left({\omega }_{c}t+\frac{\pi t}{2T_s}{a}_k+{\phi }_k\right)\left(-\sin \left({\omega }_{c}t+\psi \right)\right) \\ =\frac{1}{2}\sin \left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right)-\frac{1}{2}\sin \left(2{\omega }_{c}t+\frac{\pi t}{2T_s}{a}_k+{\phi }_k+\psi \right) \end{gathered}$$

抽取滤波，通过CIC抽取滤波器，滤除高频信号，同时对数据进行降采样：
$$\begin{gathered} Y_1=\frac{1}{2}\cos \left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right) \\ {Y}_2=\frac{1}{2}\sin \left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right) \end{gathered}$$

Y1和Y2构成的基带信号分别与f1和f2两个正交基函数在一个符号时间内作相关运算：
$$\begin{gathered} {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_1\right)dt \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_2\right)dt \end{gathered}$$

其中f1和f2为：

$$\begin{gathered} f_1=e^{j\left(\frac{\pi t}{2T_s}+\theta \right)} \\ {f}_2=e^{-j\left(\frac{\pi t}{2T_s}+\theta \right)} \end{gathered}$$

$$\begin{gathered} Y_1+j{Y}_2与f_1做相关,f_1中的\theta 是一个常量，将其归到基带信号的\psi 中，后续不显示了，没什么影响: \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_1\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right)}\cdot e^{j\frac{\pi t}{2T_s}}\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_s}\left(a_k+1\right)+{\phi }_k-\psi \right)}\right)dt \end{gathered}$$
$$\begin{gathered} Y_1+j{Y}_2与f_2做相相关,f_2中的\theta 是一个常量，将其归到基带信号的\psi 中，后续不显示了，没什么影响: \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_2\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_s}{a}_k+{\phi }_k-\psi \right)}\cdot e^{-j\frac{\pi t}{2T_s}}\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_s}\left(a_k-1\right)+{\phi }_k-\psi \right)}\right)dt \end{gathered}$$

---

$$\begin{gathered} 当a_k=1,Y_1+j{Y}_2与f_1相关得： \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_1\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_s}\left(a_k+1\right)+{\phi }_k-\psi \right)}\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{T_s}+{\phi }_k-\psi \right)}\right)dt \end{gathered}$$

$$\begin{gathered} 当a_k=1,Y_1+j{Y}_2与f_2相关得: \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_2\right)dt={\int }_0^{{\tau }_s}\left(e^{j\left(\frac{\pi t}{2{\tau }_s}\left(a_k-1\right)+{\phi }_k-\psi \right)}\right)dt={\int }_0^{T_s}\left(e^{j\left({\phi }_k-\psi \right)}\right)dt \end{gathered}$$

---

$$\begin{gathered} 当a_k=-1,Y_1+j{Y}_2与f_1相关得: \\ {\int }_0^{T_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_1\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{2T_S}\left(a_k+1\right)+{\phi }_k-\psi \right)}\right)dt={\int }_0^{{\tau }_s}\left(e^{j\left({\phi }_k-\psi \right)}\right)dt \end{gathered}$$

$$\begin{gathered} 当a_k=-1,Y_1+j{Y}_2与f_2相关得： \\ {\int }_0^{{\tau }_s}\left(\left(Y_1+j{Y}_2\right)\cdot f_2\right)dt={\int }_0^{{\tau }_s}\left(e^{j\left(\frac{\pi t}{2T_s}\left(a_k-1\right)+{\phi }_k-\psi \right)}\right)dt={\int }_0^{T_s}\left(e^{j\left(\frac{\pi t}{T_s}+{\phi }_k-\psi \right)}\right)dt \end{gathered}$$

————————————————————————————————————————————————————————
即最终都为求两个积分:

$$\begin{gathered} {\int }_0^{T_s}{e}^{j\left({\phi }_k-\psi \right)}dt=T_s\cdot e^{j\left({\phi }_k-\psi \right)} \\ {\int }_0^{T_s}{e}^{j\left(\frac{\pi t}{T_s}+{\phi }_k-\psi \right)}dt=e^{j\left({\phi }_k-\psi \right)}{\int }_0^{{\tau }_s}{e}^{j\frac{\pi t}{T_s}}dt=e^{j\left({\phi }_k-\psi \right)}\cdot \frac{T_s}{j\pi }\cdot {\int }_0^{{\tau }_s}{e}^{j\frac{\pi t}{T_s}}d\frac{j\pi }{T_s}t \\ ={e^{j\left({\phi }_k-\psi \right)}\cdot \frac{T_s}{j\pi }\cdot e^{j\frac{\pi t}{T_s}}|}_0^{T_s}=\frac{j2}{\pi }\cdot T_s\cdot e^{j\left({\phi }_k-\psi \right)} \end{gathered}$$

再进行一定的处理后，根据Y1和Y2可得

由于主要依据接收信号的包络信息进行符号判决，因而该解调算法具有较好的抗相差性能，对信号的适应性强，对载频差要求不严格。

二、突发捕获

在蓝线处抽样判决即可得到ak，但先需要通过m序列准确地检测到信号的起始时刻，才能进行后续处理。利用m序列良好的自相关特性，与本地码做相关运算

峰值检测，结果与门限相比较,即可作出判决（门限由阈值和前128个数据的平均值共同得到）

三、抽样判决

检测到m序列峰值后，在每一码元的结束时刻进行抽样判决，每组好 1 byte数据，标志位 rdy 有效（拉高），然后将8 bit 数据放入RAM中，为后续RS解码做准备。

突发捕获和抽样判决整体如下所示：

四、RS解码

当组完设定的255 byte 后，结束标志位有效，启动RS解码。

当接收数据未出错时，RS解码数据完全正确，RS解码时序如下，Latency为2606，时钟周期20ns，2606*20ns=52.12us

接收局部数据如下：

RS解码后局部数据如下：

当接收数据出现错误时（RS(255，191)，共191个数据，从191数据后为255-191个校验码，第190个数应该是190，现为142），RS解码时将数据纠正回来


当接收数据中校验数据出现错误时（RS(255，191)，共191个数据，从191数据后为255-191个校验码，第2个校验数应该是5，现为140），RS解码时也可以将数据纠正