A new time series clustering procedure, based on Functional Data Analysis techniques applied to spectral densities, is employed in this work for the detection of stationary intervals in random waves. Long records of wave data are divided into 30-minute or one-hour segments and the spectral density of each interval is estimated by one of the standard methods available. These spectra are regarded as the main characteristic of each 30-minute time series for clustering purposes. The spectra are considered as functional data and, after representation on a spline basis, they are clustered by a mixtures model method based on a truncated Karhunen-Loéve expansion as an approximation to the density function for functional data. The clustering method uses trimming techniques and restrictions on the scatter within groups to reduce the effect of outliers and to prevent the detection of spurious clusters. Simulation examples show that the procedure works well in the presence of noise and the restrictions on the scatter are effective in avoiding the detection of false clusters. Consecutive time intervals clustered together are considered as a single stationary segment of the time series. An application to real wave data is presented.