JPH042898B2

JPH042898B2 -

Info

Publication number: JPH042898B2
Application number: JP7602285A
Authority: JP
Priority date: 1985-04-10
Filing date: 1985-04-10
Publication date: 1992-01-21
Also published as: JPS61234343A

Description

【発明の詳細な説明】［産業上の利用分野］本発明は核磁気共鳴（NMR）測定方法、特に
２次元（2D）NMR法に関するものである。

［従来技術］近時、新しいNMR測定方法として2DNMRが
注目されている。この2DNMR法は、NMR信号
を２次元スペクトルとして表示することにより、
従来の方法に比べ分解能が向上しスペクトルの解
析が容易になり、該スピン間の相互作用を解明で
きる等の優れた点をもつており、今後ますます発
展するものと考えられている。

第６図にこのように2DNMR法を行うための
NMR装置の一例を示す。図において磁石１が発
生する静磁場内には試料コイル２が配置され、そ
の試料コイル２内部の空間に測定試料が挿入され
る。高周波発振器３から発生する観測核の共鳴周
波数を持つ高周波信号は0゜から360゜まで任意の位
相を選択できる可変移相回路４によつて所定の位
相が与えられた後、増幅器５及びゲート６を介し
て高周波パルスとして前記コイル２へ供給され、
試料に照射される。その高周波パルス照射後コイ
ル２に誘起された共鳴信号は、ゲート７及び受信
回路８を介して復調回路９，１０へ送られる。こ
の復調回路９，１０には前記高周波発振器からの
高周波信号が参照信号として送られるが、その内
の一方は90゜移相回路１１を介して送られるため、
２つの復調回路は90゜位相の異なる２チヤンネル
の検出系CHa，CHbを構成している。この２チ
ヤンネルの検出系から得られた自由誘導減衰信号
は、Ａ−Ｄ変換器１２，１３によつてデジタル信
号に変換されてコンピユータ１４へ送られ、付属
するメモリ１５へ格納される。１６は移相回路、
ゲート６，７及びＡ−Ｄ変換器１２，１３を制御
するパルスプログラマで、試料に照射するパルス
列の順序、パルス幅、各パルスに含まれる高周波
の位相、Ａ−Ｄ変換器１２，１３によるサンプリ
ングのタイミングが予めプログラムされており、
そのプログラムに従つて一連の測定が行われる。

このような装置による従来の2DNMR測定を、
例えば第７図ａに示すような90゜x−t₁−90゜x−t₂
のパルスシーケンスを用いて説明する。

2DNMR法における一般的測定プロセスは第７
図ａに示すように、最初の90゜パルス以前の準備
期間と、展開期間（t₁）と、検出期間（t₂）の３
つの時間領域から成る。準備期間は核の磁化を適
当な初期状態に保つために必要であり、準備パル
ス（最初の90゜パルス）によつて磁化は非平衡の
状態にされ、この状態は展開期間t₁において展開
され、そのt₁における磁化の挙動は、検出パルス
（２番目の90゜パルス）印加後の検出期間t₂におい
て検出される自由誘導減衰信号（FID信号）に位
相及び振幅情報として手渡される。そこで、t₁を
変数として段階的に例えばｎ段階に変化させ、各
段階における測定で２つの検出系CHa、CHbか
ら得られたｍ個ずつのFID信号（FIDa１〜
FIDam及びFIDb１〜FIDbm）の一方を実数部、
他方を虚数部として組合わせた複素集合データS₁
（t₁，t₂）を得れば、このデータ中にはt₂における
磁化の挙動ばかりでなく、t₁における磁化の挙動
の情報も含まれることになる。

次に、検出パルスの位相のみ90゜異ならせた第
７図ｂのパルスシーケンスを用いて上記と同一の
複数のt₁について測定を行い、それにより複素集
合データS₂（t₁，t₂）を得る。

