Using the scaling equations, quadratures involving polynomials and scaling (or wavelet) functions can be evaluated by linear algebraic equations (which are theoretically exact) instead of numerical approximations. We study two matrices which are derived from these kinds of quadratures. These particular matrices are also seen in the literature of wavelets for other purposes.