Using temperature distribution as an exteal parameter to change symmetry and measuring frequency spectrum of acoustic resonances in aluminium blocks, we investigate the statistics and dynamics of energy levels in the chaotic billiards. To extract the resonances accurately and eliminate the influence of noise, a filter-diagonalization method for harmonic inversion is used to overcome low resolution of conventional fast Fourier transform method for low quantity factor resonance systems. We present an improved and feasible simulation method to study chaotic characteristic of quantum systems experimentally.