The invention discloses a joint reconstruction method for multiple dynamic contrast enhancement nuclear magnetic resonance images. According to the method, a conjugate gradient algorithm is combined with a variable density random sampling method to reconstruct the DCE-MRIs in a high probability and efficiency mode from local k-space sampling data. The method includes the following steps of Fourier transform of the images, design of a sampling template, downsampling through the sampling template, inverse Fourier transform, design of wavelet sparse transformation matrix coefficients, design of a constrained energy function, conversion of the constrained energy function into a unconstrained problem through a Laplace operator, acquisition of a solution through the conjugate gradient algorithm based on QUOTE l2,1 l2,1 norms, quantitative evaluation of an obtained reconstruction result of the DCE-MRIs, and difference value contrast of the reconstruction result and the organic images. The reconstruction speed of the images is faster, accuracy and definition are higher, and according to the result, a good contrast mid-value sequence, the high signal to noise ratio, sufficient analysis coverage and rapid data acquisition are achieved.