Optical frequency comb is an optical spectrum with equispaced frequency lines, used in different applications, for instance, in optical communications. In this application, when an optical spectrum is propagated in single-mode fiber, temperature fluctuations, normal dispersion, and mechanical vibrations can affect the peak power and phase of the comb lines. In consequence, compensation techniques are required to correct the spectral shape, distorted by non-flatness and phase shifts of the spectrum lines. In this research, a method for analyzing optical frequency comb behavior and the spectral shape in terms of phase and intensity is used. This approach is based on fuzzy cellular automata (FCA) to catch up the dynamic of spectrum behavior, fuzzy clustering methods to classify the measured data, and intuitionistic fuzzy entropy to validate the analysis. Two settings to generate optical frequency comb are considered in the experiments: two intensity Mach- Zehnder modulators and mode- locking laser (picosecond pulsed source). Both comb line spectra are propagated through 25km of single- mode fiber. Using the pulse shaper at the optical link output, the spectrum is corrected in flatness and phase shift. The pulse shaper was controlled by a computer, where the proposed method was implemented as a control algorithm based on information obtained by the analysis results. In the experiment, inside the framework of FCA and intuitionistic fuzzy sets (IFSs), the evolution rules 27 (experiment 1) and 184 (experiment 2) were used to find the dynamic behavior of comb spectrum, chosen after a performance comparison with another evolution rules. The comparison was carried out through the calculation of mean and standard deviation of phase shift and power peak, periodicity of FCA, and computational cost. For the first experiment, the change of peak powers was reduced from 4.34dBm to 1.78dBm and the phase shift was minimized from −0.2901rad to −0.2618rad . Instead, the second experiment is observed a correction in flatness from 22.13dBm to 19.89dBm and the standard deviation of phase shift was reduced from 0.00067rad to 0.000471rad.