说明：  Algorithm optimization is very good, almost no circulation, It draws on principal component analysis algorithm (PCA), A one-dimensional transfer matrix method to calculate the phonon crystal structure.