The finite difference time domain (FDTD) has been popular to characterize wave propagation in biological tissue, which permittivity is generally described by Cole‐Cole dispersive model.The FDTD scheme by using fractional derivative to time‐domain polarization equation results in a fast and computational efficient algorithm. The accuracy obtained is, however, not generally as good as thoes of z‐transform and auxiliary differential equation methods. In this article, polarization piece‐linearity and nonlinear optimization techniques are proposed and utilized on the polarization fractional derivative appoach, which improved the accuracy while retains the high speed and efficiency advantage. © 2012 Wiley Periodicals, Inc. Microwave Opt Technol Lett 54:888–891, 2012; View this article online at wileyonlinelibrary.com. DOI 10.1002/mop.26698