关于armax的应用

yimengxueer

我要对ARMA模型进行参数估计和预测
用armax怎么实现啊？
谢谢
x=[22. 0  29. 1  34. 2  36. 0  20. 1  31. 9  19. 8  15. 2  16. 1 17. 9  29. 9  30. 8  31. 7 27. 8  35. 6  15. 2  19. 0 20. 0 20. 0 21. 0  12. 4  16. 3  21. 0  25. 0  33. 2  36. 0  28. 4 40. 0  41. 6  44. 6 37. 0  47. 5  20. 4  37. 8  24. 0 19. 5  22. 3  27. 0 30. 6  14. 4 12. 5  24. 2  28. 0  30. 8  30. 7 29. 2  26. 8  32. 5  26. 0  34. 2  39. 5  25. 0  19. 5  19. 8  34. 8  27. 0  22. 0  37. 6  21. 2  17. 4]
xx=dtrend(x)
m0=[4 5]
m= armax( xx, m0)
yp1 = predict (xx(1 :60) ,m)
总是有错误提示：
??? Error using ==> armax
Error using ==> pemdecod
   Incorrect number of orders specified.
   You should have m0 = [na nb nc nk],
   where nb and nk are row vectors of length equal to
   the number of inputs.
   For a scalar time series m0 =[na nc].

[ 本帖最后由 eight 于 2007-6-25 10:47 编辑 ]

zhlong

  Incorrect number of orders specified.
   You should have m0 = [na nb nc nk],

  m0=[4 5]

根据提示你的模型阶数设置不正确.

yimengxueer

那应该怎样设置模型参数呢？谢谢

happy

简单修改一下

1. m= armax( xx', m0)

复制代码

xjzuo

按照Happy教授的提示，稍作修改即可．
但结果似乎不好：yp1 =Empty matrix: 0-by-1
由于问题背景没讲清楚，只能你自己修改了．
%%%%%%%%%%%%%%%%%%%%%
x=[22. 0  29. 1  34. 2  36. 0  20. 1  31. 9  19. 8  15. 2  16. 1 17. 9  29. 9  30. 8  31. 7 27. 8  35. 6  15. 2  19. 0 20. 0 20. 0 21. 0  12. 4  16. 3  21. 0  25. 0  33. 2  36. 0  28. 4 40. 0  41. 6  44. 6 37. 0  47. 5  20. 4  37. 8  24. 0 19. 5  22. 3  27. 0 30. 6  14. 4 12. 5  24. 2  28. 0  30. 8  30. 7 29. 2  26. 8  32. 5  26. 0  34. 2  39. 5  25. 0  19. 5  19. 8  34. 8  27. 0  22. 0  37. 6  21. 2  17. 4];
xx=dtrend(x);
m0=[4 5];
xx1=xx';
m= armax( xx1, m0)
yp1 = predict (xx1(60:1) ,m)
%%%%%%%%%%%%%%%

kj17

这么做确实很难模拟的准确。matlab里面有个地震的振动信号，直接这么模拟也不能得到很好的效果