开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 语音压缩 > 查看源码
test_real.c
资源名称：4.tar.gz　[点击查看]
上传用户：ozl2332
上传日期：2009-12-28
资源大小：38k
文件大小：3k
源码类别：

语音压缩

开发平台：

C/C++

test_real.c：源码内容
  1. #include "kiss_fftr.h"
  2. #include "_kiss_fft_guts.h"
  3. #include <sys/times.h>
  4. #include <time.h>
  5. #include <unistd.h>
  7. static double cputime(void)
  8. {
  9.     struct tms t;
  10.     times(&t);
  11.     return (double)(t.tms_utime + t.tms_stime)/  sysconf(_SC_CLK_TCK) ;
  12. }
  14. static
  15. double snr_compare( kiss_fft_cpx * vec1,kiss_fft_cpx * vec2, int n)
  16. {
  17.     int k;
  18.     double sigpow,noisepow,err,snr,scale=0;
  19.     sigpow = noisepow = .00000000000000000001; 
  21.     for (k=0;k<n;++k) {
  22.         sigpow += (double)vec1[k].r * (double)vec1[k].r + 
  23.                   (double)vec1[k].i * (double)vec1[k].i;
  24.         err = (double)vec1[k].r - (double)vec2[k].r;
  25.         noisepow += err * err;
  26.         err = (double)vec1[k].i - (double)vec2[k].i;
  27.         noisepow += err * err;
  29.         if (vec1[k].r)
  30.             scale +=(double) vec2[k].r / (double)vec1[k].r;
  31.         /*
  32.         fprintf(stderr,"vec1=");pcpx(vec1+k);
  33.         fprintf(stderr,"vec2=");pcpx(vec2+k);
  34.         */
  35.     }
  36.     snr = 10*log10( sigpow / noisepow );
  37.     scale /= n;
  38.     if (snr<10) {
  39.         printf( "npoor snr, try a scaling factor %fn" , scale );
  40.         exit(1);
  41.     }
  42.     return snr;
  43. }
  44. #define RANDOM
  45. #ifndef RANDOM
  46. #define NFFT 8
  47. #else
  48. #define NFFT 8*3*5
  49. #endif
  51. #ifndef NUMFFTS
  52. #define NUMFFTS 1000
  53. #endif
  56. int main(void)
  57. {
  58.     double ts,tfft,trfft;
  59.     int i;
  60.     kiss_fft_cpx cin[NFFT];
  61.     kiss_fft_scalar rin[NFFT] = {0.309655,0.815653,0.768570,0.591841,0.404767,0.637617,0.007803,0.012665};
  62.     kiss_fft_cpx cout[NFFT];
  63.     kiss_fft_cpx sout[NFFT];
  64.     
  65.     kiss_fft_cfg  kiss_fft_state;
  66.     kiss_fftr_cfg  kiss_fftr_state;
  67.     
  68.     srand(time(0));
  70.     for (i=0;i<NFFT;++i) {
  71. #ifdef RANDOM        
  72.         rin[i] = (kiss_fft_scalar)(rand()-RAND_MAX/2);
  73. #endif        
  74.         cin[i].r = rin[i];
  75.         cin[i].i = 0;
  76.     }
  78.     kiss_fft_state = kiss_fft_alloc(NFFT,0,0,0);
  79.     kiss_fftr_state = kiss_fftr_alloc(NFFT,0,0,0);
  80.     kiss_fft(kiss_fft_state,cin,cout);
  81.     kiss_fftr(kiss_fftr_state,rin,sout);
  82.     printf( "nfft=%d, inverse=%d, snr=%gn",
  83.             NFFT,0, snr_compare(cout,sout,(NFFT/2)+1) );
  84. #ifdef RANDOM        
  85.     ts = cputime();
  86.     for (i=0;i<NUMFFTS;++i) {
  87.         kiss_fft(kiss_fft_state,cin,cout);
  88.     }
  89.     tfft = cputime() - ts;
  90.     
  91.     ts = cputime();
  92.     for (i=0;i<NUMFFTS;++i) {
  93.         kiss_fftr( kiss_fftr_state, rin, cout );
  94.         /* kiss_fftri(kiss_fftr_state,cout,rin); */
  95.     }
  96.     trfft = cputime() - ts;
  98.     printf("%d complex ffts took %gs, real took %gsn",NUMFFTS,tfft,trfft);
  99. #endif        
  100.     free(kiss_fft_state);
  101.     free(kiss_fftr_state);
  103.     kiss_fft_state = kiss_fft_alloc(NFFT,1,0,0);
  104.     kiss_fftr_state = kiss_fftr_alloc(NFFT,1,0,0);
  106.     kiss_fft(kiss_fft_state,cout,cin);
  107.     kiss_fftri(kiss_fftr_state,cout,rin);
  109.     for (i=0;i<NFFT;++i) {
  110.         sout[i].r = rin[i];
  111.         sout[i].i = 0;
  112.     }
  113.     
  114.     printf( "nfft=%d, inverse=%d, snr=%gn",
  115.             NFFT,1, snr_compare(cin,cin,NFFT/2) );
  116.     free(kiss_fft_state);
  117.     free(kiss_fftr_state);
  119.     return 0;
  120. }