In this paper, “least squares” and recursive methods for simultaneous identification of four nonminimum phase linear, time-invariant FIR systems are presented. The methods utilize the second- and fourth-order cumulants of outputs of the four FIR systems of which the common input is an independent, identically distributed (i.i.d.) non-Gaussian process. The new methods can be extended to the general problem of simultaneous identification of three or more FIR systems by choosing the order of the utilized cumulants to be equal to the number of systems. To illustrate the effectiveness of our methods, two simulation examples are included.