In order to identify the origin of phase coexistence at morphotropic tetragonal and rhombohedral boundary in PZT ceramics, the effect of chemical inhomogeneity on phase coexistence region was investigated. Two kinds of PZT ceramics with different chemical homogeneity were prepared by conventional solid state reaction and co-precipitation method. There was coexistence of tetragonal and rhombohedral phase over a wide composition range in PZT ceramics calcined by solid state reaction, and there was also phase coexistence of which region was reduced considerably in sintered samples. And phase coexistence region was not observed in co-precipitated PZT ceramics(within 1 mole%). Therefore compositional fluctuation is considered to be major origin of the phase coexistence at morphotropic phae boundary in PZT ceramics.