JPH067313A - Bio-magnetic field measuring instrument - Google Patents

Bio-magnetic field measuring instrument

Info

Publication number
JPH067313A
JPH067313A JP4170538A JP17053892A JPH067313A JP H067313 A JPH067313 A JP H067313A JP 4170538 A JP4170538 A JP 4170538A JP 17053892 A JP17053892 A JP 17053892A JP H067313 A JPH067313 A JP H067313A
Authority
JP
Japan
Prior art keywords
frequency
signal
refrigerator
operating frequency
noise
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.)
Granted
Application number
JP4170538A
Other languages
Japanese (ja)
Other versions
JP3106701B2 (en
Inventor
Kenichi Sata
健一 佐多
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Daikin Industries Ltd
Original Assignee
Daikin Industries Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Daikin Industries Ltd filed Critical Daikin Industries Ltd
Priority to JP04170538A priority Critical patent/JP3106701B2/en
Priority to US08/071,004 priority patent/US5343707A/en
Publication of JPH067313A publication Critical patent/JPH067313A/en
Application granted granted Critical
Publication of JP3106701B2 publication Critical patent/JP3106701B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/035Measuring direction or magnitude of magnetic fields or magnetic flux using superconductive devices
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F25REFRIGERATION OR COOLING; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS; MANUFACTURE OR STORAGE OF ICE; LIQUEFACTION SOLIDIFICATION OF GASES
    • F25BREFRIGERATION MACHINES, PLANTS OR SYSTEMS; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS
    • F25B9/00Compression machines, plants or systems, in which the refrigerant is air or other gas of low boiling point
    • F25B9/14Compression machines, plants or systems, in which the refrigerant is air or other gas of low boiling point characterised by the cycle used, e.g. Stirling cycle
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D3/00Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
    • G01D3/028Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure
    • G01D3/032Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure affecting incoming signal, e.g. by averaging; gating undesired signals
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F25REFRIGERATION OR COOLING; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS; MANUFACTURE OR STORAGE OF ICE; LIQUEFACTION SOLIDIFICATION OF GASES
    • F25BREFRIGERATION MACHINES, PLANTS OR SYSTEMS; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS
    • F25B2500/00Problems to be solved
    • F25B2500/12Sound

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Mechanical Engineering (AREA)
  • Thermal Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Measuring Magnetic Variables (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

PURPOSE:To obtain a cardiomagnetic wave signal with high accuracy by removing a noise component with high accuracy based only on the cardiomagnetic wave measurement signal. CONSTITUTION:The operation frequency of a cryogenic refrigerator 1 is changed by an operation frequency control section 4 in response to the time the frequency of heat beats equals to the operation frequency of the cryogenic refrigerator 1. The output signal from a SQUID fluxmeter 2 is supplied in this state to a signal processing section 3. The high-accuracy cardiomagnetic wave signal, from which the noise component is removed, is thereby obtd.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は生体磁場計測装置に関
し、さらに詳細にいえば、極低温冷凍機により超伝導転
移温度以下にまで冷却されたSQUID(超伝導量子干
渉素子、Superconducting Quantum Interference Devic
e )磁束計を用いて微弱な生体磁場を計測するための装
置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a biomagnetic field measuring device, and more specifically, to a SQUID (Superconducting Quantum Interference Devic) cooled to a temperature below the superconducting transition temperature by a cryogenic refrigerator.
e) A device for measuring a weak biomagnetic field using a magnetometer.

【0002】[0002]

【従来の技術】従来より、超伝導デバイスの1つとして
ジョセフソン効果を利用したSQUIDが知られてい
る。このSQUIDに超伝導ピックアップコイルを有す
る磁束入力回路を接続することにより、例えば、生体内
に流れる微小電流に伴なう磁界や生体内の微小磁性体か
らの磁界等、極めて微弱な磁束を測定できるSQUID
磁束計を得ることができる。
2. Description of the Related Art Conventionally, SQUIDs utilizing the Josephson effect have been known as one of superconducting devices. By connecting a magnetic flux input circuit having a superconducting pickup coil to this SQUID, it is possible to measure an extremely weak magnetic flux such as a magnetic field associated with a minute current flowing in the living body or a magnetic field from a minute magnetic body in the living body. SQUID
A magnetometer can be obtained.

【0003】このSQUID磁束計を極低温レベル、つ
まりSQUIDおよび超伝導コイルが超伝導状態に転移
する温度レベルにまで冷却する方法として、低温保持容
器(クライオスタット)内に液体ヘリウムを蓄え、液体
ヘリウムにSQUID磁束計を浸漬する方法がある。
尚、この方法を冷凍機を用いて採用する場合には、低温
保持容器内に寒冷発生用の冷凍機の冷却器を挿入して、
低温保持容器内で蒸発したヘリウムガスを冷凍機により
凝縮液化させることが行なわれる。
As a method for cooling the SQUID magnetometer to a cryogenic level, that is, a temperature level at which the SQUID and the superconducting coil are transformed into a superconducting state, liquid helium is stored in a cryostat and liquid helium is stored. There is a method of immersing the SQUID magnetometer.
Incidentally, when adopting this method using a refrigerator, insert the cooler of the refrigerator for cold generation in the cryogenic container,
Helium gas evaporated in the cryogenic container is condensed and liquefied by a refrigerator.

【0004】この方法を採用した場合には、SQUID
磁束計を液体ヘリウムに浸漬するので、SQUID磁束
計を全体にわたって安定してかつ短時間で冷却すること
ができる。しかし、その反面、SQUID磁束計の冷却
のために低温保持容器内にヘリウムを介在させるため、
冷却システムが大型化し、操作性も悪くなる。この点を
考慮して、SQUID磁束計を冷凍機の冷却器に直接伝
熱可能に接触させて極低温レベルまで冷却する方法が注
目されている(例えば、特開平2−302680号公報
参照)。
When this method is adopted, the SQUID
Since the magnetometer is immersed in liquid helium, the SQUID magnetometer can be cooled stably and quickly in its entirety. However, on the other hand, since helium is interposed in the cryostat for cooling the SQUID magnetometer,
The cooling system becomes large and the operability deteriorates. In consideration of this point, a method of directly contacting an SQUID magnetometer with a cooler of a refrigerator so as to be capable of heat transfer and cooling to a cryogenic level has been attracting attention (see, for example, JP-A-2-302680).

【0005】[0005]

【発明が解決しようとする課題】SQUID磁束計を冷
凍機の冷却器に直接伝熱可能に接触させて極低温レベル
まで冷却する方法を採用した場合には、冷凍機が運転時
に振動する部分を有しているのであるから、上記振動を
完全になくすことは非常に困難である。したがって、S
QUID磁束計の出力信号に冷凍機の振動に起因する雑
音が混入してしまい、SQUID磁束計本来の信号検出
が不正確になってしまうという不都合がある。
When a method of directly contacting the SQUID magnetometer with the cooler of the refrigerator so that heat can be transferred to cool it to a cryogenic level is adopted, the portion of the refrigerator that vibrates during operation is eliminated. Since it has, it is very difficult to completely eliminate the vibration. Therefore, S
There is a disadvantage that noise resulting from vibration of the refrigerator is mixed in the output signal of the QUID magnetometer, and the original signal detection of the SQUID magnetometer becomes inaccurate.

【0006】冷凍機の振動に起因する上記不都合を解消
するために、SQUID磁束計の信号系を2チャンネル
にし、一方のチャンネルで通常の測定を行なうと同時に
他方のチャンネルで冷凍機の振動に起因する雑音を測定
し、両者の差を算出することが考えられる。しかし、こ
の方法を採用すると、信号系のチャンネル数が著しく増
加してしまうという不都合がある。
In order to eliminate the above-mentioned inconvenience caused by the vibration of the refrigerator, the signal system of the SQUID magnetometer is made to have two channels, one channel is used for normal measurement and the other channel is caused by the vibration of the refrigerator. It is possible to measure the noise that occurs and calculate the difference between the two. However, when this method is adopted, there is a disadvantage that the number of channels of the signal system is significantly increased.

【0007】また、チャンネル数を増加させることな
く、先ず通常の測定を行ない、次いで冷凍機の振動に起
因する雑音の測定を行ない、両者の差を算出すれば、チ
ャンネル数の増加という不都合は解消できるが、測定時
期が大きくずれてしまうので、通常の測定時における冷
凍機の振動特性と雑音の測定時における冷凍機の振動特
性とが同じである保証がなくなり、両者の減算により冷
凍機の振動による雑音の除去が必ずしも可能ではなくな
るという不都合がある。
If the normal measurement is first performed without increasing the number of channels and then the noise caused by the vibration of the refrigerator is measured and the difference between the two is calculated, the disadvantage of an increase in the number of channels is eliminated. However, since the measurement timing will be significantly different, there is no guarantee that the vibration characteristics of the refrigerator during normal measurement will be the same as the vibration characteristics of the refrigerator during noise measurement. There is an inconvenience that it is not always possible to remove noise due to.

【0008】これらの点を考慮して本件特許出願人は、
SQUID磁束計の出力信号を冷凍機の振動周期分ずつ
加算平均して周期性雑音のテンプレートを作り、SQU
ID磁束計の出力信号がテンプレートを減算する方法を
既に特許出願している。図4は本件特許出願人が既に特
許出願した方法を説明するための概略図であり、SQU
ID磁束計から同図(A)に示す心磁計測信号が出力さ
れた場合に、冷凍機の振動周期分ずつ加算平均すること
により同図(B)に示す信号が得られ、同図(A)の信
号から同図(B)の信号を減算することにより同図
(C)に示すように殆ど雑音のみからなる信号が得られ
る。そして、同図(C)に示す信号を冷凍機の振動周期
分ずつ加算平均することにより同図(D)に示す信号
(テンプレート)が得られる。したがって、同図(A)
の心磁計測信号から同図(D)のテンプレートを減算す
ることにより同図(E)に示すように殆ど周期性雑音を
含まない信号が得られ、同図(E)に示す信号を加算平
均することにより同図(F)に示すように雑音を殆ど含
まない心磁信号が得られる。
Considering these points, the applicant of the present patent application is
The output signal of the SQUID magnetometer is added and averaged for each vibration period of the refrigerator to create a template of periodic noise.
A patent has already been applied for a method of subtracting the template from the output signal of the ID magnetometer. FIG. 4 is a schematic diagram for explaining a method that the applicant of the present patent has already applied for a patent.
When the magnetocardiogram measurement signal shown in FIG. 7A is output from the ID magnetometer, the signal shown in FIG. 7B is obtained by adding and averaging the vibration cycles of the refrigerator. By subtracting the signal of FIG. 7B from the signal of FIG. 7B, a signal consisting of almost only noise can be obtained as shown in FIG. Then, the signal (template) shown in FIG. 7D is obtained by averaging the signals shown in FIG. 7C for each vibration cycle of the refrigerator. Therefore, the same figure (A)
By subtracting the template of (D) in the figure from the magnetocardiogram measurement signal of, a signal containing almost no periodic noise is obtained as shown in (E), and the signal shown in (E) is added and averaged. By doing so, a magnetocardiographic signal containing almost no noise can be obtained as shown in FIG.

【0009】即ち、上記の方法は、冷凍機の振動周波数
と心磁信号の周波数とが異なっているか、冷凍機の振動
周波数が安定している反面、心磁信号の周波数は不安定
であることを前提としているので、この条件を充足する
限りにおいては、加算平均処理および減算処理を行なう
ことにより高精度のテンプレートを得ることができ、ひ
いては高精度の心磁信号を得ることができる。しかし、
心磁信号の周期が必ず冷凍機の振動の周期に比べて不安
定である保証はなく、被測定者の体調等によっては心磁
信号が所定の計測期間内において殆ど変化せず、しかも
心磁信号の周期が冷凍機の振動周期とほぼ一致すること
があり、このような場合には、心磁計測信号を冷凍機の
振動周期分ずつ加算平均しても周期的振動成分が十分に
は除去できず、周期的振動成分が不十分にのみ除去され
た信号を用いてテンプレートを作成しても周期的振動成
分に対してかなり誤差を含むテンプレートしか作成でき
なくなってしまい、ひいては、最終的に得られる心磁信
号に周期的振動成分がかなり残った状態になってしまう
という不都合がある。
That is, in the above method, the vibration frequency of the refrigerator is different from the frequency of the magnetocardiographic signal, or the vibration frequency of the refrigerator is stable, but the frequency of the magnetocardiographic signal is unstable. Therefore, as long as this condition is satisfied, a high-accuracy template can be obtained by performing the averaging process and the subtraction process, and thus a high-accuracy magnetocardiographic signal can be obtained. But,
There is no guarantee that the cycle of the magnetocardiographic signal will be more stable than the cycle of the vibration of the refrigerator, and the magnetocardiographic signal hardly changes within the prescribed measurement period depending on the physical condition of the person being measured. The cycle of the signal may be almost the same as the vibration cycle of the refrigerator.In such a case, even if the magnetocardiogram measurement signal is added and averaged by the vibration cycle of the refrigerator, the periodic vibration component is sufficiently removed. However, even if a template is created using a signal in which the periodic vibration component has been removed only inadequately, only a template with a considerable error in the periodic vibration component can be created. There is an inconvenience that a considerably large periodic vibration component remains in the generated magnetocardiographic signal.

【0010】尚、この不都合は心磁計測信号に限らず、
同様な傾向を有する生体磁束計測信号であれば、同様の
不都合が生じることになる。
This inconvenience is not limited to the magnetocardiographic measurement signal,
If the biological magnetic flux measurement signals have the same tendency, the same inconvenience will occur.

【0011】[0011]

【発明の目的】この発明は上記の問題点に鑑みてなされ
たものであり、生体磁束計測信号に含まれる本来の信号
の周期が冷凍機の振動周期とほぼ等しい場合にも、上記
テンプレート法等を用いることにより本来の生体磁束計
測信号を高精度に得ることができる生体磁束計測装置を
提供することを目的としている。
SUMMARY OF THE INVENTION The present invention has been made in view of the above problems, and even when the cycle of the original signal included in the biomagnetic flux measurement signal is substantially equal to the vibration cycle of the refrigerator, the template method or the like is used. It is an object of the present invention to provide a biomagnetic flux measuring device that can obtain an original biomagnetic flux measurement signal with high accuracy by using.

【0012】[0012]

【課題を解決するための手段】上記の目的を達成するた
めの、請求項1の生体磁束計測装置は、極低温冷凍機の
運転周波数を制御する運転周波数制御手段と、計測対象
となる生体信号の固有の周波数を検出する生体信号周波
数検出手段と、生体信号の固有の周波数と極低温冷凍機
の運転周波数との差が所定の閾値周波数以下になったこ
とを検出して、両周波数の差を大きくすべく運転周波数
制御手段に対して運転周波数変更指令を供給する運転周
波数変更指令手段とを含んでいる。
In order to achieve the above object, a biological magnetic flux measuring device according to claim 1 is an operating frequency control means for controlling an operating frequency of a cryogenic refrigerator, and a biological signal to be measured. Of the biological signal frequency detecting means for detecting the natural frequency of the, the difference between the natural frequency of the biological signal and the operating frequency of the cryogenic refrigerator is below a predetermined threshold frequency, the difference between the two frequencies. And an operating frequency change command means for supplying an operating frequency change command to the operating frequency control means.

【0013】[0013]

【作用】請求項1の生体磁束計測装置であれば、極低温
を発生させ得る極低温冷凍機によりSQUID磁束計を
超伝導転移温度以下にまで冷却し、SQUID磁束計に
より得られる電気信号に基づいて極低温冷凍機に起因す
る雑音成分を除去して生体磁場を計測する場合に、計測
対象となる生体信号の固有の周波数を生体信号周波数検
出手段により検出し、生体信号の固有の周波数と極低温
冷凍機の運転周波数との差が所定の閾値周波数以下にな
ったことを検出して、両周波数の差を大きくすべく運転
周波数制御手段に対して運転周波数変更指令を供給する
ことにより極低温冷凍機の運転周波数を変化させるので
あるから、生体信号の固有の周波数と極低温冷凍機の振
動周波数がほぼ一致することに起因する雑音成分除去性
能の低下を防止でき、生体信号の固有の周波数の如何に
拘らず高精度の生体磁束計測信号を得ることができる。
According to the biological magnetic flux measuring device of claim 1, the SQUID magnetometer is cooled to a temperature below the superconducting transition temperature by a cryogenic refrigerator capable of generating an extremely low temperature, and the SQUID magnetometer is based on an electric signal obtained by the SQUID magnetometer. When measuring the biomagnetic field by removing the noise component caused by the cryogenic refrigerator, the biosignal frequency detection means detects the natural frequency of the biosignal to be measured, and By detecting that the difference between the operating frequency of the low temperature refrigerator and the operating frequency is below a predetermined threshold frequency and supplying an operating frequency change command to the operating frequency control means to increase the difference between the two frequencies, the cryogenic temperature is lowered. Since the operating frequency of the refrigerator is changed, it is possible to prevent the noise component removal performance from deteriorating due to the natural frequency of the biological signal and the vibration frequency of the cryogenic refrigerator being approximately the same. , It is possible to obtain a biological flux measurement signal unique irrespective precision of the frequency of biological signals.

【0014】[0014]

【実施例】以下、実施例を示す添付図面によって詳細に
説明する。図1はこの発明の生体磁束計測装置の一実施
例を示す概略図であり、心磁波の計測に適用した場合を
示している。この生体磁束計測装置は、複数段の冷却ス
テージを含む冷凍機1の最終冷却ステージ(例えば、約
4Kの冷却ステージ)1aに、SQUID2aおよび磁
束入力回路2bを含むSQUID磁束計2を配置し、S
QUID磁束計2からの出力信号を信号処理部3に供給
している。また、運転周波数制御部4により冷凍機1の
運転周波数を制御している。但し、1段の冷却ステージ
のみでSQUID2a等の超伝導転移温度を得ることが
できる冷凍機1を採用することが可能である。
Embodiments will be described in detail below with reference to the accompanying drawings showing embodiments. FIG. 1 is a schematic view showing an embodiment of the biomagnetic flux measuring device of the present invention, and shows a case where it is applied to measurement of a magnetocardiogram. In this biological magnetic flux measurement device, a SQUID magnetometer 2 including a SQUID 2a and a magnetic flux input circuit 2b is arranged on a final cooling stage (for example, a cooling stage of about 4K) 1a of a refrigerator 1 including a plurality of cooling stages, and S
The output signal from the QUID magnetometer 2 is supplied to the signal processing unit 3. Further, the operating frequency control unit 4 controls the operating frequency of the refrigerator 1. However, it is possible to employ the refrigerator 1 that can obtain the superconducting transition temperature such as SQUID 2a with only one cooling stage.

【0015】図2は運転周波数制御部4の構成を詳細に
示すブロック図であり、現時点における冷凍機1の運転
周波数と心拍周波数検出部4dにより検出された心拍周
波数との差周波数を算出する差周波数算出部4aと、差
周波数の絶対値が所定の閾値以下であるか否かを判別す
る判別部4bと、差周波数の絶対値が所定の閾値以下で
あることを示す判別部4bの判別結果に応答して、差周
波数の絶対値を増加させるべく冷凍機1に運転周波数指
令を供給し、差周波数の絶対値が所定の閾値より大きい
ことを示す判別部4bの判別結果に応答して、冷凍機1
に従前の運転周波数指令を供給する運転周波数指令供給
部4cとを有している。尚、運転周波数指令は、例え
ば、冷凍機1のバルブモータに供給する電源周波数を所
定周波数に設定すべきことを指示するものである。
FIG. 2 is a block diagram showing the configuration of the operating frequency control unit 4 in detail. The difference for calculating the difference frequency between the operating frequency of the refrigerator 1 at the present time and the heartbeat frequency detected by the heartbeat frequency detecting unit 4d is shown. The frequency calculation unit 4a, the determination unit 4b that determines whether the absolute value of the difference frequency is less than or equal to a predetermined threshold, and the determination result of the determination unit 4b that indicates that the absolute value of the difference frequency is less than or equal to the predetermined threshold value. In response to, the operating frequency command is supplied to the refrigerator 1 to increase the absolute value of the difference frequency, and in response to the determination result of the determination unit 4b indicating that the absolute value of the difference frequency is larger than a predetermined threshold value, Refrigerator 1
The operating frequency command supply unit 4c for supplying the previous operating frequency command. The operating frequency command instructs, for example, that the power supply frequency supplied to the valve motor of the refrigerator 1 should be set to a predetermined frequency.

【0016】即ち、心磁信号の周波数は心拍周波数(通
常0.8〜2Hz)として検出可能であり、しかも心拍
周波数は心電図等により正確に検出できる。また、冷凍
機1の現時点における運転周波数(通常2〜2.4H
z)は例えばバルブモータの電源周波数により正確に検
出できる。したがって、両周波数の差周波数を差周波数
算出部4aにより算出し、差周波数の絶対値が所定の閾
値以下(通常、所定の閾値は著しく小さい値に設定され
るのであるから、差周波数の絶対値が著しく小さい)で
あると判別部4bにより判別された場合にのみ、差周波
数の絶対値を増加させるべく運転周波数指令供給部4c
により冷凍機1に新たな運転周波数指令を供給し、心拍
周波数と冷凍機1の運転周波数とがほぼ一致する状態を
排除できる。
That is, the frequency of the magnetocardiographic signal can be detected as a heartbeat frequency (usually 0.8 to 2 Hz), and the heartbeat frequency can be accurately detected by an electrocardiogram or the like. In addition, the operating frequency of the refrigerator 1 at the present time (usually 2 to 2.4H
z) can be accurately detected by the power supply frequency of the valve motor, for example. Therefore, the difference frequency between the two frequencies is calculated by the difference frequency calculation unit 4a, and the absolute value of the difference frequency is equal to or less than a predetermined threshold value (usually, the predetermined threshold value is set to a significantly small value. Only when the determination unit 4b determines that the operating frequency command supply unit 4c increases the absolute value of the difference frequency.
Thus, a new operating frequency command can be supplied to the refrigerator 1, and a state in which the heartbeat frequency and the operating frequency of the refrigerator 1 substantially match can be eliminated.

【0017】図3は信号処理部3の構成を示すブロック
図であり、心磁波計測信号を心磁波の所定のタイミング
に合せて加算平均することにより周期性雑音および非周
期性雑音が除去された心磁波信号を得る第1加算平均部
3aと、心磁波計測信号から第1加算平均部3aにより
得られた心磁波信号を減算して雑音信号を得る第1減算
部3bと、第1減算部3bにより得られる雑音信号を冷
凍機1の振動周期分ずつ加算平均して雑音テンプレート
を得る第2加算平均部3cと、心磁波計測信号から雑音
テンプレートを減算する第2減算部3dと、第2減算部
3dにより得られた信号を心磁波の所定のタイミングに
合せて加算平均して雑音成分が除去された心磁波信号を
得る第3加算平均部3eとを有している。
FIG. 3 is a block diagram showing the configuration of the signal processing unit 3. The periodic noise and the aperiodic noise are removed by adding and averaging the cardiac magnetic wave measurement signals in accordance with the predetermined timing of the cardiac magnetic wave. A first arithmetic mean unit 3a for obtaining a cardiac magnetic wave signal, a first subtractor unit 3b for subtracting a cardiac magnetic wave signal obtained by the first arithmetic mean unit 3a from a cardiac magnetic wave measurement signal to obtain a noise signal, and a first subtractor unit A second addition and averaging unit 3c that obtains a noise template by adding and averaging the noise signal obtained by 3b for each vibration period of the refrigerator 1, a second subtraction unit 3d that subtracts the noise template from the magnetocardiogram measurement signal, and a second It has a third arithmetic mean unit 3e which obtains a magnetocardiographic wave signal from which noise components are removed by arithmetically averaging the signals obtained by the subtraction unit 3d at a predetermined timing of the magnetocardiogram wave.

【0018】図4は図3の信号処理部3における処理を
説明する信号波形図であり、本件特許出願人が既に特許
出願した処理を示す信号波形と同じである。SQUID
磁束計2から同図(A)に示す心磁波計測信号が得られ
た場合に、第1加算平均部3aにより心磁波計測信号を
心磁波の所定のタイミングに合せて加算平均すれば、同
図(B)に示すように、周期性雑音および非周期性雑音
がかなり除去された心磁波信号を得ることができる。し
たがって、第1減算部3bにより同図(A)の心磁波計
測信号から同図(B)の心磁波信号を減算すれば、同図
(C)に示すように、かなり高精度の雑音信号を得るこ
とができる。次いで、第2加算平均部3cにより同図
(C)の雑音信号を冷凍機1の振動周期分ずつ加算平均
して、同図(D)に示すように、周期性雑音のみからな
る雑音テンプレートを得ることができる。その後、第2
減算部3dにより同図(A)の心磁波計測信号から同図
(D)の雑音テンプレートを減算すれば、同図(E)に
示すように、周期性雑音成分が高精度に除去された心磁
波信号を得ることができる。最後に、第3加算平均部3
eにより同図(E)の心磁波信号を心磁波の所定のタイ
ミングに合せて加算平均すれば、同図(F)に示すよう
に、非周期性雑音成分が除去された心磁波信号を得るこ
とができる。
FIG. 4 is a signal waveform diagram for explaining the processing in the signal processing unit 3 of FIG. 3, which is the same as the signal waveform showing the processing already filed by the applicant of the present patent application. SQUID
When the magnetocardiographic wave measurement signal shown in FIG. 1A is obtained from the magnetometer 2, if the electrocardiographic wave measurement signal is arithmetically averaged by the first arithmetic mean unit 3a in accordance with a predetermined timing of the magnetocardiographic wave. As shown in (B), it is possible to obtain a magnetocardiogram signal in which periodic noise and aperiodic noise are significantly removed. Therefore, when the first subtraction unit 3b subtracts the magnetocardiogram signal shown in FIG. 7B from the magnetocardiogram measurement signal shown in FIG. 9A, a noise signal with considerably high accuracy is obtained as shown in FIG. Obtainable. Next, the noise signal of FIG. 6C is added and averaged by the second averaging unit 3c for each vibration cycle of the refrigerator 1, and a noise template consisting of only periodic noise is generated as shown in FIG. Obtainable. Then the second
By subtracting the noise template shown in FIG. 7D from the magnetocardiogram measurement signal shown in FIG. 7A by the subtraction unit 3d, the heart from which the periodic noise component has been removed with high precision as shown in FIG. A magnetic wave signal can be obtained. Finally, the third averaging unit 3
By performing the averaging of the magnetocardiographic wave signals of FIG. 6E in accordance with e in accordance with e, the magnetocardiographic wave signals from which the aperiodic noise components have been removed are obtained as shown in FIG. be able to.

【0019】以上のように信号処理部3は、周期性雑音
の周波数を規定する冷凍機1の運転周波数と心拍周波数
(心磁波信号周波数と等しい)とが一致していないこと
を前提として高精度の雑音除去を達成するものであり、
この条件を充足させるために運転周波数制御部4により
冷凍機1の運転周波数を心拍周波数と十分に異なる周波
数に設定するようにしているので、信号処理部3による
高精度の雑音除去を達成でき、高精度の心磁波信号を得
ることができる。
As described above, the signal processing unit 3 is highly accurate on the assumption that the operating frequency of the refrigerator 1 defining the frequency of the periodic noise does not match the heartbeat frequency (equal to the magnetocardiogram signal frequency). To achieve the noise elimination of
In order to satisfy this condition, the operating frequency control unit 4 sets the operating frequency of the refrigerator 1 to a frequency sufficiently different from the heartbeat frequency, so that highly accurate noise removal by the signal processing unit 3 can be achieved. A highly accurate magnetocardiogram signal can be obtained.

【0020】さらに詳細に説明すると、心拍周波数は一
般的にかなり変動しているのであるから、通常は冷凍機
1の運転周波数を変化させなくても信号処理部3により
高精度の雑音除去を達成できるのであるが、約1分程度
の心磁波計測期間において心拍周波数が一定になり、し
かも一定になった心拍周波数が冷凍機1の運転周波数と
ほぼ一致する場合があり、このような場合には単に加算
平均を反復しても雑音を十分には除去できないことにな
る。しかし、この実施例においては、心拍周波数が冷凍
機1の運転周波数とほぼ一致した場合に、両周波数の差
が大きくなるように冷凍機1の運転周波数を設定するよ
うにしているのであるから、心拍周波数の状態に拘らず
信号処理部3による高精度の雑音除去を達成できる。
More specifically, since the heartbeat frequency generally fluctuates considerably, the signal processing unit 3 normally achieves highly accurate noise removal without changing the operating frequency of the refrigerator 1. Although it is possible, the heartbeat frequency becomes constant during the magnetocardiogram measurement period of about 1 minute, and the constant heartbeat frequency may be almost the same as the operating frequency of the refrigerator 1. In such a case, The noise cannot be sufficiently removed by simply repeating the averaging. However, in this embodiment, when the heartbeat frequency substantially matches the operating frequency of the refrigerator 1, the operating frequency of the refrigerator 1 is set so that the difference between the two frequencies becomes large. It is possible to achieve highly accurate noise removal by the signal processing unit 3 regardless of the state of the heartbeat frequency.

【0021】尚、この発明は上記の実施例に限定される
ものではなく、例えば、両周波数の差周波数の絶対値と
所定の閾値との大小を判別する代わりに心拍周波数と冷
凍機1の運転周波数とが一致するか否かを判別すること
が可能であるほか、心磁波以外の生体磁束であって、心
磁波と同様に周期性を有する生体磁束であっても同様に
適用することが可能であり、さらに、冷凍機1の運転周
波数を手動で変化させることが可能であるほか、この発
明の要旨を変更しない範囲内において種々の設計変更を
施すことが可能である。
The present invention is not limited to the above-described embodiment. For example, instead of discriminating the magnitude of the absolute value of the difference frequency between the two frequencies and the predetermined threshold value, the heartbeat frequency and the operation of the refrigerator 1 are changed. It is possible to determine whether or not the frequencies match, and it is also applicable to biomagnetic fluxes other than cardiac magnetic waves that have periodicity similar to cardiac magnetic waves. In addition, the operating frequency of the refrigerator 1 can be manually changed, and various design changes can be made within the range not changing the gist of the present invention.

【0022】[0022]

【発明の効果】以上のように請求項1の発明は、生体信
号の固有の周波数と極低温冷凍機の振動周波数がほぼ一
致することに起因する雑音成分除去性能の低下を防止で
き、生体信号の固有の周波数の如何に拘らず高精度の生
体磁束計測信号を得ることができるという特有の効果を
奏する。
As described above, according to the first aspect of the present invention, it is possible to prevent the deterioration of the noise component removal performance caused by the natural frequency of the biomedical signal and the vibration frequency of the cryogenic refrigerator being substantially equal to each other. A unique effect that a highly accurate biomagnetic flux measurement signal can be obtained regardless of the natural frequency of

【図面の簡単な説明】[Brief description of drawings]

【図1】この発明の生体磁束計測装置の一実施例を示す
概略図である。
FIG. 1 is a schematic view showing an embodiment of a biological magnetic flux measuring device of the present invention.

【図2】運転周波数制御部の構成を詳細に示すブロック
図である。
FIG. 2 is a block diagram showing in detail the configuration of an operating frequency control unit.

【図3】信号処理部の構成を示すブロック図である。FIG. 3 is a block diagram showing a configuration of a signal processing unit.

【図4】図3の信号処理部における処理を説明する信号
波形図である。
FIG. 4 is a signal waveform diagram illustrating processing in the signal processing unit in FIG.

【符号の説明】[Explanation of symbols]

1 冷凍機 2 SQUID磁束計 4a 差周波数算出部 4b 判別部 4c 運転周波数指令供給部 4d 心拍周波数検出
1 Refrigerator 2 SQUID magnetometer 4a Difference frequency calculation part 4b Discrimination part 4c Operating frequency command supply part 4d Heartbeat frequency detection part

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】 極低温を発生させ得る極低温冷凍機
(1)によりSQUID磁束計(2)を超伝導転移温度
以下にまで冷却し、SQUID磁束計(2)により得ら
れる電気信号に基づいて極低温冷凍機(1)に起因する
雑音成分を除去して生体磁場を計測する装置であって、
極低温冷凍機(1)の運転周波数を制御する運転周波数
制御手段(4c)と、計測対象となる生体信号の固有の
周波数を検出する生体信号周波数検出手段(4d)と、
生体信号の固有の周波数と極低温冷凍機(1)の運転周
波数との差が所定の閾値周波数以下になったことを検出
して、両周波数の差を大きくすべく運転周波数制御手段
(4c)に対して運転周波数変更指令を供給する運転周
波数変更指令手段(4a)(4b)とを含むことを特徴
とする生体磁場計測装置。
1. A SQUID magnetometer (2) is cooled to a temperature below a superconducting transition temperature by a cryogenic refrigerator (1) capable of generating an extremely low temperature, and based on an electric signal obtained by the SQUID magnetometer (2). A device for measuring a biomagnetic field by removing a noise component caused by a cryogenic refrigerator (1),
An operating frequency control means (4c) for controlling the operating frequency of the cryogenic refrigerator (1), a biological signal frequency detecting means (4d) for detecting the unique frequency of the biological signal to be measured,
The operating frequency control means (4c) detects that the difference between the natural frequency of the biological signal and the operating frequency of the cryogenic refrigerator (1) is less than or equal to a predetermined threshold frequency, and increases the difference between the two frequencies. And a driving frequency change command means (4a) (4b) for supplying a driving frequency change command to the biomagnetic field measuring apparatus.
JP04170538A 1992-06-29 1992-06-29 Biomagnetic field measurement device Expired - Fee Related JP3106701B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP04170538A JP3106701B2 (en) 1992-06-29 1992-06-29 Biomagnetic field measurement device
US08/071,004 US5343707A (en) 1992-06-29 1993-06-04 Methods and apparatus for removing cyclic noise from the output signal of a magnetic sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04170538A JP3106701B2 (en) 1992-06-29 1992-06-29 Biomagnetic field measurement device

Publications (2)

Publication Number Publication Date
JPH067313A true JPH067313A (en) 1994-01-18
JP3106701B2 JP3106701B2 (en) 2000-11-06

Family

ID=15906768

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04170538A Expired - Fee Related JP3106701B2 (en) 1992-06-29 1992-06-29 Biomagnetic field measurement device

Country Status (2)

Country Link
US (1) US5343707A (en)
JP (1) JP3106701B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10290792A (en) * 1997-04-18 1998-11-04 Toshiba Corp Biological magnetic field measuring instrument
JP2000217798A (en) * 1999-02-02 2000-08-08 Hitachi Ltd Biological magnetic field measuring instrument
JP2015183514A (en) * 2014-03-25 2015-10-22 ハム アーゲーHamm AG Method for correction of measured value chart by eliminating periodically occurring measurement artefact, in particular in soil compactor
US10918293B2 (en) 2016-03-03 2021-02-16 Ricoh Company, Ltd. Magnetic measuring apparatus

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5574369A (en) * 1994-08-19 1996-11-12 Hibbs; Andrew D. Detection and communications device employing stochastic resonance
JP2829375B2 (en) * 1995-03-13 1998-11-25 九州大学長 Weak magnetic field measuring device and measuring method
DE19520506A1 (en) * 1995-06-03 1996-12-05 Bosch Gmbh Robert Cooling device for a high frequency receiver
US5640124A (en) * 1995-10-20 1997-06-17 Massachusetts Institute Of Technology System and method of eliminating systematic noise in stimulus-response systems
US5836165A (en) * 1996-10-30 1998-11-17 Hughes Electronics Adaptive feedforward vibration control system and method
US5765378A (en) * 1996-12-31 1998-06-16 Helix Technology Corporation Method and apparatus for detecting a loss of differential pressure in a cryogenic refrigerator
JP4193382B2 (en) * 2001-07-19 2008-12-10 株式会社日立製作所 Magnetic field measuring device
KR100400755B1 (en) * 2001-11-08 2003-10-08 엘지전자 주식회사 SQUID sensor using secondary sensor
JP4972443B2 (en) * 2007-03-30 2012-07-11 株式会社日立ハイテクノロジーズ Biomagnetic field measurement device
US8282261B2 (en) 2009-06-01 2012-10-09 Apple, Inc. White point adjustment for multicolor keyboard backlight
US8451146B2 (en) 2010-06-11 2013-05-28 Apple Inc. Legend highlighting
US9275810B2 (en) 2010-07-19 2016-03-01 Apple Inc. Keyboard illumination
US8378857B2 (en) 2010-07-19 2013-02-19 Apple Inc. Illumination of input device
CN101999895B (en) * 2010-11-15 2015-09-16 中国科学院上海微系统与信息技术研究所 A kind of multi-point scanning navigation system for magnetocardiograph and localization method
CN106679217B (en) * 2016-12-16 2020-08-28 复旦大学 Mechanical vibration isolation liquid helium recondensation low-temperature refrigeration system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0361137A1 (en) * 1988-09-16 1990-04-04 Siemens Aktiengesellschaft Magnetometer device with a cryostat for the measurement of weak magnetic fields
EP0359864B1 (en) * 1988-09-23 1993-12-01 Siemens Aktiengesellschaft Apparatus and method for the measurement of weak magnetic fields dependent upon position and time
JPH02302680A (en) * 1989-05-17 1990-12-14 Daikin Ind Ltd Gradiometer cooling apparatus
JPH0443978A (en) * 1990-06-11 1992-02-13 Seiko Instr Inc High sensitivity magnetic field detecting device
FI912656A (en) * 1990-06-25 1991-12-26 Siemens Ag KYLANORDNING FOER EN SQUID-MAETANORDNING.

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10290792A (en) * 1997-04-18 1998-11-04 Toshiba Corp Biological magnetic field measuring instrument
JP2000217798A (en) * 1999-02-02 2000-08-08 Hitachi Ltd Biological magnetic field measuring instrument
JP2015183514A (en) * 2014-03-25 2015-10-22 ハム アーゲーHamm AG Method for correction of measured value chart by eliminating periodically occurring measurement artefact, in particular in soil compactor
US10338095B2 (en) 2014-03-25 2019-07-02 Hamm Ag Method for the correction of a measured value curve by eliminating periodically occurring measurement artifacts, in particular in a soil compactor
US10918293B2 (en) 2016-03-03 2021-02-16 Ricoh Company, Ltd. Magnetic measuring apparatus

Also Published As

Publication number Publication date
US5343707A (en) 1994-09-06
JP3106701B2 (en) 2000-11-06

Similar Documents

Publication Publication Date Title
JP3106701B2 (en) Biomagnetic field measurement device
Ruha et al. A real-time microprocessor QRS detector system with a 1-ms timing accuracy for the measurement of ambulatory HRV
JPH0966037A (en) Magnetism of living body measuring instrument
JP2001511021A (en) An improved method for compensating for thermal drift during cardiac output measurements
JPH07246194A (en) Mr system
Hukkinen et al. Instantaneous fetal heart rate monitoring by electromagnetic methods
JP3134499B2 (en) Method and apparatus for removing noise from SQUID magnetometer output
JP3236710B2 (en) Measurement device for RMS values
JP3230460B2 (en) Biomagnetic measurement device
JP4543327B2 (en) SQUID sensor device
JP2579280B2 (en) Calibration method for SQUID magnetometer
JP3576658B2 (en) Driving waveform generation device for gradient magnetic field, eddy current estimation device, and MRI device
JPH02107230A (en) Method and device for compensating eddy current of mr device
JP4373542B2 (en) Biomagnetic measurement device
JPH04160380A (en) Digital squid flux meter
JPS6323785B2 (en)
CA2379800A1 (en) Method and device for measuring biomagnetic and in particular cardiomagnetic fields
Lee et al. Double relaxation oscillation SQUID systems for biomagnetic multichannel measurements
Kazami et al. A Drung-type magnetometer mounted on a GM cryocooler
JP4179576B2 (en) Magnetic resonance imaging apparatus and imaging parameter calculation method
JP2752884B2 (en) Life activity current source estimation method
JP2004327874A (en) Cooling system and biomagnetism measuring apparatus equipped with coolant consumption amount monitoring function
JPH0435643A (en) Measuring method and equipment for electromagnetic waves generated from living body
CN117562543A (en) Active magnetic compensation method and system
JP4010645B2 (en) Power supply apparatus and magnetic resonance imaging apparatus using the same

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees