How can I generate random sparse stable matrices ?

I need to generate random matrices that are sparse and stable and control the sparsity level, I tried using matlab functions sprand and sprandn and use conditional "if" to check if the generated matrices have negative eigenvalues, but the generated matrices are often unstable and do not converge to the steady state.

 Accepted Answer

There is simply no way that use of sprand or sprandn will yield matrices with all positive eigenvalues with any viable frequency.
However, you might use those tools to generate the random sparse matrices, then shift the eigenvalues. Gershgorin disks will suffice to know how much of a shift is necessary. The result may not be as purely "random" as you might want, but it very possibly is as good as you will get.

1 Comment

Thanks. with your answer and another answer posted on the webpage http://www.mathworks.com/matlabcentral/newsreader/view_thread/26871 I am able to generate random sparse stable matrix.

Sign in to comment.

More Answers (0)

Categories

Find more on Sparse Matrices in Help Center and File Exchange

Asked:

on 15 Feb 2014

Commented:

on 17 Feb 2014

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!