关于浮点运算中除数为0或者NaN非数值的情况，大家有什么好的处理思路

eric2013

DSP库里面当前只有一个信噪比的函数做了简单处理，仅仅是简单处理，判断除数是否为NAN值0x7FC00000（NaN的一种情况）。

大家在使用有什么好的建议。

1. float arm_snr_f32(float *pRef, float *pTest, uint32_t buffSize)
   
2. {
   
3.   float EnergySignal = 0.0, EnergyError = 0.0;
   
4.   uint32_t i;
   
5.   float SNR;
   
6.   int temp;
   
7.   int *test;
   
8. 
   
9.   for (i = 0; i < buffSize; i++)
   
10.     {
    
11.            /* Checking for a NAN value in pRef array */
    
12.           test =   (int *)(&pRef[i]);
    
13.       temp =  *test;
    
14. 
    
15.           if (temp == 0x7FC00000)
    
16.           {
    
17.                           return(0);
    
18.           }
    
19. 
    
20.           /* Checking for a NAN value in pTest array */
    
21.           test =   (int *)(&pTest[i]);
    
22.       temp =  *test;
    
23. 
    
24.           if (temp == 0x7FC00000)
    
25.           {
    
26.                           return(0);
    
27.           }
    
28.       EnergySignal += pRef[i] * pRef[i];
    
29.       EnergyError += (pRef[i] - pTest[i]) * (pRef[i] - pTest[i]);
    
30.     }
    
31. 
    
32.         /* Checking for a NAN value in EnergyError */
    
33.         test =   (int *)(&EnergyError);
    
34.     temp =  *test;
    
35. 
    
36.     if (temp == 0x7FC00000)
    
37.     {
    
38.                   return(0);
    
39.     }
    
40. 
    
41. 
    
42.   SNR = 10 * log10 (EnergySignal / EnergyError);
    
43. 
    
44.   return (SNR);
    
45. 
    
46. }

复制代码

ps：IEEE 754-1985中，用指数部分全为1、小数部分非零表示NaN

eric2013

matlab的处理：

xiashuang

实际用MDK发现除数为0不用考虑这个问题，实测STM32F4/F3/F1没问题，F0不行

eric2013

xiashuang 发表于 2020-3-29 17:37
实际用MDK发现除数为0不用考虑这个问题，实测STM32F4/F3/F1没问题，F0不行

谢谢帮忙测试。

WZH

#define UTILS_IS_NAN(x)                ((x) != (x))
#define UTILS_NAN_ZERO(x)        (x = UTILS_IS_NAN(x) ? 0.0 : x)
可以用类似的宏定义实现，C语言里规定Nan!=Nan

eric2013

WZH 发表于 2023-10-9 11:00
#define UTILS_IS_NAN(x)                ((x) != (x))
#define UTILS_NAN_ZERO(x)        (x = UTILS_IS_NAN(x) ? 0.0 : x)
可 ...

谢谢分享。

emwin

想复杂了，value != 0.0f 就可以了。
Floating Point to Hex Converter (gregstoll.com)

eric2013

emwin 发表于 2023-10-10 09:42
想复杂了，value != 0.0f 就可以了。
Floating Point to Hex Converter (gregstoll.com)

0对应的是全0

sean4194304

不知道是不是类似的问题。 曾经 用   H7 做大量浮点运算，根据结果显示波形。极少情况下，会出现结果超级大，波形跑到显示区外，死机了。实在找不出原因，只能在显示波形之前，判断一下，是否超出范围，超出就强制赋个最大值再显示。    后来再没发现死机了。