home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 9.ddi / IDENT.DI$ / PRESENTA.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.0 KB  |  33 lines
  1. function presenta(eta)
  2. %PRESENTA
  3. %
  4. %    presenta(th)
  5. %
  6. %    Subroutine to present
  7.  
  8. %    L. Ljung 10-2-90
  9. %    Copyright (c) 1990 by the MathWorks, Inc.
  10. %    All Rights Reserved.
  11.  
  12. [nr,nc]=size(eta);
  13. T=eta(1,2);
  14. time=eta(2,2:6);
  15. cmd=eta(2,7);
  16. s=1:8;text(1,s)='ARX     ';text(2,s)='ARX2TH  '; text(3,s)='PEM     ';
  17. text(4,s)='IV4     ';text(5,s)='FIXPAR  '; text(6,s)='UNFIXPAR';
  18. text(7,s)='THINIT  ';text(8,s)='IVX     '; text(9,s)='MKETAARX';
  19. %if any(eta(2,7)==[22 25 26 27]),cmd=cmd+10;end
  20. disp(['This ARX-model was created by the command ',text(cmd-30,:),' on ',...
  21. int2str(eta(2,3)),'/',int2str(eta(2,4)),' ',int2str(eta(2,2)*100),' at ',...
  22. int2str(eta(2,5)),':',int2str(eta(2,6))])
  23. disp(['Loss fcn: ' num2str(eta(1,1)) '   The FPE criterion: ' num2str(eta(2,1)) ' Sampling interval ' num2str(eta(1,2))])
  24. disp(['The arx model parameters with their standard deviations given as imaginary parts are'])
  25. [A,B,dA,dB]=th2arx(eta);
  26. i=sqrt(-1);
  27. A=A+i*dA
  28. B=B+i*dB
  29. disp('The covariance matrix of the innovations is')
  30. [p,P,lam]=th2par(eta);
  31. lam
  32.  
  33.