#author("2019-08-26T13:48:43+09:00","default:Miyashita","Miyashita")
#author("2019-11-15T15:42:14+09:00","default:Miyashita","Miyashita")
*Julia演習問題:トレンド解析の基礎 [#l201750e]
***問題 [#me6bacc6]
+csv(コンマ区切り)データ "dat_climate.csv" を読み込みなさい.~
読み込みに使う関数は自由.~
#codeprettify(lang-julia){{
DelimitedFiles.readdlm(filename, ',', skipstart=1) # 例
}}
+不等間隔のデータを 1 日毎に補間しなさい.~
#codeprettify(lang-julia){{
using Interpolations
interpolate((torg, ), val, Gridded(Linear()))
}}
+補間データを元に,線形トレンドを求めて図にしなさい.~
#codeprettify(lang-julia){{
Polynomials.polyfit(t, val, 1)
plt = Plots.plot(t, Polynomials.polyval(...))
}}
+データのスペクトル解析を行い,固有日数 (1/f) を求めなさい.~
#codeprettify(lang-julia){{
pdg = DSP.welch_pgram(val, ...; fs=1.0)
}}

***データ [#mff0554e]
-csvファイル
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/python/dat_climate.csv,dat_climate.csv)~
1行目はヘッダー~
2行目以降は日時,データ
***解答例(参考) [#f917a104]
1年周期になるように意図的なパラメータを調整した.~
1年周期になるように意図的にパラメータを調整した.~
厳密性が求められない演習用データにつきご容赦ください.~

&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/julia/uniformed_data.png,500x375);
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/julia/trend.png,500x375);

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes