A response function method using externally applied random force is developed to identify linear or nonlinear systems. The probability distribution function of fluctuations and probability flow in phase space are used as the means for these identifications. It is also demonstrated that anomalous behavior of these quantities can be utilized as the forerunner phenomenon of the catastrophe of unstable systems.