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