In this paper, using the equatorial beta-plane approximation of the linear barotropic perturbation equations and introducing the reduced gravitational acceleration, we obtain the analytical solutions of anomalous equatorial ocean stationary wave and give the calculation results of the solutions. Then we compare the results with the modes of complex EOF analysis about abnormal circulation of the real tropical Pacific Ocean and Indian Ocean. The main conclusions are: in the first mode of anomalous equatorial ocean stationary wave, the current disturbance throughout the whole ocean is the half wave, which appears as the consistent zonal flow. The maximum disturbance appears at the middle of the tropical ocean and decays rapidly from equator to north and south, which is restricted in about 2 degree range on both sides of the equator. In the second mode, the current disturbance throughout the whole ocean is the full wave and has the opposite flow direction at east and west of the ocean. The degree of the attenuation of the current disturbance from equator to north and south is as that of the first mode. Anomalous equatorial ocean stationary wave meets the boundary conditions of the east and west coast directing along the longitude. The coefficient is inversely proportional to the square root of the product of the reduced gravity acceleration and the upper water standard depth, which determines the decay rate of anomalous equatorial ocean stationary wave on both sides of the equator. If the square root values take the same, the decay rates are the same. The oscillation frequency of anomalous stationary wave is proportional to the modal number and the square root values, which is inversely proportional to the width of tropical ocean. The modal number is lower and the width is larger, the frequency is lower and the corresponding oscillation period is longer; the first mode of the oscillation period is the longest. Taking every parameters as the typical values and the modal number as one, then taking the width of the equatorial Pacific Ocean and Indian Ocean respectively, the calculation results show that the spatial distribution and interannual variability of the first mode are the same as the corresponding mode of the real abnormal circulation obtaining from the complex EOF analysis; this means that the nature of the first mode above is the anomalous equatorial ocean stationary wave and the anomalous stationary wave is one of the generating mechanism of ENSO and IOD.