Determining competitive adsorption isotherms is an open problem in liquid chromatography. Since traditional experimental trial-and-error approaches are too complex and expensive, a modern technique of obtaining adsorption isotherms is to solve the inverse problem so that the simulated batch separation coincides with actual experimental results. This is a typical ill-posed problem. Moreover, in almost all cases the observed concentration at the outlet is the total response of all components, which makes the problem more difficult. In this work, we tackle the ill-posedness with a new regularization method, which is based on the fact that the adsorption isotherms do not depend on the injection profile. The proposed method transfers the original problem to an optimization problem with a time-dependent convection-diffusion equation constraint. Iterative algorithms for solving constraint optimization problems for both the equilibrium-dispersive and the transport-dispersive models are developed. The mass transfer resistance is also estimated by the proposed inverse method. A regularization parameter selection method and the convergence property of the proposed algorithm are discussed. Finally, numerical tests for both synthetic problems and real-world problems are given to show the efficiency and feasibility of the proposed regularization method.