そして、２つの複素集合データS₁（t₁，t₂），S₂
（t₁，t₂）を加算又は減算した後、t₁，t₂について
２重重フーリエ変換することにより所謂Ｐ型
2DNMRスペクトルデータSp（ω₁，ω₂）あるいは
Ｎ型2DNMRスペクトルデータSn（ω₁，ω₂）が得
られる。ω₁，ω₂は夫々t₁，t₂のフーリエ成分であ
る。

第８図は得られる２次元スペクトルの一例を模
式的に表わした図であり、対角ピークP₁，P₂と、
交差ピークP₃，P₄が対称的に現われる。ここ
で、・印は吸収波形ピークを、□印は分散波形ピ
ークを夫々示し、□・印は吸収波形と分散波形が混
じり合つたピークを示す。

［発明が解決しようとする問題点］ところで、スペクトルを解析する場合、ピーク
の波形が純粋な吸収波形であることが好ましい
が、従来このようにして得られる２次元スペクト
ルの対角、交差ずれのピークも吸収波形と分散波
形が混じり合つた波形をしており、解析が困難な
スペクトルになる。

その対策として従来は、絶対値スペクトルを求
め所謂パワースペクトルとして表示する方法が採
用されているが、このパワースペクトル表示には
スペクトルピークの裾が広がるという大きな欠点
がある。

そこで、States等は特殊なデータ処理を施すこ
とにより純粋な吸収波形のスペクトルを得ること
を提案している（Journal of Magnetic
Resonance48p286−，1982）が、このデータ処理
を行うには従来のNMR装置が備えているデータ
処理のためのプログラムを根本的に変更しなけれ
ばならず、装置の価格が上昇することは避けられ
ない。

本発明はこの点に鑑みてなされたものであり、
従来のNMR装置が備えているプログラムの大幅
な変更を必要とせずに、対角ピーク、交差ピーク
のいずれかが純粋な吸収波形他方が分散波形とな
つて解析が容易な2DNMRスペクトルを取得する
ことができ、しかもSN比の向上を図ることので
きる2DNMR測定方法を提供することを目的とし
ている。

［問題点を解決するための手法］この目的を達成するため、本発明にかかる
2DNMR測定方法は、 (a) 準備パルス又はパルス列の照射後展開期間t₁
をおいて検出パルス又はパルス列を照射し、こ
の検出パルス又はパルス列照射後検出期間t₂に
わたつて試料からの自由誘導減衰信号を90゜位
相の異なる２つの検出チヤンネルで検出すると
いうシーケンスを用い、異なつた複数のt₁の値
について測定した複数のFID信号から成る複素
集合データS₁（t₁，t₂）を得ること、 (b) 前記(a)におけるシーケンスと同一のシーケン
スで且つ準備パルス又はパルス列と検出パルス
又はパルス列との間の位相が前記(a)におけるシ
ーケンスの場合と90゜異なるシーケンスを用い、
(a)と同一の複数のt₁について測定した複数の
FID信号から成る複素集合データS₂（t₁，t₂）を
得ること、 (c) 前記集合データS₁（t₁，t₂），S₂（t₁，t₂）を
夫々h₁，t₂について２重複素フーリエ変換して
得られるスペクトルデータをS₁（ω₁，ω₂），S₂
（ω₁，ω₂）とし、該S₁（ω₁，ω₂），S₂（ω₁，ω₂）
を夫々ω₁軸について反転させたスペクトルデ
ータをS₁（−ω₁，ω₂），S₂（−ω₁，ω₂）とした
時、｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）
−S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，
ω₂）−S₂（ω₁，ω₂）｝なるスペクトルデータを作
成すること、より成ることを特徴としている。

［本発明の基本思想］本発明にかかる2DNMR測定方法は、第６図に
示した構成を持つNMR装置で基本的に実施する
ことができる。

以下、第１図に示す流れ図に従い本発明の基本
思想を説明する。

