(57)【特許請求の範囲】
【請求項１】
分散信号u、吸収信号vの測定量から誤差磁場δを求め、そのδ値にフィードバック・ゲインを乗算して磁場補正値とするNMR磁場安定化方法であって、
dv／dt、u、vに加え，du／dtを引数として持つ関数と、NMR磁場安定化装置で測定したuとvと前記uとvを統計的に処理して求めたdv／dtとdu／dtを用い、前記関数の値を前記δ値とすることを特徴とするNMR磁場安定化方法。
【請求項２】
請求項１において、誤差磁場δをu／（u2＋v2）と（u×dv／dt−v×du／dt）／（u2＋v2）の線形結合で求めることを特徴とするNMR磁場安定化方法。
【請求項３】
請求項１において、前記分散信号uおよび前記吸収信号vにおいて、定常状態でのuの振幅を正規化係数として保存し、uおよびv信号の値を前記正規化係数で割算することで、uおよびvの振幅を正規化するとともに、前記uおよびvのS／Nにより請求項１のフィードバック・ゲインを調整することを特徴とするNMR磁場安定化方法。
【請求項４】
ＮＭＲ用磁石に挿入されたプローブから交流信号であるNMR信号を測定し、前記磁石の磁場強度を測るとともに、前記磁場強度を予め指定された値に補正するNMR磁場安定化装置において、
前記交流信号の受信部に１チャンネルの交流信号入力を受けるディジタル検波器を有し、前記ディジタル検波器はその内部にアナログ／ディジタル変換を行うＡＤ変換部と1チャンネルのディジタル信号を２チャンネルに分岐し直交検波を行い分散信号uと吸収信号vとを得る直交検波部を持つことで、前記分散信号uと前記吸収信号vを別々にアナログ／ディジタル変換した場合に両信号のＡＤ変換部の差により生ずる磁場補正の誤差を低減し、
前記直交検波後の分散信号uおよび吸収信号vの各々のデータ量が、予め決めた帯域幅を再現するためにナイキスト定理が要求する分散信号uおよび吸収信号vの各々のデータ量の８倍以上であり、線型回帰の統計的データ処理により、前記分散信号uおよび吸収信号vの時間変化率を求め、
前記交流信号受信部のディジタル検波器において、アナログ／ディジタル変換部（ADC）およびディジタル直交検波部に同一ビット長で周波数を記述する素子を用い、前記交流信号の送信部にも前記ADCおよびディジタル直交検波部と同じ周波数記述ビット長で周波数を表すディジタル信号合成器（DDS）を用いて、前記交流信号の送信部と受信部の間で位相の同期を確保し、
前記ADCとディジタル検波器およびDDSの各部分は、各々の動作タイミングの基準となるクロック信号を入力として持ち、前記各部分のクロック信号を一つの信号発生器の出力から得ることと、前記ディジタル検波器は、誤差磁場δをu／(u2+v2)と(u×dv／dt−v×du／dt)／(u2+v2)の線形結合で求める、ことを特徴とするNMR磁場安定化装置。
【請求項５】
請求項４において、
前記DDSが出力する交流信号の周波数を時間に比例して変えながらNMR信号を測定し、NMR信号の電圧変動から磁石の磁場強度を測り、
前記磁石の磁場強度が予め指定された値から外れた場合、ロックオンからロックオンする周波数を決める周波数掃引に切り替え、DDSが出力する交流信号の周波数を時間に比例して変えながらNMR信号を測定して磁石の磁場強度を測り、該磁場強度と前記予め指定された磁場強度との差をロックオン値にオフセットとして加えてロックオンに復帰する、ことを特徴とするNMR磁場安定化装置。
【請求項６】
請求項５において、吸収信号の最大値を正規化係数として保存し、吸収信号の最大値を前期正規化係数で割算することで、uおよびvの振幅を正規化するとともに、前記uおよび前記vのS／Nによりフィードバック・ゲインを調整することを特徴とするNMR磁場安定化装置。
【請求項７】
静磁場を発生するとともに試料挿入孔を持つ磁石と、試料挿入孔に挿入されるプローブと、プローブと電気的に接続され交流信号を送受信するNMR磁場安定化装置と、NMR磁場安定化装置を制御するコントローラと、コントローラとデータ伝送手段により接続されユーザーインタフェースを提供するコンピュータとを有するNMR装置において、
前記NMR磁場安定化装置に、請求項４乃至請求項６の何れか１項に記載のNMR磁場安定化装置を用いることを特徴とするNMR装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は核磁気共鳴分光器(以下NMR分光器)におけるＮＭＲ磁場安定化装置（以下、ロックシステムと呼ぶ）に関する。
【背景技術】
【０００２】
近年NMR分光器は、有機物や固体から生体高分子に至る幅広い試料の測定に使われている。NMR分光器は、静磁場中にある試料が持つ共鳴周波数を測定する装置であり、非特許文献１に記載のように、１ppb／hr以下の時間変動率を持つ安定した静磁場が要求される。安定した静磁場を作るために、NMR分光器の殆どが超伝導電磁石を装備しているが、それによっても時間変動率が１ppb／hr以下の安定性を達成することは困難である。そのため、フィードバック制御により静磁場を安定化するロックシステムが、NMR分光器の性能を向上させる上で重要な要素となっている。
【０００３】
非特許文献２に、ロックシステムに関する従来技術（以下、従来技術１）が記載されている。ロック送信器からNMR装置用磁石中に置かれた試料にラジオ周波数を持つ交流磁場を照射し、試料中の磁化ベクトルを励起させる。励起された磁化ベクトルの運動を検出用アンテナから受信し、受信回路を用いて検波する。検波の際に使われる参照信号の位相を調整することにより、分散信号（以下，ｕ）と呼ばれる成分を得る。次に、ｕに予め決めた常数ｇを乗算し、ロック電源に渡す。ロック電源はｇ×ｕ（以下，補正信号ｓ）に比例する電流を試料の周辺に設けられたコイル（以下，ロックコイル）に追加し、上記試料の周辺にｓに比例する補正磁場を追加する。補正磁場を追加した後で得られるｕが０に近づくように、比例常数ｇを決めることにより静磁場を安定化する。
【０００４】
以上で述べたように、従来技術１は分散信号ｕに比例してロックコイルに流れる電流量を制御することで静磁場の安定化を図るので、ロック補正器が、例えばアナログ回路だけを用いても、簡単に実現できる。尚、上記ロック検波器では分散信号ｕと位相が90度異なる吸収信号ｖを求めることもできるが、補正信号ｓの算出には使われない。
【０００５】
特許文献１に、ロックシステムに使われる従来技術２が記載されている。ロック補正器から出力される補正信号ｓを算出するために、分散信号ｕと共に吸収信号ｖを用い、速い磁場変動を検出するために磁石周辺に変動磁場検出用の磁場センサを設け、補正信号ｓと磁場センサで検出した変動磁場量を演算して磁場補正を行う点で、従来技術１と異なる。
【０００６】
従来技術２に置いて、補正信号ｓはｕ／ｖと（１／ｖ）（dｕ／dｔ）に重みを付けた組み合わせである。これによれば、擾乱の周波数が低い場合に、ｕ／ｖは静磁場の擾乱による周波数オフセットｄｆと完全に比例し、その比例常数は１／Ｔ２である。Ｔ２は信号ｕおよびｖの減衰率に関わる時間であり、横緩和時間と呼ばれる。
【０００７】
従来技術２によれば、Ｔ２は制御ループの中で求めることができるため、ロック制御部でｄｆ＝（１／Ｔ２）（ｕ／ｖ）を算出し、補正信号ｓとして用いることもできる。ｄｆは、また、（１／Ｔ２）（ｕ／ｖ）＋（１／ｖ）（dｕ／dｔ）とも等しい。この等式は正確であり、どんな条件にも関連しない利点と、雑音成分を含む不利な点を併せ持つ。従って、上記式は擾乱振幅が比較的に大きい場合の組調整に適する。
【０００８】
また、従来技術２によれば、磁石の外からの速い磁場変動は磁石外に設けた磁場センサで検出する。速い磁場変動は、原理的にｄｆ＝（１／Ｔ２）（ｕ／ｖ）＋（１／ｖ）（dｕ／dｔ）により検出できるが、この量はノイズを多く含む。この磁場センサは磁石外に置かれ、磁石の中心に位置するNMR試料と、磁場センサとの間には磁性を持つ複雑な構造物が置かれる。従って、磁場センサが検出する磁場変動量とNMR試料が受ける磁場変動量との間には、大きさと分布において差が生じる。この差は数学的なデータ処理により補正される。
【０００９】
従来技術２は、ｄｆと補正信号ｓとが比例関係であり単調性が良いため、従来技術１に比べて少ない制御ループ数で磁場を安定値に戻すことができる。また、速い磁場変動の補正ができる点も優れている。
【００１０】
【特許文献１】日本特許第２５０４６６６号
【非特許文献１】荒田洋治著  「NMRの書」丸善  ２０００年
【非特許文献２】Benedict W. Bangerter、認ield/frequency lock monitor for signal averaging with high resolution NMR spectrometers煤AReview of Scientific Instruments、46巻、617項、1975年
【非特許文献３】Macro Villa, Feng Tian, Pacifico Cofrancesco, Josef Halamek and Miroslav Kasal、滴igh-resolution digital quadrature detection煤AReview of Scientific Instruments、67巻、2123項、1996年
【発明の開示】
【発明が解決しようとする課題】
【００１１】
従来技術１は、速い磁場変動の補正において課題を残している。従来技術１で用いている分散信号ｕと磁場変動量との関係は、分散信号と吸収信号の時間変化率dｕ／dｔとdｖ／dｔが０の定常状態で導かれているためである。そのため、従来技術１によるロックシステムは、超電導電流の減衰による磁場ドリフトを補正する目的には適しているものの、dｕ／dｔとdｖ／dｔが大きい値を持つ速い磁場変動を補正することは困難である。
【００１２】
従来技術２は、速い磁場変動の検出において課題を残している。従来技術２によれば、ロックシステムに新たに磁場センサを追加する必要があり、その校正作業も必要となる。また、磁場センサで検出した磁場変動量からNMR試料が受ける磁場変動量に換算するデータ処理が必要となる。更に、前記ロック補正器から出力される補正信号ｓと磁場センサによる補正量との重みを調整する操作も必要になる。
【００１３】
本発明の目的は、上記従来技術の問題点に鑑み、新たな装置の追加や校正作業が低減でき、かつ速い磁場変動を補正することのできるNMR磁場安定化装置及びそれを用いた高性能のNMR装置を提供することにある。
【課題を解決するための手段】
【００１４】
上記課題を解決するために、本発明のNMR磁場安定化方法は、分散信号u、吸収信号vの測定量から誤差磁場δを求め、そのδ値にフィードバック・ゲインを乗算して磁場補正値とするNMR磁場安定化方法であって、dv／dt、u、vに加え，du／dtを引数として持つ関数と、NMR磁場安定化装置で測定したuとvと前記uとvを統計的に処理して求めたdv／dtとdu／dtを用い、前記関数の値を前記δ値とすることを特徴とする。
【００１５】
本発明によるNMR磁場安定化装置は、前記プローブと電気的に接続され交流信号を送受信する交流信号送信部および交流信号受信部とを有し、NMR信号を測定し前記磁石の磁場強度を測る探索モードと磁石の磁場強度を予め指定された値に保つ補正モードを持つものであって、前記交流信号受信部に１チャンネルの交流信号入力を受けるディジタル検波器を有し、前記ディジタル検波器はその内部にアナログ／ディジタル変換を行うＡＤ変換部と1チャンネルのディジタル信号を２チャンネルに分岐し直交検波を行い分散信号と吸収信号とを得る直交検波部を持ち、前記分散信号と前記吸収信号を別々にアナログ／ディジタル変換した場合に両信号のＡＤ変換部の差により生ずる磁場補正の誤差を低減することを特徴とする。
【００１６】
また、NMR磁場安定化装置は、直交検波後のデータ量が、予め決めた帯域幅を再現するためにナイキスト（Nyquist）定理が要求するデータ量の８倍以上とする。
【００１７】
また、信号の時間変化率を求めるために、線型回帰の統計的データ処理を用いる。また、分散信号ｕ、吸収信号ｖの測定量から誤差磁場δを求める関数として、du／dtとdv／dt、uとvの因数を持つ関数を用いる。また、誤差磁場δをu／（u2＋v2）と（u×dv／dt−v×du／dt）／（u2＋v2）の線型結合で求める。また、振幅を正規化し、S／Nによりフィードバック・ゲインを調整する。また、受信器のADCおよびディジタル検波器と位相の同期を確保するため、送信にディジタル信号合成器（DDS）を用いる。また、単一周波数源からクロックを上記DDS、ADCおよびディジタル検波器に供給することで、位相誤差を45度／スキャン以下に抑える。また、周波数掃引で共鳴周波数探索を行う。また、ロック外れの発生時に、観測を続行しながら周波数掃引でロックオンに復帰する。
【発明の効果】
【００１８】
本発明によれば、NMR磁場安定化装置に対する装置追加と校正作業が低減でき、かつ速い磁場変動を補正することができる効果がある。また、本NMR磁場安定化装置を用いた高性能のNMR装置を提供することができる。
【発明を実施するための最良の形態】
【００１９】
本発明のNMR磁場安定化装置は、速い磁場変動を補正でき、かつ装置の追加や校正作業を最小限に抑えることが出来る。以下、本発明の実施の形態を図面を参照しながら説明する。
【００２０】
図１はNMR装置の概念的な構成図を示す。NMR装置は高周波の照射と検出を行うプローブ６、静磁場を発生させる磁石７、ロック制御に使われるロックコイル８及び試料管９を備えている。また、ユーザーインタフェースであるユーザーコンピュータ２０１、NMR装置の制御を担当するコントローラ２０２、NMR信号の観測を担当する観測システム２０３、磁場安定化を担当する磁場安定化装置であるロックシステム２０４を有している。
【００２１】
ユーザーは、プローブ６を磁石７の中に挿入する。観測対象の核とロック制御用の核とを含む物質が入った試料管９は、プローブ６中に挿入される。プローブ６を磁石７の中に入れた後に、試料管９をプローブ６中に挿入することが好ましいが、試料管９を予め挿入したプローブ６を磁石７中に挿入することも可能である。
【００２２】
続いて、ユーザーコンピュータ２０１とコントローラ２０２を通じて、ロックシステム２０４を稼動させる。ロックシステム２０４は、プローブ６に予めユーザーが決めた周波数の電気信号を送る。プローブ６は電気信号を同じ周波数の磁場に変え試料管９に照射する。照射される磁場の周波数は試料管９中のロック制御用核が持つ共鳴周波数に近いように決められているため、試料管９中のロック制御用核は照射された磁場と共鳴し回転する。ロック制御用核の回転はプローブ６で電気信号として検出され、ロックシステム２０４に送られる。ロックシステム２０４は、プローブ６で検出された前記電気信号を用いて試料管９が受ける磁場強度を測定し、磁場安定化制御を行う。
【００２３】
ロックシステム２０４により磁場が一定変動範囲内、典型的には１Hz内、に保たれる状態になる。ユーザーは、磁場安定性が保たれたことを確認した後に、ユーザーコンピュータ２０１とコントローラ２０２を通じて、観測システム２０３を稼動させる。
【００２４】
図２は本発明の一実施例によるロックシステムの詳細な構成図を示す。コントローラ１はロックシステム２０４の制御を行うと共に、ユーザーコンピュータ２０１とロックシステム２０４とのインタフェースの機能も持つ。すなわち、コントローラ１は、ユーザーコンピュータ２０１からの入力により、ロックシステム２０４のロック送受信器２とディジタル検波器３を設定し、ディジタル検波器３から出力されるロック信号データをユーザーコンピュータ２０１に送信する。
【００２５】
コントローラ１は、図１記載のコントローラ２０２と同じであっても良いし、異なってもよい。コントローラ１とコントローラ２０２が一致する場合は、装置構成が単純になる反面、コントローラの処理負荷が大きくなる。
【００２６】
ロック送受信器２とコントローラ１とのインタフェースはCPU１１である。CPU１１は、また、ディジタル信号合成器（DDS：Direct  Digital  Synthesizer）１４の設定、送信ゲート１９と送受信切替回路２１とディジタル直交検波を行うディジタル直交検波部（FPGA：Field  Programmable  Gate  Array）４２の動作トリガーの役割も行う。ロックシステム全体の同期は周波数合成器（FRQSYN：Frequency  Synthesizer）１３により保たれる。
【００２７】
FRQSYN１３は外部の周波数源１２から種クロックを受け、ロックシステムの動作に必要な周波数の信号を生成する。DDS１４とFPGA４２には素子の動作クロックを供給する。また、送信部ミキサ１６と１６'、受信部ミキサ３２に、局部発振信号を供給する。局部発振信号の周波数はfrf−fifであり、frfはプローブ６から照射される高周波信号の周波数、fifはDDS１４が出力する中間周波数信号の周波数である。
【００２８】
DDS１４は、90度異なる位相を持つ2つの中間周波数信号、IF−QとIF−Iを出力する。IF−Q信号とIF−Iは送信部ミキサ１６と１６'で局部信号LOと混合され、合成器(combiner)１７で合成される。合成器１７の出力信号はfrfの周波数を持つ。局部信号LOは2つに分けられ、その1つは90度位相変換器１５を経由し送信部ミキサ１６に入る。
【００２９】
位相変換器１５とミキサ１６、１６'、合成器１７の組合せはイメージ・リジェクション・ミキサ（Image  Rejection  Mixer）と呼ばれ、フィルタを使わない信号分離に使われている。イメージ・リジェクション・ミキサの代わりに、ミキサとフィルタを用いる単純な構成も可能である。
【００３０】
合成器１７の出力信号は、増幅器１８と送信ゲート１９と送受切替回路２１を経由し、プローブ６から試料管９に照射される。試料管９中の核が高周波照射により示す挙動は、プローブ６で中心周波数frfを持つ電気信号の形で検出され送受切替回路２１を経由し低ノイズ増幅器２２で増幅される。増幅された信号はミキサ３１で再び増幅される。ミキサ３１の前か後にはフィルタを入れてもよい。受信信号は、受信部ミキサ３２で局部信号LOと混合され、フィルタ３３で濾過される。フィルタ３３から出力される信号IFは中心周波数fifを持ち、ディジタル検波器３に入力される。
【００３１】
ディジタル検波器３は、入力されるアナログ信号IFをディジタルデータに変換するADC（Analog  Ditial  Converter）４１と、ADC４１からディジタルデータを受けディジタル直交検波を行うFPGA４２と、ディジタル直交検波で得られる分散信号uと吸収信号vを受けて補正信号sを出力するDSP（Digital  Signal  Processor）４３を含む。FPGA４２の代わりに、ディジタル検波機能を実装できる他の電子素子、例えばASIC（Application  Specific  Integrated  Circuit）などを用いてもよい。
【００３２】
ディジタル検波器３は、この他に電源回路、増幅回路、フィルタ回路、データ・バスとの接続を担当するブリッジ回路、メモリなどを含むが、本発明の内容と直接関連がないため、これらの回路は説明しない。
【００３３】
直交検波を行うディジタル検波器は非特許文献３に開示されている。非特許文献３によれば、１９７８年Hewlett  Packard社の技術者により提案されたディジタル直交検波は、１９８４年Dieter  ZiessowによりNMR装置に最初に適用された。
【００３４】
ディジタル直交検波では、分岐されてない単一の中間周波数（IF：Intermediate  Frequency）のアナログ入力信号が単一のADCによりディジタル信号に変換され、uとvの信号への分離はディジタルデータの状態でディジタル信号処理器により行われる。一方、アナログ直交検波では、例えば従来技術２に開示されたように、アナログ入力信号を位相の９０度異なるuとvの２つに直交検波された後にADCに入れる。ディジタル直交検波は、アナログ直交検波に比べて、uとv間の相対位相精度と振幅バランスが良好である。
【００３５】
ディジタルデータの分離は数学的な２進演算によりコンピュータと同様に安定的に実施できる反面、アナログ信号の分岐と位相遅延は素子のばらつきや環境変化による変動に弱い。従来技術では、素子の性質がばらつくため装置毎に異なる調整が必要である。これは、装置温度が変わると分岐したuとv信号の相対位相や振幅が変動するためである。u信号とv信号の相対位相および振幅の精度と安定性は、後述する本発明による補正信号sの演算効果を高めるために重要である。
【００３６】
以下、本発明におけるディジタル検波器３の構成と動作を説明する。アナログ入力信号は分岐されてない単一の入力部からADC４１に入る。アナログ入力信号は、中間周波数IFを中心周波数として持ち、典型的には１kHz以下の帯域幅を持つ。ADC４１はアナログ入力信号の変換で得たディジタルデータをFPGA４２へ出力する。
【００３７】
ADC４１の動作周波数、即ち１秒の間にADC４１がアナログ信号からサンプルを採取し変換するディジタルデータの数は、入力信号の帯域幅の８倍若しくは８倍以上である。ADCの動作周波数が入力信号の帯域幅の2倍を超える点は、非特許文献３記載のオーバーサンプリングに類似しているが、オーバーサンプリングとは目的及び後処理の面で異なる。
【００３８】
オーバーサンプリングの目的はディジタル化に伴う変換誤差を低減することであり、後処理ではナイキスト定理により規定された数のデータだけが演算に用いられる。ナイキスト定理によれば、ADCの動作周波数が入力信号の帯域幅より２倍以上であれば、ディジタルデータから本来のアナログ信号を再現することができる。オーバーサンプリングにおいて、ナイキスト定理で規定される入力信号の帯域幅以上に変換したディジタルデータは、間引き(decimation)により捨てられる。
【００３９】
しかし、本発明によるADC４１の動作周波数は、後述する補正信号sの演算方法の効果を高めるために、ナイキスト定理の最低値より４倍以上と設定されている。また、変換されたデータに対する処理も以下で述べるように、非特許文献３記載のオーバーサンプリングと異なる。
【００４０】
図３は本発明におけるディジタル検波器のデータ処理を示す機能構成図である。この機能構成は、公知であるADC４１からディジタルフィルタ（DF：Digital  Filter）105までのディジタル直交検波の基本操作と、線形フィルタ（LF：Linear  Filter）１０６から出力OUT１１７までの本発明に固有の誤差磁場の制御操作からなる。
【００４１】
ADC４１で変換されたディジタルデータsadは、以下の操作からなるディジタル直交検波を受ける。ディジタルデータsadは２つに複製されディジタル乗算器（Mult）１０２に入力される。ディジタル乗算器１０２は、複製されたディジタルデータsadに、NCO（Numerically  Controlled  Oscillator）１０３が入力クロック信号CLKに同期して出力する正弦波データsncoを乗算する。正弦波データsncoの一方は位相をπ／２遅延させる遅延回路１０４を経由する。乗算されたデータumとvmはディジタルフィルタDF（Digital  Filter）105で濾過され、必要な周波数帯域の情報のみを持つがudfとvdfとなる。以上がディジタル直交検波の基本操作であり、非特許文献３の記載と類似する。
【００４２】
udfとvdfに対する処理は、従来のディジタル直交検波と本発明とでは異なる。従来のディジタル直交検波では、間引き処理により、ユーザーが予め決めた周波数帯域幅を再現するために必要なデータだけを残す。例えば、１００Hzの帯域幅を再現するために必要なデータは、ナイキスト定理によれば、２００点／秒である。従来のディジタル直交検波において、間引き処理で残されるデータは多くてもナイキスト定理で定められた量の８倍未満である。
【００４３】
一方、本発明では線形フィルタ（ＬＦ）１０６によりudfとvdfを重みつき平均する。平均は移動平均でもよい。LF１０６から出力されるulfとvlfは、ナイキスト定理で定められる量の８倍以上のデータ量を有する。例えば、１００Hzの帯域幅に対し、ulfとvlfは８００点／秒以上のデータ量を持つ。ulfとvlfは位相補正器（PhCorr：Phase  Corrector）１０７により位相補正されたuphとvphになる。θラジアンを補正するため、PhCorr１０７では式（１）の処理を行う。
【００４４】
【数１】
0004811995.tif 000002
【００４５】
位相補正を、アナログ回路に比べ安定度の高いディジタル信号処理で行うことで、位相の安定性が改善できる。また、アナログ回路の複雑度を低減できるメリットもある。
【００４６】
正規化器（Nrm：Normalizer）１０８は、uphとvphに式（２）の正規化処理を加え、UとVを出力する。
【００４７】
【数２】
0004811995.tif 000003
【００４８】
vph0は、当該設定における吸収信号の最大値であり、ロックシステム２０４を用いた実測により容易に決定できる。式（２）の如く、分散信号と吸収信号を吸収信号の最大値で正規化することにより、ロック用溶媒と図３の１０９、１１０、１１４、１１５の演算に含まれるフィード・バック制御パラメータとの対応を保つことができる。
【００４９】
前記フィード・バック制御パラメータに対応するフィード・バック制御系は、NMR信号から誤差磁場を検出し誤差磁場量に対応する電流をロック電源から試料管周囲のロックコイルへ流すことにより磁場安定化を行うロックシステム全体である。
【００５０】
正規化しなければ、ロックシステム２０４の設定により分散信号および吸収信号の強度が変わった場合に、フィード・バック制御パラメータを設定し直さなければならない。例えば、図２の送信増幅器１８のゲインを上げ送信信号パワーを大きくし、取得される分散信号および吸収信号の強度が大きくなった場合に、正規化なしではフィード・バックのゲインを下げなければ制御が不安定になる。
【００５１】
正規化は、図３の定常関数演算器１０９以下の演算に入る前に、その入力の大きさを式（２）により整えることにより、不安定の問題を低減する。なお、正規化により雑音成分がより大きくなる場合に備え、最終ゲインを調整するGAIN１１６を設ける。GAIN１１６は、Nrm１０８から受けるnorm＝ vph0を基に、補正信号の大きさを調整する。調整は、例えば単純な比例調整、閾値カットオフ調整から、より高次の関数を用いた調整まで可能である。また、正規化の代わりに、フィード・バック制御パラメータを変えることも可能であり、その場合、図３に示す１１０、１１５を調整すればよい。
【００５２】
正規化されたUとVは定常関数演算器１０９により、分散信号と吸収信号から誤差磁場値を求める、誤差磁場マッピング処理に入る。以下では、表記を簡単にするために、UとVをuとvと書く。
【００５３】
誤差磁場マッピングのために、本発明ではuとvとを用いる定常関数処理１０９と、uとvと共にdu／dt、dv／dtを用いる過渡関数処理１１４を実施する。du／dtおよびdv／dtは、時間変化率計算処理１１２と１１３で求められ、過渡関数処理１１４に入力される。時間変化率計算処理１１２と１１３は、線型フィルタ１０６でナイキスト定理による必要量の８倍以上の量で残したデータを対象に統計的なデータ処理を行うことで、du／dtおよびdv／dtの誤差を著しく軽減できる。
【００５４】
前記データ処理は、例えば、線型回帰がある。ある時点に時間変化率計算処理１１２に保持されるu値がN個であるとし、{u1， u2，…，uN}と表す。{u1，u2，…，uN}が1つの直線上にあると仮定し、線型回帰によりその直線の傾斜率を求めれば、その値はdu／dtと等しい。線型回帰の手法としては、最小2乗法や中間値法など種々の手法がある。du／dtを計算し出力する周期は、図２の送受切替スイッチ回路２１の切替周期と等しくすればよいが、その他にも種々の変形が可能である。
【００５５】
定常関数処理１０９および過渡関数処理１１４の出力は各々スカラー量になり、定常関数重みWs（Weight−Static）１１０と過渡関数重み（Wt：Weight-Transient）１１５を乗算し、ディジタル合成器１１１で足し合わされる。ディジタル合成器１１１の出力はGAIN１１６でフィード・バック・ゲインの調整を受け、出力回路１１７で図２に示したLock電源４に適合する形式で出力される。
【００５６】
次に、誤差磁場マッピングにおける定常関数１０９と過渡関数１１４の具体例を示す。そのために、以下では分散信号および吸収信号と誤差磁場との関係をモデル化したBloch方程式を用いる。
【００５７】
Bloch方程式は、磁化ベクトルの成分{u0、v0、Mz}と共鳴周波数オフセットδとの方程式であり、式（３）〜（５）の連立微分方程式で表す。
【００５８】
【数３】
0004811995.tif 000004
【００５９】
但し、α＝１／T1、β＝１／T2、δ＝γＢ０−ω、ω１＝γＢ１であり、T1とT2は各々縦と横緩和時間、γは核の磁気回転率（magnetogyric ratio）、Ｂ０はマグネットの磁場強度、Ｂ１はプローブから照射される高周波磁場の振幅である。またＭ０は磁化ベクトルの大きさである。分散信号uと吸収信号vは、各々、上記u0とv0に比例する。また、誤差磁場はδと等しい。u0、v0と分散信号u、吸収信号v間の比例係数は、以下の議論で不要なため、今後１とし、u0とv0の代わりにuとvと書く。実際においても、ディジタル検波器３の内部で正規化し、比例係数の影響をなくす。
【００６０】
誤差磁場マッピング関数は、uおよびvからδを求める関数である。特に、du／dt、dv／dt、dMz／dtが小さい場合のマッピング関数は、式（３）〜（５）の時間微分項を０とした定常解から得ることができる。その定常解は、微分項を含まない簡単な連立方程式を解くことで、以下の式（６）〜（８）の如く得られる。
【００６１】
【数４】
0004811995.tif 000005
【００６２】
従来技術１は式（６）のみを用いる。以下の式（９）を満たす場合、
【００６３】
【数５】
0004811995.tif 000006
【００６４】
式（６）は式（１０）に書き換えることができる。
【００６５】
【数６】
0004811995.tif 000007
【００６６】
誤差磁場予測量Ｓ０は、単純に式（１１）でよい。
【００６７】
【数７】
0004811995.tif 000008
【００６８】
適切なゲインc＝β２／Ｍoω１を選択すれば、S０＝δとなる。一方、従来技術２の誤差磁場予測量Ｓ０は、式（１２）である。
【００６９】
【数８】
0004811995.tif 000009
【００７０】
適切なゲインcを選択することで、S０＝δとなる。
【００７１】
式(１２)は、式（６）を式（７）で除算することで直ちに得られる。式（９）のような近似条件がないことは、式（１２）の長所である。一方、式（１２）は、式（１１）に比べノイズと位相誤差に対するロバスト性が低い短所を持つ。また、係数が物理量β＝１／T2に依存する点も短所である。βは、原理的には横緩和時間を測定し、その逆数を用いればよいが、横緩和時間を簡単に測定することは難しい。横緩和時間は、溶媒の物性だけでなく、磁場均一度や温度などの測定環境にも依存するためである。測定環境が変わる度に横緩和時間を測定し、係数を補正することは実用的ではない。
【００７２】
上記の式（１１）と（１２）は、式（３）〜（５）の時間微分項を全て０にして得られたため、以下では定常型関数と呼ぶことにする。定常型関数は、その基本条件のため、定常解周辺の領域を離れると誤差磁場の予測量S０と真の誤差磁場量δとの差が大きくなる。この制約は、しかし、ロックシステムで補正しようとする誤差磁場が超電導磁石内部の永久電流の減衰による磁場ドリフトなら問題にならない。磁場ドリフトによる磁場変動率は１ｐｐｂ／時間に過ぎないため、式（３）〜（５）の時間微分項を０にする近似は適切である。
【００７３】
磁場変動率が大きい場合、定常型関数の精度は低下する。即ち、急激な磁場変動が加わり定常解から遠く離れた状態になると、上記定常型関数ではS０と誤差磁場量との差が増加する。この点を改善するため、時間微分項を０にしない新たな誤差磁場マッピング関数を開発した。新たなマッピング関数の特徴を、従来技術２と比較して説明する。従来技術２は、δが大きい場合に式（３）の変形である。
【００７４】
【数９】
0004811995.tif 000010
【００７５】
式（１３）による誤差磁場マッピング関数は式（１４）となる。
【００７６】
【数１０】
0004811995.tif 000011
【００７７】
式（１４）の右辺第１項は式（１２）と同じ定常状態を表す項である。第２項はdu／dtを含み過度状態を表す。定常状態に到達すればdu／dtが小さくなり、式（１４）は式（１２）に近くなる。係数はc2とcω／c2で定義する。c2はフィード・バックの安定度などにより選択し、cω／c2はβと等しく選択すれば、S０＝δとなる。
【００７８】
式（１４）は磁場変動率の大きい場合でも適用できる長所を持つ。しかし、吸収信号vによる除算を含むため、ノイズと位相誤差に対する低いロバスト性と、困難な係数調整の面では式（１２）と同じ短所を持つ。また、du／dtの計算方法によっては、ノイズが増大する新たな問題もある。
【００７９】
過度状態にも対応する新しい誤差磁場マッピング関数は、式（３）と（４）を連立して得た式（１５）の関係を用いる。
【００８０】
【数１１】
0004811995.tif 000012
【００８１】
式（１５）は近似を含まないδとu、v間の厳密な関係である。しかし、式（１５）は観測できない磁化ベクトルのz成分Mｚを含む問題がある。この問題は式（５）を用いることで解決する。式（５）を変形すると式（１６）になる。
【００８２】
【数１２】
0004811995.tif 000013
【００８３】
右辺第３項はdMz／dtを含む。Z方向に大きな静磁場B0が存在する条件において、Mzの変化はエネルギーの逸散を必要とする。そのため、dMz／dtはdu／dt、dv／dtに比べ著しく小さいと考えられる。そして、ロックオン状態の時分割送受信は、Mzを常に一定な値に戻す。
【００８４】
従って、本発明のマッピング関数では式（１６）の第３項は無視する。式（１６）の第1項と第２項は何れも観測可能な量である。本発明のマッピング関数では第1項を用いた誤差磁場マッピング関数を実装しテストした。ロックオン送信では核の飽和を避けることが重要であり、ω１＝γＢ１を小さくすることが一般的になっているためである。
【００８５】
式（１５）に式（１６）の第1項を代入すると、式（１７）の誤差磁場マッピング関数を得る。
【００８６】
【数１３】
0004811995.tif 000014
【００８７】
定常型関数の場合と同じく、式（１７）は式（１４）に比べ、ノイズと位相誤差に対するロバスト性の面で有利である。また、振幅により係数を決めるため、横緩和時間から係数を決める式（１４）に比べ、係数調整の自動化の面でも有利である。
【００８８】
図３に示したディジタル検波器の機能はソフトウェアでも実現できる。ソフトウェアは、例えば、コントローラ１に実装すればよい。
【００８９】
図４はディジタル検波器の機能をソフトウェアで実現した例を示すフローチャートである。検波処理の最初にデータ読込み処理ｓ１０１がA／D変換されたデータｄ１０２を読み込む。ｓ１０１に読み込まれたデータは、IとQの2つに分岐され乗算処理ｓ１０３とｓ１０３'に入る。
【００９０】
乗算処理ｓ１０３とｓ１０３'は、図３に示した１０２ブロックに相当する処理であり、角度とその角度での正弦波値を収納したテーブルｄ１０４および余弦波値を収納したテーブルｄ１０４'から正弦波値および余弦波値を読み出してｓ１０１から入力されたデータに乗算する。ｄ１０４およびｄ１０４'は、図３に示した１０３、１０４ブロックに相当する。
【００９１】
乗算されたIとQのデータは、ディジタル・フィルタリングｓ１０５、線形フィルタリングｓ１０６、位相補正ｓ１０７、正規化ｓ１０８を受ける。ｓ１０５からｓ１０８までの処理は、各々、図３の１０５、１０６、１０７、１０８ブロックに対応する。
【００９２】
続いて、ｓ１０９、ｓ１１０、ｓ１１１、ｓ１１２、ｓ１１３、ｓ１１４の処理は、ｓ１０８処理で正規化された分散信号および吸収信号を用い、式（１７）の計算を行う。処理ｓ１０９とｓ１１０は式（１７）の第１項を求める処理であり、処理ｓ１１２とｓ１１３は式（１７）の第２項を求める処理である。処理ｓ１１１は、du／dtとdv／dtを求める処理である。処理ｓ１１４は式（１７）の第1項と第2項とを足し合せる処理である。
【００９３】
ｓ１１５処理は、正規化処理ｓ１０８の結果に応じてフィード・バック・ゲインをｓ１１４処理の出力値に乗算するもので、図３に示した１１６ブロックに対応する。ｓ１１６処理は、磁場補正値データｄ１１７をロック電源４へ出力する処理である。図３の１１７ブロックに相当する。
【００９４】
図２から明らかなように、本発明のロックシステム構成は、中間周波数fifに1つの局部信号LOを混ぜ最終的な高周波周波数frfを得る、１段周波数変換方式である。送信部においてはイメージ・リジェクション・ミキサ１５〜１７が、受信部においてはミキサ３２とフィルタ３３が、各々周波数変換部である。
【００９５】
１段周波数変換方式は、温度などの環境変動に弱くまた素子毎のばらつきの大きいアナログ回路を減らすことができるため、多段周波数変換方式よりシステムの安定性と保守性、コストの面でメリットを持つ。
【００９６】
一方、１段周波数変換方式では、周波数変換における周波数選択性に限界があるため、fifはfrfに比べ一定程度以上の高い周波数でなければならない制約がある。例えば、frfが92MHzの時、fifが1MHz、LOの周波数fLOが93MHzなら、周波数変換部では92MHz、93MHz、94MHzが混ざった信号の中から92MHz信号だけを選び出さなければならない。
【００９７】
しかし、中心周波数に対する選択帯域の比が１%になるような周波数変換は、不可能ではないにしても実現が非常に困難である。回路として容易に実現するためには、前記比が10%以上であることが好ましい。回路の実現が容易であることは、安定性や保守性、コストの面で有利である。回路を容易に実現でき、かつ１段周波数変換方式を用いるためには、中間周波数fifを高くしなければいけない。
【００９８】
高い周波数の信号を高速でアナログ／ディジタル変換を行う本発明の構成において、DDS１４を用いて中間周波数信号を発生させることは位相の長時間安定性に重要な利点を持つ。
【００９９】
ところで、中間周波数fifがFPGA４２中のNCO１０３が作る信号周波数と１mHzのオフセットを持ったら、10秒に3．6度の位相差が生じることになる。位相差は分散信号uと吸収信号vを式（１）により混ぜることになり、ロックシステムの動作に致命的な障害となる。fifが高い周波数、例えば10MHzなら、１mHzのオフセットは10の−10乗の小さい量であり制御が極めて困難である。
【０１００】
本発明では、NCO１０３と同じ周波数表示ビット数とクロック周波数を持つDDS１４を採用し中間周波数信号を出力することで、この問題を解決する。また、DDS１４とADC４１およびNCO１０３のクロックを同じ周波数源FRQSYN１３から供給することで、位相誤差は如何なる場合でも45度／スキャン以下に抑えることができる。
【０１０１】
中間周波数fifの出力にDDS１４を用いる本発明の構成は、更に、ロックシステムの掃引動作において大きな利点を持つ。掃引とは、ロックオンに入る前に、ロックオンする周波数を決めるための操作である。
【０１０２】
掃引は、磁場もしくは送受信周波数を次々と変えながら、ユーザーが「掃引幅」として指定した周波数帯域全体に渡ってNMR信号を取得する操作である。従来、掃引は磁場を変えて行う磁場掃引が一般的だった。磁場掃引はロック電源の指令値を一定に変えるだけで簡単に実現できるためである。一方、周波数を変える周波数掃引は、送信と受信の周波数を短い時間で変える必要があるため、周波数と位相の精度を実現することは複雑である。しかし、DDSは周波数を１μ秒未満の時間で変えることができるため、DDS１４を用いる本発明の構成では周波数掃引をシステムの構成を複雑にすることなく実現できる。
【０１０３】
周波数掃引はロックシステム２０４の送受信周波数を変えるだけであり、その影響はロックシステム内に限定される。一方、磁場掃引はマグネットの磁場を変えるため、NMRシステムの全体に影響を及ぼす。磁場を変えると、ロック信号だけでなく、観測核からのNMR信号も周波数が変わるためである。
【０１０４】
本発明の構成によれば、NMR信号の観測に影響を与えることなく、掃引で共鳴周波数を探すことができる。本発明のシステムは、ロック外れの発生時に、観測を続行しながら周波数掃引でロックオンに復帰する動作を行うことができる。
【０１０５】
本発明の磁場安定化装置は、NMRの他にも、MRIやEPR（Electron  Paramagnetic  Resonance）など、磁場の安定性を必要とする分野で利用できる。
【図面の簡単な説明】
【０１０６】
【図１】NMR装置の概念図。
【図２】本発明の一実施例によるロックシステムを示す構成図。
【図３】一実施例によるロックシステムのディジタルデータ処理を示す機能構成図。
【図４】一実施例によるロックシステムのディジタルデータ処理を示すフローチャート。
【符号の説明】
【０１０７】
１…コントローラ、２…ロック送受信器、３…ディジタル検波器、４…ロック電源、５…プリアンプボックス、６…プローブ、７…磁石、８…ロックコイル、９…試料管、１１…CPU、１２…      周波数源、１３…周波数合成器、１４…DDS、１５…90度位相変換器、１６，１６’…    送信部ミキサ、１７…合成器、１８…増幅器、１９…送信ゲート、２１…送受切替回路、２２…低ノイズ増幅器、３１…増幅器、３２…受信部ミキサ、３３…フィルタ、３４…増幅器、４１…ADC、４２…FPGA、４３…DSP、１０２…ディジタル乗算器、１０３…NCO、１０４…π／２遅延回路、１０５…ディジタルフィルタ、１０６…線形フィルタ、１０７…位相補正器、１０８…正規化器、１０９…定常関数演算器、１１０…重み付き乗算器、１１１…ディジタル合成器、１１２，１１３…時間変化率演算器、１１４…過渡関数演算器、１１５…重み付き乗算器、１１６…フィード・バック・ゲイン、１１７…出力回路、２０１…ユーザーコンピュータ、２０２…コントローラ、２０３…観測システム、２０４…ロックシステム。
【図１】
0004811995.tif 000015
【図２】
0004811995.tif 000016
【図３】
0004811995.tif 000017
【図４】
0004811995.tif 000018