请教懂音乐频谱的，烦请解释一下面的各条语句的含义。

rjx007 · 发表于 2017-10-14 09:31:54

请教懂音乐频谱的，解释一下面的各条语句的含义。主要是这些语句中除256，乘32768，乘65536等是什么作用？最好能给每句加个注释，本人智商低下，看不懂，只能求助了。谢谢！！！

void powerMag()
{       int32_t lX,lY;
            uint32_t i;
            for (i=0;i<256;i++)
            {
                     lX=(lBUFOUT<<16)>>16; /* sine_cosine --> cos */
                     lY=(lBUFOUT>>16); /* sine_cosine --> sin */
                     {
                              float X=256*((float)lX)/32768;
                              float Y=256*((float)lY)/32768;
                              float Mag = sqrt(X*X+ Y*Y)/256;  // 先平方和,再开方
                              lBUFMAG =(long)(Mag*65536);
                     }
            }
}

rjx007 · 发表于 2017-10-14 09:31:55

本帖最后由 rjx007 于 2017-10-17 13:44 编辑

   这是在网上找到的最详细的解释，贴上来共享：
   之所以语句这么写是基于芯片的计算习惯，DSP最擅长的就是浮点数的乘机运算。语句分析如下：
首先定义的的一个16位的有符号的数据IX 和IY 这两个只是中间变量，然后定义的i，是32位的无符号型。
语句的目的的确是你理解的Mag = sqrt(X*X + Y*Y)/NPT。但直接这么写不符合DSP的计算习惯，也就是不符合浮点运算的习惯。因此语句在for函数i写道
lX  = (lBUFOUT << 16) >> 16 //就是取32位的i的低16位数据，
lY  = (lBUFOUT >> 16);       //是取高16位数据。下面的两句
float X = NPT * ((float)lX) /32768;
float Y = NPT * ((float)lY) /32768
目的就是把数据浮点化，至于为什么是除以32768 。可以这么说，浮点化就好像10进制里面的科学计数法。32768=2的15次方。除以32768也就是去除了浮点数后面的那个基数，只剩下前面的。比如1991 改写成1.991*10的三次幂，再除以10的三次方，只剩下1.991，便于余下的运算。
之后的语句就不用解释了吧？至于最后一句要乘以65536，是因为我们定义的数据和我们需要求得的数据都是无符号32位的，
之前已经把32位的数据拆开又分别浮点化了，又开了个根号，所以再把它变回来只需要乘以2的16次，也就是65536.
打个比方吧。比如说问你什么时候生日，你说是19911030，然而DSP是不习惯这么干的，他需要把它拆开为1991和1030。再写成1.991x10的3次方和1.030x10的3次方。然后才能进行其他的运算。