[1] 測定測定によるデータの取得は先に説明した従来の
場合と同様である。第７図ａのシーケンスを用い
た測定により複素集合データS₁（t₁，t₂）が得ら
れ、第７図ｂのシーケンスを用いた測定により複
素集合データS₂（t₁，t₂）が得られる。この２つ
の複素集合データは、メモリ１５内に夫々格納さ
れる。

得られた複素集合データS₁（t₁，t₂），S₂（t₁，
t₂）は、 S₁（t₁，t₂）＝exp（iω₁t₁）ｅ×ｐ（iω₂t₂） (1) S₂（t₁，t₂）＝exp｛ｉ（ω₁t₁＋90）｝ ×exp｛ｉ（ω₂t₂＋90）｝ (2) と表わすことができる。今、説明を簡単にするた
めt₁に関する項の実数部分のみに着目すれば、S₁
（t₁，t₂），S₂（t₁，t₂）は夫々cosω₁t₁exp（iω₂t₂
）
(1′) S₂（t₁，t₂）＝sinω₁t₁exp｛ｉ（ω₂t₂＋90）｝ (2′) と表わされる。

尚、この２つのシーケンスによる測定は、測定
条件が異ならないようできるだけ接近した時刻に
行う必要があり、できれば２つのシーケンスによ
る測定が時分割的に並行して行われることが好ま
しい。

[2] 和データ及び差データの作成上記２つの複素集合データの加算及び減算を行
う。これにより、和データSp（t₁，t₂）及び差デ
ータSn（t₁，t₂）が下式の通り得られる。

Sp（t₁，t₂）＝S₁（t₁，t₂）＋S₂（t₁，t₂）＝cos（ω₂t₂＋ω₁t₁）＋ｉ sin（ω₂t₂＋ω₁t₁） (3) Sn（t₁，t₂）＝S₂（t₁，t₂）−S₁（t₁，t₂）＝cos（ω₂t₂−ω₁t₁）＋ｉ sin（ω₂t₂−ω₁t₁） (4) 尚、和データと差データの作成は、上述の如く
［１］の測定終了後におこなつても良いが、測定
中に並行しておこなつても良い。そうすれば、測
定終了時点で既に和データと差データの作成が終
了していることになる。

[3] 和データ及び差データをt₂，t₁について２重
フーリエ変換 Sp（t₁，t₂）及びSn（t₁，t₂）を先ずt₂について
夫々複素フーリエ変換し、Sp（t₁，ω₂）及びSn
（t₁ω₂）を得る。

Sp（t₁，ω₂）＝exp（iω₁t₁）×｛ａ（ω₂）＋id（ω₂）｝ (5) Sn（t₁，ω₂）＝exp（−iω₁t₁）×｛ａ（ω₂）＋id（ω₂）｝(6) ここでａ（ω₂）は、T⁽²⁾ ₂をω₂軸に関する横緩和
時間、ω⁽²⁾ ₀をω₂軸の観測周波数に対応した基準の
角周波数とした時下式で表わされ、ａ（ω₂）＝T⁽²⁾ ₂／｛１＋T⁽²⁾ ₂ ²（ω₂＋ω⁽²⁾ ₀）²｝吸収波形を与える偶数関数である。

又、ｄ（ω₂）は下式で表わされ、ｄ（ω₂）＝T⁽²⁾ ₂ ²（ω₂＋ω⁽²⁾ ₀）／｛１＋T⁽²⁾ ₂ ²（ω₂＋ω⁽²⁾ ₀）²｝分散波形を与える奇関数である。

Sp（t₁，ω₂）及びSn（t₁，ω₂）を更にt₁について
夫々複素フーリエ変換し、2Dスペクトルデータ
Sp（ω₁，ω₂）及びSn（ω₁，ω₂）を得る。

