creating random dense symmetric matrix with desired condition number

Question

othello nag 2017-5-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/339173-creating-random-dense-symmetric-matrix-with-desired-condition-number

编辑： John D'Errico 2017-5-10

I want to test my code on a large size matrix which should a dense symmetric with specific condition number. I wonder if there are built-in functions in Matlab for doing that?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Goodmanson 2017-5-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/339173-creating-random-dense-symmetric-matrix-with-desired-condition-number#answer_266008

编辑：David Goodmanson 2017-5-8

在 MATLAB Online 中打开

[modified to address the first comment below]

Hi othello, I don't know about built-in functions but there are ways to make your own matrices, e.g.

a = 2*rand(10,10)-1 + 2*rand(10,10)-1;     
a = a+a';              % symmetric with random entries beween -2 and 2
C = 2345;              % desired condition number
[u s v] = svd(a);
s = diag(s);           % s is vector
% ===== linear stretch of existing s
s = s(1)*( 1-((C-1)/C)*(s(1)-s)/(s(1)-s(end))) ;
% =====
s = diag(s);           % back to matrix
b = u*s*v';
cond(b)
ans =  2.3450e+03
%
%
% ===== log spacing example, where also norm(b) = 1  (induced matrix norm)
% s = exp(linspace(0,-log(C),length(s)));
% =====

The condition number is the ratio between the first and last diagonal entries in matrix s (which are sorted), so s is modified to make that ratio equal to C. This also works if the matrix has complex values and is Hermitian. You can also use normal random variables instead of uniform, etc. Seems to work, although something bad is going to happen eventually if you make C too large.

5 个评论
显示 3更早的评论隐藏 3更早的评论

John D'Errico 2017-5-8

在 MATLAB Online 中打开

Nope. A decent idea, but this actually does NOT work. It MAY work some of the time, depending on the value of C. But it is relatively easy to choose a value of C that will make this scheme fail.

So, why do I claim that this MAY fail?

The condition number of a matrix is NOT the ratio of the first and last singular values. It is the ratio of the largest and SMALLEST singular values.

a = 2*rand(10,10)-1;     
a = a+a';
svd(a)
ans =
          4.58080178736462
          4.53709793776632
          3.97533050963259
          3.73031220692586
          2.53840739105403
          2.50126338658107
          2.16399284277328
         0.873311342704311
         0.676947566203514
         0.126414987305412

So quick inspection of the singular values of a tells me that if C=5, your scheme will fail.

C = 5;
[u s v] = svd(a);
s(end,end) = s(1,1)/C;
b = u*s*v';
cond(b)
ans =
      6.76684874288694

As you can see, setting C=5 did not yield a condition number of 5.

The point is, re-setting the singular values CAN be made to work. But you need to be slightly more careful about the algorithm you choose.

David Goodmanson 2017-5-10

编辑：David Goodmanson 2017-5-10

I forgot to mention in comments that I modified the code so now nothing is tbd. Works in every case I tried.

John D'Errico 2017-5-10

编辑：John D'Errico 2017-5-10

Yes. A complete linear or log transformation of the singular values will suffice. Anything like that will work, where all of the singular values move together. That way it does not allow the singular values to effectively change their sequence. Of course, if the small singular value was exactly zero, then there would be a problem. That event would be unlikely, starting from a random matrix. But even that can be solved. +1 now.

请先登录，再进行评论。

creating random dense symmetric matrix with desired condition number

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

5 个评论
显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

creating random dense symmetric matrix with desired condition number

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

5 个评论 显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论