Abstarct: In this research, an analytical solution is presented for piezoelectric heterogeneous variable wall cylindrical shells that are subjected to mechanical, thermal, and electrical loads in their inner and outer laxyers. The non-homogeneity variations of properties are considered as a power function. Using the first-order electric potential theory, the first-order temperature theory, and the first-order shear deformation theory, and applying the energy method, the governing equations for different boundary conditions at the two ends of the cylinder are derived. The governing differential equation system is a system of ordinary differential equations with variable coefficients, which is solved by using the analytical method of the matched asymptotic expansion method of the perturbation technique, for a cylinder with two clamped heads and with specific potential and temperature. The results of the electric potential distribution, temperature, displacement field, and stress distribution for variable wall cylinders are extracted and the non-homogeneity effects of the properties on the electro-thermo-mechanical behavior of the cylinder are investigated and compared with the results predicted by the finite element method. The results show that the non-homogeneity of properties has a significant effect on the electro-thermo-mechanical behavior of the structure and the present analytical solution can be used as an objective function to optimize and, design a structure with optimal behavior.