This work is devoted to the study of the Fokker--Planck equation for a stochastic heat equation with an additive $Q$-Wiener noise and non-homogeneous boundary conditions. We explicitly construct the probability density function and establish the associated Fokker--Planck equation by applying the eigenfunction expansion technique. Moreover, the Feynman--Kac formula is used to obtain the probabilistic representation of the solution. The analysis is further extended to cases with multiplicative noise involving nonlocal diffusion operators under homogeneous boundary conditions, as well as the corresponding Kardar--Parisi--Zhang (KPZ) equation. Notably, the evolution of the probability density function for the stochastic heat equation depends critically on the spatial location.