This paper offers conditions ensuring the existence of solutions of linear boundary value problems for systems of dynamic equations on time scales. Utilizing a method of Moore-Penrose pseudo-inverse matrices leads to an analytical form of a criterion for the existence of solutions in a relevant space and, moreover, to the construction of a family of linearly independent solutions of such problems in a general case with the number of boundary conditions (defined by a linear vector functional) not coinciding with the number of unknowns of a system of dynamic equations. As an example of an application of the presented results, the problem of bifurcation of solutions of boundary value problems for systems of dynamic equations on time scales with a small parameter is considered.