Cell migration through 3D extracellular matrices (ECMs) is crucial to the normal development of tissues and organs and in disease processes, yet adequate analytical tools to characterize 3D migration are lacking. The motility of eukaryotic cells on 2D substrates in the absence of gradients has long been described using persistent random walks (PRWs). Recent work shows that 3D migration is anisotropic and features an exponential mean cell velocity distribution, rendering the PRW model invalid. Here we present a protocol for the analysis of 3D cell motility using the anisotropic PRW model. The software, which is implemented in MATLAB, enables statistical profiling of experimentally observed 2D and 3D cell trajectories, and it extracts the persistence and speed of cells along primary and nonprimary directions and an anisotropic index of migration. Basic computer skills and experience with MATLAB software are recommended for successful use of the protocol. This protocol is highly automated and fast, taking <30 min to analyze trajectory data per biological condition.