An algorithm for output spectrum analysis of nonlinear system with correlative sources is proposed. First, the output terms of nonlinear system are analyzed based on the power series. Then, the algorithm is proposed. In the proposed method, different terms with the same frequency are picked out and added together by considering different phases, and the amplitudes of the output spectrum are obtained from the summarized results. Two numerical examples are used to verify the proposed algorithm. It is found that the amplitudes of output spectrum vary with different combinations of phases. Finally, the proposed algorithm is applied to control the output spectrum. Results show that the output amplitude of a specific frequency can be suppressed through changing both the amplitudes and the phases of the input signals.