Pythonの全くの初心者で、分かりにくい質問かと思いますが、投稿させていただきます。
時系列モデルのシミュレーションを行うためのコード作成に関する質問です。
具体的には、
X(t)=0.5*X(t-1)+e(t) for t<=50
X(t)=X(t-1)+e(t) for t>50
(ただし、e(t)は、平均0、標準偏差1の正規分布の乱数)
について、
(1): t=100のときのX(100)の値をシミュレートする関数のコードを作成する
(2): 適当な値X(0)を入力した場合のX(1),X(2),X(3),...X(100)の値をプロットしたグラフ(横軸にt=1,2,3,...100、縦軸にX(1),X(2),X(3),...X(100))を作成する
が目標です。
(1)については、以下のコードを作成しましたが、正しいでしょうか。
import numpy as np
x = 0
for i in range(50):
x = 0.5*x + np.random.randn()
print(x)
for i in range(50):
x += np.random.randn()
print(x)
(2)については、例えば、numpyを利用する方法も試しましたが、関数(x = 0.5*x + np.random.randn())から導かれる計算結果を行列やベクトルに変換したり組み込む方法が分かりませんでした。
なお、
Python 3
Jupyter Notebook (anaconda 3)
を使用しています。
かなり色々なサイトで探したのですが、類似のケースが見つからず、こちらに投稿させていただいた次第です。
分かりにくい質問で申し訳ありませんが、何卒よろしくお願いいたします。