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

Я экспериментирую с 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) результат следующий:

введите описание изображения здесь

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

Где я взял ошибку?

1

Решение

Задача ещё не решена.

Другие решения

Других решений пока нет …

По вопросам рекламы [email protected]