Sp（ω₁，ω₂）＝｛ａ（ω₁）＋id（ω₁）｝ ×｛ａ（ω₂）＋id（ω₂）｝ (7) Sn（ω₁，ω₂）＝｛ａ（−ω₁）＋id（−ω₁）｝× ｛ａ｛ω₂）＋id（ω₂）｝ (8) ここでａ（ω₁）は、T⁽¹⁾ ₂をω₁軸に関する横緩和
時間、ω⁽¹⁾ ₀をω₁軸の観測周波数に対応した基準の
角周波数とした時下式で表わされ、ａ（ω₁）＝T⁽¹⁾ ₂／｛１＋T⁽¹⁾ ₂ ² （ω₁＋ω⁽¹⁾ ₀）²｝吸収波形を与える偶関数である。

又、ｄ（ω₁）は下式で表わされ、ｄ（ω₁）＝T⁽¹⁾ ₂ ²（ω₁＋ω⁽¹⁾ ₀）／｛１＋T⁽¹⁾ ₂ ²（ω₁＋ω⁽¹⁾ ₀）²｝分散波形を与える奇関数である。

ここまでの過程は先に説明した従来方法と変わ
りがなく、従来はこのようにして得られた2Dス
ペクトルデータSp（ω₁，ω₂），Sn（ω₁，ω₂）に基
づいて2DNMRスペクトルを表示していた。例え
ば、Sp（ω₁，ω₂）の実数部ａ（ω₁）ａ（ω₂）−ｄ
（ω₁）ｄ（ω₂）に基づくスペクトルは先に説明し
た第８図における対角ピークP₁，P₂となつて現
われる。

尚、説明を簡単にするため先に(1)，(2)式の実数
部のみに着目したが、実際には説明が省略された
虚数部についても全く同様な演算処理が施され、
それによつて−ａ（ω₁）ａ（ω₂）＋ｄ（ω₁）ｄ（ω₂）
と表わされる交差ピークP₃，P₄が対角ピークP₁，
P₂と対称的に現われる。

同様にSn（ω₁，ω₂）に基づくスペクトルは第３
図ａのようになる。第８図及び第３図ａ中のいず
れのピークもａ（ω）とｄ（ω）を含むことから、
吸収波形と分散波形が混合された波形になつてし
まうことは避けられない。

本発明は更に以下の処理を付加することを特徴
としている。

[4] 一方のスペクトルデータをω₁に関して周波
数反転例えばスペクトルデータSn（ω₁，ω₂）をω₁に
関して周波数反転し、Sn（−ω₁，ω₂）を求める
と、ａ（ω）が偶関数、ｄ（ω）が奇関数であるこ
とから、Sn（−ω₁，ω₂）は下式のようになる。

Sn（−ω₁，ω₂）＝｛ａ（ω₁）−id（ω₁）｝ ×｛ａ（ω₂）＋id（ω₂）｝ (9) 具体的な周波数反転処理は、例えば以下のよう
に行えば良い。

即ち、スペクトルデータSn（ω₁，ω₂）が第２図
に示すようなＭ行Ｎ列のメモリ領域に格納されて
いるとした時、１行目のスペクトルとＭ行目のデ
ータを交換し、次に２行目のデータとＭ−１行目
のデータを交換し、以下全データについて同様な
交換処理を施す。このようにして反転したデータ
に基づくスペクトルを表示すれば第３図ｂのよう
になる。

尚、２のベキを基数とする高速フーリエ変換
（FFT）用の２次元データの場合には、例えば以
下のような反転操作を施すことが好ましい。

即ち、第２図に示すようなＭ行Ｎ列のメモリ領
域に格納されたスペクトルデータに対し、１行目
のデータは固定しておき、２行目のデータとＭ行
目のデータを交換し、以下、３行目のデータとＭ
−１行目のデータ、４行目のデータとＭ−２行目
のデータ、……、Ｍ／２行目のデータとＭ／２＋
２行目のデータを順次交換する。尚、Ｍ／２＋１
行目のデータは反転の中心になるため１行目のデ
ータと同じく移動しない。

