システム同定の線形近​似(ARXモデル)で​ゲイン線図を見る

行いたいこと：一次遅れ{1/(s+1)}に対して-1,1のランダム2値信号を入力した際の出力データ(sample.csv)をARXで近似，ゲイン線図を見たい．

行ったこと：データからARX関数を用いて算出したパラメータで伝達関数を作り，ゲイン線図を見た．

問題点：明らかに一次遅れのゲイン線図からかけ離れていること(sample.jpg)がわかるため行っていることがおかしいと考える．

考えられること：おそらく伝達関数の作成が間違えていると思います．どなたかご教授お願いします．

補足

sample.csvはlsimを用いて作成しました．

clear
data = csvread("hogehoge\sample.csv");
time = data(:,1);%周期
input = data(:,2);%入力
output = data(:,3);%出力

data_2 = iddata(input,output,0.01);
System = arx(data_2,[1 1 1])
%%以下はARX関数が算出したパラメータを代入
Gs = tf([0.02634],[1 0.001301],0.01)
bodemag(Gs)

Answer 1

[milk]

解決できました．inputとoutputの順番が逆でした．．．．