为何我的复数求模结果不对为什么？？求助

水中的天空

q15_t pSrc2[10] = {200, 200, 200, 200,300, 400, 400, 300, 300, 400};
q15_t pDst2[5];
q15_t pDst3[5];
void QUMO1(void)
{
   arm_cmplx_mag_q15(pSrc2, pDst2, 5);
   arm_sqrt_q15(pSrc2[0],&pDst3[0]);
}

eric2013

这个是源码的计算问题，如果实部和虚部的平方和小于2^17 = 131072 将按照0算。

水中的天空

但是为何浮点运算却是正确的
而且也能正常求FFT变换
q15_t TestIn[32]={100,200,300,400,100,200,300,400,100,200,300,400,100,200,300,400,100,200,300,400,100,200,300,400,100,200,300,400,100,200,300,400};
uint8_t i;
  //q15格式的32点复数FFT
  for(i=0; i<32; i++)
        {
                testInput[i*2+1] = 0;
                testInput[i*2] = TestIn;
        }
  arm_cfft_q15(&arm_cfft_sR_q15_len32,testInput,0, 1);

水中的天空

eric2013 发表于 2018-9-6 15:56
这个是源码的计算问题，如果实部和虚部的平方和小于2^17 = 131072 将按照0算。

大神，那你的意思就是库自己有问题或者是我有问题，我要把他放大2^15倍后再计算，还是直接转到f32型算法计算，我只想计算32点的FFT，代码效率越高越好，若浮点比定点快我宁愿浮点

eric2013

水中的天空 发表于 2018-9-6 15:57
但是为何浮点运算却是正确的
而且也能正常求FFT变换
q15_t TestIn[32]={100,200,300,400,100,200,300,400 ...

这就所谓的定点数。超出的他能够表示的范围。都将被刷掉。

水中的天空

发现好像也有个规律
对q15格式1024求开方得到5793
(1024/32768)开方得到（5793/32768）

eric2013

水中的天空 发表于 2018-9-6 16:30
发现好像也有个规律
对q15格式1024求开方得到5793
(1024/32768)开方得到（5793/32768）

这就是定点，你是不是没有看我们的教程，抽空看看定点数。
https://forum.anfulai.cn/forum.ph ... 3886&extra=page%3D1

定点数就是用的整数来表示小数，看看下面的数据表就更明白了