Soil heterogeneity in an infinite unsaturated slope increases the difficulty of analysing its stability. A common method of modelling the spatial variability in the slope stability analysis is to use random fields together with Monte Carlo simulation, but this requires an impressive amount of computational efforts. To reduce the computing time, this paper develops an analytical approach to study the stability of infinite heterogeneous unsaturated slopes. The analytical method provides a newly time-efficient way of analysing such kinds of slopes. It has been used to investigate the stability of four different kinds of unsaturated soil slopes. Meanwhile, the influences of the correlation between the strength parameters and the autocorrelation length of the strength parameters are investigated. The analysed results show that the suction has influences on the distribution of the FOS for both homogeneous slopes and heterogeneous slopes. In addition, the correlation coefficient between the strength parameters and the autocorrelation length of the strength parameters affects the mean and standard deviation of FOS for homogeneous slopes and the probability of failure for heterogeneous slopes.