Сигнал в частотной области с OpenCV dft - PullRequest
Новая
       47

Сигнал в частотной области с OpenCV dft

0 голосов
/ 13 октября 2018

Я экспериментирую с cv::dft: синусоидальный сигнал 1 Гц генерируется и отображается в частотной области.Но по какой-то причине он не имеет максимальной составляющей на 1 Гц.Мой код следующий: 

       const int FRAME_RATE = 20;             //!< sampling rate in [Hz]
       const int WINDOW_SIZE = 256;         

       double len = double(WINDOW_SIZE)/double(FRAME_RATE);   // signal length in seconds
       double Fb = 1./len;                                     // frequency bin in Hz

       // Constructing frequency vector
       std::vector<double> f;
       double freq_step = 0;
       for (int i = 0; i < WINDOW_SIZE; ++i)
       {
          f.push_back(freq_step);
          freq_step += Fb;
       }


   // Create time vector
   std::vector<double> t;
   double time_step = 0;
   for(int i = 0; i<WINDOW_SIZE; ++i)
   {
       t.push_back(time_step);
       time_step += 1./double(FRAME_RATE);
   }

   // Creating sin signal with 1Hz period
   std::vector<double> y;
   for(auto val : t)
   {
       y.push_back(sin(1*FRAME_RATE*val));
   }

   // Compute DFT
   cv::Mat fd;
   cv::dft(y, fd, cv::DFT_REAL_OUTPUT);
   fd = cv::abs(fd);

Если я нанесу сигнал во временной и частотной области: plot(t, y); plot(f, fd) результат будет следующим:

enter image description here

Сигнал времени хороший, но частота сигнала максимальная около 6 Гц вместо 1 Гц.

Где я допустил ошибку?

...