In this paper, we consider the problem of Hardy space decomposition on multiangular domain. By using rational approximation, we achieve that a function f in L p ð∂D σ k Þ can be decomposed into a sum ∑ 2 n k¼1 f ∂D σ k in the sense of L p ð∂D σ k Þ, where f ∂D σ k are the boundary limits of functions in H p ðD σ k Þ.