A sequence of special functions in Hardy space [Formula: see text] are constructed from Cauchy kernel on unit disk 𝔻. Applying projection operator of the sequence of functions leads to an analytic sampling approximation to f, any given function in [Formula: see text]. That is, f can be approximated by its analytic samples in 𝔻s. Under a mild condition, f is approximated exponentially by its analytic samples. By the analytic sampling approximation, a signal in [Formula: see text] can be approximately decomposed into components of positive instantaneous frequency. Using circular Hilbert transform, we apply the approximation scheme in [Formula: see text] to Ls(𝕋2) such that a signal in Ls(𝕋2) can be approximated by its analytic samples on ℂs. A numerical experiment is carried out to illustrate our results.