From 83d24ce79e7a2b68e7def296e6546d72d46e8b69 Mon Sep 17 00:00:00 2001
From: mcrucifix <michel.crucifix@uclouvain.be>
Date: Thu, 19 Jun 2025 16:44:53 +0200
Subject: [PATCH] periodogram
---
R/periodogram.R | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/R/periodogram.R b/R/periodogram.R
index faf7d44..3bab598 100644
--- a/R/periodogram.R
+++ b/R/periodogram.R
@@ -16,16 +16,16 @@ periodogram <- function(xdata){
N2 <- ceiling(N/2)
freqs <- ((seq(N)-1)/dt/N)[0:N2]
f <- fft(xdata)[0:N2]
- Power <- Mod(f)^2
+ Power <- Mod(f)^2 / (dt*N2)^2
Phase <- Arg(f)
- out <- list(Freq=freqs, Power=Power, Phase=Phase)
+ OUT <- list(Freq=freqs, Power=Power, Phase=Phase)
attr(OUT, "class") = "periodogram"
return(OUT)
}
#' @rdname periodogram
#' @export
-plot.periodogram <- function(X,...){
- plot(X$Freq, X$Mod, type='l', log='xy',...)
+plot.periodogram <- function(X,log='xy', xlabel= "Frequency", ylabel="Power density", ...){
+ plot(X$Freq, X$Power, type='l', log=log,xlabel=xlabel, ylabel=ylabel, ...)
}
--
GitLab