home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / ORTC.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  497 b   |  20 lines
  1. function [a1cmp] = ortc(a1)
  2. %
  3. % [A1CMP] = ORTC(A1) produces an orthgonal complement of a column
  4. %         type matrix via QR decomposition.
  5. %
  6.  
  7. % R. Y. Chiang & M. G. Safonov 7/85
  8. % Copyright (c) 1988 by the MathWorks, Inc.
  9. % All Rights Reserved.
  10. % ----------------------------------------------------------------
  11. %
  12. [ra1,ca1] = size(a1);
  13. if ra1 == ca1
  14.      a1cmp = zeros(ra1);
  15. else
  16.      [u,r] = qr(a1);
  17.      a1cmp = u(:,ca1+1:ra1);
  18. end
  19. %
  20. % ----- End of ORTC.M ---- RYC/MGS 7/85 %