JPS60105979A - 信号処理装置 - Google Patents
信号処理装置Info
- Publication number
- JPS60105979A JPS60105979A JP58213102A JP21310283A JPS60105979A JP S60105979 A JPS60105979 A JP S60105979A JP 58213102 A JP58213102 A JP 58213102A JP 21310283 A JP21310283 A JP 21310283A JP S60105979 A JPS60105979 A JP S60105979A
- Authority
- JP
- Japan
- Prior art keywords
- data
- signal
- spectrum
- processing device
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の技術分野〕
この発明は、たとえば核磁気共鳴現象などを利用して周
波数スペクトルを作るための信号処理装置、特に自動的
に誤差補正を行なって周波数スペクトルを作ることがで
きる信号処理装置に関するものである。
波数スペクトルを作るための信号処理装置、特に自動的
に誤差補正を行なって周波数スペクトルを作ることがで
きる信号処理装置に関するものである。
核磁気共鳴現象などを利用して誤差補正関数を用いずに
周波数スペクトルを作るためには次の様に装置を調整し
なければならない。すなわち、９０゜の位相差を保った
ままで位相調整された２つの参照信号をそれぞれ出力信
号の大きさを揃えられｔコ２つの位相敏感検波器に供給
し同期検波を行なし）吸収成分と分散成分とを出力する
ようにする。この時、増幅器に固有の不感時間があるた
めに信号の欠落を生じる可能性があるが、この信号の欠
落を防止するためにエコー信号を発生させて同期検波せ
ねばならない。この検波出力から周波数スペクトルを作
るためにコンピュータで信号処理を行なうのであるが、
コンピュータで信号処理を行なうために検波出力をディ
ジタル化するが、エコー信号の吸収成分と分散成分との
中心がそれぞれディジタル化されていなければならない
。この時に、増幅器や位相敏感検波器の経時変化による
エコー信号の変化が生じないように計器の制御を行ない
、エコー信号の中心以前の前半分の信号又は中心以後の
後半分の信号をスペクトル解析に用いる。
周波数スペクトルを作るためには次の様に装置を調整し
なければならない。すなわち、９０゜の位相差を保った
ままで位相調整された２つの参照信号をそれぞれ出力信
号の大きさを揃えられｔコ２つの位相敏感検波器に供給
し同期検波を行なし）吸収成分と分散成分とを出力する
ようにする。この時、増幅器に固有の不感時間があるた
めに信号の欠落を生じる可能性があるが、この信号の欠
落を防止するためにエコー信号を発生させて同期検波せ
ねばならない。この検波出力から周波数スペクトルを作
るためにコンピュータで信号処理を行なうのであるが、
コンピュータで信号処理を行なうために検波出力をディ
ジタル化するが、エコー信号の吸収成分と分散成分との
中心がそれぞれディジタル化されていなければならない
。この時に、増幅器や位相敏感検波器の経時変化による
エコー信号の変化が生じないように計器の制御を行ない
、エコー信号の中心以前の前半分の信号又は中心以後の
後半分の信号をスペクトル解析に用いる。
しかしながら、位相敏感検波器の出力が吸収成分と分散
成分とになるように参照信号の位相調整を行なうことは
困難であり、計器の特性の経時変化ヲ無くすことも困難
である。更に、ディジタル化する時にエコー信号の中心
をディジタル化することが困難であることから、スペク
トル解析にエコー信号の中心から半分の信号を使うこと
が困難となり、作られた周波数スペクトルに誤差が生じ
ることになる。このために誤差補正が必要となる。
成分とになるように参照信号の位相調整を行なうことは
困難であり、計器の特性の経時変化ヲ無くすことも困難
である。更に、ディジタル化する時にエコー信号の中心
をディジタル化することが困難であることから、スペク
トル解析にエコー信号の中心から半分の信号を使うこと
が困難となり、作られた周波数スペクトルに誤差が生じ
ることになる。このために誤差補正が必要となる。
次に従来法による誤差補正の実施例の動作を第１図、第
２図を参照しながら説明する。第１図は実施例のシステ
ムの構成を示すブロック包である。
２図を参照しながら説明する。第１図は実施例のシステ
ムの構成を示すブロック包である。
第１図から明らかなように、受信器（１）が受信した信
号を増幅する増幅器（２）と、増幅器（２）の出力を２
つに分配する第１の分配器（３）と、発振器（４）から
出力される参照信号の位相を変える移相器（５）の出力
を２つに分配する第２の分配器（６）と、第２の分配器
（６）の出力で第１の分配器（３）の出力を同期検波す
る第１の位相敏感検波器（８）と、第２の分配器（６）
の出力を第１の位相敏感検波器（８）に出力された参照
信号の位相と９０°位相を変化させる９０°移相器（７
）の出力で第１の分配器（３）の出力を同期検波する第
２の位相敏感検波器（９）とを設け、ｍ１１の位相敏感
検波器（８）と第２の位相敏感検波器（９）との出力を
それぞれディジタル化する第１のＡ／Ｄ変換器ＯＱと第
２のＡ／Ｄ変換器αｐと、第１のＡ／Ｄ変換器００と第
２のＡ／Ｄ変換器の出力をコンピュータ（６）の内部に
ある入力回路α谷からｃ ｐ Ｕ （１３を経由してメ
モリａ→にストアし、ＣＰＵＱ■でメモリαゆの内容に
対して演算を行なって出力回路ＯＱから周波数スペクト
ルとしてＣＲＴ （＋７１に出力し、ＣＲＴ（ｌηが周
波数スペクトルを表示する。
号を増幅する増幅器（２）と、増幅器（２）の出力を２
つに分配する第１の分配器（３）と、発振器（４）から
出力される参照信号の位相を変える移相器（５）の出力
を２つに分配する第２の分配器（６）と、第２の分配器
（６）の出力で第１の分配器（３）の出力を同期検波す
る第１の位相敏感検波器（８）と、第２の分配器（６）
の出力を第１の位相敏感検波器（８）に出力された参照
信号の位相と９０°位相を変化させる９０°移相器（７
）の出力で第１の分配器（３）の出力を同期検波する第
２の位相敏感検波器（９）とを設け、ｍ１１の位相敏感
検波器（８）と第２の位相敏感検波器（９）との出力を
それぞれディジタル化する第１のＡ／Ｄ変換器ＯＱと第
２のＡ／Ｄ変換器αｐと、第１のＡ／Ｄ変換器００と第
２のＡ／Ｄ変換器の出力をコンピュータ（６）の内部に
ある入力回路α谷からｃ ｐ Ｕ （１３を経由してメ
モリａ→にストアし、ＣＰＵＱ■でメモリαゆの内容に
対して演算を行なって出力回路ＯＱから周波数スペクト
ルとしてＣＲＴ （＋７１に出力し、ＣＲＴ（ｌηが周
波数スペクトルを表示する。
第２図は第１図のコンピュータ（６）のメモリａ→に記
憶された周波数スペクトルを作るプログラムを示すフロ
ーチャートである。
憶された周波数スペクトルを作るプログラムを示すフロ
ーチャートである。
まず入力回路ａＧにエコー信号が第１のＡ／Ｄ変換器０
０と第２のＡ／Ｄ変換器αυとでディジタル化されて入
力される入力過程Ｑメがスタートする。入力回路ＯＱか
ら信号が入力されるとデータ作成過程（イ）によって被
処理データ５ｄｔ）と５ｓ（ｔ）とが作られる。
０と第２のＡ／Ｄ変換器αυとでディジタル化されて入
力される入力過程Ｑメがスタートする。入力回路ＯＱか
ら信号が入力されるとデータ作成過程（イ）によって被
処理データ５ｄｔ）と５ｓ（ｔ）とが作られる。
５ｃ（ｔ） ””（第１のＡ／Ｄ変換器α０の出カフ０
≦（≦ｔｆＳｓ （ｔ） ＝ ’（第２のＡ／Ｄ変換ｈ
’ｓ （１１） ノ出力Ｊ Ｏ≦ｔ≦ｔｆここで１＝０
はエコー信号をＡ／Ｄ変換し始めた時間、ｔ ＝ ｔｆ
はエコー信号をＡ／Ｄ変換し終えた時間であり、１＝０
ではエコー信号の中心付近がＡ／Ｄ変換され、ｔ ＝
ｔｆではエコー信号はＡ／Ｄできない程十分小さくなっ
ているとする。
≦（≦ｔｆＳｓ （ｔ） ＝ ’（第２のＡ／Ｄ変換ｈ
’ｓ （１１） ノ出力Ｊ Ｏ≦ｔ≦ｔｆここで１＝０
はエコー信号をＡ／Ｄ変換し始めた時間、ｔ ＝ ｔｆ
はエコー信号をＡ／Ｄ変換し終えた時間であり、１＝０
ではエコー信号の中心付近がＡ／Ｄ変換され、ｔ ＝
ｔｆではエコー信号はＡ／Ｄできない程十分小さくなっ
ているとする。
第１および第２の位相敏感検波器（８） 、 （９）の
参照信号の位相が９０°ずれており、ゲインと特性とが
同じであれば５ｃ（ｔ）と５ｓ（ｔ）とは次の様に書け
る。
参照信号の位相が９０°ずれており、ゲインと特性とが
同じであれば５ｃ（ｔ）と５ｓ（ｔ）とは次の様に書け
る。
５ｃ（ｔ） ＝ ｆ：：Ａｆ（Ω）ｃｏｓ（Ω（ｔ−ｔ
ｏ ） −（ｔ’Ω十〇））ｄΩ５ｓ（ｔ） ＝ Ｊ：
：ＡｆｃΩ）ｓｉｎ （Ω（ｔ−ｔｏ） −’（ｔ’Ω
＋０））ｄΩここでＡは第１および第２の位相敏感検波
器（８）。
ｏ ） −（ｔ’Ω十〇））ｄΩ５ｓ（ｔ） ＝ Ｊ：
：ＡｆｃΩ）ｓｉｎ （Ω（ｔ−ｔｏ） −’（ｔ’Ω
＋０））ｄΩここでＡは第１および第２の位相敏感検波
器（８）。
（９）の出力の大きさ、ｆｔＡ＞はめたい角周波数Ωの
スペクトル、ｔｏはエコー信号の中心から実際にＡ／Ｄ
変換された最初の点までの時間のずれ、ｔ′Ωとして、
誤差の無い場合を考える。ｔｏ＝ ｏ、 ｔ’Ω＋０＝
０であるから、 ５ｃ（ｔ）十１Ｓｓ（ｔ） ＝ Ａ、／：ｌ： ｆ帆）
ｃｏｓｎｔｄΩ十ｌへｆ：ニｆ（Ω）ｓｉｎΩｔｄΩ＝
５（ｔ） と表わせることから、ｆ（Ｑ）はＳ （ｔ）を逆フーリ
エ変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）をそれ
ぞれ実部と虚部としてデータ変換過程翰によって逆フー
リエ変換を行なってｈ［有］）が作られる。
スペクトル、ｔｏはエコー信号の中心から実際にＡ／Ｄ
変換された最初の点までの時間のずれ、ｔ′Ωとして、
誤差の無い場合を考える。ｔｏ＝ ｏ、 ｔ’Ω＋０＝
０であるから、 ５ｃ（ｔ）十１Ｓｓ（ｔ） ＝ Ａ、／：ｌ： ｆ帆）
ｃｏｓｎｔｄΩ十ｌへｆ：ニｆ（Ω）ｓｉｎΩｔｄΩ＝
５（ｔ） と表わせることから、ｆ（Ｑ）はＳ （ｔ）を逆フーリ
エ変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）をそれ
ぞれ実部と虚部としてデータ変換過程翰によって逆フー
リエ変換を行なってｈ［有］）が作られる。
ｈ（Ω）＝（が１〜）ｅ−曹ｔ
２０
＝ 、、７”、”（Ｓｄｔ）＋１ｓｓ（ｔ） ） ｅ
ｄｔ＝ 古”、 ”（Ｓｃ（ｔ）ｃｏｓΩｔ＋５ｓ（ｔ
）ｓｉｎΩｔ＋ ｉ （Ｓｓ（を沖ＳΩｔ−５ｃ（ｔ）
ｓｉｎΩｔ））ｄｔｈ（Ω）の実部のみをとるとｆ（Ω
）になるので、ｆ（Ω）＝（ｈ（Ω））嘱部 従来の方法では、この関係を使って周波数スペクトルを
めている。
ｄｔ＝ 古”、 ”（Ｓｃ（ｔ）ｃｏｓΩｔ＋５ｓ（ｔ
）ｓｉｎΩｔ＋ ｉ （Ｓｓ（を沖ＳΩｔ−５ｃ（ｔ）
ｓｉｎΩｔ））ｄｔｈ（Ω）の実部のみをとるとｆ（Ω
）になるので、ｆ（Ω）＝（ｈ（Ω））嘱部 従来の方法では、この関係を使って周波数スペクトルを
めている。
誤差のある場合を考える。
５ｃ（ｔ）＋１Ｓｓ（ｔ）＝ ｋｆ：： ：Ｆ（ｉｌ）
ＣＯ３（Ω（ｔ ｔｏ） （ビΩ＋０））ｄΩ＋ｌＡｆ
：：ｆ（Ω）ｓｉｎ （Ω（ｔ ｔｏ） （ｔ’Ω十〇
））ｄｆｌ＝ ｅ−”５（ｔ−ｔｏ−ｔ／） と表わせることから、ｆ（Ω）はＳ （ｔ）を逆フーリ
エ変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）とをそ
れぞれ実部と虚部としてデータ変換過程翰によって逆フ
ーリエ変換を行なってｇｃＱ）が作られる。
ＣＯ３（Ω（ｔ ｔｏ） （ビΩ＋０））ｄΩ＋ｌＡｆ
：：ｆ（Ω）ｓｉｎ （Ω（ｔ ｔｏ） （ｔ’Ω十〇
））ｄｆｌ＝ ｅ−”５（ｔ−ｔｏ−ｔ／） と表わせることから、ｆ（Ω）はＳ （ｔ）を逆フーリ
エ変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）とをそ
れぞれ実部と虚部としてデータ変換過程翰によって逆フ
ーリエ変換を行なってｇｃＱ）が作られる。
ｈ（Ω） ＝ −’−５１（ｔ） ｅ −”’ｄｉ２π
Ａ。
Ａ。
２二Ｊ″〜
＝ ｆ（ｎ）＋ ｉ −（Ｓｍｔ）ωＳΩｔ ５ｃ（ｔ
）ｓｉｐΩｔ ）ｄｔとおくと、 一土猥−”ＳＣｔ ＋６ ｔ’）ｅ−”’ｄｔ、２πＡ ＝２ｙｒＪ：：−ｊ７ｅ−”Ｏ±Ω”ｏ＋ｔ’）＋５（
ｔ ｔ、ｌｒｔ’）ｅ−’Ω（ｔ−１４−ビ）ｄ（ｔ
ｔｃｒ−ｔ’）＝＝′−ｅ−１（０１ｆｌ（＋６＋ｊ’
））（Ｊ’；”＋、／ｆｆｉ、−，）ｓ（ｔｔｏ−１ｚ
）ｅ−＋４Ｘｔ−＋６−ｔ’ｋｔ−、ｔｚ）２πＡ 註泉−イ ＝６−”叶Ω（ｔｏ＋ビ）ｌ（ｈ（ｎ）＋ ５（ｔ−ｔ
ｏイ）８−１Ω（ｔ−ｊ、、−ビ）ｄ（ｔ−ｔｏ−ｔ’
））データ変換過程（ハ）で作られたｇ（Ω）はそのま
までは誤差を含んでいるので、データ処理過程（ハ）で
ｇ （Ｑ）に含まれている係数ｅ−４［ｏ＋Ｑ（ｔｏ＋
ｆｆ）ｌを消去する。データ処理過程（ハ）では０＋Ω
（ｔｏ＋ｔ’）の値を推定してｇ［有］）とｅｉ （０
＋ｎ（ｔｏ＋ｔ’）］との積を作るとｈ（ｎ） ＋ｆ、
−ｔＳ（ｔ ｔｏｔ’ ）ｅ”Ω３１鳴イ’ｄ（ｔ−ｔ
ｏ−ｔつ２πＡ が得られる。データ処理過程（財）の結果は周波数スペ
クトルになっていないので判定過程に）によって補正過
程（ハ）にデータ処理過程（財）の結果が送られ、仏几
−，ｓ（ｔ−１ｏ−ｔ’）ｅ−””−’−’ｄ（ｔ −
恥−ｔ’）の項をゼロにする補正を行なう。この補正過
程（イ）で行なわれる補正は解析的でなく経験的関数で
ある。補正過程（イ）によって周波数スペクトルが得ら
れ、判定過程（ハ）からスペクトル出力過程（イ）によ
っテｈ（Ω）の実部を取り出し、スペクトルｆ（Ω）が
作５られる。スペクトル出力過程ｅカで作られたスペク
トルは出力回路ＯＱからＣＲＴ （＋７１に出力されス
ペクトル表示過程（ハ）が行なわれる。
）ｓｉｐΩｔ ）ｄｔとおくと、 一土猥−”ＳＣｔ ＋６ ｔ’）ｅ−”’ｄｔ、２πＡ ＝２ｙｒＪ：：−ｊ７ｅ−”Ｏ±Ω”ｏ＋ｔ’）＋５（
ｔ ｔ、ｌｒｔ’）ｅ−’Ω（ｔ−１４−ビ）ｄ（ｔ
ｔｃｒ−ｔ’）＝＝′−ｅ−１（０１ｆｌ（＋６＋ｊ’
））（Ｊ’；”＋、／ｆｆｉ、−，）ｓ（ｔｔｏ−１ｚ
）ｅ−＋４Ｘｔ−＋６−ｔ’ｋｔ−、ｔｚ）２πＡ 註泉−イ ＝６−”叶Ω（ｔｏ＋ビ）ｌ（ｈ（ｎ）＋ ５（ｔ−ｔ
ｏイ）８−１Ω（ｔ−ｊ、、−ビ）ｄ（ｔ−ｔｏ−ｔ’
））データ変換過程（ハ）で作られたｇ（Ω）はそのま
までは誤差を含んでいるので、データ処理過程（ハ）で
ｇ （Ｑ）に含まれている係数ｅ−４［ｏ＋Ｑ（ｔｏ＋
ｆｆ）ｌを消去する。データ処理過程（ハ）では０＋Ω
（ｔｏ＋ｔ’）の値を推定してｇ［有］）とｅｉ （０
＋ｎ（ｔｏ＋ｔ’）］との積を作るとｈ（ｎ） ＋ｆ、
−ｔＳ（ｔ ｔｏｔ’ ）ｅ”Ω３１鳴イ’ｄ（ｔ−ｔ
ｏ−ｔつ２πＡ が得られる。データ処理過程（財）の結果は周波数スペ
クトルになっていないので判定過程に）によって補正過
程（ハ）にデータ処理過程（財）の結果が送られ、仏几
−，ｓ（ｔ−１ｏ−ｔ’）ｅ−””−’−’ｄ（ｔ −
恥−ｔ’）の項をゼロにする補正を行なう。この補正過
程（イ）で行なわれる補正は解析的でなく経験的関数で
ある。補正過程（イ）によって周波数スペクトルが得ら
れ、判定過程（ハ）からスペクトル出力過程（イ）によ
っテｈ（Ω）の実部を取り出し、スペクトルｆ（Ω）が
作５られる。スペクトル出力過程ｅカで作られたスペク
トルは出力回路ＯＱからＣＲＴ （＋７１に出力されス
ペクトル表示過程（ハ）が行なわれる。
従来の誤差補正方法は以上の様に構成されているので、
０とｔｏ＋ｔ’との２つの値をＪＩＵ定したうえに、最
後に残る誤差の補正が経験的に行なわれることから、信
号処理装置を設計する場合に誤差補正方法作りに任意性
があり、エコー信号が半分しか利用できないという欠点
があった。
０とｔｏ＋ｔ’との２つの値をＪＩＵ定したうえに、最
後に残る誤差の補正が経験的に行なわれることから、信
号処理装置を設計する場合に誤差補正方法作りに任意性
があり、エコー信号が半分しか利用できないという欠点
があった。
この発明は、かかる問題を改善するためになされたもの
であり、データ作成手段に基づいてエコー信号全体を被
処理データにすることによって、誤差補正関数を作るこ
となく自動的に誤差補正を行なうスペクトル作成手段を
備え、周波数スペクトルを作るようにした信号処理装置
を提案するものである。
であり、データ作成手段に基づいてエコー信号全体を被
処理データにすることによって、誤差補正関数を作るこ
となく自動的に誤差補正を行なうスペクトル作成手段を
備え、周波数スペクトルを作るようにした信号処理装置
を提案するものである。
第８図はこの発明による信号処理装置の一実施例の全体
構成図である。この実施例は第８図から明らかなように
、（１）から（＋１）は第１図に示したものと同じであ
り、第１および第２のＡ／Ｄ変換器ｌ１ｌＯ。
構成図である。この実施例は第８図から明らかなように
、（１）から（＋１）は第１図に示したものと同じであ
り、第１および第２のＡ／Ｄ変換器ｌ１ｌＯ。
（＋１）の出力から被処理データを作るデータ作成手段
（８１）と、このデータ作成手段（８υによって作られ
た被処理データに対してデータ変換を行なうデータ変換
手段（８２）と、このデータ変換手段（８のによって作
られたデータから周波数スペクトルを作るスペクトル作
成手段（財）と、このスペクトル作成手段（８３）によ
って作られた周波数スペクトルを表示するＣＲＴ（１７
）とから構成されている。第８図は第１および第２の位
相敏感検波器（８）、＋９）を用いた場合を示している
。
（８１）と、このデータ作成手段（８υによって作られ
た被処理データに対してデータ変換を行なうデータ変換
手段（８２）と、このデータ変換手段（８のによって作
られたデータから周波数スペクトルを作るスペクトル作
成手段（財）と、このスペクトル作成手段（８３）によ
って作られた周波数スペクトルを表示するＣＲＴ（１７
）とから構成されている。第８図は第１および第２の位
相敏感検波器（８）、＋９）を用いた場合を示している
。
次に上記実施例の動作を第４図を参照しながら説明する
。第４図はコンピュータ（２）のメモリ０４１に記憶さ
れた周波数スペクトルを作るプログラムを示すフローチ
ャートである。
。第４図はコンピュータ（２）のメモリ０４１に記憶さ
れた周波数スペクトルを作るプログラムを示すフローチ
ャートである。
まず入力回路ＯＱにエコー信号が第１および第２のＡ／
Ｄ変換器αす、αηでディジタル化されて入力される信
号入力過程（４１）がスタートする。信号が入力される
とデータ作成過程（４２）によって被処理データ５ｃ（
ｔ）と５ｓ（ｔ）とが作られる。この５ｃ（ｔ）と５ｓ
（ｔ）とは誤差を含んでおり次の様に鶴ける。
Ｄ変換器αす、αηでディジタル化されて入力される信
号入力過程（４１）がスタートする。信号が入力される
とデータ作成過程（４２）によって被処理データ５ｃ（
ｔ）と５ｓ（ｔ）とが作られる。この５ｃ（ｔ）と５ｓ
（ｔ）とは誤差を含んでおり次の様に鶴ける。
５ｃ（ｔ） ＝ （第１のＡ／Ｄ変換器（１０の出力）
０≦ｔ≦ｔｆＳｓ（ｔ） ＝ （第２のＡ／Ｄ変換器θ
υの出力）０≦ｔ≦１ｆここで１ ＝ ０はエコー信号
をＡ／Ｄ変換し始めた時間であり、ｔ ＝ ｔｆはエコ
ー信号をＡ／Ｄ変換し終えた時間であり、１＝０とｔ
＝ ｔｆとはエコー信号はＡ／Ｄ変換できない程十分小
さくなっており、ｔ−０とｔ ＝ ｔｆとの間にエコー
信号がある。
０≦ｔ≦ｔｆＳｓ（ｔ） ＝ （第２のＡ／Ｄ変換器θ
υの出力）０≦ｔ≦１ｆここで１ ＝ ０はエコー信号
をＡ／Ｄ変換し始めた時間であり、ｔ ＝ ｔｆはエコ
ー信号をＡ／Ｄ変換し終えた時間であり、１＝０とｔ
＝ ｔｆとはエコー信号はＡ／Ｄ変換できない程十分小
さくなっており、ｔ−０とｔ ＝ ｔｆとの間にエコー
信号がある。
第１および第２の位相敏感検波器（８）、（９）の参照
信号の位相が９０°ずれているため５ｃ（ｔ）と５ｓ（
ｔ）とは次の様に書ける。
信号の位相が９０°ずれているため５ｃ（ｔ）と５ｓ（
ｔ）とは次の様に書ける。
５ｃ（ｔ） ＝Ｊ：：＃（Ｑ）ＣＯＳ （Ω（ｔ ｔｏ
） （ビΩ十〇））ｄΩ５ｓ（ｔ）−Ｊ↑ニＭ＠ｓｉｎ
（Ω（ｔ−ｔｏ）−（ｔ’Ω十〇））ｄΩここでＡは
第１および第２の位相敏感検波器（８）。
） （ビΩ十〇））ｄΩ５ｓ（ｔ）−Ｊ↑ニＭ＠ｓｉｎ
（Ω（ｔ−ｔｏ）−（ｔ’Ω十〇））ｄΩここでＡは
第１および第２の位相敏感検波器（８）。
（９）の大きさであり、ｔｏはエコー信号の中心までの
時間であ°す、ｆ（Ω）、ｔ′Ω＋０は従来法の説明で
述べたものと同じであり、周波数スペクトルと誤差とを
示している。
時間であ°す、ｆ（Ω）、ｔ′Ω＋０は従来法の説明で
述べたものと同じであり、周波数スペクトルと誤差とを
示している。
５（ｔ） ＝ ＡＪ：：ｆ（ｎ）ｅ”’ｄｔ＝ ｅ”５
（ｔ−ｔ、ｔ’） と表わせることから、ｆ（Ｑ）は５（ｔ）を逆フーリエ
変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）とをそれ
ぞれ実部と虚部としてデータ変換過程（財）で逆フーリ
エ変換を行なって敢Ω）が作られる。
（ｔ−ｔ、ｔ’） と表わせることから、ｆ（Ｑ）は５（ｔ）を逆フーリエ
変換するとめられる。５ｃ（ｔ）と５ｓ（ｔ）とをそれ
ぞれ実部と虚部としてデータ変換過程（財）で逆フーリ
エ変換を行なって敢Ω）が作られる。
ｋ（Ω）＝−−ｆ’（５ｃ（ｔ）＋１Ｓｓ（ｔ））ｅ
”’ｄｔ２πＡ 盲シ”ｅ−”５（ｔ−ｔｏ−ｔ’）ｅ−”’ｄｔ： ｅ
−Ｎｏ＋Ω（ｔｏ＋ｔ’） ）ｆθ２゜データ変換過Ｎ
（４Ｂ）によって作られたｋ［有］）を使ってスペク
トル作成過程（ロ）によって周波数スペクトルｆ（Ω）
が作られる。
”’ｄｔ２πＡ 盲シ”ｅ−”５（ｔ−ｔｏ−ｔ’）ｅ−”’ｄｔ： ｅ
−Ｎｏ＋Ω（ｔｏ＋ｔ’） ）ｆθ２゜データ変換過Ｎ
（４Ｂ）によって作られたｋ［有］）を使ってスペク
トル作成過程（ロ）によって周波数スペクトルｆ（Ω）
が作られる。
：Ｆ（ｆｌ） ＝ ｌ ｋ（Ω）Ｉ
スペクトル作成過程■にょってデータ変換過程（財）の
出力の絶対値をとることでスペクトルｆ（Ω）が作られ
るのでスペクトル出力過程（４５）によって出力回路α
ＱからＣＲＴ（ｌηにスペクトルｆ（Ω）が出力され、
スペクトル表示過ｉ （４６）が行なわれ、周波数スペ
クトルがＣＲ′ｒ Ｑ７１上に表示される。
出力の絶対値をとることでスペクトルｆ（Ω）が作られ
るのでスペクトル出力過程（４５）によって出力回路α
ＱからＣＲＴ（ｌηにスペクトルｆ（Ω）が出力され、
スペクトル表示過ｉ （４６）が行なわれ、周波数スペ
クトルがＣＲ′ｒ Ｑ７１上に表示される。
なお、上記実施例は周波数スペクトルを作るのにコンピ
ュータを用いているが、他の計算装置で実行することも
可能である。
ュータを用いているが、他の計算装置で実行することも
可能である。
以上のようにこの発明によれば、誤差補正のパラメータ
を推定することなく、エコー信号を十分長い時間計測す
ることによって、データ作成手段。
を推定することなく、エコー信号を十分長い時間計測す
ることによって、データ作成手段。
データ変換手段、スペクトル作成手段に基づいて自動的
に誤差補正された周波数スペクトルが得られ、経験的な
誤差補正関数を作らなくてもよく、信号処理装置の設計
が容易になるという効果がある。
に誤差補正された周波数スペクトルが得られ、経験的な
誤差補正関数を作らなくてもよく、信号処理装置の設計
が容易になるという効果がある。
第１図は従来の信号処理装置の構成を示すブロック図、
第２図は周波数スペクトルを作成する過程を示すフロー
チャート、第８図はこの発明による信号処理装置の一実
施例の全体構成図、第４図は第８図における動作を示す
フローチャートである。 図中、（７）は９０’移相器、（３１）はデータ作成手
段、（８町よデータ変換手段、（８Ｂ）はスペクトル作
成手段、０″４はコンピュータ、■はスペクトル作成過
程である。なお図中同一符号は同一部分を示す。 代理人大岩増雄
第２図は周波数スペクトルを作成する過程を示すフロー
チャート、第８図はこの発明による信号処理装置の一実
施例の全体構成図、第４図は第８図における動作を示す
フローチャートである。 図中、（７）は９０’移相器、（３１）はデータ作成手
段、（８町よデータ変換手段、（８Ｂ）はスペクトル作
成手段、０″４はコンピュータ、■はスペクトル作成過
程である。なお図中同一符号は同一部分を示す。 代理人大岩増雄
Claims (3)
- （１）信号を受信する受信器、この受信器の受信した信
号を増幅する増幅器、この増幅器の増幅した信号を２つ
に分配する分配器、この分配器からの信号を入力として
発振器からの基準信号を参照して同期検波を行なう２つ
の位相敏感検波器、この２つの位相敏感検波器の出力を
ディジタル化する２つのＡ／Ｄ変換器、この２つのＡ／
Ｄ変換器の出力から被処理データを作るデータ作成手段
、このデータ作成手段によって作られた被処理データに
対してデータ変換を行なうデータ変換手段、このデータ
変換手段によって作られたデータから自動的に誤差補正
を行なって周波数スペクトルを作るスペクトル作成手段
を備えた信号処理装置。 - （２） データ変換手段は予め定められたデータ変換式
によって被処理データを変換するようになっている特許
請求の範囲第（１）項に記載の信号処理装置。 - （３） スペクトル作成手段は予め定められたスペクト
ル作成式によって周波数スペクトルを作成するようにな
っている特許請求の範囲第（１）項をこ記載の信号処理
装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58213102A JPS60105979A (ja) | 1983-11-11 | 1983-11-11 | 信号処理装置 |
| GB08428336A GB2149518B (en) | 1983-11-11 | 1984-11-09 | Signal processing unit |
| DE19843441290 DE3441290A1 (de) | 1983-11-11 | 1984-11-12 | Signalverarbeitungsanlage |
| US06/670,579 US4719581A (en) | 1983-11-11 | 1984-11-13 | Signal processing unit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58213102A JPS60105979A (ja) | 1983-11-11 | 1983-11-11 | 信号処理装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS60105979A true JPS60105979A (ja) | 1985-06-11 |
Family
ID=16633600
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58213102A Pending JPS60105979A (ja) | 1983-11-11 | 1983-11-11 | 信号処理装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4719581A (ja) |
| JP (1) | JPS60105979A (ja) |
| DE (1) | DE3441290A1 (ja) |
| GB (1) | GB2149518B (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4931977A (en) * | 1987-10-30 | 1990-06-05 | Canadian Marconi Company | Vectorial adaptive filtering apparatus with convergence rate independent of signal parameters |
| US4857846A (en) * | 1988-03-31 | 1989-08-15 | The Regents Of The University Of California | Rapid MRI using multiple receivers producing multiply phase-encoded data derived from a single NMR response |
| DE3825160A1 (de) * | 1988-07-23 | 1990-01-25 | Rohde & Schwarz | Messgeraet fuer hochfrequenzsignale |
| US5323332A (en) * | 1992-02-27 | 1994-06-21 | Infotec Development, Inc. | Instrument landing system calibrator |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5520495A (en) * | 1978-07-20 | 1980-02-13 | Univ California | Method and apparatus for mapping atomic nuceus density within object by using nucear magnetic resonance |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3609317A (en) * | 1970-01-29 | 1971-09-28 | Varian Associates | Variable frequency audio frequency modulator for r.f. spectrometer |
| US3824451A (en) * | 1972-11-15 | 1974-07-16 | Varian Associates | Pulse modulated fourier transform magnetic resonance spectrometer utilizing quadrature phase detection and high pass filter to eliminate undesired strong resonance lines |
| US3903401A (en) * | 1974-06-27 | 1975-09-02 | Bell Telephone Labor Inc | Spectrum analyzer using delta modulation encoding |
| US4070708A (en) * | 1976-09-01 | 1978-01-24 | Varian Associates, Inc. | Analytical instrument optimum data acquisition method and apparatus |
| JPS53129079A (en) * | 1977-04-18 | 1978-11-10 | Hitachi Ltd | Spectrum change monitoring system |
| JPS5516229A (en) * | 1978-07-21 | 1980-02-04 | Hitachi Ltd | Fourier transformation type nuclear magnetic resonance device |
| US4340862A (en) * | 1978-12-13 | 1982-07-20 | Picker International Limited | Imaging systems |
| GB2076541B (en) * | 1980-05-21 | 1984-05-31 | Emi Ltd | Nmr imaging apparatus and methods |
| GB2076540B (en) * | 1980-05-21 | 1984-02-01 | Emi Ltd | Sampling arrangement for nmr imaging apparatus |
| US4468621A (en) * | 1982-01-20 | 1984-08-28 | National Research Development Corporation | Investigation of samples by N.M.R. techniques |
-
1983
- 1983-11-11 JP JP58213102A patent/JPS60105979A/ja active Pending
-
1984
- 1984-11-09 GB GB08428336A patent/GB2149518B/en not_active Expired
- 1984-11-12 DE DE19843441290 patent/DE3441290A1/de active Granted
- 1984-11-13 US US06/670,579 patent/US4719581A/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5520495A (en) * | 1978-07-20 | 1980-02-13 | Univ California | Method and apparatus for mapping atomic nuceus density within object by using nucear magnetic resonance |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2149518B (en) | 1987-06-10 |
| DE3441290A1 (de) | 1985-05-23 |
| GB2149518A (en) | 1985-06-12 |
| DE3441290C2 (ja) | 1988-10-20 |
| GB8428336D0 (en) | 1984-12-19 |
| US4719581A (en) | 1988-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH07135612A (ja) | 画像信号変換装置、ｓｓｔｖ信号復調装置及び映像信号復調装置 | |
| JPS60105979A (ja) | 信号処理装置 | |
| JPS5938636A (ja) | 核磁気共鳴装置 | |
| Spangler | Radio observations of the binary stars UX Arietis and HR 1099. | |
| JPH0380497B2 (ja) | ||
| JPH0746074B2 (ja) | 真空計 | |
| Lomtev et al. | Improvement of signal compression-restoration algorithms for remote measurement systems | |
| US12359917B2 (en) | Precision gyroscope mode-matching insensitive to rate input | |
| JPH08128825A (ja) | 傾斜角検出方法及び傾斜計 | |
| JP3602078B2 (ja) | 人工衛星の姿勢検出装置 | |
| CN117686003A (zh) | 数字输出陀螺仪相频特性测试方法及系统 | |
| EP1535028A1 (en) | Methods and apparatus for generating a sinusoidal motor drive signal for a mems gyroscope | |
| US6463364B2 (en) | Offline apparatus for detecting attitude of artificial satellite in wide-band and method thereof | |
| JPH0779431B2 (ja) | ゴ−スト除去装置 | |
| JP4417042B2 (ja) | 測定する電子的対象の雑音レベルを測定する方法と装置 | |
| JP2904425B2 (ja) | スペクトラム・アナライザ | |
| Sugimoto et al. | Nuclear-spin conversion analysis of ν2+ ν4 combination band of crystalline methane in phase II | |
| JPH01191055A (ja) | クロマトグラフ用データ処理装置 | |
| JPH07218345A (ja) | 高温物体の測温装置 | |
| JPS63293466A (ja) | クロマトグラフ用デ−タ処理装置 | |
| CN113418536B (zh) | 一种基于相关信号对消的陀螺仪在轨精度评价方法和系统 | |
| Fessler et al. | Use of an infrared spectroscopic method for isotopic analysis of gaseous uranium hexafluoride | |
| Donini et al. | Low-frequency photoacoustic spectroscopy of solids | |
| JPS63302721A (ja) | デジタル型方向地絡継電器 | |
| JPH09101276A (ja) | デ−タ処理装置 |