home *** CD-ROM | disk | FTP | other *** search
/ The Datafile PD-CD 5 / DATAFILE_PDCD5.iso / utilities / r / rlab / CTB / lyaptest < prev    next >
Encoding:
Text File  |  1995-11-15  |  805 b   |  56 lines
  1.  
  2. rfile lyap2
  3. rfile lyap
  4. rfile lyapmsf
  5.  
  6. A=rand(4,4);
  7.  
  8. E=eig(A);
  9. Eval=[-4.0,-3.0,-2.0,-5.0];
  10. Anew=conj(E.vec')*diag(Eval)*E.vec;
  11. A=real(Anew);
  12.  
  13. B=rand(4,4);
  14. E=eig(B);
  15. Eval=[-12,-3.4,-5.8,-1.3];
  16. Bnew=conj(E.vec')*diag(Eval)*E.vec;
  17. B=real(Bnew);
  18. B=eye(size(B));
  19.  
  20. C=rand(4,4);
  21. E=eig(C);
  22. Eval=[-0.4,-3.4,-9.5,-6.3];
  23. Cnew=conj(E.vec')*diag(Eval)*E.vec;
  24. C=real(Cnew);
  25. C=eye(size(C));
  26.  
  27. //--------------------------------------------------------
  28.  
  29. printf("%s"," \n");
  30. printf("%s","A,B Solution.\n");
  31.  
  32. X=lyap(A,B);
  33.  
  34. X2=lyap2(A,B);
  35.  
  36. X3=lyapmsf(A,B);
  37.  
  38. Xnorm=norm(X-X2)
  39. X2norm=norm(X-X3)
  40. A*X3+X3*A'+B;
  41.  
  42. //--------------------------------------------------------
  43. printf("%s"," \n");
  44. printf("%s","A,B,C Solution.\n");
  45.  
  46. X=lyap(A,B,C);
  47.  
  48. X2=lyap2(A,B,C);
  49.  
  50. X3=lyapmsf(A,B,C);
  51.  
  52. Xnorm=norm(X-X2)
  53. X2norm=norm(X-X3)
  54. A*X3+X3*B+C;
  55.  
  56.