This paper is a study on model method for estimating system inputs from vibration responses, which is one of indirect input identification methods in frequency domain. The method has advantages over direct inverse method especially when points of operational inputs are inaccessible so that artificial excitation forces cannot be applied to obtain frequency response functions of the complete system. Procedures of extended modal model method are proposed and checked by numerical experiment. Mechanisms of error propagation, i.e., how errors in modal parameters such as poles nad mode shape vectors affect estimation of the input forces, are illustrated. Then, in order to counteract the error propagation, discrete modal filter approach is taken in this paper to compute the inversion of modal matrix in which the most serious errors seem to be generated. Further, a Reduced form of Modified Reciprocal Modal Vector(RMRMV) is proposed for estimating multiple inputs. It is shown to have smaller orthogonality error than MRMV.