Polarized differential phase and modulation ratios were obtained for 1,6-diphenyl-1,3,5-hexatriene (DPH) in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles by using multifrequency phase fluorometry. Data were analyzed in terms of both empirical sums-of-exponentials modeling and directly in terms of the orientational distribution functions. The orientational analysis model was used to recover the angular distribution of DPH and the rotational diffusion coefficient in the various membrane model systems throughout the phase transition. A global analysis methodology was utilized to obtain an internally consistent set of parameters that fit all of the data simultaneously. The rank order parameters (P2) and (P4) were extracted from the experimental data, and the angular distribution functions of DPH were calculated. When the time-zero anisotropy (r0) of several sets of data taken at various temperatures were linked in a single global analysis, better recovery of the fourth rank order parameter (P4), diffusion constant D, and r0 was obtained with respect to the unlinked analysis. From these recovered values, a detailed picture concerning the orientational distribution of DPH in membranes as a function of temperature was obtained. The results suggest that a single population of DPH molecules was present in the bilayers with their orientational distributions dependent upon the physical state of the membranes in the pure phases. During the phase transition, a superposition of two populations corresponding to the population of the pure phases was present. As the temperature increased in the transition region, one population was increasing at the expense of the other.