The present invention particularly relates to a method for analyzing at least one sample comprising at least one characterization species and a reference species using nuclear magnetic resonance (NMR), comprising one or more antennas, Obtaining (1330) one or more complex free induction decay (FID) signals S (t) and obtaining a FID spectrum S (ω) by applying a Fourier transform to the signal S (t) (1323); The obtained FID spectrum S (ω) includes two portions (UFR, DFR) each extending from the resonance frequency (F0Ref) of the reference species to each side of F0Ref, and the frequency of the species to be characterized Is located in a part of the spectrum taken from the two parts (UFR, DFR) and a signal Sref (t) of the reference species Modeling from the real and complex parts of at least one complex FID signal S (t) (1321) and including only the reference species by applying a Fourier transform to the modeling of the reference species signal Sref (t) Step (1322) of obtaining a spectrum Sref (ω) of a reference species and correcting by replacing the corresponding part of the spectrum Sref (ω) with the part of the FID spectrum S (ω) that does not contain the species to be characterized Obtaining the FID spectrum {tilde over (1)} (1324), applying the inverse Fourier transform to {circumflex over (2)} (1325) to obtain the modified FID signal {circumflex over (3)}, and Calculating a module (1326).本発明は、特に、少なくとも1つの特性決定対象の種および参照種を含む少なくとも1つの試料を、核磁気共鳴(NMR)を使用して分析する方法であって、 1つまたは複数のアンテナによって、1つまたは複数の複素自由誘導減衰(FID)信号S(t)を取得するステップ(1330)と、 信号S(t)にフーリエ変換を適用することによってFIDスペクトルS(ω)を得るステップ(1323)であり、この得られるFIDスペクトルS(ω)が、それぞれ参照種の共鳴周波数(F0Ref)からF0Refのそれぞれの側に延びる2つの部分(UFR、DFR)を含み、特性決定対象の種の周波数が、前記2つの部分(UFR、DFR)から取られたスペクトルの一部分に位置するステップと、 参照種の信号Sref(t)を、前記少なくとも1つの複素FID信号S(t)の実数部および複素部からモデル化するステップ(1321)と、 参照種の信号Sref(t)のモデル化にフーリエ変換を適用することによって参照種のみを含む参照種のスペクトルSref(ω)を得るステップ(1322)と、 スペクトルSref(ω)の対応する部分で、特性決定対象の種を含まないFIDスペクトルS(ω)の部分を置換することによって、修正FIDスペクトル【数1】を得るステップ(1324)と、 逆フーリエ変換を【数2】に適用して(1325)、修