今回は，Ｒを使った微分方程式の数値計算についての説明というよりも，数値計算結果を3次元プロットする方法の紹介です．私が以下で使っている方法は，ベストではないかもしれませんが，とりあえず，3次元空間内の軌道を立体的に描くことができます．結果は…

今日のセミナーの課題は， 問題 数値的に生成した1次自己回帰過程のサンプル時系列の，自己共分散関数とパワースペクトル（ピリオドグラム）を推定し，理論値と比較せよ．推定は，複数のサンプル時系列を使って平均として求めよ．解答としては，以下のような…