The method of “auxiliary variables” (collective coordinates) used by Bogolyubov and Zubarev is formulated in a somewhat different manner. Improvements are made in two points: 1) The Hamiltonian is given in an Hermitian form and 2) the number of the auxiliary variables can be limited, if necessary. We calculate the ground state energy and the excitation energy spectrum of the spinless bosons taking account of phonon-phonon interactions. The expression for the ground state energy is shown to be equivalent to the one derived by Bogolyubov and Zubarev.