The primary aim of the research is to provide a comprehensive design process of the piezoelectric device. This simulation is based on Mason equivalent circuit model and transfer matrix method. Besides, the Smith chart, Return Loss and Transducer Loss have been carried out. The theory is applied to a piezoelectric device which excite 1.064μm acoustic wavelength to propagate in the substrate, three different materials of substrate have been chosen in the discussion. By using the Matlab simulation, several important factors including the influence of the bonding layer that affects the performance of piezoelectric transducers have been studied. Furthermore, the designs of matching network will be discussed. Finally, a piezoelectric device’s design procedure to offer a systematic way for the optimum designs of piezoelectric devices is presented.