JP2012141282A - 気圧推定方法及び気圧推定装置 - Google Patents
気圧推定方法及び気圧推定装置 Download PDFInfo
- Publication number
- JP2012141282A JP2012141282A JP2011194655A JP2011194655A JP2012141282A JP 2012141282 A JP2012141282 A JP 2012141282A JP 2011194655 A JP2011194655 A JP 2011194655A JP 2011194655 A JP2011194655 A JP 2011194655A JP 2012141282 A JP2012141282 A JP 2012141282A
- Authority
- JP
- Japan
- Prior art keywords
- atmospheric pressure
- oscillator
- oscillation frequency
- phase difference
- frequency
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000010355 oscillation Effects 0.000 claims abstract description 177
- 238000009530 blood pressure measurement Methods 0.000 claims abstract description 22
- 230000007613 environmental effect Effects 0.000 claims description 8
- 239000013078 crystal Substances 0.000 abstract description 101
- 238000012545 processing Methods 0.000 abstract description 46
- 230000008569 process Effects 0.000 description 28
- 238000003860 storage Methods 0.000 description 22
- 230000008859 change Effects 0.000 description 16
- 238000001514 detection method Methods 0.000 description 9
- 239000010453 quartz Substances 0.000 description 8
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 8
- 238000012935 Averaging Methods 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000003416 augmentation Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L9/00—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means
- G01L9/0001—Transmitting or indicating the displacement of elastically deformable gauges by electric, electro-mechanical, magnetic or electro-magnetic means
- G01L9/0008—Transmitting or indicating the displacement of elastically deformable gauges by electric, electro-mechanical, magnetic or electro-magnetic means using vibrations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L9/00—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means
- G01L9/08—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means by making use of piezoelectric devices, i.e. electric circuits therefor
- G01L9/085—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means by making use of piezoelectric devices, i.e. electric circuits therefor with temperature compensating means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
【課題】気圧測定用発振器を用いた新たな気圧推定手法の提案。
【解決手段】水晶振動子型気圧センサー1において、処理部10は、基準発振器30の発振周波数を気圧測定用の水晶発振器20の発振周波数に合わせるように制御する。そして、処理部10は、基準発振器30の基準発振信号と水晶発振器20の水晶発振信号との位相差をコスタスループを用いて検出する。そして、処理部10は、基準発振器30の発振周波数と位相差とを用いて気圧を推定する。
【選択図】図1
【解決手段】水晶振動子型気圧センサー1において、処理部10は、基準発振器30の発振周波数を気圧測定用の水晶発振器20の発振周波数に合わせるように制御する。そして、処理部10は、基準発振器30の基準発振信号と水晶発振器20の水晶発振信号との位相差をコスタスループを用いて検出する。そして、処理部10は、基準発振器30の発振周波数と位相差とを用いて気圧を推定する。
【選択図】図1
Description
本発明は、気圧推定方法及び気圧推定装置に関する。
大気の圧力を測定するセンサーとして、水晶振動子を利用した気圧センサー(水晶振動子型気圧センサー)が考案されている。この水晶振動子型気圧センサーは、水晶振動子(発振器)の発振周波数が圧力によって変化する現象を利用して、気圧を推定するものである(例えば、特許文献1)。
圧力によって発振周波数が変化する発振器を用いて気圧を推定する場合には、その発振器の発振周波数を正しく測定する必要がある。しかし、発振周波数を正しく測定するためには、高精度のクロックや内部タイムベース等、周波数を高精度に測定するための機構が必要になるという問題があった。
本発明は上述した課題に鑑みて為されたものであり、その目的とするところは、気圧測定用発振器を用いた新たな気圧推定手法を提案することにある。
以上の課題を解決するための第1の形態は、所定の発振器の発振周波数を気圧測定用発振器の発振周波数に合わせるように制御することと、前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出することと、前記発振器の発振周波数と前記位相差とを用いて気圧を推定することと、を含む気圧推定方法である。
また、他の形態として、気圧に応じて発振周波数が変化する気圧測定用発振器と、発振周波数を変更可能な発振器と、前記発振器の発振周波数を前記気圧測定用発振器の発振周波数に合わせるように制御する周波数制御部と、前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出する位相差検出部と、前記発振器の発振周波数と前記位相差とを用いて気圧を推定する気圧推定部と、を備えた気圧推定装置を構成してもよい。
この第1の形態等によれば、所定の発振器の発振周波数を気圧測定用発振器の発振周波数に合わせるように制御する。そして、これらの発振器の発振信号の位相差を検出し、発振器の発振周波数と位相差とを用いて気圧を推定する。これにより、気圧測定用発振器の発振周波数を直接正確に測定することなく、気圧を推定することができる。
また、第2の形態として、第1の形態の気圧推定方法であって、前記発振器の発振周波数を制御することは、前記位相差に基づいて制御することを含む、気圧推定方法を構成してもよい。
この第2の形態によれば、例えば、発振器の発振信号と気圧測定用発振器の発振信号との位相差に基づく位相同期手法により、発振器の発振周波数を適切に制御することができる。
また、第3の形態として、第1又は第2の形態の気圧推定方法であって、前記位相差を検出することは、コスタスループを用いて前記位相差を検出することを含む、気圧推定方法を構成してもよい。
この第3の形態によれば、コスタスループを用いて位相差の検出を簡単に行うことができる。
また、第4の形態として、第1〜第3の何れかの形態の気圧推定方法であって、前記気圧を推定することは、前記発振器の発振周波数を用いて第1の分解能で気圧を推定することと、前記位相差を用いて前記発振器の1Hz以下に相当する気圧を前記第1の分解能より高い第2の分解能で推定することと、を含む、気圧推定方法を構成してもよい。
この第4の形態によれば、発振器の発振周波数を用いて第1の分解能で気圧を推定するとともに、検出した位相差を用いて発振器の1Hz以下に相当する気圧を第1の分解能より高い第2の分解能で推定する。かかる手法により、精細な気圧推定を実現することができる。
また、第5の形態として、第1〜第4の何れかの形態の気圧推定方法であって、前記気圧を推定することは、所定時間の間に検出された前記位相差を用いて前記気圧を推定することを含み、前記気圧測定用発振器の発振周波数の安定度を判定することと、前記安定度に基づいて前記所定時間を設定することと、を更に含む、気圧推定方法を構成してもよい。
この第5の形態によれば、所定時間の間に検出された位相差を用いて気圧を推定する。そして、気圧測定用発振器の発振周波数の安定度を判定し、当該安定度に基づいて上記所定時間を設定する。気圧測定用発振器の発振周波数の安定度が高くなる時間を所定時間として設定することで、気圧推定のタイミングを適正化することができる。
また、第6の形態として、第5の形態の気圧推定方法であって、前記安定度を判定することは、前記気圧測定用発振器の発振周波数のアラン分散を求めることを含む、気圧推定方法を構成してもよい。
この第6の形態によれば、気圧測定用発振器の発振周波数のアラン分散を求めることで、気圧測定用発振器の発振周波数の安定度を適切に判定することができる。
また、第7の形態として、第1〜第6の何れかの形態の気圧推定方法であって、前記気圧を推定することは、環境温度に基づいて前記気圧を温度補償することを含む、気圧推定方法を構成してもよい。
この第7の形態によれば、環境温度に基づいて気圧を温度補償することで、環境温度に依らずに気圧推定を適切に行うことができる。
以下、図面を参照して、本発明の好適な実施形態の一例について説明する。本実施形態は、気圧推定装置を具備した水晶振動子型気圧センサーの実施形態である。但し、本発明を適用可能な実施形態が以下説明する実施形態に限定されるわけでないことは勿論である。
1.機能構成
図1は、本実施形態における水晶振動子型気圧センサー1の機能構成の一例を示す図である。水晶振動子型気圧センサー1は、処理部10と、水晶発振器20と、基準発振器30と、第1の乗算器40と、遅延回路50と、第2の乗算器60と、記憶部80とを備えて構成される。
図1は、本実施形態における水晶振動子型気圧センサー1の機能構成の一例を示す図である。水晶振動子型気圧センサー1は、処理部10と、水晶発振器20と、基準発振器30と、第1の乗算器40と、遅延回路50と、第2の乗算器60と、記憶部80とを備えて構成される。
処理部10は、水晶振動子型気圧センサー1の各部を統括的に制御する制御装置及び演算装置であり、CPU(Central Processing Unit)等のプロセッサーを有して構成される。処理部10は、主要な機能部として、位相比較部11と、ループフィルター処理部13と、気圧推定部15と、推定時間間隔較正部17と、カウンタ部19とを有する。
基準発振器30から出力される基準発振信号は、そのままの信号をI相の基準発振信号とし、遅延回路50を通過させることによってQ相の基準発振信号とされる。そして、水晶発振器20から出力される水晶発振信号と、I相及びQ相の基準発振信号とが第1の乗算器40及び第2の乗算器60で乗算されてI相の乗算結果信号及びQ相の乗算結果信号となり、位相比較部11にそれぞれ入力される。
位相比較部11は、水晶発振器20の水晶発振信号の位相と、基準発振器30の基準発振信号の位相とを比較する。具体的には、第1の乗算器40から出力されるI相の乗算結果信号と、第2の乗算器60から出力されるQ相の乗算結果信号とを乗算し、その乗算結果をループフィルター処理部13に出力する。位相比較部11は、位相同期回路として知られるPLL(Phase Locked Loop)回路における位相比較器に相当する機能ブロックである。
ループフィルター処理部13は、位相比較部11の比較結果に対する平均化処理を行い、比較結果を安定的な値(信号)に変換する。位相比較部11とループフィルター処理部13とで、基準発振器30の基準発振信号と水晶発振器20の水晶発振信号との位相差を検出する位相差検出部が構成される。
また、検出された位相差でもって基準発振器30の発振周波数が制御される。具体的には、基準発振器30の発振周波数と水晶発振器20の発振周波数との差が位相差として検出され、水晶発振器20の発振周波数に合わせるように基準発振器30の発振周波数が制御される。つまり、位相比較部11及びループフィルター処理部13は周波数制御部でもある。以下の説明では、水晶発振器20の発振周波数を「水晶発振周波数」と称し、基準発振器30の発振周波数を「基準発振周波数」と称する。
気圧推定部15は、記憶部80に記憶された気圧推定プログラム81に従って気圧を推定する。具体的には、位相比較部11の比較結果、或いは、ループフィルター処理部13により平均化された後の値に基づいて、水晶発振周波数の基準発振周波数からのずれの大きさを、1Hz以下の所定の精度で算出する。随時検出される位相差の所定の単位時間分の時間変化に基づいて、水晶発振周波数と基準発振周波数との周波数ずれの大きさを見積もるのである。
そして、所定の推定タイミングにおいて、基準発振周波数と算出した周波数ずれとを用いて水晶発振周波数を推定する。水晶発振周波数を推定したならば、水晶発振周波数と気圧とを対応付ける関係式或いは参照テーブルを用いて、水晶発振周波数の推定値を気圧に換算する。
本実施形態では、周波数ずれを算出する時間間隔のことを「単位時間間隔」と定義し、単位時間間隔に相当する時間のことを「単位時間」と定義する。また、水晶発振周波数及び気圧の推定を行う時間間隔のことを「推定時間間隔」と定義し、推定時間間隔に相当する時間のことを「推定時間」と定義する。単位時間間隔は、例えば「1ミリ秒」とされる。この場合は、1ミリ秒毎に周波数ずれが算出される。また、推定時間間隔は、推定時間間隔較正部17により、単位時間間隔よりも長い時間が設定される。
推定時間間隔較正部17は、記憶部80に記憶された推定時間間隔較正プログラム811に従って推定時間間隔を較正する。より具体的には、水晶発振周波数の安定度をアラン分散を利用して判定し、その判定結果に基づいて推定時間間隔を較正・設定する。
カウンタ部19は、外部入力されるクロック信号に基づいて水晶発振周波数の初期値を測定する。具体的には、電源投入後の初期設定時に、水晶発振器20の水晶発振信号を入力して、水晶発振周波数の概測値を測定する。概測値は、水晶発振信号の周波数の大まかな値(例えば周波数の整数部分)であればよい。そして、測定した概測値で基準発振器30の発振周波数を初期設定する。
水晶発振器20は、気圧に応じて発振周波数が変化する気圧測定用発振器であり、例えば水晶振動子と水晶発振回路とを搭載した水晶デバイスとして構成される。水晶振動子と水晶発振回路とをパッケージ化して1チップとして製造することも可能である。水晶振動子としては、例えば双音叉型水晶振動子を用いることができる。
基準発振器30は、水晶発振信号に対する同期用の発振器であり、水晶発振周波数に合わせるように基準発振周波数が制御される。基準発振器30は、発振周波数を変更可能に構成された可変周波数発振器であり、例えばNCO(Numerical Controlled Oscillator)により構成される。
第1の乗算器40は、水晶発振器20の水晶発振信号と、基準発振器30の基準発振信号とを乗算する乗算器である。水晶発振信号と基準発振信号とが乗算されることで、水晶発振周波数と基準発振周波数との周波数差の信号(I相乗算結果信号)に変換される。I相乗算結果信号は位相比較部11に出力される。
遅延回路50は、基準発振器30の基準発振信号を90度(=π/2)だけ時間的に遅延させる遅延回路であり、遅延させた基準発振信号を直交基準発振信号として第2の乗算器60に出力する。
第2の乗算器60は、水晶発振器20の水晶発振信号と直交基準発振信号とを乗算する乗算器である。水晶発振信号と直交基準発振信号とが乗算されることで、水晶発振周波数と直交基準発振周波数との周波数差の信号(Q相乗算結果信号)に変換される。Q相乗算結果信号は位相比較部11に出力される。
第1の乗算器40→位相比較部11→ループフィルター処理部13→基準発振器30→第1の乗算器40の第1のループと、第2の乗算器60→位相比較部11→ループフィルター処理部13→基準発振器30→遅延回路50→第2の乗算器60の第2のループとでコスタスループが形成される。本実施形態では、このコスタスループにより位相差を検出する。
記憶部80は、処理部10が水晶振動子型気圧センサー1の各部を統括的に制御するためのシステムプログラムや、気圧推定処理、推定時間間隔較正処理といった各種の処理を実行するための各種プログラムやデータ等を記憶した記憶装置である。記憶部80は、ROM(Read Only Memory)やフラッシュROM、RAM(Random Access Memory)等のメモリーを有して構成される。
2.原理
図2は、本実施形態における気圧推定の原理の説明図である。図2において、横軸は気圧“P”を示し、縦軸は水晶発振周波数“f”を示す。また、温度“T1〜T4”(T1<T2<T3<T4)それぞれにおける気圧対水晶発振周波数を、黒四角形、黒三角形、バツ印及び黒丸でそれぞれプロットしている。
図2は、本実施形態における気圧推定の原理の説明図である。図2において、横軸は気圧“P”を示し、縦軸は水晶発振周波数“f”を示す。また、温度“T1〜T4”(T1<T2<T3<T4)それぞれにおける気圧対水晶発振周波数を、黒四角形、黒三角形、バツ印及び黒丸でそれぞれプロットしている。
このグラフを見ると、温度に関わらず、気圧“P”と水晶発振周波数“f”とは、ほぼ線形の関係にあることがわかる。また、温度が高くなるほど、水晶発振周波数“f”の総体的な大きさは小さくなる傾向がある。かかる気圧対水晶発振周波数の関係から、環境温度と水晶発振周波数の推定値とを用いて気圧を推定することができる。
処理部10は、コスタスループを用いて検出した位相差“Δθ”の単位時間分の時間変化に基づいて、動作条件の変化等に起因する水晶発振周波数の基準発振周波数からの“ずれ”を算出する。この水晶発振周波数の基準発振周波数からのずれ“Δf”を「周波数ずれ」と定義する。また、周波数ずれ“Δf”を基準発振周波数“fb”で除算した値“Δf/fb”のことを「周波数偏差」と定義する。
式(1)から、1周期(360°)に相当する位相差“Δθ”の時間変化が“1[Hz]”の周波数ずれ“Δf”に相当するものとして計算する。そして、単位時間分の位相差“Δθ”の時間変化量“dΔθ/dt”から、周波数ずれ“Δf”を算出する。この周波数ずれ“Δf”の算出を単位時間毎に行う。そして、所定の推定時間が経過する毎に、当該推定時間の間に算出された周波数ずれ“Δf”を用いて水晶発振周波数“f”を推定する。
例えば、単位時間が“1ミリ秒”であり、推定時間が“10ミリ秒”であるとすると、推定タイミングが到来するまでの間に、10個の周波数ずれ“Δf”の算出データが得られる。これらの算出データのうち、例えば、周波数ずれ“Δf”の最大値や中央値、平均値といった値を代表値として用いて水晶発振周波数“f”を推定する。つまり、基準発振周波数“fb”に周波数ずれ“Δf”の代表値を加算した値を水晶発振周波数“f”と推定する(f=fb+Δf)。水晶発振周波数“f”を推定したならば、図2の気圧対水晶発振周波数の関係から気圧を推定する。
基準発振周波数“fb”を用いて気圧を推定した場合、気圧の値は大まかにしか求まらない。しかし、上記のように位相差“Δθ”を用いて1Hz以下の周波数ずれ“Δf”を算出し、この周波数ずれ“Δf”を加味して気圧を推定することで、高い分解能で気圧を推定することができる。つまり、本実施形態の手法は、基準発振周波数“fb”を用いて第1の分解能で気圧を推定するとともに、位相差“Δθ”を用いて基準発振器30の1Hz以下に相当する気圧を第1の分解能より高い第2の分解能で推定することに相当する。
本願発明者が行った実験によれば、位相差“Δθ”が1周期分(360°)変化したとした場合の周波数ずれ“Δf”は57パスカル[Pa]程度である。従って、位相差“Δθ”が1°変化したとすれば、周波数ずれ“Δf”は約0.16[Pa]となる。100[m]の高度変化で気圧が1000[Pa]変化するとして換算すると、位相差“Δθ”の1°の変化は、約1.6[cm]の高度変化となる。誤差を見て位相差“Δθ”の変化を30°単位で求めるとしても、約48[cm]の高度変化を検出することができる。このことから、本実施形態の気圧推定方法では、1Hz以下の周波数ずれを用いて高精度な気圧推定を実現可能であることがわかる。
図3及び図4は、本実施形態における推定時間間隔較正の原理の説明図である。本実施形態では、水晶発振周波数の安定度(周波数安定度)を判定し、この周波数安定度に基づいて推定時間間隔を較正する。周波数安定度の判定は、水晶発振周波数のアラン分散を求めることで行う。
アラン分散は、ある発振器が、安定した周波数の信号をどの程度の期間に亘って発振することができるかを示す指標値である。時間領域における周波数安定度の尺度であると言える。アラン分散は、周波数偏差を所定の平均化時間に亘って平均化し、サンプル数を2個として分散を計算する2標本分散として定義される。
今、推定時間間隔を“τ”と表記する。このとき、上記の周波数偏差“y(t)”を、ある推定時間間隔“τn”の区間毎に平均化することで周波数偏差平均値“yk(τn)”を算出する。推定時間間隔“τ”は、周波数偏差“y(t)”を平均化する平均化時間に相当する。
具体的には、周波数偏差平均値“yk(τn)”を、次式(3)に従って算出する。
但し、添え字の“k”は周波数偏差平均値の番号を示し、“yk(τn)”は、推定時間間隔“τn”の区間毎で平均した周波数偏差平均値のうちのk番目の値であることを示す。
このとき、推定時間間隔“τn”での水晶発振周波数“f”のアラン分散“σy(τn)”を、次式(4)に従って算出する。
なお、“σy(τn)”は正確にはアラン標準偏差であるが、本明細書ではアラン標準偏差はアラン分散と同義であるものとして説明する。
上記のアラン分散“σy(τn)”を、推定時間間隔“τn”を変化させながら算出する。例えば、図3に示すように、N種類の推定時間間隔“τn={τ1,τ2,τ3,・・・, τN}”を設定する。図3において、横軸は時間軸であり、最上段のラインの下向きの矢印が周波数ずれ“Δf”の算出タイミングを示している。また、2段目以降の横方向の帯は、各推定時間間隔“τn”での周波数偏差平均値“yk(τn)”のデータを示しており、1つの矩形が1つのデータに相当する。推定時間間隔“τn”は、例えば10ミリ秒〜100秒までの時間範囲に含まれる離散的な値として設定する。
図4は、推定時間間隔とアラン分散との対応関係を示すグラフの一例である。横軸は推定時間間隔“τn”を示し、縦軸はアラン分散“σy(τn)”を示す。このグラフを見ると、アラン分散“σy(τn)”は、推定時間間隔“τn”が増加するにつれて徐々に減少することがわかる。そして、ある推定時間間隔において最小となった後、再び増加する傾向があることがわかる。
アラン分散“σy(τn)”は、推定時間間隔“τn”に相当する期間において、水晶発振周波数がどの程度の誤差の広がりを持つかを示す値である。そのため、アラン分散“σy(τn)”が小さいほど、水晶発振周波数“f”は安定であると言える。
例えば、図4のグラフでは、推定時間間隔“τ5”でのアラン分散“σy(τ5)”が最小となっている。つまり、水晶発振周波数が最も安定するのは推定時間間隔“τ5”の場合であり、推定時間間隔“τ5”に相当する期間であれば、水晶発振周波数“f”の誤差が限りなく小さくなると判断できる。そこで、本実施形態では、アラン分散“σy(τn)”が最小となる推定時間間隔を推定時間間隔の適正値として設定する。これにより、気圧推定のタイミングを適正化することができる。
3.データ構成
図1に示すように、記憶部80には、処理部10により読み出され、気圧推定処理(図5参照)として実行される気圧推定プログラム81が記憶されている。また、気圧推定プログラム81には、推定時間間隔較正処理(図6参照)として実行される推定時間間隔較正プログラム811がサブルーチンとして含まれる。これらの処理については、フローチャートを用いて詳細に後述する。
図1に示すように、記憶部80には、処理部10により読み出され、気圧推定処理(図5参照)として実行される気圧推定プログラム81が記憶されている。また、気圧推定プログラム81には、推定時間間隔較正処理(図6参照)として実行される推定時間間隔較正プログラム811がサブルーチンとして含まれる。これらの処理については、フローチャートを用いて詳細に後述する。
また、記憶部80には、温度依存オフセット値テーブル82と、基準発振周波数83と、推定時間間隔較正用データ84と、推定時間間隔適正値85と、位相差検出データ86と、周波数ずれ算出データ87と、周波数推定値88と、気圧推定値89とが記憶される。
温度依存オフセット値テーブル82は、水晶発振周波数のオフセット値が温度と対応付けて記憶されたテーブルである。温度依存オフセット値テーブル82は、環境温度に基づいて気圧を温度補償するために用いられる。
基準発振周波数83は、基準発振器30の発振周波数である。電源投入後の初期設定時には、水晶発振周波数の概測値が測定され、基準発振周波数83として初期設定される。その後、ループフィルター処理により位相差が無くなるように基準発振周波数が制御され、基準発振周波数83は随時更新される。
推定時間間隔較正用データ84は、推定時間間隔の較正に用いられるデータである。原理で説明した周波数偏差や、周波数偏差平均値、水晶発振周波数のアラン分散といったデータがこれに含まれる。
推定時間間隔適正値85は、推定時間間隔の適正値のデータである。推定時間間隔較正処理が行われる毎に、推定時間間隔適正値85が設定・更新される。
位相差検出データ86は、ループフィルター処理により随時検出される位相差のデータである。また、周波数ずれ算出データ87は、単位時間間隔で算出される周波数ずれのデータである。
4.処理の流れ
図5は、記憶部80に記憶されている気圧推定プログラム81が処理部10により読み出されて実行されることで、水晶振動子型気圧センサー1において実行される気圧推定処理の流れを示すフローチャートである。
図5は、記憶部80に記憶されている気圧推定プログラム81が処理部10により読み出されて実行されることで、水晶振動子型気圧センサー1において実行される気圧推定処理の流れを示すフローチャートである。
先ず、カウンタ部19は、水晶発振周波数の概測値を測定する(ステップA1)。つまり、水晶発振器20から出力される水晶発振信号の周波数の整数部分を、所定のクロック信号に基づいて測定する。そして、測定した概測値を基準発振周波数83として、基準発振器30を初期設定する(ステップA3)。
次いで、ループフィルター処理部13はループフィルター処理を開始し、検出した位相差を位相差検出データ86として記憶部80に随時記憶させる(ステップA5)。そして、処理部10は、電源投入後の初回推定であるか否かを判定し(ステップA7)、初回推定であると判定した場合は(ステップA7;Yes)、記憶部80に記憶されている推定時間間隔較正プログラム811に従って推定時間間隔較正処理を行う(ステップA9)。
図6は、推定時間間隔較正処理の流れを示すフローチャートである。
先ず、推定時間間隔較正部17は、推定時間間隔の複数の候補値それぞれについて、水晶発振周波数のアラン分散を算出する(ステップB1)。
先ず、推定時間間隔較正部17は、推定時間間隔の複数の候補値それぞれについて、水晶発振周波数のアラン分散を算出する(ステップB1)。
次いで、推定時間間隔較正部17は、ステップB1で算出したアラン分散が最小となる推定時間間隔の候補値を選択する(ステップB3)。そして、推定時間間隔較正部17は、選択した候補値を推定時間間隔適正値85として記憶部80に記憶させた後(ステップB5)、推定時間間隔較正処理を終了する。
図5の気圧推定処理に戻って、ステップA7において初回推定ではないと判定した場合は(ステップA7;No)、気圧推定部15は、推定時間間隔の較正タイミングであるか否かを判定する(ステップA11)。較正タイミングとしては、種々のタイミングを設定可能である。例えば、所定時間が経過したタイミングとしてもよいし、環境温度が所定温度以上変化したタイミングとしてもよい。また、ユーザーにより較正が指示されたタイミングとしてもよい。
較正タイミングであると判定した場合は(ステップA11;Yes)、気圧推定部15は、記憶部80に記憶されている推定時間間隔適正値85をリセットする(ステップA13)。そして、ステップA9へと処理を移行し、推定時間間隔較正部17が再び推定時間間隔較正処理を実行する。
ステップA9において推定時間間隔較正処理を行った後、又は、ステップA11において較正タイミングではないと判定した場合は(ステップA11;No)、気圧推定部15は、位相差検出データ86に記憶されている位相差の単位時間分の時間変化に基づいて周波数ずれを算出し、周波数ずれ算出データ87として記憶部80に記憶させる(ステップA15)。
気圧推定部15は、推定タイミングが到来するまでの間(ステップA17;No)、ステップA15の処理を繰り返し行う。そして、推定タイミングが到来したならば(ステップA17;Yes)、気圧推定部15は、記憶部80に記憶されている基準発振周波数83と、周波数ずれ算出データ87に格納されている周波数ずれとを用いて水晶発振周波数を推定し、周波数推定値88として記憶部80に記憶させる(ステップA19)。
次いで、気圧推定部15は、環境温度に基づいて周波数推定値88を補正する(ステップA21)。具体的には、記憶部80に記憶されている温度依存オフセット値テーブル82を参照し、温度センサー等から取得した環境温度に対応する発振周波数のオフセット値を読み出す。そして、周波数推定値88からオフセット値を減算する。
次いで、気圧推定部15は、周波数推定値88を気圧に換算し、気圧推定値89として記憶部80に記憶させる(ステップA23)。そして、気圧推定部15は、処理を終了するか否かを判定し(ステップA25)、まだ終了しないと判定した場合は(ステップA25;No)、ステップA7に戻る。また、処理を終了すると判定した場合は(ステップA25;Yes)、気圧推定処理を終了する。
5.作用効果
水晶振動子型気圧センサー1において、処理部10は、基準発振器30の発振周波数を気圧測定用の水晶発振器20の発振周波数に合わせるように制御する。そして、基準発振器30の発振信号と水晶発振器20の発振信号との位相差を検出し、基準発振器30の発振周波数と位相差とを用いて気圧を推定する。
水晶振動子型気圧センサー1において、処理部10は、基準発振器30の発振周波数を気圧測定用の水晶発振器20の発振周波数に合わせるように制御する。そして、基準発振器30の発振信号と水晶発振器20の発振信号との位相差を検出し、基準発振器30の発振周波数と位相差とを用いて気圧を推定する。
具体的には、基準発振信号を用いて水晶発振信号をIQ分離し、コスタスループを用いて基準発振信号と水晶発振信号との位相差を検出する。そして、単位時間分の位相差の時間変化に基づいて周波数ずれを算出し、基準発振周波数に周波数ずれを加算することで、水晶発振周波数を推定する。かかる構成により、水晶発振器20の発振信号から直接正確に発振周波数を測定せずとも、水晶発振信号と基準発振信号との位相差に基づいて、水晶発振器20の発振周波数を推定することができる。
また、推定時間間隔較正部17は、初回推定時及び所定の較正タイミングにおいて、水晶発振周波数及び気圧の推定を行う時間間隔を較正する。具体的には、水晶発振周波数のアラン分散を求めることで水晶発振器20の周波数安定度を判定し、周波数安定度が最も高くなる時間間隔を推定時間間隔の適正値として設定する。これにより、周波数安定度に基づく適切なタイミングで気圧を推定することが可能となる。
6.変形例
6−1.適用例
上記の実施形態の水晶振動子型気圧センサー1は、例えば高度計に搭載して利用することができる。具体的には、水晶振動子型気圧センサー1を搭載した高度計において、処理部は、水晶振動子型気圧センサー1から出力される気圧推定値を高度に換算・推定する。そして、推定した高度を表示部に表示させる。
6−1.適用例
上記の実施形態の水晶振動子型気圧センサー1は、例えば高度計に搭載して利用することができる。具体的には、水晶振動子型気圧センサー1を搭載した高度計において、処理部は、水晶振動子型気圧センサー1から出力される気圧推定値を高度に換算・推定する。そして、推定した高度を表示部に表示させる。
6−2.位相差の検出
上記の実施形態では、水晶発振信号と基準発振信号との位相差の検出を、処理部がデジタル信号処理としてソフトウェア的に行うものとして説明した。しかし、位相比較器、ループフィルター及びVCO(Voltage Controlled Oscillator)を有して構成されるPLL回路によりコスタスループを形成し、位相差の検出をハードウェア的に行うこととしてもよいのは勿論である。
上記の実施形態では、水晶発振信号と基準発振信号との位相差の検出を、処理部がデジタル信号処理としてソフトウェア的に行うものとして説明した。しかし、位相比較器、ループフィルター及びVCO(Voltage Controlled Oscillator)を有して構成されるPLL回路によりコスタスループを形成し、位相差の検出をハードウェア的に行うこととしてもよいのは勿論である。
6−3.気圧の推定
水晶発振周波数の推定を次のように行ってもよい。すなわち、基準発振器の発振周波数(基準発振周波数)を気圧に換算し、気圧基準値とする。また、位相差に基づいて算出した周波数ずれを用いて、気圧基準値からの気圧変化分を推定する。そして、気圧基準値に気圧変化分を加算することで気圧を推定する。
水晶発振周波数の推定を次のように行ってもよい。すなわち、基準発振器の発振周波数(基準発振周波数)を気圧に換算し、気圧基準値とする。また、位相差に基づいて算出した周波数ずれを用いて、気圧基準値からの気圧変化分を推定する。そして、気圧基準値に気圧変化分を加算することで気圧を推定する。
6−4.周波数安定度の判定
上記の実施形態では、アラン分散を用いて水晶発振周波数の安定度を判定したが、周波数安定度の判定方法はこれに限られない。アラン分散は2標本分散であるが、例えば、周波数偏差のサンプル数を2個よりも多くして分散値を計算し、その分散値に基づいて周波数安定度を判定することとしてもよい。推定時間間隔の較正に適用可能な周波数安定度の判定方法であれば、任意の手法を適用可能である。
上記の実施形態では、アラン分散を用いて水晶発振周波数の安定度を判定したが、周波数安定度の判定方法はこれに限られない。アラン分散は2標本分散であるが、例えば、周波数偏差のサンプル数を2個よりも多くして分散値を計算し、その分散値に基づいて周波数安定度を判定することとしてもよい。推定時間間隔の較正に適用可能な周波数安定度の判定方法であれば、任意の手法を適用可能である。
6−5.他の適用例
上記の実施形態では、第1の乗算器40→位相比較部11→ループフィルター処理部13→基準発振器30→第1の乗算器40の第1のループと、第2の乗算器60→位相比較部11→ループフィルター処理部13→基準発振器30→遅延回路50→第2の乗算器60の第2のループとで形成されるコスタスループを、気圧推定装置の一種である水晶振動子型気圧センサー1に適用する場合の実施形態について説明した。
上記の実施形態では、第1の乗算器40→位相比較部11→ループフィルター処理部13→基準発振器30→第1の乗算器40の第1のループと、第2の乗算器60→位相比較部11→ループフィルター処理部13→基準発振器30→遅延回路50→第2の乗算器60の第2のループとで形成されるコスタスループを、気圧推定装置の一種である水晶振動子型気圧センサー1に適用する場合の実施形態について説明した。
しかし、上記のコスタスループを適用可能な装置は何も気圧推定装置に限られるわけでないことは勿論である。例えば、衛星測位システムを利用して位置算出を行う位置算出装置に適用することも可能である。そこで、衛星測位システムの一種であるGPS(Global Positioning System)を利用して位置算出を行うGPS位置算出装置にコスタスループを適用する場合の適用例について説明する。なお、図1で説明した水晶振動子型気圧センサー1と同一の構成要素については同一の符号を付して、再度の説明を省略する。
図7は、GPS位置算出装置3の機能構成の一例を示す図である。GPS位置算出装置3は、RF(Radio Frequency)受信回路部210と、ベースバンド処理回路部220とを備えて構成される。なお、RF受信回路部210とベースバンド処理回路部220とは、それぞれ別のLSI(Large Scale Integration)として製造することも、1チップとして製造することも可能である。
RF受信回路部210は、不図示のGPSアンテナで受信されたRF信号を処理する受信回路である。RF受信回路部210の構成としては、受信したRF信号をA/D変換器でデジタル信号に変換し、デジタル信号を処理する受信回路を構成してもよい。また、GPSアンテナで受信されたRF信号をアナログ信号のまま信号処理し、最終的にA/D変換することで、デジタル信号をベースバンド処理回路部220に出力する構成としてもよい。
ベースバンド処理回路部220は、RF受信回路部210から出力された受信信号に基づいて、GPS衛星から送出されているGPS衛星信号を捕捉する回路部である。GPS衛星信号は、拡散符号の一種であるC/A(Coarse and Acquisition)コードによって、スペクトラム拡散方式として知られるCDMA(Code Division Multiple Access)方式によって変調された1.57542[GHz]の信号である。C/Aコードは、コード長1023チップを1PNフレームとする繰返し周期1msの擬似ランダム雑音符号であり、各GPS衛星に固有のコードである。
ベースバンド処理回路部220は、受信信号に対してキャリア(搬送波)の除去や相関演算を、専用の回路によってハードウェア的に、或いは、デジタル信号処理としてソフトウェア的に行うことで、GPS衛星信号を捕捉する。そして、捕捉したGPS衛星信号から抽出した衛星軌道情報や時刻情報等を利用して、GPS位置算出装置3の位置(位置座標)や時計誤差(クロックバイアス)を算出する。
ベースバンド処理回路部220は、例えば、水晶発振器20と、基準発振器30と、第1の乗算器40と、遅延回路50と、第2の乗算器60と、処理部100と、記憶部800とを備えて構成される。
処理部100は、ベースバンド処理回路部220の各機能部を統括的に制御する制御装置及び演算装置であり、CPU等のプロセッサーを有して構成される。処理部100は、例えば、気圧推定部15と、推定時間間隔較正部17と、カウンタ部19と、衛星捕捉部110と、位置算出部120とを機能部として有する。
衛星捕捉部110は、RF受信回路部210から出力されるデジタル化された受信信号に対して、キャリア除去や相関演算等のデジタル信号処理を実行する。そして、そのデジタル信号処理の結果に基づいて、捕捉対象とするGPS衛星に係るメジャメント情報830(コード位相、ドップラー周波数、擬似距離、擬似距離変化率等)を演算する。
位置算出部120は、衛星捕捉部110によって演算されたメジャメント情報830と、気圧推定部15によって推定された気圧推定値89とを用いて、所定の位置算出計算を行ってGPS位置算出装置3の位置及び時計誤差を算出する。
具体的には、位置算出部120は、例えば、メジャメント情報830を利用した3次元の位置算出計算を行って、緯度、経度及び高度で表される3次元の位置を算出する。そして、気圧推定値89を用いて高度方向の位置成分を補正することで、最終的な位置を求める。或いは、メジャメント情報830を利用した2次元の位置算出計算を行って、緯度及び経度で表される2次元の位置を算出する。そして、気圧推定値89から求まる高度を含めた3次元の位置をGPS位置算出装置3の位置として算出することとしてもよい。
記憶部800には、プログラムとして、例えば、位置算出部120によって位置算出処理として実行される位置算出プログラム820が記憶される。位置算出プログラム820は、気圧推定処理(図5参照)として実行される気圧推定プログラム81をサブルーチンとして含む。また、気圧推定プログラム81は、推定時間間隔較正処理(図6参照)として実行される推定時間間隔較正プログラム811をサブルーチンとして含む。
また、記憶部800には、データとして、例えば、温度依存オフセット値テーブル82と、基準発振周波数83と、推定時間間隔較正用データ84と、推定時間間隔適正値85と、位相差検出データ86と、周波数ずれ算出データ87と、周波数推定値88と、気圧推定値89と、メジャメント情報830と、算出位置データ840とが記憶される。
図7のGPS位置算出装置3は、携帯型電話機やカーナビゲーション装置、携帯型ナビゲーション装置、パソコン、PDA(Personal Digital Assistant)、デジタルカメラ、腕時計といった各種の電子機器に搭載して利用することが可能である。
また、この場合に適用可能な衛星測位システムはGPSに限らず、WAAS(Wide Area Augmentation System)やQZSS(Quasi Zenith Satellite System)、GLONASS(GLObal NAvigation Satellite System)、GALILEO等の衛星測位システムであってもよい。
1 水晶振動子型気圧センサー、 3 GPS位置算出装置、 10 処理部、 11 位相比較部、 13 ループフィルター処理部、 15 気圧推定部、 17 推定時間間隔較正部、 19 カウンタ部、 20 水晶発振器、 30 基準発振器、 40 第1の乗算器、 50 遅延回路、 60 第2の乗算器、 80 記憶部
Claims (8)
- 所定の発振器の発振周波数を気圧測定用発振器の発振周波数に合わせるように制御することと、
前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出することと、
前記発振器の発振周波数と前記位相差とを用いて気圧を推定することと、
を含む気圧推定方法。 - 前記発振器の発振周波数を制御することは、前記位相差に基づいて制御することを含む、
請求項1に記載の気圧推定方法。 - 前記位相差を検出することは、コスタスループを用いて前記位相差を検出することを含む、
請求項1又は2に記載の気圧推定方法。 - 前記気圧を推定することは、
前記発振器の発振周波数を用いて第1の分解能で気圧を推定することと、
前記位相差を用いて前記発振器の1Hz以下に相当する気圧を前記第1の分解能より高い第2の分解能で推定することと、
を含む、
請求項1〜3の何れか一項に記載の気圧推定方法。 - 前記気圧を推定することは、所定時間の間に検出された前記位相差を用いて前記気圧を推定することを含み、
前記気圧測定用発振器の発振周波数の安定度を判定することと、
前記安定度に基づいて前記所定時間を設定することと、
を更に含む、
請求項1〜4の何れか一項に記載の気圧推定方法。 - 前記安定度を判定することは、前記気圧測定用発振器の発振周波数のアラン分散を求めることを含む、
請求項5に記載の気圧推定方法。 - 前記気圧を推定することは、環境温度に基づいて前記気圧を温度補償することを含む、
請求項1〜6の何れか一項に記載の気圧推定方法。 - 気圧に応じて発振周波数が変化する気圧測定用発振器と、
発振周波数を変更可能な発振器と、
前記発振器の発振周波数を前記気圧測定用発振器の発振周波数に合わせるように制御する周波数制御部と、
前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出する位相差検出部と、
前記発振器の発振周波数と前記位相差とを用いて気圧を推定する気圧推定部と、
を備えた気圧推定装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011194655A JP2012141282A (ja) | 2010-12-17 | 2011-09-07 | 気圧推定方法及び気圧推定装置 |
US13/329,159 US20120152018A1 (en) | 2010-12-17 | 2011-12-16 | Atmospheric pressure estimation method and atmospheric pressure estimation device |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010281339 | 2010-12-17 | ||
JP2010281339 | 2010-12-17 | ||
JP2011194655A JP2012141282A (ja) | 2010-12-17 | 2011-09-07 | 気圧推定方法及び気圧推定装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2012141282A true JP2012141282A (ja) | 2012-07-26 |
Family
ID=46232610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011194655A Withdrawn JP2012141282A (ja) | 2010-12-17 | 2011-09-07 | 気圧推定方法及び気圧推定装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20120152018A1 (ja) |
JP (1) | JP2012141282A (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9435767B2 (en) * | 2012-09-04 | 2016-09-06 | Syracuse University | Sensor for sensing substances in an environment |
JP2016519505A (ja) * | 2013-04-04 | 2016-06-30 | テレフオンアクチーボラゲット エルエム エリクソン(パブル) | 無線ノードにより実行される測定に環境条件が与える影響の考慮 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4638644A (en) * | 1978-04-07 | 1987-01-27 | Edward Gidseg | Circulating air refrigerator with removable divider shelf |
JPH01313728A (ja) * | 1988-06-13 | 1989-12-19 | Seiko Electronic Components Ltd | 水晶式真空計 |
SE510713C2 (sv) * | 1997-10-10 | 1999-06-14 | Ericsson Telefon Ab L M | Faslåsningskrets samt metod för reglering av spänningsstyrd oscillator |
US8502455B2 (en) * | 2009-05-29 | 2013-08-06 | Agilent Technologies, Inc. | Atmospheric inductively coupled plasma generator |
-
2011
- 2011-09-07 JP JP2011194655A patent/JP2012141282A/ja not_active Withdrawn
- 2011-12-16 US US13/329,159 patent/US20120152018A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
US20120152018A1 (en) | 2012-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5034935B2 (ja) | 測位方法、プログラム、測位装置及び電子機器 | |
JP5186873B2 (ja) | 測位方法、プログラム、測位装置及び電子機器 | |
US20090262018A1 (en) | High Accuracy Satellite Receiving Controller and Associated Method | |
JP5423036B2 (ja) | 位置算出方法及び位置算出装置 | |
WO2016050055A1 (zh) | 一种信号捕获方法、装置及计算机存储介质 | |
US20120319898A1 (en) | Method of detecting multipath, gnss receiving apparatus , and mobile terminal | |
JP2011128071A (ja) | 位置算出方法及び位置算出システム | |
US8358241B2 (en) | Position calculating method and position calculating device | |
JP2006189320A (ja) | 測位演算機、測位装置および測位演算方法 | |
JP5526492B2 (ja) | 擬似距離算出方法、測位方法、プログラム及び測位装置 | |
US8779973B2 (en) | Satellite signal tracking method, position calculating method, and position calculating device | |
EP2006706B1 (en) | Coherent integration enhancement method, positioning method, storage medium, coherent integration enhancement circuit, positioning circuit, and electronic instrument | |
US20130141279A1 (en) | Positioning satellite signal receiver, positioning satellite signal receiving method, and computer readable storage medium | |
JP5228433B2 (ja) | 測位方法、プログラム、測位装置及び電子機器 | |
JP4848146B2 (ja) | 測位信号を送信するための装置、その装置を備える測位システムおよび測位信号を送信するシステム | |
US8482459B2 (en) | Positioning device, positioning method and storage medium storing program | |
JP2012141282A (ja) | 気圧推定方法及び気圧推定装置 | |
US20200120632A1 (en) | Mobile terminal, information processor, information processing method, and program | |
JP6047944B2 (ja) | 受信装置及び相関積算処理方法 | |
JP5375773B2 (ja) | 位置算出方法及び位置算出装置 | |
JP2014228540A (ja) | 擬似距離の誤差を推定するためのシステムおよび方法 | |
JP2011164088A (ja) | 衛星信号追尾方法、位置算出方法、衛星信号追尾装置及び位置算出装置 | |
JP5267478B2 (ja) | 衛星信号追尾方法及び位置算出装置 | |
US11747487B2 (en) | GNSS receiver clock frequency drift detection | |
JP2013253814A (ja) | 測位方法、測位プログラム、測位装置、および情報機器端末 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20141202 |