In this paper, we propose a statistical approach to reconstruct the brain neuronal activity based only on recorded EEG data. The brain zones with the strongest activity are expressed at a macro level by a few number of active brain dipoles. Normally, for solving the EEG inverse problem, fixed dipole locations are assumed, independently of the different stimuli that excite the brain. The proposed particle filter (PF) framework presents a shift in the current paradigm by estimating dynamic brain dipoles, which may vary from one location to another in the brain depending on internal/external stimuli that may affect the brain. Also, in contrast to previous solutions, the proposed PF algorithm estimates simultaneously, the number of the active dipoles, their moving locations and their respective oscillations in the three dimensional head geometry.