Nonlinear intermodulation distortion inside duplexers is a limitation for the sensitivity of mobile devices receivers. A perturbation approach for non linear harmonic acoustic generation along the propagation of surface acoustic wave is derived. It is used to extend the P matrix model to non linear generation of harmonics and mixing products. Interface of the new model to harmonic balance simulators is discussed. Results are presented showing the validity of the approach.