We present a novel method for modal decomposition of a composite beam guided by a large-mode-area fiber by means of direct far-field pattern measurements with a multi-variable optimization algorithm. For reconstructing far-field patterns, we use finite-number bases of Hermite Gaussian modes that can be converted from all the guided modes in the given fiber and exploit a stochastic parallel gradient descent (SPGD)-based multi-variable optimization algorithm equipped with the D4σ technique in order for completing the modal decomposition with compensating the centroid mismatch between the measured and reconstructed beams. We measure the beam intensity profiles at two different distances, which justifies the uniqueness of the solution obtained by the SPGD algorithm. We verify the feasibility and effectiveness of the proposed method both numerically and experimentally. We have found that the fractional error tolerance in terms of the beam intensity overlap could be maintained below 1 × 10-7 and 3.5 × 10-3 in the numerical and experimental demonstrations, respectively. As the modal decomposition is made uniquely and reliably, such a level of the error tolerance could be maintained even for a beam intensity profile measured at a farther distance.