[5] 周波数反転したスペクトルデータを他方の
スペクトルデータと加算又は減算得られた反転データSn（−ω₁，ω₂）を他方のデ
ータSp（ω₁，ω₂）と例えば加算し、結合データＳ
（ω₁，ω₂）を求める。結合データＳ（ω₁，ω₂）は
下式で表わされる。

Ｓ（ω₁，ω₂）＝Sp（ω₁，ω₂）＋Sn（−ω₁，ω₂）＝2a（ω₁）｛１（ω₂）＋id（ω₂）｝ (9) (9)式の実数部が2a（ω₁）ａ（ω₂）と表わされる
ことから分るように、このデータの実数部に基づ
いて表示されるスペクトルピークは第３図ｃにお
いてP₁，P₂で示すように純粋な吸収波形になる。

尚、計算過程は省略するが、(1)、(2)式の虚数部
についても全く同様な演算処理が施され、それに
より2d（ω₁）ｄ（ω₂）と表わされる交差ピークP₃，
P₄が対角ピークP₁，P₂と対称的に現われる。

従つて、本発明によつて得られるスペクトル中
のピークは純粋な吸収波形と純粋な分散波形とな
り、解析が極めて容易となる。しかも、そのため
に新たに加わる処理はデータの反転と加算又は減
算のみであり、従来のNMR装置が備えているプ
ログラムの大幅な変更を必要としない。

更に重要なことは、本発明では(9)式の実数部が
2a（ω₁）ａ（ω₂）となつていることから分るよう
に、スペクトル強度が従来求めていた(7)式で表わ
される2Dスペクトルデータの２倍になつており、
そのためSN比が√２倍改善されるという大きな
効果が伴せて得られることである。従つて、同じ
SN比であれば積算回数を1/2にすることができ、
測定時間を1/2に短縮することが可能となる。

尚、上述した説明では［５］において周波数反
転したスペクトルデータSn（−ω₁，ω₂）を他方の
スペクトルデータSp（ω₁，ω₂）と加算したが、両
者の差を求めても良い。ただし、その場合には対
角ピークが純粋な分散波形となり、交差ピークが
純粋な吸収波形になる。

又、上記説明では［４］の説明においてSn
（ω₁，ω₂）の方を周波数反転した例を示したが、、
Sp（ω₁，ω₂）の方を周波数反転しSn（ω₁，ω₂）と
加算又は減算しても全く同様な結果が得られる。

これらのことをまとめると、本発明において
は、２種の測定により得られる集合データS₁（t₁，
t₂），S₂（t₁，t₂）に基づいて｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝ ±｛S₁（−ω₁，ω₂）−S₂（−ω₁，ω₂）｝ (10) 又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝ ±｛S₁（ω₁，ω₂）−S₂（ω₁，ω₂）｝(11) なるスペクトルデータが最終的に作成されれば良
く、作成にあたつてのデータの加算あるいは減算
処理、フーリエ変換処理及びデータの反転処理の
順序は任意に選択できる。

例えば、上記説明では時間域のデータである集
合データS₁（t₁，t₂），S₂（t₁，t₂）を２重フーリエ
変換により周波数域のスペクトルデータに変換し
た後、ω₁軸に関して周波数反転したが、時間域
のデータについてt₁軸に関するデータ反転をおこ
なつてS₁（−t₁，t₂），S₂（−t₁，t₂）を求め、次に
前記（11）又は（12）式に対応する｛S₁（t₁，t₂）＋S₂（t₁，t₂）｝ ±｛S₁（−t₁，t₂）−S₂（−t₁，t₂）｝(12) 又は｛S₁（−t₁，t₂）＋S₂（−t₁，t₂）｝ ±｛S₁（t₁，t₂）−S₂（t₁，t₂）｝ (13) なる時間域の集合データを作成し、該時間域の集
合データをt₁，t₂について２重複素フーリエ変換
することにより前記(10)、(11)式で表わされる周波数
域のスペクトルデータを求めるようにしても全く
同じ結果が得られる。この場合、比較的時間を要
するフーリエ変換処理が１回ですむので、処理時
間を短くすることができる。

