1. A method implemented by a processing unit of a medical imaging analysis system for producing an estimate of a haemodynamic parameter in a dynamic artery/tissue/vein system, on the basis of experimental intensity signals S(t), wherein the method comprises: obtaining a signal S(t), wherein said signal S(t) represents the evolution of a voxel of interest of a plurality of voxels over time t, wherein the plurality of voxels constitute elementary volumes of an organ; deriving, from the signal S(t), a concentration curve c for the voxel of interest; performing, by the processing unit, a step for breaking down a convolution matrix A canonically into singular values in the form A=U·S·V T where S≡diag(σ 1 ,. .. , σ L ) is the diagonal matrix of singular values classified by increasing order, V T is the transpose of a matrix V, and V=(v ij ) and U=(u ij ) are two unitary and real square matrices of dimension L×L, L≧N, N being a defined number of samples, A and c respectively having dimensions L×L and L×1; performing, by the processing unit, at least one iterative step for: producing a pseudo-inverse Ã 1 −1 of A, in the form Ã 1 −1 =V·W l ·U T , where U T designates the transpose of U, and where W l = diag ( 0 , … ⁢ , 0 ︸ l , 1 / σ l + 1 , … ⁢ , 1 / σ L ) ⁢ 0 < l ≤ L ; ⁢ producing d l =Ã −1 ·c; performing, by the processing unit, a step for producing the estimate of the haemodynamic parameter on the basis of an estimate of a quantity of interest {circumflex over (d)}=d l F produced with the iteration l F , where {circumflex over (d)} has dimensions L×1, and l F is positive and less than or equal to L; producing, by the processing unit, the estimate of the haemodynamic parameter in a format suitable for display on a man-machine interface; and displaying, on the man-machine interface, an illustration of the estimate.