Test matrices with special properties are important tools for testing
numerical algorithms and software. In this paper we describe how to ea
sily generate random test matrices that are particularly suited for te
sting regularization algorithms. Our algorithm is based on the singula
r value decomposition (SVD) of bidiagonal matrices with positive eleme
nts.