Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function B = generateQR(n,s)
- function A = integermatrix(n)
- A = rand(n);
- for i=1:n
- for j=1:n
- A(i,j) = round(5*A(i,j));
- end
- end
- end
- A = triu(integermatrix(n));
- function Pi = permute(i,pi)
- Pi = eye(n);
- Pi(i,i) = 0;
- Pi(pi,pi) = 0;
- Pi(pi,i) = 1;
- Pi(i,pi) = 1;
- end
- for j=1:s
- i = round((n-1)*rand(1)) + 1;
- pi = round((n-1)*rand(1)) + 1;
- Pi = permute(i,pi);
- A = Pi*A;
- end
- B = A;
- end
Add Comment
Please, Sign In to add comment