This article shows that the bidiagonal decomposition of many important matrices of q-integers can be constructed to high relative accuracy (HRA). This fact can be used to compute with HRA the eigenvalues, singular values, and inverses of these matrices. These results can be applied to collocation matrices of q-Laguerre polynomials, q-Pascal matrices, and matrices formed by q-Stirling numbers. Numerical examples illustrate the theoretical results.