using TySystemIdentification
using TyControlSystems

输入输出数据

u = [1.4601,0.8849,1.1854,1.0887,1.413,1.3096,1.0651,0.7148,1.3571,1.0557,1.1923,1.3335,
1.4374,1.2905,0.841,1.0245,1.4483,1.4335,1.0282,1.4149,0.7463,0.9822,1.3505,0.7078,0.8111,0.8622,0.8589,1.183,0.9177,0.859,0.7122,1.2974,1.056,1.4454,1.0727,1.0349,1.3769,1.1201,0.8621,1.2377,1.3704,0.7157,1.245,1.0035,1.3654,1.1022,1.2675,1.0431]

y = [0,0,8.7606,13.1939,17.4100,17.6361,18.7627,18.5296,17.0414,13.4154,14.4539,14.5900,
16.1104,17.6853,19.4981,19.5935,16.4106,14.3359,15.7463,18.1179,17.7840,18.8104,15.3086,13.7004,14.8178,13.2354,12.2993,11.6001,11.6074,13.7662,14.1950,13.7630,11.8713,13.8566,14.6944,17.8659,17.6543,16.6386,17.1071,16.5373,14.643,15.0862,16.8058,14.7641,15.4976,14.6790,16.6552,16.6301]

na, nb = 4, 4
U = iddata(y, u, 0.1)
T = arx(U, na, nb)

H=tf(T)
t=0:0.1:4.7
y_lsim=lsim(H,reshape(u,1,length(u)),t)
y_lsim_vec=vec(y_lsim.y)
plot(t,y_lsim_vec)
hold("on")
plot(t,y,"o")
plot(t,u)

ERROR: KeyError: key :y not found