［実施例］第４図は本発明を2DNMR測定法の一種である
COSY法に適用して得られた2DNMRスペクトル
を示す。第５図は前述した［３］のフーリエ変換
までの処理で得られた従来の2DNMRスペクトル
である。第５図の従来スペクトルでは吸収波形と
分散波形が混じり合つて解析が困難なスペクトル
になつていたものが、本発明の適用により第４図
のように単純化され極めて解析のし易いスペクト
ルになつていることが分る。

そのため、第５図において円で囲んだ部分のス
ペクトルでは交差ピークの同定は困難であるのに
対し、第４図では同じ部分で交差ピークが分離し
て見えており、同定を容易におこなうことができ
る。

以上、COSY法を例にとり本発明を説明した
が、本発明はNOESY法をはじめ各種2DNMR法
に全く同様に適用することができる。

［発明の効果］以上詳述した如く、本発明によればスペクトル
のピークを純粋な吸収波形及び純粋な分散波形で
表示することができ、しかもSN比を改善するこ
とのできる２次元核磁気共鳴測定方法が実現され
る。

【図面の簡単な説明】

第１図は本発明の基本思想を説明するための流
れ図、第２図はスペクトルデータの反転処理を説
明するための図、第３図は本発明により得られる
スペクトルを模式的に示す図、第４図は本発明を
2DNMR測定法の一種であるCOSY法に適用して
得られた2DNMRスペクトルを示す図、第５図は
前記［３］のフーリエ変換までの処理で得られた
従来の2DNMRスペクトルを示す図、第６図は
2DNMR法を行うためのNMR装置の一例を示す
図、第７図は2DNMR測定に用いられるパルスシ
ーケンスを示す図、第８図は従来法で得られるス
ペクトルを模式的に示す図である。１：磁石、２：試料コイル、３：高周波発振
器、４：可変移相回路、６，７：ゲート、９，１
０：復調回路、１１：90゜移相回路、１２，１
３：Ａ−Ｄ変換器、１４：コンピユータ、１５：
メモリ、１６：パルスプログラマ。

Claims

