A channel estimation method, applied to an orthogonal frequency comprises performing a first number of sampling points inverse fast Fourier transformation (IFFT) operation on a preliminary frequency-domain channel response having a second number of response values to generate a first time-domain channel impulse response (CIR), the second number being greater than the first number; performing a time-domain windowing operation on the first time-domain CIR to generate a second time-domain CIR; performing a smoothing operation on a plurality of second time-domain CIRs of successive time points to generate a smooth time-domain CIR; and performing FFT operation on the smooth time-domain CIR to generate a frequency-domain channel response.