Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- n = 100;
- m = 50;
- A = 1/sqrt(m) * randn(m, n);
- K = 10;
- rmin = ones(K, 1);
- rmax = ones(K, 1);
- trials = 1e4;
- for k = 1:K
- for i = 1:trials
- s = randperm(n);
- s = s(1:k);
- sv = svd(A(:,s));
- rmin(k) = min(rmin(k), min(sv));
- rmax(k) = max(rmax(k), max(sv));
- end
- end
- alpha = rmax ./ rmin;
- eps = (alpha - 1) ./ (alpha + 1);
- bar(1:K, eps);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement