"所定の発振器の発振周波数を気圧測定用発振器の発振周波数に合わせるように制御することと、前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出することと、前記発振器の発振周波数と前記位相差とを用いて気圧を推定することと、を含む気圧推定方法。前記発振器の発振周波数を制御することは、前記位相差に基づいて制御することを含む、請求項1に記載の気圧推定方法。前記位相差を検出することは、コスタスループを用いて前記位相差を検出することを含む、請求項1又は2に記載の気圧推定方法。前記気圧を推定することは、前記発振器の発振周波数を用いて第1の分解能で気圧を推定することと、前記位相差を用いて前記発振器の1Hz以下に相当する気圧を前記第1の分解能より高い第2の分解能で推定することと、を含む、請求項1~3の何れか一項に記載の気圧推定方法。前記気圧を推定することは、所定時間の間に検出された前記位相差を用いて前記気圧を推定することを含み、前記気圧測定用発振器の発振周波数の安定度を判定することと、前記安定度に基づいて前記所定時間を設定することと、を更に含む、請求項1~4の何れか一項に記載の気圧推定方法。前記安定度を判定することは、前記気圧測定用発振器の発振周波数のアラン分散を求めることを含む、請求項5に記載の気圧推定方法。前記気圧を推定することは、環境温度に基づいて前記気圧を温度補償することを含む、請求項1~6の何れか一項に記載の気圧推定方法。気圧に応じて発振周波数が変化する気圧測定用発振器と、発振周波数を変更可能な発振器と、前記発振器の発振周波数を前記気圧測定用発振器の発振周波数に合わせるように制御する周波数制御部と、前記発振器の発振信号と前記気圧測定用発振器の発振信号との位相差を検出する位相差検出部と、前記発振器の発振周波数と前記位相差とを用いて気圧を推定する気圧推定部と、を備えた気圧推定装置。圧力によって発振周波数が変化する発振器を用いて気圧を推定する場合には、その発振器の発振周波数を正しく測定する必要がある。しかし、発振周波数を正しく測定するためには、高精度のクロックや内部タイムベース等、周波数を高精度に測定するための機構が必要になるという問題があった。本発明は上述した課題に鑑みて為されたものであり、その目的とするところは、気圧測定用発振器を用いた新たな気圧推定手法を提案することにある。以上の課題を解決するための第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とを備えて構成される。処理部10は、水晶振動子型気圧センサー1の各部を統括的に制御する制御装置及び演算装置であり、CPUCentralProcessingUnit等のプロセッサーを有して構成される。処理部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は、位相同期回路として知られるPLLPhaseLockedLoop回路における位相比較器に相当する機能ブロックである。ループフィルター処理部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は、発振周波数を変更可能に構成された可変周波数発振器であり、例えばNCONumericalControlledOscillatorにより構成される。第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は、ROMReadOnlyMemoryやフラッシュROM、RAMRandomAccessMemory等のメモリーを有して構成される。2.原理図2は、本実施形態における気圧推定の原理の説明図である。図2において、横軸は気圧Pを示し、縦軸は水晶発振周波数fを示す。また、温度T1~T4T1<T2<T3<T4それぞれにおける気圧対水晶発振周波数を、黒四角形、黒三角形、バツ印及び黒丸でそれぞれプロットしている。このグラフを見ると、温度に関わらず、気圧Pと水晶発振周波数fとは、ほぼ線形の関係にあることがわかる。また、温度が高くなるほど、水晶発振周波数fの総体的な大きさは小さくなる傾向がある。かかる気圧対水晶発振周波数の関係から、環境温度と水晶発振周波数の推定値とを用いて気圧を推定することができる。処理部10は、コスタスループを用いて検出した位相差Δθの単位時間分の時間変化に基づいて、動作条件の変化等に起因する水晶発振周波数の基準発振周波数からのずれを算出する。この水晶発振周波数の基準発振周波数からのずれΔfを周波数ずれと定義する。また、周波数ずれΔfを基準発振周波数fbで除算した値Δf/fbのことを周波数偏差と定義する。位相差Δθと周波数ずれΔfとの間には、次式1の関係が成立する。式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標本分散として定義される。最初に、周波数偏差ytを、次式2で定義する。但し、周波数ずれΔftは、単位時間間隔で算出された時系列データであり、時間の関数で表記している。今、推定時間間隔をτと表記する。このとき、上記の周波数偏差ytを、ある推定時間間隔τnの区間毎に平均化することで周波数偏差平均値ykτnを算出する。推定時間間隔τは、周波数偏差ytを平均化する平均化時間に相当する。具体的には、周波数偏差平均値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がサブルーチンとして含まれる。これらの処理については、フローチャートを用いて詳細に後述する。また、記憶部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において実行される気圧推定処理の流れを示すフローチャートである。先ず、カウンタ部19は、水晶発振周波数の概測値を測定するステップA1。つまり、水晶発振器20から出力される水晶発振信号の周波数の整数部分を、所定のクロック信号に基づいて測定する。そして、測定した概測値を基準発振周波数83として、基準発振器30を初期設定するステップA3。次いで、ループフィルター処理部13はループフィルター処理を開始し、検出した位相差を位相差検出データ86として記憶部80に随時記憶させるステップA5。そして、処理部10は、電源投入後の初回推定であるか否かを判定しステップA7、初回推定であると判定した場合はステップA7;Yes、記憶部80に記憶されている推定時間間隔較正プログラム811に従って推定時間間隔較正処理を行うステップA9。図6は、推定時間間隔較正処理の流れを示すフローチャートである。先ず、推定時間間隔較正部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の発振周波数と位相差とを用いて気圧を推定する。具体的には、基準発振信号を用いて水晶発振信号をIQ分離し、コスタスループを用いて基準発振信号と水晶発振信号との位相差を検出する。そして、単位時間分の位相差の時間変化に基づいて周波数ずれを算出し、基準発振周波数に周波数ずれを加算することで、水晶発振周波数を推定する。かかる構成により、水晶発振器20の発振信号から直接正確に発振周波数を測定せずとも、水晶発振信号と基準発振信号との位相差に基づいて、水晶発振器20の発振周波数を推定することができる。また、推定時間間隔較正部17は、初回推定時及び所定の較正タイミングにおいて、水晶発振周波数及び気圧の推定を行う時間間隔を較正する。具体的には、水晶発振周波数のアラン分散を求めることで水晶発振器20の周波数安定度を判定し、周波数安定度が最も高くなる時間間隔を推定時間間隔の適正値として設定する。これにより、周波数安定度に基づく適切なタイミングで気圧を推定することが可能となる。6.変形例6-1.適用例上記の実施形態の水晶振動子型気圧センサー1は、例えば高度計に搭載して利用することができる。具体的には、水晶振動子型気圧センサー1を搭載した高度計において、処理部は、水晶振動子型気圧センサー1から出力される気圧推定値を高度に換算推定する。そして、推定した高度を表示部に表示させる。6-2.位相差の検出上記の実施形態では、水晶発振信号と基準発振信号との位相差の検出を、処理部がデジタル信号処理としてソフトウェア的に行うものとして説明した。しかし、位相比較器、ループフィルター及びVCOVoltageControlledOscillatorを有して構成されるPLL回路によりコスタスループを形成し、位相差の検出をハードウェア的に行うこととしてもよいのは勿論である。6-3.気圧の推定水晶発振周波数の推定を次のように行ってもよい。すなわち、基準発振器の発振周波数基準発振周波数を気圧に換算し、気圧基準値とする。また、位相差に基づいて算出した周波数ずれを用いて、気圧基準値からの気圧変化分を推定する。そして、気圧基準値に気圧変化分を加算することで気圧を推定する。6-4.周波数安定度の判定上記の実施形態では、アラン分散を用いて水晶発振周波数の安定度を判定したが、周波数安定度の判定方法はこれに限られない。アラン分散は2標本分散であるが、例えば、周波数偏差のサンプル数を2個よりも多くして分散値を計算し、その分散値に基づいて周波数安定度を判定することとしてもよい。推定時間間隔の較正に適用可能な周波数安定度の判定方法であれば、任意の手法を適用可能である。6-5.他の適用例上記の実施形態では、第1の乗算器40→位相比較部11→ループフィルター処理部13→基準発振器30→第1の乗算器40の第1のループと、第2の乗算器60→位相比較部11→ループフィルター処理部13→基準発振器30→遅延回路50→第2の乗算器60の第2のループとで形成されるコスタスループを、気圧推定装置の一種である水晶振動子型気圧センサー1に適用する場合の実施形態について説明した。しかし、上記のコスタスループを適用可能な装置は何も気圧推定装置に限られるわけでないことは勿論である。例えば、衛星測位システムを利用して位置算出を行う位置算出装置に適用することも可能である。そこで、衛星測位システムの一種であるGPSGlobalPositioningSystemを利用して位置算出を行うGPS位置算出装置にコスタスループを適用する場合の適用例について説明する。なお、図1で説明した水晶振動子型気圧センサー1と同一の構成要素については同一の符号を付して、再度の説明を省略する。図7は、GPS位置算出装置3の機能構成の一例を示す図である。GPS位置算出装置3は、RFRadioFrequency受信回路部210と、ベースバンド処理回路部220とを備えて構成される。なお、RF受信回路部210とベースバンド処理回路部220とは、それぞれ別のLSILargeScaleIntegrationとして製造することも、1チップとして製造することも可能である。RF受信回路部210は、不図示のGPSアンテナで受信されたRF信号を処理する受信回路である。RF受信回路部210の構成としては、受信したRF信号をA/D変換器でデジタル信号に変換し、デジタル信号を処理する受信回路を構成してもよい。また、GPSアンテナで受信されたRF信号をアナログ信号のまま信号処理し、最終的にA/D変換することで、デジタル信号をベースバンド処理回路部220に出力する構成としてもよい。ベースバンド処理回路部220は、RF受信回路部210から出力された受信信号に基づいて、GPS衛星から送出されているGPS衛星信号を捕捉する回路部である。GPS衛星信号は、拡散符号の一種であるC/ACoarseandAcquisitionコードによって、スペクトラム拡散方式として知られるCDMACodeDivisionMultipleAccess方式によって変調された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は、携帯型電話機やカーナビゲーション装置、携帯型ナビゲーション装置、パソコン、PDAPersonalDigitalAssistant、デジタルカメラ、腕時計といった各種の電子機器に搭載して利用することが可能である。また、この場合に適用可能な衛星測位システムはGPSに限らず、WAASWideAreaAugmentationSystemやQZSSQuasiZenithSatelliteSystem、GLONASSGLObalNAvigationSatelliteSystem、GALILEO等の衛星測位システムであってもよい。"