We present an analytic formalism to compute the fluctuating component of the \ion{H}{1} signal and extend it to take into account the effects of partial Lyman-$\alpha$ coupling during the era of cosmic dawn. We use excursion set formalism to calculate the size distribution of randomly distributed self-ionized regions. These ionization bubbles are surrounded by partially heated and Lyman-$\alpha$ coupled regions, which create spin temperature $T_S$ fluctuations. We use the ratio of number of Lyman-$\alpha$ to ionizing photon ($f_L$) and number of X-ray photons emitted per stellar baryons ($N_{\rm heat}$) as modeling parameters. Using our formalism, we compute the global \ion{H}{1} signal, its auto-correlation and power spectrum in the redshift range $10 \le z \le 30$ for the $\Lambda$CDM model. We check the validity of this formalism for various limits and simplified cases. Our results agree reasonably well with existing results from N-body simulations, in spite of following a different approach and requiring orders of magnitude less computation power and time. We further apply our formalism to study the fluctuating component corresponding to the recent EDGES observation that shows an unexpectedly deep absorption trough in global \ion{H}{1} signal in the redshift range $15 <z< 19$. We show that, generically, the EDGES observation predicts larger signal in this redshift range but smaller signal at higher redshifts. We also explore the possibility of negative real-space auto-correlation of spin temperature and show it can be achieved for partial Lyman-$\alpha$ coupling in many cases corresponding to simplified models and complete model without density perturbations.

We present an analytic formulation to model the fluctuating component of the HI signal from the epoch of reionization during the phase of partial heating. During this phase, we assume self-ionized regions, whose size distribution can be computed using excursion set formalism, to be surrounded by heated regions. We model the evolution of heating profile around these regions (near zone) and their merger into the time-dependent background (far zone). We develop a formalism to compute the two-point correlation function for this topology, taking into account the heating auto-correlation and heating-ionization cross-correlation. We model the ionization and X-ray heating using four parameters: efficiency of ionization, $\zeta$, number of X-ray photons per stellar baryon, $N_{\rm heat}$, the spectral index of X-ray photons, $\alpha$, and the minimum frequency of X-ray photons, $\nu_{\rm min}$. We compute the HI signal in the redshift range $10 < z < 20$ for the $\Lambda$CDM model for a set of these parameters. We show that the HI signal for a range of scales $1\hbox{-}8 \, \rm Mpc$ show a peak strength $100\hbox{-}1000 \, \rm (mK)^2$ during the partially heated era. The redshift at which the signal makes a transition to uniformly heated universe depends on modelling parameters, e.g. if $\nu_{\rm min}$ is changed from $100 \, \rm eV$ to $1 \, \rm keV$, this transition moves from $z \simeq 15$ to $z \simeq 12$. This result, along with the dependence of the HI signal on modelling parameters, is in reasonable agreement with existing results from N-body simulations.