[請益]系統辨識問題

[JasonScott (傑生史考特)]

小弟目前在解系統辨識的問題，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);

謝謝

--
推 jj5110:解剖學老師：大體上來說
推 AesopChen:大體老師：上去講台說嗎? 好等我一下
推 Aragorn:                                    我先把自己縫起來
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.202.124
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1443541319.A.257.html