This paper describes perturbation analysis of nonlinear signal generation in apodized BAW devices using a multi-branch LCR equivalent circuit where all necessary simulation parameters are derived by the two-dimensional (2D) scalar wave model. Nonlinear signals are evaluated by adding two voltage sources representing nonlinear strain and electric flux, and their amplitude is determined by currents in the circuit for linear analysis. Simulated second and third harmonic responses agree well with experimental ones. The magnification mechanism of spurious resonances is also discussed.