In this paper, we develop the Littman-Stampacchia-Weinberger duality approach to obtain global W^1,p estimates for a class of elliptic problems involving Leray-Hardy operators and measure sources in a distributional framework associated with a dual formulation with a specific weight function.