MATLAB函数详解系列一、qammod与qamdemode函数.

目录

正文开始前的一些絮叨

一、qammod函数详解

二、不用qammod函数来写出代码

三、qamdemod的函数详解

四、不用qamdemod函数的代码实现

---

正文开始前的一些絮叨

本文是不使用任何matlab自带函数的方法来实现4QAM,以及更高阶MQAM的调制的调理，不使用MATLAB的qammod和qamdemod函数。

写这篇文章的原因是我在面试通信物理层算法岗的时候，说我自己做了整个OFDM模块，面试官很吃惊，表示我全都是自己写的吗？我说是的，之后我才明白他的意思是不套用网上存在的任何代码，只按照公式或者算法来写代码。这个我承认我真做不到。一般我们做科研，是在csdn或者pudn上下载别人的代码，根据代码不断地去改参数，改模块，从而写代码的。

面试官语重心长的跟我饿即使原因（此处划重点）：算法研究人员应该做到不依赖库里的函数，自己去写物理层的代码。我想了想，是呀一般做物理层都是用C++，不可能会像matlab一下子有个berawgn或者pskmod函数调用的，于是这三天的时间内，我单独去学习了下QAM调制，并且想办法不依赖任何matlab函数，并用matlab自带的这些函数进行验证我写的算法是否正确。

一、qammod函数详解

我直接给出qammod函数的代码用法和结论

M=4;%M是高阶调制的阶数 M=4就是 4QAM M=64就是64QAM
A=[0 1 2 3 ];
qammod(A,M);

 这里发现MATLAB输出结果为

我直接先告诉大家结论，matlab的qammod函数输出结果是已经经过格雷码调制的，所以大家使用的时候不需要再转换为格雷了！

为什么这么说呢，我们直接画出这四个点的图，其中A矩阵的0 1 2 3二进制结果分别代表是'00','01','10','11'。 很明显的看出每个星座只相差一个1，16QAM大家可以自己去验证。

二、不用qammod函数来写出代码

这个我们就需要两个步骤了，我们以4QAM为例，

1.首先生成二进制码 A=[0 1 2 3];

2.进行格雷码转换 Gray_A=[0 1 3 2]； %这步其实也就是格雷码编码规则

M=4;
A=[0 1 2 3];
Gray_A=[0 1 3 2];%转化规则的代码我就不放出来了网上都有
sqrt(2)*exp(j*(2*pi/M*Gray_A+pi/4)) %

运行之后的结果如下， 

大家发现其实这个公式的结果与qammod函数的结果相差了 一个pi/2的相位，原因就是qammod的00是从-1+i 开始的，而我们的函数的00位置是从1+i开始的。但是它们都属于格雷码只不过是位置不同而已。符合格雷码要求相隔的码只相差一个比特。

三、qamdemod的函数详解

四、不用qamdemod函数的代码实现

第三四部分涉及到LLR进行硬判决，计算高阶QAM的LLR比较复杂，以后有时间在更新