(57)【特許請求の範囲】
【請求項１】
アナログの正弦波信号Ｖｉｎの周波数情報を検出する周波数情報検出装置であって、
アナログの正弦波信号Ｖｉｎを、予め定めたサンプリング周期（サンプリング周波数ｆｓの逆数：１／ｆｓ）でデジタル信号に変換することによって得られる時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を順次出力することにより、複数のサンプリングデータで構成される正弦波信号に変換するＡ／Ｄ変換部と、
時刻［ｔ］における入力信号Ｖｉｎの位相情報の推定値を第１推定位相情報として出力する第１推定位相情報出力部と、
時刻［ｔ］における入力信号Ｖｉｎの周波数情報の推定値を第１推定周波数情報として出力する第１推定周波数情報出力部と、
入力信号Ｖｉｎの設定周波数情報に補正周波数情報を加算した周波数情報を同期周波数情報としたときに、前記同期周波数情報に基づいて正弦波信号を生成する正弦波信号生成部と、
時刻［ｔ］における前記正弦波信号生成部から出力される正弦波信号の位相情報の推定値を第２推定位相情報として出力する第２推定位相情報出力部と、
時刻［ｔ］における前記正弦波信号生成部から出力される正弦波信号の周波数情報の推定値を第２推定周波数情報として出力する第２推定周波数情報出力部と、
時刻［ｔ］における前記第１推定周波数情報に対する前記第２推定周波数情報の誤差（誤差周波数情報）の移動平均値を演算し、誤差周波数移動平均値情報として出力する誤差周波数情報出力部と、
時刻［ｔ］における前記第１推定位相情報に対する前記第２推定位相情報の誤差（誤差位相情報）の移動平均値を演算し、誤差位相移動平均値情報として出力する誤差位相情報出力部と、
前記誤差周波数移動平均値情報が予め定めた範囲外の場合は、前記誤差周波数移動平均値情報が小さくなるように前記補正周波数情報を定めて出力し、前記誤差周波数移動平均値情報が予め定めた範囲内の場合は、前記誤差位相移動平均値情報が小さくなるように前記補正周波数情報を定めて出力する周波数補正部と、
を備えた周波数情報検出装置。
【請求項２】
前記第１推定位相情報出力部は、
時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を微分した余弦値（ｃｏｓ（α［ｔ］））の推定値を、時刻［ｔ−１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］及び既知の値を用いて演算する余弦値推定部と、
時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を正弦要素（ｓｉｎ要素）とし、前記余弦値推定部で推定された時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を余弦要素（ｃｏｓ要素）として、逆正接関数（ｔａｎ−１）を用いて時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））の位相α［ｔ］を演算し、演算した位相α［ｔ］を第１推定位相情報として出力する位相推定部と、を有し、
前記第１推定周波数情報出力部は、
前記位相推定部で演算された時刻［ｔ−１］における位相α［ｔ−１］及び時刻［ｔ］における位相α［ｔ］に基づいて、時刻［ｔ］における位相変位量Δα［ｔ］を演算する位相変位量演算部と、
前記位相変位量演算部で演算された時刻［ｔ］における位相変位量Δα［ｔ］に基づいて、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報を演算し、演算した推定周波数情報を第１推定周波数情報として出力する周波数推定部と、
を有していることを特徴とする請求項１に記載の周波数情報検出装置。
【請求項３】
前記余弦値推定部は、前記入力信号Ｖｉｎの設定周波数情報を用いて演算される定数をＫとしたときに、余弦値（ｃｏｓ（α［ｔ］））を式（１）によって演算することを特徴とする請求項２に記載の周波数情報検出装置。
ｃｏｓ（α［ｔ］）＝Ｋ・｛Ｖｉｎ［ｔ＋１］−Ｖｉｎ［ｔ−１］｝  ・・・（１）
【請求項４】
前記定数Ｋを式（２）によって演算することを特徴とする請求項３に記載の周波数情報検出装置。ただし、設定周波数情報を「Ｆｓｅｔ」とする。
Ｋ＝１／｛２ｓｉｎ（２π・（Ｆｓｅｔ／ｆｓ）［ｔ］）｝  ・・・（２）
【請求項５】
請求項１〜４のいずれかに記載の正弦波信号生成部から出力された正弦波信号に基づいて進行波電圧を出力する高周波電源装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、アナログの正弦波信号の周波数情報を検出する周波数情報検出装置及びこの周波数情報検出装置を用いた高周波電源装置に関するものである。
【背景技術】
【０００２】
周波数情報検出装置は、例えばプラズマエッチング、プラズマＣＶＤを行うプラズマ処理装置等の負荷に電力を供給する高周波電源等で用いることが可能である。
【０００３】
図６は、高周波電力供給システムの概略構成例である。
第１高周波電源１１０及び第２高周波電源１２０は、共通の発振器１００から出力される高周波信号を進行波電力出力部（図略）で増幅して高周波電力を出力する。第１高周波電源１１０及び第２高周波電源１２０から出力された高周波電力は、負荷となるプラズマ処理装置１３０内の電極（１３１，１３２）の２箇所に供給される。この図６では、電極１３１及び電極１３２のぞれぞれに高周波電力を供給する例を示したが、例えば、図７に示すように、上側の電極１３１の２箇所に高周波電力を供給することも可能である。
【０００４】
なお、第１高周波電源１１０及び第２高周波電源１２０からプラズマ処理装置１３０に向かう高周波電力を進行波電力といい、プラズマ処理装置１３０で反射されて高周波電源側に戻ってくる高周波電力を反射波電力という。
【０００５】
プラズマ処理装置１３０では、第１高周波電源１１０及び第２高周波電源１２０から供給された高周波電力を用いてプラズマ１３３を発生させて、エッチング等を行う。
【０００６】
なお、第１高周波電源１１０とプラズマ処理装置１３０との間には、第１高周波電源１１０とプラズマ処理装置１３０とのインピーダンスを整合させるインピーダンス整合器が用いられることがある。第２高周波電源１２０とプラズマ処理装置１３０との間にも同様にインピーダンス整合器が用いられることがある。
【０００７】
図６及び図７の高周波電力供給システムでは、第１高周波電源１１０及び第２高周波電源１２０の制御が、それぞれ別の制御回路で行われることになる。制御回路が異なると、たとえ共通の発振器１００から高周波信号が与えられても、第１高周波電源１１０及び第２高周波電源１２０の出力周波数が僅かに異なる場合がある。この原因としては、例えば、制御クロックが異なることによる誤差等が考えられる。特に、第１高周波電源１１０と第２高周波電源１２０とが別のメーカのものであれば、若干の違いが生じる可能性がある。
【０００８】
図６及び図７のような高周波電力供給システムにおいて、２つの高周波電源１１０及び１２０の出力周波数が僅かに異なる場合、電極（１３１，１３２）には、周波数が僅かに異なる高周波電力が供給されて混合される。すなわち、電極（１３１，１３２）に印加される２種類の高周波電圧の周波数が僅かに異なっているので、２種類の高周波電圧が重なり合って生じる電圧の波形は変調されたように振幅が周期的に変動する。例えば、２種類の高周波電圧の周波数が、１３．５６ＭＨｚ近傍で１Ｈｚ異なっていると、約１ｓの周期で振幅が変動する。その結果、プラズマ処理装置１３０で発生するプラズマの電位も周期的に変動するため、悪影響を及ぼす。
例えば、一方の高周波電圧波形の山と他方の高周波電圧波形の山が同じ位相になると重なり合って生じる電圧の波形は最大値となるが、一方の高周波電圧波形の山と他方の高周波電圧波形の谷が同じ位相になると重なり合って生じる電圧の波形は最小値となる。この最大値と最小値との差は大きいので、２種類の高周波電圧が重なり合って生じる電圧の波形は大きな振幅の変化を伴って周期的に変動する。
【０００９】
したがって、第１高周波電源１１０及び第２高周波電源１２０の出力周波数を同一にする必要がある。そのためには、一方の高周波電源の出力周波数を検出し、他方の高周波電源の出力周波数を検出した周波数に合わせれば、上記の問題を解決できる。そのため、高周波電力供給システムの分野においても、周波数を検出する周波数情報検出装置が必要である。また、高周波電源の出力周波数をモニタする等の用途が考えられる。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開昭６３−６６４７２
【発明の概要】
【発明が解決しようとする課題】
【００１１】
従来、高周波領域における周波数情報検出装置としては、例えば、特許文献１（特開昭６３−６６４７２）のようなものが提案されていた。しかし、この特許文献１では、周波数を検出することはできるが、位相情報を検出することができない。そのため、従来の周波数情報検出装置を用いても、第１高周波電源１１０から出力される進行波電圧及び第２高周波電源１２０（他の高周波電源）から出力される進行波電圧の周波数及び位相を合わせることができない。
【００１２】
本発明は、２つの高周波電源から出力される進行波電圧の周波数及び位相を合わせることを可能とする周波数情報検出装置、及びこの周波数情報検出装置を備えている高周波電源装置を提供することを目的としている。
【課題を解決するための手段】
【００１３】
第１の発明によって提供される周波数情報検出装置は、
アナログの正弦波信号Ｖｉｎの周波数情報を検出する周波数情報検出装置であって、
アナログの正弦波信号Ｖｉｎを、予め定めたサンプリング周期（サンプリング周波数ｆｓの逆数：１／ｆｓ）でデジタル信号に変換することによって得られる時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を順次出力することにより、複数のサンプリングデータで構成される正弦波信号に変換するＡ／Ｄ変換部と、
時刻［ｔ］における入力信号Ｖｉｎの位相情報の推定値を第１推定位相情報として出力する第１推定位相情報出力部と、
時刻［ｔ］における入力信号Ｖｉｎの周波数情報の推定値を第１推定周波数情報として出力する第１推定周波数情報出力部と、
入力信号Ｖｉｎの設定周波数情報に補正周波数情報を加算した周波数情報を同期周波数情報としたときに、前記同期周波数情報に基づいて正弦波信号を生成する正弦波信号生成部と、
時刻［ｔ］における前記正弦波信号生成部から出力される正弦波信号の位相情報の推定値を第２推定位相情報として出力する第２推定位相情報出力部と、
時刻［ｔ］における前記正弦波信号生成部から出力される正弦波信号の周波数情報の推定値を第２推定周波数情報として出力する第２推定周波数情報出力部と、
時刻［ｔ］における前記第１推定周波数情報に対する前記第２推定周波数情報の誤差（誤差周波数情報）の移動平均値を演算し、誤差周波数移動平均値情報として出力する誤差周波数情報出力部と、
時刻［ｔ］における前記第１推定位相情報に対する前記第２推定位相情報の誤差（誤差位相情報）の移動平均値を演算し、誤差位相移動平均値情報として出力する誤差位相情報出力部と、
前記誤差周波数移動平均値情報が予め定めた範囲外の場合は、前記誤差周波数移動平均値情報が小さくなるように前記補正周波数情報を定めて出力し、前記誤差周波数移動平均値情報が予め定めた範囲内の場合は、前記誤差位相移動平均値情報が小さくなるように前記補正周波数情報を定めて出力する周波数補正部と、
を備えている。
【００１４】
第２の発明によって提供される周波数情報検出装置は、
前記第１推定位相情報出力部が、
時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を微分した余弦値（ｃｏｓ（α［ｔ］））の推定値を、時刻［ｔ−１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］及び既知の値を用いて演算する余弦値推定部と、
時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を正弦要素（ｓｉｎ要素）とし、前記余弦値推定部で推定された時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を余弦要素（ｃｏｓ要素）として、逆正接関数（ｔａｎ−１）を用いて時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））の位相α［ｔ］を演算し、演算した位相α［ｔ］を第１推定位相情報として出力する位相推定部と、を有し、
前記第１推定周波数情報出力部が、
前記位相推定部で演算された時刻［ｔ−１］における位相α［ｔ−１］及び時刻［ｔ］における位相α［ｔ］に基づいて、時刻［ｔ］における位相変位量Δα［ｔ］を演算する位相変位量演算部と、
前記位相変位量演算部で演算された時刻［ｔ］における位相変位量Δα［ｔ］に基づいて、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報を演算し、演算した推定周波数情報を第１推定周波数情報として出力する周波数推定部と、
を有していることを特徴としている。
【００１５】
第３の発明によって提供される周波数情報検出装置は、
前記余弦値推定部は、前記入力信号Ｖｉｎの設定周波数情報を用いて演算される定数をＫとしたときに、余弦値（ｃｏｓ（α［ｔ］））を下式によって演算することを特徴としている。
ｃｏｓ（α［ｔ］）＝Ｋ・｛Ｖｉｎ［ｔ＋１］−Ｖｉｎ［ｔ−１］｝
【００１６】
第４の発明によって提供される周波数情報検出装置は、
前記定数Ｋを下式によって演算することを特徴としている。ただし、設定周波数情報を「Ｆｓｅｔ」とする。
Ｋ＝１／｛２ｓｉｎ（２π・（Ｆｓｅｔ／ｆｓ）［ｔ］）｝
ことを特徴としている。
【００１７】
第５の発明によって提供される高周波電源装置は、
上記第１〜４のいずれかの発明に記載の正弦波信号生成部から出力された正弦波信号に基づいて進行波電圧を出力することを特徴としている。
【発明の効果】
【００１８】
本発明では、周波数情報検出装置に備わっている正弦波信号生成部から出力される正弦波信号の周波数及び位相が、入力信号の周波数及び位相に近づくように収束していくので、正弦波信号生成部から出力される正弦波信号を用いることにより、他の高周波電源の出力周波数及び位相を入力信号Ｖｉｎの周波数及び位相に合わせることが可能となる。
【図面の簡単な説明】
【００１９】
【図１】高周波電圧の検出例を示すブロック図である。
【図２】周波数情報検出装置１の構成例である。
【図３】正弦波信号生成部４０の構成例を示す図である。
【図４】入力信号Ｖｉｎの周波数ｆｉｎが設定周波数ｆｓｅｔからずれていると想定した場合に、移動平均部７１から出力される差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］のシミュレーション結果である。
【図５】入力信号Ｖｉｎの周波数ｆｉｎが設定周波数ｆｓｅｔからずれていると想定した場合に、移動平均部７１から出力される差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］のシミュレーション結果である。
【図６】高周波電力供給システムの概略構成例である。
【図７】高周波電力供給システムの他の概略構成例である。
【発明を実施するための形態】
【００２０】
以下、図面を参照して本発明の実施形態を詳細に説明する。なお、従来と同一又は同様の構成には、同一符号を付している。
【００２１】
図１は、高周波電圧の検出例を示すブロック図である。
第１高周波電源１１０は、発振器１００から出力される高周波信号を増幅し、無線周波数帯域の出力周波数を有する進行波電力ＰＦを出力して負荷となるプラズマ処理装置１３０に供給するための装置である。高周波電源１１０から出力された進行波電力ＰＦは、伝送線路１４０を介してプラズマ処理装置１３０に供給される。なお、一般にこの種の高周波電源では、数百ｋＨｚ以上の無線周波数帯域の周波数（例えば、１３ＭＨｚ，４０ＭＨｚ等の周波数）を有する進行波電力ＰＦを出力している。
【００２２】
負荷となるプラズマ処理装置１３０は、加工部を備え、その加工部の内部に搬入したウエハ、液晶基板等の被加工物を加工（エッチング、ＣＶＤ等）するための装置である。なお、加工部には、電極（図６及び図７参照）が設けられている。この加工部にプラズマ放電用ガスが導入され、電極に高周波電源１１０から出力された進行波電力ＰＦが供給されると、放電が生じてプラズマが発生する。プラズマ処理装置１３０は、このプラズマを利用して被加工物を加工している。
【００２３】
方向性結合器１５０は、第１高周波電源１１０とプラズマ処理装置１３０との間に挿入されて、高周波電源１１０からプラズマ処理装置１３０に向かう進行波電圧ＶＦを検出し、その検出信号を進行波検出信号Ｖｆとして出力する。なお、方向性結合器１５０は、負荷で反射された反射波電圧を検出する機能も有するが、この例では必要ないので省略している。
【００２４】
フィルタ１６０は、ローパスフィルタ又はバンドパスフィルタであり、方向性結合器１５０から出力される進行波検出信号Ｖｆから高調波成分を除去し、基本周波数成分を通過させる。これにより、フィルタ１６０から出力される信号は正弦波信号となり、周波数情報検出装置１へ送られる。本明細書では、フィルタ１６０から出力される信号を周波数情報検出装置１への入力信号Ｖｉｎとする。
なお、入力信号Ｖｉｎは、アナログの電圧信号であるので、出力周波数をｆ、時間をｔ、位相オフセットをθ、角周波数をω（＝２π・ｆ）とすると、入力信号Ｖｉｎは、式（１）のように表すことができる。なお、ここでは、入力信号Ｖｉｎの振幅を「１」としている。
Ｖｉｎ＝ｓｉｎ（２π・ｆ・ｔ＋θ）
＝ｓｉｎ（ω・ｔ＋θ）  ・・・（１）
【００２５】
また、フィルタ１６０の後段には、後述するＡ／Ｄ変換部１０の入力範囲に適するように信号のレベルを変換するレベル変換回路を設けてもよいが、この図１では省略している。
【００２６】
図２は、周波数情報検出装置１の構成例である。周波数情報検出装置１は、図２に示すように、Ａ／Ｄ変換部１０、第１推定位相情報出力部２０、第１推定周波数情報出力部３０、正弦波信号生成部４０、第２推定位相情報出力部５０、第２推定周波数情報出力部６０、誤差周波数情報出力部７０、誤差位相情報出力部８０及び周波数補正部９０を備えている。
【００２７】
＜Ａ／Ｄ変換部１０＞
Ａ／Ｄ変換部１０は、フィルタ１６０から出力されるアナログの入力信号Ｖｉｎを予め定めたサンプリング周期（サンプリング周波数ｆｓの逆数：１／ｆｓ）でデジタル信号に変換する。Ａ／Ｄ変換されたサンプリングデータは、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］として順次出力される。これにより、交流のアナログ信号波形が、複数のサンプリングデータで構成されるデジタル信号波形に変換される。なお、フィルタ１６０から出力される入力信号Ｖｉｎが正弦波信号であれば、Ａ／Ｄ変換部１０から出力されるデジタル信号波形も正弦波信号となる。
なお、Ａ／Ｄ変換部１０には、Ａ／Ｄコンバータが備わっていて、上記のようなＡ／Ｄ変換を行う。また、Ａ／Ｄコンバータの後段にローパスフィルタ又はバンドパスフィルタを設けて、基本周波数成分を通過させるようにしてもよい。
【００２８】
ここで、入力信号Ｖｉｎの周波数をｆｉｎ、サンプリング周波数をｆｓ、位相オフセットをθ’、相対角周波数をωｉｎ＝２π・（ｆｉｎ／ｆｓ）とすると、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］は、式（２）のように表すことができる。なお、ここでは、入力信号Ｖｉｎの振幅を「１」としている。また時間データ「ｔ」はサンプリング周期毎にインクリメントされる変数である。
Ｖｉｎ［ｔ］＝ｓｉｎ（２π・（ｆｉｎ／ｆｓ）［ｔ］＋θ’）
＝ｓｉｎ（ωｉｎ［ｔ］＋θ’）            ・・・・・（２）
【００２９】
＜第１推定位相情報出力部２０＞
第１推定位相情報出力部２０は、時刻［ｔ］における入力信号Ｖｉｎの位相情報の推定値α［ｔ］を第１推定位相情報として出力する。この第１推定位相情報出力部２０は、余弦値推定部２１、位相推定部２２を備えている。
【００３０】
［余弦値推定部２１］
余弦値推定部２１は、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を微分した余弦値「ｃｏｓ（ωｉｎ［ｔ］＋θ’）」の推定値を、時刻［ｔ−１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］及び既知の値を用いて演算する。推定した余弦値「ｃｏｓ（ωｉｎ［ｔ］＋θ’）」は位相推定部２２に送られる。以下、具体的に説明する。
【００３１】
「α［ｔ］＝ωｉｎ［ｔ］＋θ’＝２π・（ｆｉｎ／ｆｓ）［ｔ］＋θ’」とすると、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））は、式（３）で表すことができる。
ｃｏｓ（α［ｔ］）
＝｛（２ｓｉｎ（ωｉｎ［ｔ］）・ｃｏｓ（α［ｔ］）｝／（２ｓｉｎ（ωｉｎ［ｔ］））
＝｛ｓｉｎ（α［ｔ］＋ωｉｎ［ｔ］）−ｓｉｎ（α［ｔ］−ωｉｎ［ｔ］）｝／（２ｓｉｎ（ωｉｎ［ｔ］））
＝｛ｓｉｎ（２ωｉｎ［ｔ］＋θ’）−ｓｉｎ（θ’）｝／（２ｓｉｎ（ωｉｎ［ｔ］））  ・・・（３）
【００３２】
式（２）を参照すると、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］は「ｓｉｎ（ωｉｎ［ｔ］＋θ’）」であるので、入力信号Ｖｉｎの周波数ｆｉｎが変化しなければ、サンプリング周期（１／ｆｓ）毎の相対角周波数ωｉｎの変位量は「ωｉｎ［ｔ］」で一定である。したがって、式（３）の分子は、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］から、時刻［ｔ−１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］を減算することを表している。
【００３３】
また、式（３）の分母の「ｓｉｎ（ωｉｎ［ｔ］）」は、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］から位相オフセットθ’を省略したものとなっている。
【００３４】
したがって、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））は、時刻［ｔ−１］、［ｔ］及び［ｔ＋１］における３つの入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］、Ｖｉｎ［ｔ］及びＶｉｎ［ｔ＋１］を用いて推定することができる。
【００３５】
ただし、上記のように、「ｓｉｎ（ωｉｎ［ｔ］）」は、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］から位相オフセットθ’を省略したものとなっているので、式（３）の分母には誤差が生じる。ここで、正弦関数の値（ｓｉｎの値）は、０を中心として±１の範囲で変化するものであるから、誤差は正になることもあれば負になることもある。
【００３６】
後述するように、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を推定した後、位相α［ｔ］を演算し、さらに位相α［ｔ］のサンプリング周期毎の位相変位量Δα［ｔ］に基づいて入力信号Ｖｉｎの周波数の推定値ｆｅｓｔ［ｔ］又は角周波数の推定値ωｅｓｔ［ｔ］を演算する。その後、演算した周波数の推定値ｆｅｓｔ［ｔ］又は角周波数の推定値ωｅｓｔ［ｔ］の移動平均値を演算するので、式（３）の分母の誤差は、殆ど相殺される。
【００３７】
そこで、もともと誤差のある式（３）の分母を構成する入力信号Ｖｉｎの周波数ｆｉｎを、式（４）のように、設定周波数ｆｓｅｔに置き換えても殆ど影響はない。ここで、設定周波数ｆｓｅｔとは、入力信号Ｖｉｎの設定周波数であるので、誤差も小さい。例えば、図６及び図７の例では、発振器１００から出力される高周波信号の周波数である。この設定周波数ｆｓｅｔは、予め分かっているので、例えば余弦値推定部２１に入力しておけばよい。もちろん、設定角周波数ωｓｅｔ（＝２π・ｆｉｎ）を用いてもよい。
２ｓｉｎ（ωｉｎ［ｔ］）＝２ｓｉｎ（２π・（ｆｉｎ／ｆｓ）［ｔ］）
→  ２ｓｉｎ（２π・（ｆｓｅｔ／ｆｓ）［ｔ］）  ・・・（４）
【００３８】
なお、本明細書では、設定周波数ｆｓｅｔ及び設定角周波数ωｓｅｔを総称して「設定周波数情報Ｆｓｅｔ」という。
【００３９】
ここで、式（４）を構成する各要素は、全て既知の値である。そのため、「１／｛２ｓｉｎ（２π・（ｆｓｅｔ／ｆｓ）［ｔ］）｝」を定数Ｋで表すと、式（３）は式（５）のように変形できる。したがって、本来であれば、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を推定する際には、式（３）を用いる必要があるので複雑な演算が必要であるが、式（５）を用いることによって、演算式を簡略化でき、演算負荷を低減できる。
ｃｏｓ（α［ｔ］）＝Ｋ・｛ｓｉｎ（２ωｉｎ［ｔ］＋θ’）−ｓｉｎ（θ’）｝
＝Ｋ・｛Ｖｉｎ［ｔ＋１］−Ｖｉｎ［ｔ−１］｝  ・・・（５）
【００４０】
上記のように、余弦値推定部２１は、Ａ／Ｄ変換部１０によってデジタル信号となった複数の入力信号Ｖｉｎのデータを用いて時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を演算する。そのため、余弦値推定部２１は図示しないメモリを有し、そのメモリにＡ／Ｄ変換部１０から出力された入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を順次記憶していく。
【００４１】
なお、上記のように、余弦値推定部２１では、時刻［ｔ−１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ−１］と、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］とを用いるが、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を用いない。しかし、後述するように位相推定部２２で行う演算には時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を用いるので、メモリには、少なくとも瞬時値Ｖｉｎ［ｔ−１］、瞬時値Ｖｉｎ［ｔ］及び瞬時値Ｖｉｎ［ｔ＋１］の連続する３つのデータを記憶しておく。
【００４２】
また、上記の例では、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を演算することができる。
【００４３】
［位相推定部２２］
位相推定部２２は、時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］を正弦要素（ｓｉｎ要素）とし、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））を余弦要素（ｃｏｓ要素）として、式（６）に示すように、逆正接関数（ｔａｎ−１）を用いて時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））の位相α［ｔ］を演算する。この位相α［ｔ］は、±π［単位：ｒａｄ］の範囲で演算される。演算された位相α［ｔ］は、位相変位量演算部３１に送られる。
α［ｔ］＝ｔａｎ−１（Ｖｉｎ［ｔ］／ｃｏｓ（α［ｔ］））  ・・・（６）
【００４４】
なお、上記のように、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））が演算されるので、位相α［ｔ］も、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算される。また、位相α［ｔ］の演算に必要な時刻［ｔ］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ］は、余弦値推定部２１に設けたメモリから読み出せばよい。
【００４５】
＜第１推定周波数情報出力部３０＞
第１推定周波数情報出力部３０は、時刻［ｔ］における入力信号Ｖｉｎの周波数情報の推定値を第１推定周波数情報Ｆｅｓｔ［ｔ］として出力する。この第１推定周波数情報出力部３０は、位相変位量演算部３１、周波数推定部３２を備えている。
【００４６】
［位相変位量演算部３１］
位相変位量演算部３１は、位相推定部２２で演算された時刻［ｔ−１］における位相α［ｔ−１］及び時刻［ｔ］における位相α［ｔ］に基づいて、時刻［ｔ］における位相変位量Δα［ｔ］を演算する。演算された位相変位量Δα［ｔ］は、周波数推定部３２に送られる。
【００４７】
ここで、位相変位量演算部３１は、時刻［ｔ−１］から時刻［ｔ］のサンプリング周期の間に生じた位相変位量Δα［ｔ］を演算するのであるが、位相推定部２２から出力される時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））の位相α［ｔ］が、上記のように±π［単位：ｒａｄ］の範囲で演算されるため、｛「α［ｔ］−α［ｔ−１］｝が正の場合は、位相変位量Δα［ｔ］を式（７）を用いて演算し、｛「α［ｔ］−α［ｔ−１］｝が負の場合は、位相変位量Δα［ｔ］を式（８）を用いて演算する。
Δα［ｔ］＝α［ｔ］−α［ｔ−１］            ・・・（７）
Δα［ｔ］＝（α［ｔ］−α［ｔ−１］）＋２π  ・・・（８）
【００４８】
なお、上記のように、位相変位量演算部３１は、位相推定部２２によって演算された時刻［ｔ−１］における位相α［ｔ−１］及び時刻［ｔ］における位相α［ｔ］を用いて時刻［ｔ］における位相変位量Δα［ｔ］を演算するので、位相変位量演算部３１は、図示しないメモリを有し、そのメモリに位相推定部２２から出力された位相α［ｔ］を順次記憶していく。このメモリには、少なくとも時刻［ｔ−１］における位相α［ｔ−１］及び時刻［ｔ］における位相α［ｔ］を記憶しておく。
【００４９】
また、上記のように、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に、時刻［ｔ］における位相α［ｔ］が演算されるので、位相変位量Δα［ｔ］も、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算される。
【００５０】
［周波数推定部３２］
周波数推定部３２は、位相変位量演算部３１で演算された時刻［ｔ］における位相変位量Δα［ｔ］に基づいて、時刻［ｔ］における入力信号Ｖｉｎの角周波数の推定値ωｅｓｔ［ｔ］又は周波数の推定値ｆｅｓｔ［ｔ］を演算する。演算された入力信号Ｖｉｎの角周波数の推定値ωｅｓｔ［ｔ］又は周波数の推定値ｆｅｓｔ［ｔ］は、誤差周波数情報出力部７０に送られる。
なお、入力信号Ｖｉｎの周波数の推定値（検出値）は、周波数推定部３２から出力される推定周波数情報Ｆｅｓｔ［ｔ］を移動平均することによって、誤差を小さくし、精度のよい周波数情報を得ることができる。そこで、周波数推定部３２の後段に移動平均部（図略）を設けて、移動平均部（図略）の出力を入力信号Ｖｉｎの周波数の推定値（検出値）とすることができる。
【００５１】
ここで、位相変位量Δα［ｔ］は、サンプリング周期の間に生じた位相変位量であるから、式（９）に示すように、位相変位量Δα［ｔ］にサンプリング周波数ｆｓを乗算すると、１秒間に生じる位相変位量を求めることができる。すなわち、時刻［ｔ］における入力信号Ｖｉｎの角周波数の推定値ωｅｓｔ［ｔ］［単位：ｒａｄ／ｓ］を演算することができる。また、式（１０）に示すように、時刻［ｔ］における位相変位量Δα［ｔ］にサンプリング周波数ｆｓを乗算するとともに「２π」で除算することにより、時刻［ｔ］における周波数の推定値ｆｅｓｔ［ｔ］［単位：Ｈｚ］を演算することができる。
ωｅｓｔ［ｔ］＝Δα・ｆｓ              ・・・（９）
ｆｅｓｔ［ｔ］＝Δα・ｆｓ／（２π）    ・・・（１０）
【００５２】
また、角周波数の推定値ωｅｓｔ［ｔ］と設定角周波数ωｓｅｔとの差角周波数の推定値ωｄｉｆ［ｔ］は、式（１１）のように演算できる。なお、設定角周波数ωｓｅｔは、「２π・ｆｓｅｔ」である。また周波数の推定値ｆｅｓｔ［ｔ］と設定周波数ｆｓｅｔとの差周波数の推定値ｆｄｉｆ［ｔ］は、式（１２）のように演算できる。そのため、周波数推定部３２は、差角周波数の推定値ωｄｉｆ［ｔ］、差周波数の推定値ｆｄｉｆ［ｔ］を出力することもできる。この場合、例えば、周波数推定部３２に設定周波数ｆｓｅｔを入力すればよい。もちろん、周波数推定部３２に設定角周波数ωｓｅｔを入力してもよい。また設定周波数ｆｓｅｔから設定角周波数ωｓｅｔを演算してもよいし、その逆に、設定角周波数ωｓｅｔから設定周波数ｆｓｅｔを演算してもよい。
ωｄｉｆ［ｔ］    ＝ωｅｓｔ［ｔ］−ωｓｅｔ［ｔ］    ・・・（１１）
ｆｄｉｆ［ｔ］    ＝ｆｅｓｔ［ｔ］−ｆｓｅｔ［ｔ］    ・・・（１２）
【００５３】
なお、本明細書では、入力信号Ｖｉｎの角周波数の推定値ωｅｓｔ［ｔ］、周波数の推定値ｆｅｓｔ［ｔ］、差角周波数の推定値ωｄｉｆ［ｔ］及び差周波数の推定値ｆｄｉｆ［ｔ］を総称して時刻［ｔ］における入力信号Ｖｉｎの「推定周波数情報Ｆｅｓｔ［ｔ］」という。
【００５４】
また、入力信号Ｖｉｎの角周波数の推定値ωｅｓｔ［ｔ］、周波数の推定値ｆｅｓｔ［ｔ］、差角周波数の推定値ωｄｉｆ［ｔ］及び差周波数の推定値ｆｄｉｆ［ｔ］は、いずれか１つを出力すればよいが、複数を出力するようにしてもよい。
【００５５】
なお、上記のように、時刻［ｔ］における位相変位量Δα［ｔ］は、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算されるので、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報Ｆｅｓｔ［ｔ］も、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算される。
【００５６】
＜正弦波信号生成部４０＞
正弦波信号生成部４０は、設定周波数情報Ｆｓｅｔに後述する補正周波数情報Ｆｒｅｖ[ｔ]を加算した同期周波数情報Ｆｓｙｃ[ｔ]を指令値として、同期周波数情報Ｆｓｙｃ[ｔ]と同じ周波数（角周波数で表現してもよい）の正弦波信号Ｖｏｕｔ［ｔ］を生成して出力する。
なお、内部の発振器４１（図３参照）から出力する正弦波信号Ｖｏｕｔの時刻［ｔ］における瞬時値をＶｏｕｔ［ｔ］としている。
【００５７】
この正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔ［ｔ］は、入力信号Ｖｉｎ［ｔ］に同期する信号となる。すなわち、後述するように内部の発振器４１から出力する正弦波信号Ｖｏｕｔの周波数及び位相が、入力信号Ｖｉｎの周波数及び位相に近づくように収束していく。その結果、内部の発振器４１から出力する正弦波信号Ｖｏｕｔは、入力信号Ｖｉｎに同期する信号となる。
【００５８】
図３は、正弦波信号生成部４０の構成例を示す図である。この図３に示すように、正弦波信号生成部４０は、例えば図３（ａ）〜（ｂ）のような構成例が考えられるが、いずれの場合も、正弦波信号生成部４０には、例えば、ダイレクト・デジタル・シンセサイザー（ＤＤＳ））で構成される発振器４１が備わっている。
発振器４１は、他の高周波電源（例えば図６及び図７に示した第２高周波電源１２０）の発振器として機能する。なお、「他の高周波電源」の「他」とは、入力信号Ｖｉｎの発生源である高周波電源（この実施形態では図１に示した第１高周波電源１１０）に対して他の高周波電源という意味である。この図３では、他の高周波電源４３と表記している。
【００５９】
（１）図３（ａ）の正弦波信号生成部４０は、発振器４１とＡ／Ｄ変換部４２とを備えている。この例では、発振器４１の出力を他の高周波電源４３に送り、他の高周波電源４３において無線周波数帯域の周波数を有する進行波電力ＰＦを出力させるとともに、発振器４１から出力される正弦波信号ＶｏｕｔをＡ／Ｄ変換部４２を介して第２推定位相情報出力部５０に送る。
【００６０】
ここで、発振器４１から出力される正弦波信号Ｖｏｕｔの周波数をｆｏｕｔ、時間をｔ、位相オフセットをθ”、角周波数をω（＝２π・ｆ）とすると、発振器４１から出力される正弦波信号Ｖｏｕｔは、式（１３）のように表すことができる。なお、ここでは、正弦波信号Ｖｏｕｔの振幅を「１」としている。
Ｖｏｕｔ＝ｓｉｎ（２π・ｆｏｕｔ・ｔ＋θ”）
＝ｓｉｎ（ω・ｔ＋θ”）  ・・・（１３）
【００６１】
Ａ／Ｄ変換部４２は、Ａ／Ｄ変換部１０と同様に、発振器４１から出力される正弦波信号Ｖｏｕｔを予め定めたサンプリング周期（サンプリング周波数ｆｓの逆数：１／ｆｓ）でデジタル信号に変換する。Ａ／Ｄ変換されたサンプリングデータは、時刻［ｔ］における正弦波信号Ｖｏｕｔの瞬時値Ｖｏｕｔ［ｔ］として順次出力される。これにより、交流のアナログ信号波形が、複数のサンプリングデータで構成されるデジタル信号波形に変換される。なお、発振器４１から出力される正弦波信号Ｖｏｕｔが正弦波信号であれば、Ａ／Ｄ変換部４２から出力されるデジタル信号波形も正弦波信号となる。
なお、Ａ／Ｄ変換部４２には、Ａ／Ｄコンバータが備わっていて、上記のようなＡ／Ｄ変換を行う。また、Ａ／Ｄコンバータの後段にローパスフィルタ又はバンドパスフィルタを設けて、基本周波数成分を通過させるようにしてもよい。
【００６２】
ここで、正弦波信号Ｖｏｕｔの周波数をｆｏｕｔ、サンプリング周波数をｆｓ、位相オフセットをθ”、相対角周波数をωｏｕｔ＝２π・（ｆｏｕｔ／ｆｓ）とすると、時刻［ｔ］における正弦波信号Ｖｏｕｔの瞬時値Ｖｏｕｔ［ｔ］は、式（１４）のように表すことができる。なお、ここでは、正弦波信号Ｖｏｕｔの振幅を「１」としている。また時間データ「ｔ」はサンプリング周期毎にインクリメントされる変数である。
Ｖｏｕｔ［ｔ］＝ｓｉｎ（２π・（ｆｏｕｔ／ｆｓ）［ｔ］＋θ”）
＝ｓｉｎ（ωｏｕｔ［ｔ］＋θ”）            ・・・・・（１４）
【００６３】
（２）図３（ｂ）の正弦波信号生成部４０は、発振器４１の出力を他の高周波電源４３で増幅し、その出力である進行波電圧ＶＦ２を方向性結合器４４で検出する。そして、その検出信号である進行波検出信号Ｖｆ２をフィルタ４５及びＡ／Ｄ変換部４６を介して第２推定位相情報出力部５０に送る。
なお、方向性結合器４４、フィルタ４５及びＡ／Ｄ変換部４６は、それぞれ図１に示した方向性結合器１５０、フィルタ１６０及びＡ／Ｄ変換部１０と同様の機能を有する。そのため詳細説明は省略する。また、図３（ｂ）では、フィルタ４５の出力を正弦波信号Ｖｏｕｔ’、Ａ／Ｄ変換部４６の出力を正弦波信号Ｖｏｕｔ’［ｔ］としている。また、説明を簡略化するために、図３（ａ）及び図３（ｂ）では、発振器４１の符号を同じにしている。
【００６４】
図３（ｂ）の場合でも、発振器４１から出力される正弦波信号Ｖｏｕｔに基づいて他の高周波電源４３から進行波電圧ＶＦ２が出力されるので、正弦波信号Ｖｏｕｔと進行波電圧ＶＦ２とは、周波数及び位相が同期する。そのため、その後段にある方向性結合器４４、フィルタ４５及びＡ／Ｄ変換部４６を介して第２推定位相情報出力部５０に送られる正弦波信号Ｖｏｕｔ’［ｔ］も入力信号Ｖｉｎに同期する信号となる。
【００６５】
図２及び図３に示す正弦波信号生成部４０から出力する正弦波信号Ｖｏｕｔ[ｔ]の周波数は、同期周波数情報Ｆｓｙｃ[ｔ]である。そして同期周波数情報Ｆｓｙｃ[ｔ]は、第１高周波電源１１０からプラズマ処理装置１３０に向かう進行波電圧ＶＦの周波数を表すので、他の高周波電源用の発振器４１の指令信号として同期周波数情報Ｆｓｙｃ[ｔ]を用いることにより、第１高周波電源１１０の出力周波数と他の高周波電源の出力周波数とを合わせることができる（同期させることができる）。
【００６６】
＜第２推定位相情報出力部５０＞
第２推定位相情報出力部５０は、第１推定位相情報出力部２０と同様の構成であるが、第１推定位相情報出力部２０が、時刻［ｔ］における入力信号Ｖｉｎを入力とし、入力信号Ｖｉｎの位相α［ｔ］を出力するのに対して、第２推定位相情報出力部５０は、正弦波信号生成部４０から出力された時刻［ｔ］における正弦波信号Ｖｏｕｔ［ｔ］を入力とし、正弦波信号Ｖｏｕｔ［ｔ］の位相α２［ｔ］を出力する点が異なる。その他は、第１推定位相情報出力部２０と同様なので、説明を省略する。
【００６７】
＜第２推定周波数情報出力部６０＞
第２推定周波数情報出力部６０は、第１推定周波数情報出力部３０と同様の構成であるが、第１推定周波数情報出力部３０が、第１推定位相情報出力部２０から出力された入力信号Ｖｉｎの位相α［ｔ］を入力し、時刻［ｔ］における入力信号Ｖｉｎの周波数情報の推定値を第１推定周波数情報Ｆｅｓｔ［ｔ］として出力するのに対して、第２推定周波数情報出力部６０は、第２推定位相情報出力部５０から出力された時刻［ｔ］における正弦波信号Ｖｏｕｔ［ｔ］の位相α２［ｔ］を入力し、時刻［ｔ］における正弦波信号Ｖｏｕｔ［ｔ］の周波数情報の推定値を第２推定周波数情報Ｆｅｓｔ２［ｔ］として出力する点が異なる。その他は、第１推定周波数情報出力部３０と同様なので、説明を省略する。
【００６８】
＜誤差周波数情報出力部７０＞
誤差周波数情報出力部７０は、時刻［ｔ］における第１推定周波数情報Ｆｅｓｔ［ｔ］に対する第２推定周波数情報Ｆｅｓｔ２［ｔ］の誤差を、時刻［ｔ］における誤差周波数情報Ｆｅｒｒ［ｔ］として出力する。この誤差周波数情報出力部７０は、誤差周波数情報Ｆｅｒｒ［ｔ］の移動平均値を演算する移動平均部７１が備わっている。
【００６９】
なお、第１推定周波数情報Ｆｅｓｔ［ｔ］と第２推定周波数情報Ｆｅｓｔ２［ｔ］とは同じ種類の情報で比較する必要がある、例えば、第１推定周波数情報Ｆｅｓｔ［ｔ］が周波数で表されているのであれば、第２推定周波数情報Ｆｅｓｔ２［ｔ］も周波数で表されている必要がある。角周波数、差周波数又は差角周波数の場合も同様である。
【００７０】
［移動平均部７１］
移動平均部７１は、誤差周波数情報Ｆｅｒｒ［ｔ］の移動平均値を演算し、誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］として出力する（周波数情報検出装置１の外部に出力することもできる）。移動平均値は、予め定めたデータ数を用いて演算される。
【００７１】
上記のように、移動平均部７１は、誤差演算された時刻［ｔ］における誤差周波数情報Ｆｅｒｒ［ｔ］に基づいて、予め定めたデータ数の移動平均値を演算するので、移動平均部７１は、図示しないメモリを有し、そのメモリに誤差周波数情報Ｆｅｒｒ［ｔ］を順次記憶していく。このメモリには、少なくとも移動平均値の演算に必要な予め定めた数のデータを記憶しておく。
【００７２】
また、上記のように、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報Ｆｅｓｔ［ｔ］が演算されるので、角周波数の移動平均値ωａｖｅ［ｔ］、周波数の移動平均値ｆａｖｅ［ｔ］も、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算される。そのため、対応する誤差周波数情報Ｆｅｒｒ［ｔ］も、時刻［ｔ＋１］における入力信号Ｖｉｎの瞬時値Ｖｉｎ［ｔ＋１］が余弦値推定部２１に入力された後に演算される。
【００７３】
また、データのサンプリングを開始した直後は、移動平均値を演算するために必要な
予め定めた数のデータがメモリに保存されていない。そのため、予め定めた数のデータがメモリに保存されていない間は、予め定めた数よりも少ないデータを用いて移動平均値を演算してもよい。または、予め定めた数のデータがメモリに保存された後に、移動平均値を演算してもよい。どのように移動平均値を演算するかは、予め定めておけばよい。
【００７４】
＜移動平均値を演算する理由＞
上記のように、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））には、誤差が含まれているので、位相α［ｔ］にも誤差が含まれる。もちろん、時刻［ｔ−１］における余弦値（ｃｏｓ（α）［ｔ−１］）、位相α［ｔ−１］にも誤差が含まれる。その結果、位相変位量Δα［ｔ］にも誤差が含まれる。そのため、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報Ｆｅｓｔ［ｔ］にも誤差が含まれる。
【００７５】
すなわち、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報Ｆｅｓｔ［ｔ］に含まれる誤差は、時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））に起因する。
この時刻［ｔ］における余弦値（ｃｏｓ（α［ｔ］））の誤差は、上記のように正になることもあれば負になることもあるので、時刻［ｔ］における入力信号Ｖｉｎの推定周波数情報Ｆｅｓｔ［ｔ］の誤差も正になることもあれば負になることもある。そこで、移動平均部７１によって、複数のデータの移動平均を行えば、誤差が相殺するので、誤差を小さくし、精度のよい周波数情報を得ることができる。
【００７６】
誤差周波数情報出力部７０での移動平均は、第１推定周波数情報Ｆｅｓｔ［ｔ］に対して行われるわけではなく、時刻［ｔ］における誤差周波数情報Ｆｅｒｒ［ｔ］に対して行われる。この際、第２推定周波数情報Ｆｅｓｔ２［ｔ］にも誤差が含まれるため、上記と同様に移動平均の効果を得ることができる。
【００７７】
このように、移動平均部７１を設けて移動平均を行うことが好ましいが、第１推定周波数情報Ｆｅｓｔ［ｔ］及び第２推定周波数情報Ｆｅｓｔ２［ｔ］の誤差が小さければ省略可能な場合もある。
【００７８】
＜誤差位相情報出力部８０＞
誤差位相情報出力部８０は、時刻［ｔ］における第１推定位相情報位相α［ｔ］に対する第２推定位相情報α２［ｔ］の誤差を、時刻［ｔ］における誤差位相情報αｅｒｒ［ｔ］として出力する。この誤差位相情報出力部８０は、誤差位相情報αｅｒｒ［ｔ］の移動平均値を演算する移動平均部８１が備わっている。
【００７９】
［移動平均部８１］
移動平均部８１は、誤差位相情報αｅｒｒ［ｔ］の移動平均値を演算し、誤差位相移動平均値情報αｅｒｒ＿ａｖｅ［ｔ］として出力する（周波数情報検出装置１の外部に出力することもできる）。移動平均値は、予め定めたデータ数を用いて演算される。
移動平均部７１が誤差周波数情報Ｆｅｒｒ［ｔ］の移動平均値を演算するのに対して、移動平均部８１は、誤差位相情報αｅｒｒ［ｔ］の移動平均値を演算する点が異なるがその他は、同様なので説明を省略する。
【００８０】
＜周波数補正部９０＞
周波数補正部９０は、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］）が予め定めた範囲外の場合は、誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が小さくなるように時刻［ｔ］における補正周波数情報Ｆｒｅｖ[ｔ]を定めて出力する。また、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲内の場合は、誤差位相移動平均値情報αｅｒｒ＿ａｖｅ［ｔ］が小さくなるように補正周波数情報Ｆｒｅｖ[ｔ]を定めて出力する。以下、周波数補正部９０について、更に説明する。
【００８１】
上述したように、設定周波数情報Ｆｓｅｔと第１高周波電源１１０からプラズマ処理装置１３０に向かう進行波電圧ＶＦの周波数とに誤差が生じる場合がある。そのために、実際の進行波電圧ＶＦの周波数を検出し、その周波数の情報を他の高周波電源（例えば図６及び図７に示した第２高周波電源１２０）の発振器の指令信号にすることにより、第１高周波電源１１０の出力周波数と他の高周波電源の出力周波数とを合わせることができる（同期させることができる）。
【００８２】
なお、進行波電圧ＶＦの周波数は、周波数推定部３２から出力される推定周波数情報Ｆｅｓｔ［ｔ］によって得られる。また、誤差を低減させて精度のよい周波数情報の推定値を得るために、周波数推定部３２の後段に移動平均部を設けている。
【００８３】
しかし、周波数推定部３２から出力される推定周波数情報Ｆｅｓｔ［ｔ］だけでは、進行波電圧ＶＦの位相情報を得ることができない。
そこで、進行波電圧ＶＦの検出信号である進行波検出信号Ｖｆ（より具体的には、フィルタ１６０を通過した入力信号Ｖｉｎ）と同期させるための正弦波信号Ｖｏｕｔを正弦波信号生成部４０で発生させている。また、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの周波数及び位相を入力信号Ｖｉｎの周波数及び位相に近づける制御を行っている。
そして、正弦波信号Ｖｏｕｔの周波数及び位相が入力信号Ｖｉｎの周波数及び位相に近づき、両者の誤差が十分に小さい状態にすることで、第１高周波電源１１０と他の高周波電源の出力周波数及び位相を合わせることができる（同期させることができる）。
【００８４】
そのために、誤差周波数情報出力部７０によって、時刻［ｔ］における第１推定周波数情報Ｆｅｓｔ［ｔ］に対する第２推定周波数情報Ｆｅｓｔ２［ｔ］の誤差を求めている。この誤差が誤差周波数情報Ｆｅｒｒ［ｔ］として出力され、移動平均部７１を介して誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が周波数補正部９０に入力される。
【００８５】
一方、位相情報に関しては、誤差位相情報出力部８０によって、時刻［ｔ］における第１推定位相情報位相α［ｔ］に対する第２推定位相情報α２［ｔ］の誤差を求めている。この誤差が誤差位相情報αｅｒｒ［ｔ］として出力され、移動平均部８１を介して誤差位相移動平均値情報αｅｒｒ＿ａｖｅ［ｔ］が周波数補正部９０に入力される。
【００８６】
そして、上述したように、周波数補正部９０では、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲外の場合は、誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が小さくなるように時刻［ｔ］における補正周波数情報Ｆｒｅｖ[ｔ]を定めて出力する。
すなわち、第１推定周波数情報Ｆｅｓｔ［ｔ］の方が、第２推定周波数情報Ｆｅｓｔ２［ｔ］よりも大きい（周波数が高い）場合は、正弦波信号生成部４０への指令値である同期周波数情報Ｆｓｙｃ[ｔ]が大きく（周波数が高く）なるように、正（プラス）の補正周波数情報Ｆｒｅｖ[ｔ]を出力する。反対に、第１推定周波数情報Ｆｅｓｔ［ｔ］の方が、第２推定周波数情報Ｆｅｓｔ２［ｔ］よりも小さい（周波数が低い）場合は、負（マイナス）の補正周波数情報Ｆｒｅｖ[ｔ]を出力する。
【００８７】
なお、第１推定周波数情報Ｆｅｓｔ［ｔ］と第２推定周波数情報Ｆｅｓｔ２［ｔ］との誤差の大きさ（より詳細には誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］の大きさ）に応じて補正周波数情報Ｆｒｅｖ[ｔ]が定まるが、誤差に乗じる係数（Ｋｆ）を調整できるようにしておくと、制御がし易くなるので好ましい。
【００８８】
上記のような制御を行うことによって、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの周波数を入力信号Ｖｉｎの周波数に近づけることができる。
【００８９】
そして、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲内になると、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの周波数が入力信号Ｖｉｎの周波数に近づいたと判断して、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの位相を入力信号Ｖｉｎの位相に近づける制御に移行する。
【００９０】
このとき、周波数補正部９０では、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲内の場合は、誤差位相情報αｅｒｒ［ｔ］が小さくなるように補正周波数情報Ｆｒｅｖ[ｔ]を定めて出力する。
すなわち、第１推定位相情報位相α［ｔ］の方が、第２推定位相情報α２［ｔ］よりも大きい場合は、正弦波信号生成部４０への指令値である同期周波数情報Ｆｓｙｃ[ｔ]が大きく（周波数が高く）なるように、正（プラス）の補正周波数情報Ｆｒｅｖ[ｔ]を出力する。反対に、第１推定位相情報位相α［ｔ］の方が、第２推定位相情報α２［ｔ］よりも小さい場合は、負（マイナス）の補正周波数情報Ｆｒｅｖ[ｔ]を出力する。
【００９１】
なお、第１推定位相情報位相α［ｔ］と第２推定位相情報α２［ｔ］との誤差の大きさ（より詳細には誤差位相移動平均値情報αｅｒｒ＿ａｖｅ［ｔ］の大きさ）に応じて補正周波数情報Ｆｒｅｖ[ｔ]が定まるが、誤差に乗じる係数（Ｋα）を調整できるようにしておくと、制御がし易くなるので好ましい。
【００９２】
周波数補正部９０から補正周波数情報Ｆｒｅｖ[ｔ]が出力されると、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの周波数が変化するので、第１推定位相情報位相α［ｔ］に対する第２推定位相情報α２［ｔ］の関係が変化して、誤差位相情報αｅｒｒ［ｔ］が小さくっていく。
【００９３】
そして、誤差位相情報αｅｒｒ［ｔ］が予め定めた範囲内になると、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの位相が入力信号Ｖｉｎの位相に近づいたと判断する。
上記のような制御を行うことによって、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの位相を入力信号Ｖｉｎの位相に近づけることができる。
【００９４】
なお、正弦波信号生成部４０から出力される正弦波信号Ｖｏｕｔの位相を入力信号Ｖｉｎの位相に近づける制御を行う際に、正弦波信号Ｖｏｕｔの周波数を変更しているので、誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲外になっていることが考えられる。そのため、周波数補正部９０は、再度、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲外になっているか否かを確認し、上記のような制御を繰り返して実行する。
【００９５】
これにより、正弦波信号Ｖｏｕｔの周波数及び位相が、入力信号Ｖｉｎの周波数及び位相に近づくように収束していく。そして、周波数補正部９０は、時刻［ｔ］における誤差周波数移動平均値情報Ｆｅｒｒ＿ａｖｅ［ｔ］）が予め定めた範囲内であり、且つ、誤差位相移動平均値情報αｅｒｒ＿ａｖｅ［ｔ］が予め定めた範囲内になったときに、収束したと判定する。なお、収束したと判定した際には、その旨を示す収束信号Ｓｃｏｎを外部に出力してもよい。
【００９６】
＜シミュレーション結果＞
上記のように、進行波電圧ＶＦの周波数は、周波数推定部３２から出力される推定周波数情報Ｆｅｓｔ［ｔ］を移動平均することによって得ることができる。そこで、周波数推定部３２の後段に移動平均部を設けた場合のシミュレーション結果を示す。
【００９７】
図４は、入力信号Ｖｉｎの周波数ｆｉｎが設定周波数ｆｓｅｔからずれていると想定した場合に、移動平均部７１から出力される差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］のシミュレーション結果である。図４（ａ）は、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］のシミュレーション結果であり、図４（ｂ）は、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］のシミュレーション結果であり、図４（ｃ）は、図４（ｂ）の一部拡大図である。また、図４（ａ）及び図４（ｂ）は、サンプリング開始後、約１５０〜１６０μｓ間のデータであり、図４（ｃ）は、サンプリング開始後、約１５９〜１６０μｓ間のデータである。
【００９８】
なお、シミュレーション条件は次のとおりである。
（１）サンプリング周波数ｆｓ      ：５０ＭＨｚ
（２）設定周波数ｆｓｅｔ          ：１３．５６ＭＨｚ
（３）入力信号Ｖｉｎの周波数ｆｉｎ：１３．５６３ＭＨｚ
（設定周波数ｆｓｅｔと３，０００Ｈｚだけずれていると想定）
（４）移動平均値の演算に用いるデータ数：５００個（１０μｓの移動平均値）
【００９９】
図５は、入力信号Ｖｉｎの周波数ｆｉｎが設定周波数ｆｓｅｔからずれていると想定した場合に、移動平均部７１から出力される差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］の他のシミュレーション結果である。図５（ａ）は、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］のシミュレーション結果であり、図５（ｂ）は、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］のシミュレーション結果であり、図５（ｃ）は、図５（ｂ）の一部拡大図である。
なお、この図５は、入力信号Ｖｉｎの周波数ｆｉｎが１３．５６０３ＭＨｚ（設定周波数ｆｓｅｔと３００Ｈｚだけずれていると想定）であることを除き、図４と同じ条件でのシミュレーション結果である。
【０１００】
図４（ａ）のように、入力信号Ｖｉｎの周波数が設定周波数ｆｓｅｔ（１３．５６ＭＨｚ）に対して３，０００Ｈｚずれていると想定した場合、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］は、検出すべき周波数（３，０００Ｈｚ）に対して約±４００Ｈｚの範囲でばらついているが、図４（ｂ）及び図４（ｃ）に示すように、差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］は、検出すべき周波数（３，０００Ｈｚ）に対して約±０．６Ｈｚの範囲に収まっていることが分かる。
このように、入力信号Ｖｉｎの周波数が設定周波数ｆｓｅｔに対して大幅にずれている場合であっても、移動平均部７１から出力される差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］は、検出すべき周波数に対して誤差が小さい。
【０１０１】
また、図５（ａ）のように、入力信号Ｖｉｎの周波数ｆｉｎが設定周波数ｆｓｅｔ（１３．５６ＭＨｚ）に対して３００Ｈｚずれていると想定した場合、入力信号Ｖｉｎの差周波数の推定値ｆｄｉｆ［ｔ］は、検出すべき周波数（３００Ｈｚ）に対して約±４０Ｈｚの範囲でばらついているが、図５（ｂ）及び図５（ｃ）に示すように、差周波数の推定値ｆｄｉｆ［ｔ］の移動平均値ｆｄｉｆ＿ａｖｅ［ｔ］は、検出すべき周波数（３００Ｈｚ）に対して約±０．０５Ｈｚの範囲に収まっていることが分かる。したがって、精度良く入力信号Ｖｉｎの周波数を検出できていることが分かる。
【０１０２】
すなわち、入力信号Ｖｉｎの周波数ｆｉｎと設定周波数ｆｓｅｔとの誤差が小さい程、より精度良く入力信号Ｖｉｎの周波数を検出できていることが分かる。
上記の実施形態では、正弦波信号Ｖｏｕｔの周波数が、入力信号Ｖｉｎの周波数に近づくように制御するので、制御を行うにつれて、より精度良く入力信号Ｖｉｎの周波数を検出することができることを示している。
【０１０３】
なお、本実施例の周波数情報検出装置は、アナログの正弦波信号Ｖｉｎを、予め定めたサンプリング周期（サンプリング周波数ｆｓの逆数：１／ｆｓ）でデジタル信号にＡ／Ｄ変換することによって得られるデジタルの正弦波信号の周波数情報を検出する周波数情報検出装置であり、特に用途は限定されない。そのため、上記では、周波数情報検出装置を、例えば図６及び図７のような高周波電力供給システムに用いる例を示したが、これに限定されるものではない。例えば、商用周波数帯域（５０Ｈｚ〜６０Ｈｚ）でも適用できる。
【符号の説明】
【０１０４】
１      周波数情報検出装置
１０    Ａ／Ｄ変換部
２０    第１推定位相情報出力部
２１    余弦値推定部
２２    位相推定部
３０    第１推定周波数情報出力部
３１    位相変位量演算部
３２    周波数推定部
４０    移動平均部
５０    第２推定位相情報出力部５０
６０    第２推定周波数情報出力部６０
７０    誤差周波数情報出力部
８０    誤差位相情報出力部８０
９０    周波数補正部
１００    発振器
１１０    第１高周波電源
１２０    第２高周波電源
１３０    プラズマ処理装置（負荷）
１３１    電極
１３２    電極
１３３    プラズマ
１４０    伝送線路
１５０    方向性結合器
１６０    フィルタ
【図１】
0006396157.tif 000002
【図２】
0006396157.tif 000003
【図３】
0006396157.tif 000004
【図４】
0006396157.tif 000005
【図５】
0006396157.tif 000006
【図６】
0006396157.tif 000007
【図７】
0006396157.tif 000008