C语言函数库atan2(a, b)与atan(a/b)的区别@项目简介：基于PSS序列（时域）估计小数倍频偏

atan2和atan函数都在#include<math.h>函数库中，区别如下：

 （1）atan2(a,b)是4象限反正切，它的取值不仅取决于正切值a/b，还取决于点 (b, a) 落入哪个象限：

   当点(b, a) 落入第一象限时，atan2(a,b)的范围是 0 ~ pi/2;
　当点(b, a) 落入第二象限时，atan2(a,b)的范围是 pi/2 ~ pi;
   当点(b, a) 落入第三象限时，atan2(a,b)的范围是 －pi～－pi/2;
　当点(b, a) 落入第四象限时，atan2(a,b)的范围是 -pi/2～0.


  （2） atan(a/b) 仅仅根据正切值为a/b求出对应的角度 （可以看作仅仅是2象限反正切）：
   当 a/b > 0 时，atan(a/b)取值范围是 0 ~ pi/2；
   当 a/b < 0 时，atan(a/b)取值范围是 -pi/2～0


故 atan2(a,b) = atan(a/b) 仅仅发生在 点 (b, a) 落入第一象限 （b>0, a>0）或 第四象限（b>0, a<0）。当点 (b, a) 落入第二、三象限时，很显然atan2(a,b) 不等于 atan(a/b) ，并且atan2(a,b)也不可能等于 2*atan(a/b) 。这是因为，假如点 (b, a) 落入第二象限，则 a/b<0,  故atan(a/b)取值范围始终是 -pi/2～０，2*atan(a/b) 的取值范围是－pi～0，然而，atan2(a,b)的范围是 pi/2 ~ pi，故不可能有atan2(a,b) = 2*atan(a/b) 。假如点(b, a) 落入第三象限，则则 a/b>0 , 故 atan(a/b) 取值范围是 0 ~ pi/2，2*atan(a/b) 的取值范围是 0 ~ pi，而此时atan2(a,b)的范围是 －pi～－pi/2，很显然，atan2(a,b) = 2*atan(a/b) 

举个最简单的例子，a = 1, b = -1，则 atan(a/b) = atan(-1) = -pi/4, 而 atan2(a,b) = 3*pi/4

项目中用来求复数sum的相位，即FractionalFreOffset =  atan2(sum.imag,sum.real) / PI，所以只能使用atan2函数，如果使用atan函数时会造成第二象限和第三象限的值不对。

附录：求小数倍频偏代码。

#include "complexMath.h"
#include <math.h>
/*
 * @Function FractionalFreOffsetDetect
 * @Descrption 估计小数倍频偏
 * @Input complexRecPssData        接收序列中的PSS序列（时域）
 * @Output float       FractionalFreOffset 小数倍频偏值
 * @Author 高世皓
 * @Date 2016-08-03
 */
float FractionalFreOffsetDetect(complex* RecPssData)
{
    int cell_id0 = 0;
    int k;
    float FractionalFreOffset = 0;//小数倍频偏值
    complex PssSymbol0[fftSIZE];//本地PSS序列（时域）
    complex data1;//存储前一项运算的结果
    complex data2;//存储后一项运算的结果
    complex sum;//放置 fftSIZE/2 项相加之和
    sum.real = 0;
    sum.imag = 0;
    PssSymbolGen(cell_id0, fftSIZE, PssSymbol0);//产生本地PSS序列（时域）
    for(k=0; k<fftSIZE/2; k++)
    {
        data1 = complex_conj( complex_conj_mul(RecPssData[k], PssSymbol0[k]));//公式前一项
        data2 = complex_conj_mul(RecPssData[k + fftSIZE/2], PssSymbol0[k + fftSIZE/2]);//公式前一项
        sum = complex_add( sum, complex_mul(data1, data2) );
    }
   FractionalFreOffset =  atan2(sum.imag,sum.real) / PI;
   return FractionalFreOffset;


}


int main()
{
     int i=0,j;
     float FractionalFreOffset = 0;
     FILE *fp;//文件指针
     dataType realInformation[fftSIZE];//接收数据中PSS序列的实部
     dataType imagInformation[fftSIZE];//接收数据中PSS序列的虚部
     complex testData[fftSIZE];//接收数据中PSS序列


   /******读取Matlab生成的实部******/


    fp = fopen("C:\\Users\\Bug\\Desktop\\MFreOffsetDetect\\FractionalFreOffsetDetect\\DataReal.txt", "r");//以文本方式打开文件。
    if(fp == NULL) //打开文件出错。
        return -1;
    while(fscanf(fp, "%lf", &realInformation[i]) != EOF) //读取数据到数组，直到文件结尾(返回EOF)
        i++;
    fclose(fp);//关闭文件
//    for(j = 0; j < i; j ++)//循环输出数组元素。
//    {
//        printf("%lf ", realInformation[j]);
//    }




    /******读取Matlab生成的虚部******/


    fp = fopen("C:\\Users\\Bug\\Desktop\\MFreOffsetDetect\\FractionalFreOffsetDetect\\DataImag.txt", "r");//以文本方式打开文件。
    if(fp == NULL) //打开文件出错。
        return -1;
    i=0;


    while(fscanf(fp, "%lf", &imagInformation[i]) != EOF) //读取数据到数组，直到文件结尾(返回EOF)
        i++;
    fclose(fp);//关闭文件


//    for(j = 0; j < i; j ++)//循环输出数组元素。
//    {
//        printf("%lf ", imagInformation[j]);
//    }


for (i = 0; i < fftSIZE; i++){
testData[i] = dcomplex(realInformation[i], imagInformation[i]);
}
FractionalFreOffset = FractionalFreOffsetDetect(testData);
printf("\n\FractionalFreOffset=%f\n",FractionalFreOffset);
    return 0;
}