小弟目前在解系统辨识的问题，u w q theta分别是输入，delta是输出
以下是小弟的code
uu=[delta];
yy=[u w q theta];
dry=iddata(yy,uu,0.1);
ze=dry(1:400);
ze=detrend(ze);
na=[1 1 0 1;1 1 1 0;1 1 1 0;0 0 1 1];
nb=[1;1;1;0];
nk=[1;1;1;0];
LS=arx(ze,[na nb nk]);
PEM=pem(ze,LS);
compare(ze,LS,PEM);
出现错误
Error using pem (line 102)
Index exceeds matrix dimensions.
Error in test (line 10)
PEM=pem(ze,LS);
谢谢