The analytical solutions of light beams under two-dimensional harmonic oscillator modulation were obtained with self-similar technique, and the change rules with the quantum number of transmission properties of light beams were discussed by numerical simulation when their transmission distance were taken fixed. It is found that, under two-dimensional harmonic oscillator modulation, the solutions of light beams are Hermite polynomials, whose energy peaks are the highest in the middle of the wave packet which is similar to the bright solitons, and the energy distribution forms a matrix wave packets or a square matrix wave packets in geometry. The total number of wave packets is influenced by the quantum numbers of x and y direction. The amplitude of wave packets energy increases gradually along the positive and negative direction of x and y axis and is diagonal symmetrical on the matrix, x axis or y axis.