Efficient algorithms for estimating LSFCs with no aid of SSFCs by taking advantage of the channel hardening effect and large spatial samples available to a massive MIMO base station (BS) are proposed. The LSFC estimates are of low computational complexity and require relatively small training overhead. In the uplink direction, mobile stations (MSs) transmit orthogonal uplink pilots for the serving BS to estimate LSFCs. In the downlink direction, the BS transmits either pilot signal or data signal intended to the MSs that have already established time and frequency synchronization. The proposed uplink and downlink LSFC estimators are unbiased and asymptotically optimal as the number of BS antennas tends to infinity.