Theory of the circular harmonics expansion method for multiple-optical-fiber system