In terms of split-step Fourier transform method, the numerical solution of nonlinear Schr?dinger equation is obtained. The interaction of first-order self-Fourier optical soliton in optical fibers is studied numerically. It is proved that there is not any theoretical error in the inner part of algorithm. The results of numerical simulation show that: 1) The interaction between a couple of first-order self-Fourier solitons are different from first-order standard soliton. It is similar to the interaction between a couple of second-order solitons or a couple of second-order quasi-solitons. The solitons repel each other even more strongly after an initial attraction stage. 2) It is not enough to recognize and classify the characteristics of interaction between two solitons according by the order of soliton. The different pulse width have the different influences to the interaction of solitons or in the same order. 3) A little third-order dispersion may restrict the interaction between the neighboring first-order self-Fourier solitons.