Я экспериментирую с 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 Гц.
Где я взял ошибку?
Задача ещё не решена.
Других решений пока нет …