A new time-dependent wave operator approach to the internal eigenstate problems for large matrices