A method for estimating angle of arrival of signals in a wireless communication system include a transmitting side and a receiving side, the method comprising: receiving, by the receiving side, a plurality of signals transmitted by the transmitting side, characterised in that the plurality of signals transmitted by the transmitting side further comprises a plurality of training blocks and each one of the plurality of training block comprises a plurality of consecutive Orthogonal Frequency Division Multiplexing (OFDM) symbols, wherein the OFDM symbols comprise pilot symbols at certain subcarriers, the plurality of training blocks are transmitted by the transmitting side using different transmitting beamforming vectors, and the plurality of consecutive OFDM symbols of each one of the plurality of training blocks are received by the receiving side using different receiving beamforming vectors; extracting, by the receiving side, a plurality of pilot symbols from the received signals; estimating, by the receiving side, time-domain channel responses using a compressive sensing algorithm based on the extracted pilot symbols; recovering, by the receiving side, spatial channel responses based on the estimated time-domain channel responses; obtaining, by the receiving side, a plurality of channel taps based on the recovered spatial channel responses; calculating, by the receiving side, a correlation matrix, for each one of the plurality of channel taps, based on the recovered spatial channel responses; performing, by the receiving side, singular value decomposition on the correlation matrix for each one of the plurality of channel taps to obtain a singular value representation of the correlation matrix; determining, by the receiving side, a number of responses caused by different paths for each one of the plurality of channel taps; and estimating, by the receiving side, angle of arrival, for each one of the plurality of channel taps, based on the determined number of responses and the singular value representation of the correlation matrix. receiving, by the receiving side, a plurality of signals transmitted by the transmitting side, characterised in that the plurality of signals transmitted by the transmitting side further comprises a plurality of training blocks and each one of the plurality of training block comprises a plurality of consecutive Orthogonal Frequency Division Multiplexing (OFDM) symbols, wherein the OFDM symbols comprise pilot symbols at certain subcarriers, the plurality of training blocks are transmitted by the transmitting side using different transmitting beamforming vectors, and the plurality of consecutive OFDM symbols of each one of the plurality of training blocks are received by the receiving side using different receiving beamforming vectors; extracting, by the receiving side, a plurality of pilot symbols from the received signals; estimating, by the receiving side, time-domain channel responses using a compressive sensing algorithm based on the extracted pilot symbols; recovering, by the receiving side, spatial channel responses based on the estimated time-domain channel responses; obtaining, by the receiving side, a plurality of channel taps based on the recovered spatial channel responses; calculating, by the receiving side, a correlation matrix, for each one of the plurality of channel taps, based on the recovered spatial channel responses; performing, by the receiving side, singular value decomposition on the correlation matrix for each one of the plurality of channel taps to obtain a singular value representation of the correlation matrix; determining, by the receiving side, a number of responses caused by different paths for each one of the plurality of channel taps; and estimating, by the receiving side, angle of arrival, for each one of the plurality of channel taps, based on the determined number of responses and the singular value representation of the correlation matrix.