Guest User

Untitled

a guest
May 26th, 2018
92
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.48 KB | None | 0 0
  1. function B = generateQR(n,s)
  2.  
  3.  
  4. function A =  integermatrix(n)
  5.     A = rand(n);
  6.     for i=1:n
  7.         for j=1:n
  8.             A(i,j) = round(5*A(i,j));
  9.         end
  10.     end
  11. end
  12.  
  13. A = triu(integermatrix(n));
  14.  
  15. function Pi = permute(i,pi)
  16.     Pi = eye(n);
  17.     Pi(i,i) = 0;  
  18.     Pi(pi,pi) = 0;
  19.     Pi(pi,i) = 1;
  20.     Pi(i,pi) = 1;
  21. end
  22.  
  23.  
  24. for j=1:s
  25.     i = round((n-1)*rand(1)) + 1;
  26.     pi = round((n-1)*rand(1)) + 1;
  27.     Pi = permute(i,pi);
  28.     A = Pi*A;
  29. end
  30.  
  31. B = A;
  32.  
  33.  
  34. end
Add Comment
Please, Sign In to add comment