Я збираю дані про температуру з холодильника. Дані виглядають як хвиля. Я хотів би визначити період та частоту хвилі (щоб я міг виміряти, чи можуть зміни холодильника мати якийсь ефект).
Я використовую R, і я думаю, що мені потрібно використовувати FFT на даних, але я не впевнений, куди звернутися звідти. Я дуже новачок у аналізі R та сигналів, тому будь-яка допомога буде дуже вдячна!
Ось хвиля, яку я виробляю:
Ось мій R-код поки що:
require(graphics)
library(DBI)
library(RSQLite)
drv <- dbDriver("SQLite")
conn <- dbConnect(drv, dbname = "s.sqlite3")
query <- function(con, query) {
rs <- dbSendQuery(con, query)
data <- fetch(rs, n = -1)
dbClearResult(rs)
data
}
box <- query(conn, "
SELECT id,
humidity / 10.0 as humidity,
temp / 10.0 as temp,
ambient_temp / 10.0 as ambient_temp,
ambient_humidity / 10.0 as ambient_humidity,
created_at
FROM measurements ORDER BY id DESC LIMIT 3600
")
box$x <- as.POSIXct(box$created_at, tz = "UTC")
box$x_n <- box$temp - mean(box$temp)
png(filename = "normalized.png", height = 750, width = 1000, bg = "white")
plot(box$x, box$x_n, type="l")
# Pad the de-meaned signal so the length is 10 * 3600
N_fft <- 3600 * 10
padded <- c(box$x_n, seq(0, 0, length= (N_fft - length(box$x_n))))
X_f <- fft(padded)
PSD <- 10 * log10(abs(X_f) ** 2)
png(filename = "PSD.png", height = 750, width = 1000, bg = "white")
plot(PSD, type="line")
zoom <- PSD[1:300]
png(filename = "zoom.png", height = 750, width = 1000, bg = "white")
plot(zoom, type="l")
# Find the index with the highest point on the left half
index <- which(PSD == max(PSD[1:length(PSD) / 2]))
# Mark it in green on the zoomed in graph
abline(v = index, col="green")
f_s <- 0.5 # sample rate in Hz
wave_hz <- index * (f_s / N_fft)
print(1 / (wave_hz * 60))
Я відправив код R разом з SQLite бази даних тут .
Ось графік нормованого графіка (із середнім вилученим):
Все йде нормально. Ось графік спектральної щільності:
Потім ми збільшуємо масштаб лівої частини сюжету і зеленою лінією відзначаємо найвищий показник (що становить 70):
Нарешті обчислюємо частоту хвилі. Ця хвиля дуже повільна, тому ми перетворюємо її в хвилини за цикл і роздруковуємо це значення, яке становить 17,14286.
Ось мої дані у форматі з обмеженими вкладками, якщо хтось ще хоче спробувати.
Дякую за допомогу! Ця проблема була важкою (для мене), але я чудово провів час!
