近在做MSK的正交调制解调仿真，有一些心得体会，也踩了一些坑，再次记录一下。 我的目标是，自己写一个正交调制解调，其性能和matlab自带的mskmod和mskdemod函数性能相近，功能相同，输出相同。

发现了一些mskmod函数的一些秘密，在此记录一下

mskmod - Minimum shift keying modulation This MATLAB function outputs the complex envelope y of the modulation of the message signal x using differentially encoded minimum shift keying (MSK) modulation.

从help mskmod函数我们可以得到mskmod函数的使用方法

最基础的是直接输入需要调制的信息序列和输入上采样倍数

信息序列：输入需要是01序列，不能是+1/-1这种 上采样：即一个波形采样多少个点，比如一个波形时间为1s，采样率为100Hz，那么上采样倍数就是100，如果一个波形时间为2s，采样率为100Hz，上采样倍数就是50

在实际的仿真代码中，我是这样设置的：

mskmod函数第三个参数是是否做数据差分编码，默认情况下，mskmod是’diff’方式进行调制，也就是默认情况下，mskmod函数认为你的输入是已经经过差分编码的

’nondiff’ 表示输入的数据是未经差分编码 ’diff’ 表示输入的数据是经过差分编码的

但是要注意，这里有一个大坑！！！！ 起初我认为，我在mskmod函数里输入diff，是让mskmod函数帮我做一个差分编码，但是仿真结果总是出错，后来我扣了扣mskmod代码注解的字眼

dataenc can be either ‘diff’ for differentially encoded MSK or ‘nondiff’ for nondifferentially encoded MSK.

这里他用的encoded，意思是已经编码过的，非常坑！！ 言归正传，由于msk有相位模糊的问题，因此一般是加上差分编码，这里附上差分编码的代码：

这里也有一个误区，有两种常见的差分编码形式，需要好好区分，用对才可以

而另一种差分编码（查分曼彻斯特码）的编码方式通常用在bpsk里：

初相这个我没有太多研究，一般都是设置为0，默认也是0

经过对mskmod函数的探索，我将mskmod配置到我需要的形式，然后进行正交调制，我的目标是，正交调制的输出和mskmod函数输出相同，参考的是书上这段表格 mskmod函数代码为：

正交调制的代码为：

其正交调制框图参考书上的结构：

最终得到的波形图为 从上到下依次是，i路数据，q路数据，两路合并的数据。 通过对比可以发现，与书上是一样的，证明我们做的结果是正确的

解调部分就比较简单，就是正常的正交解调

这里需要对i路进行一个补偿和延迟，这是因为： 仔细观察书上的那个表格，其实要传输的序列应该是从0时刻开始，而-2和-1时刻是人为补偿的数据，为了能够进行差分编码，因此i路传输的数据其实是6个数据，前面包含了一个自己设置的1，而这个1与我们的数据没有关系，因此在解调是应该去除半个周期的头部，但是这样操作会使得i路数据变短了，无法解调出最后一个数据，因为缺了半个周期，因此我们这里进行补零，补零不影响数据的解调，只是补足了缺少的数据。

解调时，先进行一次本地相关，然后积分判决就可以得到iq数据，代码如下：

至此，我们完成了所有的代码，可以比对数据得到误码率了 仿真得到的误码率曲线为： 看到，仿真曲线与理论曲线基本一致，因此实现了我们的目标

都看到这里了，我就告诉你一个小秘密，不知道你注意到了没有，代码中没有低通滤波！！

是的，我故意的，没想到吧，但是别急着骂我，上面代码仿真出来的结果肯定和理论值相同的，代码没问题，仿真没问题，结果也没问题，那是我有问题 ？

从理论上来说，与本地载波自相关后肯定要过一个低通滤波来滤除高频项，而且我与本地载波进行了两次自相关，应该要用两个低通滤波才行，这就涉及到下面我说的两个问题了 . .

这里有我本人也没搞清楚的两个点，希望大家一起讨论：