home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 3.ddi / MMLE3.DI$ / ML_DEMO2.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.9 KB  |  51 lines
  1. %----------------------------------------------------------------------------
  2. %                          DEMONSTRATION PROGRAM-2
  3. %           FOR  MMLE.M  MAXIMUM LIKELIHOOD IDENTIFICATION FUNCTION
  4. %           IDENTIFIES 2*2 TRANSITION MATRIX WITH SINGLE DUMMY STATE
  5. %    - A QUADRATIC COST FUNCTION EXAMPLE TO CHECK NEWTON ALGO CONVERGENCE -
  6. %----------------------------------------------------------------------------
  7.  
  8. format compact,clc
  9. diary ml_demo2.log
  10.  
  11. '-----------------------------------------------------------------------------';
  12. rule=ans;disp(ans)
  13. disp('                          MMLE  DEMONSTRATION  No 2')
  14. disp('                      No dynamics, 2 by 2 mixing matrix')
  15. disp('                       True parameters are [ 1 1 1 1 ]')
  16. %--------------------------------------------------------------------------
  17.  
  18. p2snam='ml_p2ss2';%                   FOR DEMO_2 WE GENERATE SIMULATED DATA
  19. ptrue=[1 1 1 1];
  20. [a,phi,gam,c,d]=eval([p2snam,'(ptrue)']);%             GET SIMULATED SYSTEM
  21.  
  22. uydata=ones(100,4);%                                        CREATE AN INPUT
  23. uydata(1:50,1)=zeros(50,1);
  24.  
  25. randn('seed',0);%                 INITIALIZE RANDOM GENERATOR
  26.  
  27. uydata(:,3:4)=dlsim(phi,gam,c,d,uydata(:,1:2))+randn(100,2)/100;%   SIMULATE
  28.  
  29.  
  30. % Note: IF WE DON'T ADD SOME NOISE OR MODELING ERROR, THE RESPONSE FIT WILL
  31. %       BECOME PERFECT AND THE INNOVATIONS COVARIANCE WILL BECOME ZERO TO
  32. %       CAUSE NUMERICAL PROBLEMS IN THE gg ESTIMATION PHASE.  EITHER ADD
  33. %       NOISE OR DON'T ESTIMATE gg ON SIMULATED DATA. THIS PROBLEM ONLY
  34. %       OCCURS WITH SIMULATED DATA.
  35.  
  36. disp(rule)
  37.  
  38. % IN ADDITION TO uydata AND p2snam (DEFINED ABOVE), MUST CHOOSE VALUES FOR :
  39.  
  40. gg0=eye(2);
  41. p0=[0 0 0 0];
  42. pert=.001;
  43. pidq=1;
  44. pidm=[1:4];
  45. pidf=[1:4];
  46. opt=[0 5 5 5 .02 .01 .001 1];
  47. mmle
  48. diary off
  49. %------------------------------------------------------------- end ml_demo2.m
  50.  
  51.