Matematické Fórum

↑ Michal 64:
Ahoj,
pokud si dobře vzpomínám, tak se to dělá tak, že nejprve si určíš trend (zde by mohl být lineární) a data od trendu očistíš. Tedy v tomto případě si sestavíš lineární model (prostě metoda nejmenších čtverců...), jeho směrnice bude trend. V dalším kroku od každé naměřené hodnoty odečteš hodnotu, která vyjde po dosatení do modelu. Výsledkem by měly být hezčí data, u kterých můžeš vypočítat autokorelační funkci nebo výkonové spektrum (technicky vzato je to jedno, obojí se bude stejně v počítači dělat přes FFT) - a tam, kde bude mít autokorelační funkce/výkonové spektrum maximum, máš významnou periodicitu (pokud je maxim více, znamená to, že se v řadě ukrývá několik cyklů).


Tady konkrétně by šlo zkusit i jiný přístup, protože ta data jsou docela hezká sinusovka. Data očištěná od trendu zkus aproximovat funkcí:


, kde A, B, C jsou parametry.  Hledaná perioda pak bude:

↑ Michal 64:
Volba těch parametrů (A, B, C) je záležitostí numerické aproximace. Myslím, že R to samo o sobě neumí, aproximaci harmonickou funkcí by mohl umět např. gnuplot, případně si můžeš stáhnout udělátko do matlabu.

Takže nakonec bude asi nejjednodušší, když si uděláš fft a podíváš se kde má maximum amplitudové spektrum. To by ti mělo dát dostatečný odhad.

↑ Formol:

Ahoj tak sem postupoval následujícím způsobem pro FFT jsem použil Matlab fukci Y=fft(x), kde vektor x byly data(původní data - hodnoty lineárního trendu).No potom mi matlab vyhodil vektor čísel, který když vykreslim do grafu vypadá takto.



takže v obrázku vydim pouze jedno maximum ale moudrej z něho teda nejsem.Dokážeš z toho poznat periodu sezonní složky (postup by měl být shodný s postupem předepsaným v prvním přísppěvku) ?

Děkuji

Ahoj, ano termín FFT a fourierova transformace sem slyšel poprvé až tady dosud sem to nikdy nepoužíval, takže sem se to sanžil narychlo pochopit z internetu.No v tom případě je perioda okolo 7,4.

Každopádně díky za pomoc a rady

Zdravím,
zkusil jsem aproximaci podle vzorce:

vyšel mi tento výsledek:
y = −0.1651238651643678−1.769676299097698*x−8.912554336244837*sin(0.7403812988506794*x+1.385902605832021)
Perioda pak vychází: 8.486

↑ mák:

Aha no díky moc a mám dotaz tu aproximaci si jsi zvolil sám takhle nebo je to náký standartní model ?