arm_recip_q31 stm32 DSP 倒数使用方法

最新推荐文章于 2024-10-22 21:27:43 发布

原创最新推荐文章于 2024-10-22 21:27:43 发布 · 994 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#stm32 #dsp开发

SMT32 同时被 2 个专栏收录

2 篇文章

订阅专栏

DSP

1 篇文章

订阅专栏

文章探讨了Q31定点数表示下的倒数计算问题，通过程序验证发现对于正数结果正确，负数处理有误。作者使用了Arm的RecipTableQ31进行计算，并提取出整数部分，最终得到近似值3.001，与预期相符。

网上查了一下这个函数的使用方法，貌似没有，简单记录一下

q31是定点小数，把32位有符号整数的第一位定义为符号位，后31位为小数位，因此q31的范围是-0.9999 ~ 0.9999之间，小数的倒数必定大于1或小于-1，它的结果又是q31。所以这个函数会把大于1的整数部分存储在结果中，那么这个函数就会返回一个存储了几个bit作为整数部分的值。

简单写了个程序验证下，正数结果正确，负数不对，暂时没搞明白，如有错误欢迎指正。

extern q31_t armRecipTableQ31[64];//这个声明应该在arm_common_tables.h，但是我的STM32CubeMX没有给我这个头文件，就自己声明了下
q31_t orig_val=0.33333f*0x7FFFFFFF;//0.33333f是浮点数，下边的代码就是求它的倒数
    
uint32_t sb=arm_recip_q31(orig_val,&recip_val,armRecipTableQ31);
uint32_t sign=recip_val&0x80000000u;
recip_val&=0x7FFFFFFF;
uint32_t intMask=0;
uint32_t _sb=sb;
while(_sb){
    intMask=intMask<<1;
    intMask=intMask|1;
    _sb--;
}
uint32_t intVal=(recip_val&(intMask << (31-sb)))>>(31-sb);

float recip_val_f=1.0f*((recip_val<<(sb+1)) & 0x7FFFFFFF)/0x7FFFFFFF + intVal;
recip_val_f*= sign?-1:1;

运行结果recip_val_f的值是3.001，与预期基本吻合