ARM DSP库中有一个专门的信噪比计数函数arm_snr_f32

eric2013

函数：

1. /**
   
2. * @brief  Caluclation of SNR
   
3. * @param[in]  pRef         Pointer to the reference buffer
   
4. * @param[in]  pTest        Pointer to the test buffer
   
5. * @param[in]  buffSize        total number of samples
   
6. * @return     SNR
   
7. * The function Caluclates signal to noise ratio for the reference output
   
8. * and test output
   
9. */
   
10. 
    
11. float arm_snr_f32(float *pRef, float *pTest, uint32_t buffSize)
    
12. {
    
13.   float EnergySignal = 0.0, EnergyError = 0.0;
    
14.   uint32_t i;
    
15.   float SNR;
    
16.   int temp;
    
17.   int *test;
    
18. 
    
19.   for (i = 0; i < buffSize; i++)
    
20.     {
    
21.           /* Checking for a NAN value in pRef array */
    
22.           test =   (int *)(&pRef[i]);
    
23.       temp =  *test;
    
24. 
    
25.           if (temp == 0x7FC00000)
    
26.           {
    
27.                           return(0);
    
28.           }
    
29. 
    
30.           /* Checking for a NAN value in pTest array */
    
31.           test =   (int *)(&pTest[i]);
    
32.       temp =  *test;
    
33. 
    
34.           if (temp == 0x7FC00000)
    
35.           {
    
36.                           return(0);
    
37.           }
    
38.       EnergySignal += pRef[i] * pRef[i];
    
39.       EnergyError += (pRef[i] - pTest[i]) * (pRef[i] - pTest[i]);
    
40.     }
    
41. 
    
42.         /* Checking for a NAN value in EnergyError */
    
43.         test =   (int *)(&EnergyError);
    
44.     temp =  *test;
    
45. 
    
46.     if (temp == 0x7FC00000)
    
47.     {
    
48.                   return(0);
    
49.     }
    
50. 
    
51. 
    
52.   SNR = 10 * log10 (EnergySignal / EnergyError);
    
53. 
    
54.   return (SNR);
    
55. 
    
56. }

复制代码

说明：

wu32698318

非常好，期待使用