【特許請求の範囲】１ (a) 準備パルス又はパルス列の照射後展開期
間t₁をおいて検出パルス又はパルス列を照射
し、この検出パルス又はパルス列照射後検出期
間t₂にわたつて試料からの自由誘導減衰信号を
90゜位相の異なる２つの検出チヤンネルで検出
するというシーケンスを用い、異なつた複数の
t₁の値について測定した複数のFID信号から成
る複素集合データS₁（t₁，t₂）を得ること、 (b) 前記(a)におけるシーケンスと同一のシーケン
スで且つ準備パルス又はパルス列と検出パルス
又はパルス列との間の位相が前記(a)におけるシ
ーケンスの場合と90゜異なるシーケンスを用い、
(a)と同一の複数のt₁について測定した複数の
FID信号から成る複素集合データS₂（t₁，t₂）を
得ること、 (c) 前記集合データS₁（t₁，t₂），S₂（t₁，t₂）を
夫々t₁，t₂について２重複素フーリエ変換して
得られるスペクトルデータをS₁（ω₁，ω₂），S₂
（ω₁，ω₂）とし、該S₁（ω₁，ω₂），S₂（ω₁，ω₂）
を夫々ω₁軸について反転させたスペクトルデ
ータをS₁（−ω₁，ω₂），S₂（ω₁，ω₂）とした時、｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）
−S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，
ω₂）−S₂（ω₁，ω₂）｝なるスペクトルデータを作
成すること、より成る２次元核磁気共鳴測定方法。２前記(a)の測定と(b)の測定は時分割的に並行し
て行われる特許請求の範囲第１項記載の２次元核
磁気共鳴測定方法。３前記(c)において、前記集合データS₁（t₁，
t₂），S₂（t₁，t₂）を夫々t₁軸について反転させた
集合データをS₁（−t₁，t₂），S₂（−t₁，t₂）とした
時、先ず｛S₁（t₁，t₂）＋S₂（t₁，t₂）｝±｛S₁（−t₁，
t₂）−S₂（−t₁，t₂）｝又は｛S₁（−t₁，t₂）＋S₂（−t₁，t₂）｝±｛S₁（t₁，t₂）
−S₂
（t₁，t₂）｝なる集合データを作成し、該集合デー
タをt₁，t₂について２重複素フーリエ変換するこ
とにより前記｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）−
S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，ω₂）
−S₂（ω₁，ω₂）｝なるスペクトルデータを作成す
ることを特徴とする特許請求の範囲第１項乃至第
２項のいずれかに記載の２次元核磁気共鳴測定方
法。４前記(c)において、前記集合データS₁（t₁，
t₂），S₂（t₁，t₂）の和データと差データを求め、
該和データと差データを夫々t₁，t₂について２重
複素フーリエ変換した後前記｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）−
S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，ω₂）
−S₂（ω₁，ω₂）｝なるスペクトルデータを作成す
ることを特徴とする特許請求の範囲第１項乃至第
２項のいずれかに記載の２次元核磁気共鳴測定方
法。５前記(a)の測定と(b)の測定を時分割的に並行し
て行いつつ得られるデータの和と差を求めること
により、測定終了時点で和データSp（t₁，t₂）＝S₁
（t₁，t₂）＋S₂（t₁，t₂）及び差データSn（t₁，t₂）＝
S₁（t₁，t₂）−S₂（t₁，t₂）を得、その後該和データ
及び差データに基づいて前記｛S₁（ω₁，ω₂）＋S₂
（ω₁，ω₂）｝±｛S₁（−ω₁，ω₂）−S₂（−ω₁，ω₂）｝又
は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，ω₂）
−S₂（ω₁，ω₂）｝なるスペクトルデータを作成す
ることを特徴とする特許請求の範囲第１項記載の
２次元核磁気共鳴測定方法。６前記和データSp（t₁，t₂）及び差データSn
（t₁，t₂）の内いずれか一方のデータと、他方の
データをt₁軸について反転させた集合データとの
和又は差を求めることにより、｛S₁（t₁，t₂）＋S₂（t₁，t₂）｝±｛S₁（−t₁，t₂）−
S₂
（−t₁，t₂）｝又は｛S₁（−t₁，t₂）＋S₂（−t₁，t₂）｝±｛S₁（t₁，t₂）
−S₂
（t₁，t₂）｝なる集合データを作成し、該集合デー
タをt₁，t₂について２重複素フーリエ変換するこ
とにより前記｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）−
S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（ω₁，
ω₂）−
S₂（ω₁，ω₂）｝なるスペクトルデータを作成する
ことを特徴とする特許請求の範囲第５項記載の２
次元核磁気共鳴測定方法。７前記和データSp（t₁，t₂）と差データSn（t₁，
t₂）を夫々t₁，t₂について２重複素フーリエ変換
した後前記｛S₁（ω₁，ω₂）＋S₂（ω₁，ω₂）｝±｛S₁（−ω₁，
ω₂）−
S₂（−ω₁，ω₂）｝又は｛S₁（−ω₁，ω₂）＋S₂（−ω₁，ω₂）｝±｛S₁（ω₁
，ω₂）
−S₂（ω₁，ω₂）｝なるスペクトルデータを作成す
ることを特徴とする特許請求の範囲第５項記載の
２次元核磁気共鳴測定方法。