**4.2.4 Channel estimation**

Spilot = TrDataMod(Ip); Ypilot = TxDataMod(Ip); G = (Ep \* length(Ip))"-1 ... \* ctranspose(sqrt(Ep)\*diag(Spilot)\*ctranspose(F(l:L,Ip))); hHat = G\*Ypilot; TxDataBit = qamdemod(TxDataMod./ (fft(hHat,N)),M);
