home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 9.ddi / IDENT.DI$ / PRESENTS.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.5 KB  |  45 lines
  1. function presents(eta)
  2. %PRESENTS  presents a parametric model on the screen
  3. %
  4. %    presents(TH)
  5. %
  6. %    This function displays the model TH together with its standard 
  7. %    deviations, loss function and Akaike's Final Prediction Error 
  8. %    criterion (FPE).
  9. %    
  10. %    TH is defined in format for general state-space models, described
  11. %    by HELP TH.
  12.  
  13. %    L. Ljung 10-2-90
  14. %    Copyright (c) 1990 by the MathWorks, Inc.
  15. %    All Rights Reserved.
  16.  
  17. [nr,nc]=size(eta);
  18. T=gett(eta);
  19. nd=eta(1,5);
  20. time=eta(2,2:6);
  21. cmd=eta(2,7);
  22. s=1:8;text(1,s)='MS2TH   ';text(2,s)='MF2TH   '; text(3,s)='PEM     ';
  23. text(4,s)='CANSTART';text(5,s)='FIXPAR  '; text(6,s)='UNFIXPAR';
  24. text(7,s)='THINIT  ';
  25.  
  26. disp(['This matrix was created by the command ' text(eta(2,7)-20,:) ' on '...
  27.  int2str(eta(2,3)) '/' int2str(eta(2,4)) ' ' int2str(eta(2,2)*100) ' at '...
  28.  int2str(eta(2,5)) ':' int2str(eta(2,6))])
  29. if T>0,disp(['Loss fcn: ' num2str(eta(1,1)) '   Akaike`s FPE: ' num2str(eta(2,1)) ...
  30. ' Discrete time model with sampling interval ' num2str(T)])
  31. else disp(['Loss fcn: ' num2str(eta(1,1)) '   Akaike`s FPE: ' num2str(eta(2,1))...
  32.  ' Continuous time model estimated using sampling interval ' num2str(-T)]),end
  33. disp(['The state space matrices with their standard deviations given as imaginary parts are'])
  34. [a,b,c,d,k,x0,da,db,dc,dd,dk,dx0]=eta2ss(eta);
  35. [par,covp,lambda]=th2par(eta);
  36. i=sqrt(-1);
  37. a=a+i*da
  38. b=b+i*db
  39. c=c+i*dc
  40. disp('Strike a key to continue'),pause
  41. d=d+i*dd
  42. k=k+i*dk
  43. x0=x0+i*dx0
  44. lambda
  45.