I have the following model (which is in this form is not identifiable if the $y$'s are indeed endogenous):
(1) $y_1 = a_0 + a_1y_2 + a_2y_3 + \boldsymbol{Xa} + \boldsymbol{u}$
(2) $y_2 = b_0 + b_1y_1 + b_2y_3 + \boldsymbol{Xb} + \boldsymbol{v}$
(3) $y_3 = c_0 + c_1y_1 + c_2y_2 + \boldsymbol{Xc} + \boldsymbol{w}$
where $\boldsymbol{X}$ is a $n\times k$ matrix of exogenous variables $\boldsymbol{a, b, c}$ which are each $k\times 1$ vectors and the $\boldsymbol{u,v,w}$ are the random noise.
Since every endogenous variable (i.e., $y_i$ for $i=1,2,3$) is a function of the other $y$'s and each equation contains the whole set of exogenous variables (i.e., $X$) one cannot identify the equations given above.
So I checked my empirical data and the theory I have in mind, and I cannot relax the dependency on which of the endogenous variables are related to each other.
My next guess would be that I could check the exogenous variables to find a way to exclude at least 2 of them in each equation and each of these 2 should in the best case not be excluded in any of the other equations.
This would be just amazing, but I won't get out 6 different exogenous variables anyway.
How should I deal with this problem?