ARMAモデルの実験

ARMAモデル
G(z^{-1})=\dfrac{z^{-1}+0.5 z^{-2}}{1+1.5z^{-1}+0.7z^{-2}}とする。


\dfrac{1}{z}での伝達関数を作成するにはfilt()関数を用いる

>> filt([0,1,0.5],[1,1.5,0.7],0.1)

ans =
 
      z^-1 + 0.5 z^-2
  -----------------------
  1 + 1.5 z^-1 + 0.7 z^-2
 
Sample time: 0.1 seconds
Discrete-time transfer function.

ステップ応答
f:id:seinzumtode:20200517175238p:plain

周波数応答
f:id:seinzumtode:20200517175232p:plain

インパルス応答
f:id:seinzumtode:20200517175241p:plain

v=randn()に対する時間応答

t=0:99
v = randn(1,100)
lsim(G,v,t)

f:id:seinzumtode:20200517175734p:plain