Функция FDA () FDA R - Как это работает? - PullRequest
Новая
       8

Функция FDA () FDA R - Как это работает?

0 голосов
/ 27 января 2019

Я искал простой способ выполнить анализ Фурье временных рядов.Кажется, мне нужна функция fourier() в пакете FDA от R, но я не могу понять, как она работает.Вот простой пример: 

maxTime <- 10
omega <- (2*pi)/1
time <- seq(0, maxTime, maxTime/200)
signal <- 5*cos(omega*time)
# 
library(fda)
fit <- fourier(signal, 2)
#
summary(fit)

   const             sin1              cos1        
 Min.   :0.3162   Min.   :-0.4303   Min.   :-0.4472  
 1st Qu.:0.3162   1st Qu.:-0.2525   1st Qu.:-0.4419  
 Median :0.3162   Median : 0.0000   Median :-0.3691  
 Mean   :0.3162   Mean   : 0.0000   Mean   :-0.1376  
 3rd Qu.:0.3162   3rd Qu.: 0.2525   3rd Qu.: 0.2525  
 Max.   :0.3162   Max.   : 0.4303   Max.   : 0.4472

Я хотел бы знать, как вычисляются значения выше?Документация пакета не объясняет это;кто-нибудь знает?

1 Ответ

0 голосов
/ 27 января 2019

Из описания функции: Обратите внимание, что это было последнее изменение в 2007 году (https://github.com/cran/fda/blob/master/R/fourier.R), поэтому более новая документация может быть более точной или лучше изучение самой функции может дать больше информации.

Вычисляет производную NDERIV базиса ряда Фурье для функций NBASIS с периодом PERIOD, которые оцениваются по значениям в векторе X

Возвращает матрицу значений функций N по NBASIS

Примечание:Число базисных функций всегда нечетное. Если аргумент NBASIS четный, он увеличивается на единицу.

...