Fluorescence correlation spectroscopy (FCS) is being applied increasingly to study diffusion and interactions of fluorescently labeled macromolecules in complex biological systems. Fluctuations in detected fluorescence, delta F(t), are expressed as time-correlation functions, G(tau), and photon-count histograms, P(k;Delta T). Here, we developed a generalized simulation approach to compute G(tau) and P(k;Delta T) for complex systems with arbitrary geometry, photophysics, diffusion, and macromolecular interactions. G(tau) and P(k;Delta T) were computed from delta F(t) generated by a Brownian dynamics simulation of single-molecule trajectories followed by a Monte Carlo simulation of fluorophore excitation and detection statistics. Simulations were validated by comparing analytical and simulated G(tau) and P(k;Delta T) for diffusion of noninteracting fluorophores in a three-dimensional Gaussian excitation and detection volume. Inclusion of photobleaching and triplet-state relaxation produced significant changes in G(tau) and P(k;Delta T). Simulations of macromolecular interactions and complex diffusion were done, including transient fluorophore binding to an immobile matrix, cross-correlation analysis of interacting fluorophores, and anomalous sub- and superdiffusion. The computational method developed here is generally applicable for simulating FCS measurements on systems complicated by fluorophore interactions or molecular crowding, and experimental protocols for which G(tau) and P(k;Delta T) cannot be computed analytically.