(57)【要約】【課題】本発明は、デ−タを処理するデータ処理装置を備えるプラズマ処理装置、デ−タを処理するデータ処理装置およびデータ処理方法において、高Ｓ／Ｎ比とデ−タ処理応答性向上を両立でき、該データ処理に最適なパラメータを自動で設定できるプラズマ処理装置、データ処理装置およびデータ処理方法を提供する。【解決手段】本発明は、試料がプラズマ処理される処理室と、前記処理室内にプラズマを生成するための高周波電力を供給する高周波電源と、前記プラズマの発光データが処理されるデータ処理装置とを備えるプラズマ処理装置において、前記データ処理装置により入力データと平滑化されたデータの予測値との誤差に基づいて平滑化パラメータを変化させる応答型２段指数平滑化法を用いて前記発光データが処理され、前記平滑化パラメータの応答係数は、前記誤差をパラメータとする確率密度関数を用いて求められることを特徴とする。【選択図】図２
2020061569.tif 000002
【特許請求の範囲】
【請求項１】
２段指数平滑化法を用いてデータが処理されるデータ処理装置において、
入力データと平滑化されたデータの予測値との誤差に基づいて平滑化パラメータを変化させる応答型２段指数平滑化法を用いて前記データが処理され、
前記入力データの最新時における多項式適合法により平滑化された前記入力データが前記平滑化されたデータの予測値を算出する式の入力データとして用いられ、
前記平滑化されたデータの予測値の最新時における多項式適合法により一次微分された前記平滑化されたデータの予測値が前記平滑化されたデータの傾きの予測値を算出する式の平滑化されたデータの予測値の傾きとして用いられることを特徴とするデータ処理装置。
【請求項２】
請求項１に記載のデータ処理装置において、
前記平滑化パラメータの応答係数は、前記誤差をパラメータとする確率密度関数を用いて求められることを特徴とするデータ処理装置。
【請求項３】
請求項１に記載のデータ処理装置において、
Ｎを０以上の整数とした場合、前記平滑化パラメータの応答係数は、前記誤差の相対値の予測値を前記誤差の絶対値の予測値で除して定数が加えられた値をＮ乗して求められることを特徴とするデータ処理装置。
【請求項４】
２段指数平滑化法を用いてデータが処理されるデータ処理装置において、
入力データと平滑化されたデータの予測値との誤差に基づいて平滑化パラメータを変化させる応答型２段指数平滑化法を用いて前記データが処理され、
前記誤差が指数荷重移動平均され、
前記指数荷重移動平均の係数は、平滑化されたデータの傾きの予測値に基づいて変化させられることを特徴とするデータ処理装置。
【請求項５】
請求項１または請求項２に記載のデータ処理装置において、
前記誤差が指数荷重移動平均され、
前記指数荷重移動平均の係数は、平滑化されたデータの傾きの予測値に基づいて変化させられることを特徴とするデータ処理装置。
【請求項６】
２段指数平滑化法を用いてデータを処理するデータ処理方法において、
入力データと平滑化されたデータの予測値との誤差に基づいて平滑化パラメータを変化させる応答型２段指数平滑化法を用いて前記データを処理し、
前記入力データの最新時における多項式適合法により平滑化された前記入力データを前記平滑化されたデータの予測値を算出する式の入力データとして用い、
前記平滑化されたデータの予測値の最新時における多項式適合法により一次微分された前記平滑化されたデータの予測値を前記平滑化されたデータの傾きの予測値を算出する式の平滑化されたデータの予測値の傾きとして用いることを特徴とするデータ処理方法。
【請求項７】
請求項６に記載のデータ処理方法において、
前記誤差をパラメータとする確率密度関数を用いて前記平滑化パラメータの応答係数を求めることを特徴とするデータ処理方法。
【請求項８】
請求項６に記載のデータ処理方法において、
Ｎを０以上の整数とした場合、前記誤差の相対値の予測値を前記誤差の絶対値の予測値で除して定数が加えられた値をＮ乗して前記平滑化パラメータの応答係数を求めることを特徴とするデータ処理方法。
【請求項９】
２段指数平滑化法を用いてデータを処理するデータ処理方法において、
入力データと平滑化されたデータの予測値との誤差に基づいて平滑化パラメータを変化させる応答型２段指数平滑化法を用いて前記データを処理し、
前記誤差を指数荷重移動平均し、
平滑化されたデータの傾きの予測値に基づいて前記指数荷重移動平均の係数を変化させることを特徴とするデータ処理方法。
【請求項１０】
請求項６または請求項７に記載のデータ処理方法において、
前記誤差を指数荷重移動平均し、
平滑化されたデータの傾きの予測値に基づいて前記指数荷重移動平均の係数を変化させることを特徴とするデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、対象とされる系（装置、分析データ等）の一連の数値データにおけるデータ処理装置およびデータ処理方法に関する。
【背景技術】
【０００２】
対象とされる系において、一連の数値データをデータスムージング、データ予測する技術として、種々のモデルが提案されている。また、対象とされる系の状態を把握する目的で、データの変化点である極値（最大値、最小値）あるいは変曲点を検出するため、データの1次微分、２次微分等の微分処理が多用される。特に、ノイズ等を含む計測データ等の時系列データの場合は、高精度にデータスムージング処理、微分処理を行うことにより対象とされる系での変化点を検出し、対象とされる系の制御を行うことが重要な技術課題となっている。
【０００３】
従来のデータスムージング、データ予測する技術は、非特許文献１に記載されたようにカ−ブフィッティング法、移動平均法等がある。カーブフィッティング法には、特許文献１に記載されたように多項式適合法（Savitzky-Golay法）がある。また、デジタルフィルタとしては、バターワース型ローパスフィルタ等がある。さらに、移動平均法には、非特許文献２に記載されたように、指数平滑化法等がある。
【０００４】
非特許文献２には、1段指数平滑法（平滑化パラメータが１ケ）、2段指数平滑化法（平滑化パラメ−タが２ケ）が開示されており、需要供給予測等の経済関係の分野で利用されている。また、非特許文献３には、入力データ（特に逐次入力データ）に応じて平滑化係数が変化する、応答型１段指数平滑化法が開示されている。平滑化係数は、｜相対誤差/絶対誤差｜に応じて変化するケースとロジスティック関数に応じて変化するケースが記載されている。一般には、式(1)〜(4)に示す、平滑化係数を｜相対誤差/絶対誤差｜に応じて変化するケースを使用するケースが多い。
【０００５】
データの平滑化：S1t+1 = α1t Y1t + (1 - α1t) S1t                      式(1)
平滑化係数：  α1t =  ｜δαt/Δαt｜                                  式(2)
相対誤差：    δαt = A1（Y1t - S1t） + (1 - A1) δαt-1               式(3)
絶対誤差：    Δαt = A1 |Y1t - S1t| + (1 - A1) Δαt-1                式(4)
ここで、入力データは、例えば、時系列データY1t：t=1、2、・・・とし、1期先の出力のデータの平滑化の予測値をSt+1とする。またA1は、任意の定数である。尚、各記号は、非特許文献３と異なるが、なるべく以下の実施例での記号と同一となるように記載した。
【０００６】
従来のデータ微分処理する技術は、差分法が多用されることが多い。また、特許文献１、非特許文献４に記載されたように、多項式適合法（Savitzky-Golay法）が使用される場合もある。通常、多項式適合法(Savitzky-Golay法）は、複数の連続したデータを使用し、この連続データの期間の中心点での微分処理結果を導出する場合が多い。しかし、非特許文献４には、使用するデータの任意の点で、多項式適合法(Savitzky-Golay法）により、データスムージング処理および1次微分処理を行う方法が開示されている。
【０００７】
ノイズ等を含む計測データ等の時系列データをデータスムージング処理し、１次微分処理、２次微分処理によりデータの変化点を検出し、対象とされる系を制御する例としては、特許文献２に「プラズマエッチング処理装置において、プラズマ発光からの分光強度信号データを移動平均処理によりデータスムージング処理し、１次微分値、２次微分値により、エッチング終点を判定する方法」が開示されている。また、非特許文献３には、1段指数平滑化法において、1期先予測値の誤差の総和を最小とすることにより、最適な平滑化パラメータを推定する方法が開示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０００−２２８３９７号公報
【特許文献２】特開昭６１-５３７２８号公報
【非特許文献】
【０００９】
【非特許文献１】「データ処理の裏側」高橋和裕著、Journal of Surface Analysis Vol7.No 1(2000) p68-p77.
【非特許文献２】「時系列モデル入門」A.C.ハーベイ著、国友直人/山本拓訳、(1985)、p173、東京大学出版会、(原著：TIME SERIES MODELS by A.C.Harvey)
【非特許文献３】James W Taylor, Journal of Forecasting, 2004, (23), pp385-394.
【非特許文献４】Peter A. Gorry, Anal. Chem. 1990, (62), pp570-573.
【発明の概要】
【発明が解決しようとする課題】
【００１０】
計測データ等のように真の信号（シグナル）にノイズが加算されたデータよりノイズを除去するデータスムージング処理の場合、Ｓ/Ｎ比（シグナル/ノイズ比）を大きくすることが技術課題の１つである。データ微分処理の場合においてもＳ/Ｎ比（シグナル/ノイズ比）を大きくすることが技術課題の１つである。
【００１１】
例えば、特許文献１に記載されたように差分法によりデータ微分処理を行うとデータスムージング処理が十分でない場合、１次微分処理の出力結果は、ノイズが多く含まれ、スムーズ（平滑な）なデータでなく、S/N比（シグナル/ノイズ比）が小さくなる。上記データを使用して再度、差分法によりデータ微分処理を行うと、２次微分処理の出力結果は、更にノイズが増大し、スムーズなデータでなく、更にS/N比が小さくなる。即ち、データスムージング処理よりもデータ1次微分処理、更にデータ高次微分処理になる程、Ｓ/Ｎ比（シグナル/ノイズ比、以下、S/N比と称する）を大きくすることが難しい技術課題となる。
【００１２】
データ処理には、全てのデータを取得した後、処理を行う方法（オフライン処理）とデータを取得しながら、処理を行う方法（オンライン処理、あるいは逐次データ処理）がある。オンライン処理の場合、データ処理結果をリアルタイムで利用して、例えば、プロセス処理装置の制御等に活用する場合が多い。このため、データ処理結果の遅延時間（真の処理結果との時間的なズレ量）が小さいことが技術課題の１つである。
【００１３】
差分法によるデータ微分処理で差分間隔を増加させると、1次微分処理の出力結果および２次微分処理の出力結果は、スムーズなデータとなり、S/N比は増大するが（但し、差分間隔には最適値がある）、時系列データ等のオンライン処理の場合、データ微分処理によるデータ遅延時間が増加するという課題がある。また、ローパスフィルタ等では、カットオフ周波数を減少させると、1次微分処理の出力データおよび２次微分処理の出力データは、スムーズなデータとなり、S/N比は増大するが（但し、カットオフ周波数には最適値がある）、上記と同様に、データ微分処理によるデータ遅延時間が増加するという課題がある。
【００１４】
また、多項式適合法（Savitzky-Golay法）によりデータ微分処理を行う場合、連続した複数のデータを必要とし、そのデータ期間の中心の時間での微分値を導出するのが、一般的である。このため、逐次データ処理の場合、少なくとも、最新データの時間から上記の中心のデータまでの時間差だけ原理的に時間遅延が発生するという課題がある。非特許文献４に記載の多項式適合法（Savitzky-Golay法）により最新のデータを取得した時点で微分値を導出した場合、通常の多項式適合法（Savitzky-Golay法）によりデータ中心点で微分値を導出した場合よりもデータ遅延時間は減少するが、その一方で、Ｓ/Ｎ比は減少する。
【００１５】
このため、使用するデータ数を増加すると、Ｓ/Ｎ比は増加して改善する。しかし、その一方で、多くのデータ数を使用することにより急激な変化に対する追随性（応答性）が劣化して使用するデータ数により応答性に限界が生じるという課題がある。具体的には、真のシグナルに含まれていた高周波成分が、データ処理後のシグナルでにおいて欠落するという課題である。
【００１６】
Savitzky-Golay法は、複数のデータを２次曲線あるいは３次曲線にフィッティングするものである。このため、データ数を増加させると、フィッティング精度は向上する。しかし、データ数の期間により２次あるい３次曲線以上の変化には対応できないことが上記の高周波成分が欠落する理由である。
【００１７】
さらにプラズマエッチング処理装置でのプラズマエッチングの終点判定方法として特許文献２に開示された方法があるが、この特許文献２に開示されたプラズマエッチングの終点判定方法に対して主に以下のような２つの課題がある。
【００１８】
現状の半導体デバイスは、高性能化および高集積化によりＦｉｎ  Ｆｉｅｌｄ  Ｅｆｆｅｃｔ  Ｔｒａｎｓｉｓｔｏｒ(Ｆｉｎ  ＦＥＴ)構造等の高段差構造を有している。また、通常のプラズマエッチングでは、疎部のパターンと密部のパターンとのエッチング性能差であるマイクロローディングが発生する。さらに被エッチング膜がウエハ面内で不均一の場合がある。
【００１９】
これらのことから例えば、プラズマエッチングの終点判定に用いるプラズマ発光からの分光強度信号データの時系列データが2段で変化する場合がある。このようにプラズマ発光からの分光強度信号データの時系列データが２段で変化して１番目あるいは２番目の変化によりプラズマエッチングの終点判定を行う場合、１番目の変化と２番目の変化が短い期間で発生するため、データ処理が２回の変化に追従できずにプラズマエッチングの終点を検出できない。尚、ここでは、プラズマ発光からの分光強度信号データの時系列データが変化した時点をプラズマエッチングの終点としている。
【００２０】
また、例えば、１次微分値あるいは２次微分値の遅延時間に対して短時間のプラズマエッチングにおいても１次微分値あるいは２次微分値を算出するデータ処理がプラズマ発光からの分光強度信号データの時系列データの変化に追従できない。つまり、プラズマエッチングの終点検出の応答性が不十分という課題がある。
【００２１】
次にプラズマエッチングのマスクパターンを大別すると、溝パターンと孔パターンがある。また、通常、孔パターンのウエハの開口率は、溝パターンのウエハの開口率より小さく、１％未満の開口率の場合もある。また、開口率が小さくなるに従ってプラズマ発光の分光強度の変化は小さくなる。
【００２２】
このため、例えば、１％未満の開口率のウエハの場合、プラズマ発光からの分光強度信号データの時系列データの変化が小さすぎるため、プラズマエッチングの終点検出が困難となる。つまり、低S/N比に対応できないという課題がある。尚、ここで開口率とはウエハ全体の面積に対する被エッチング面積の割合のことである。
【００２３】
非特許文献３には、入力データ（特に逐次入力データ）に応じて平滑化係数が変化する応答型１段指数平滑化法が開示されている。平滑化係数は、｜相対誤差/絶対誤差｜に応じて変化するケースとロジスティック関数に応じて変化するケースが記載されている。上記２ケースともに平滑化係数はデータ応答型だが、ベースは1段指数平滑化法なので、S/N比とデータ応答性を両立することが難しいという課題がある。
【００２４】
また、ロジスティック関数を使用するケースは、平滑化係数の可動範囲に適した単に出力が０〜１となることが主な理由と推定され、平滑化係数と誤差との関係が十分には配慮されていなかった。このため、データスムージング処理のＳ/Ｎ比性能が必ずしも十分ではないという課題がある。さらにロジスティック関数には、複数のパラメータがあり、パラメータの設定が難しくて使い勝手が良くないという課題がある。
【００２５】
上述のようにデータスムージング処理、データ微分処理において、Ｓ/Ｎ比性能向上とデータ応答性能向上（遅延時間短縮）は、トレードオフ関係がある。データスム-ジング処理方法およびデータ微分処理方法では、データ処理の対象である入力データに応じて各データ処理方法のパラメータを最適化する必要がある。
【００２６】
従来は、データ処理のパラメータを順次変更してデータ処理を行い、各データスムージング処理波形、データ微分処理波形等のデータ処理波形、およびＳ/Ｎ比、遅延時間等の数値データを総合的に俯瞰して最適なパラメータを見出していた。しかしながら、上述の方法では、最適なパラメータを見出すのに長時間を要し、また、最適なパラメータを見出す時間短縮等には、データ処理の知識や経験が必要等の課題があった。
【００２７】
また、非特許文献３には、1段指数平滑化法において、1期先予測値の誤差の総和を最小とすることにより最適な平滑化パラメータを推定する方法が開示されている。しかし、この方法では、予測値の曲線のなめらかさについて配慮されていないため、変化点を検出する１次微分処理、２次微分処理において、ノイズが大きく、かつ、S/N比が良くないという課題があった。
【００２８】
このため、本発明は、デ−タを処理するデータ処理装置を備えるプラズマ処理装置、デ−タを処理するデータ処理装置およびデータ処理方法において、高Ｓ／Ｎ比とデ−タ遅延時間の短縮やデータ応答性向上を両立できるプラズマ処理装置、データ処理装置およびデータ処理方法を提供する。また、入力データをデータスムージング処理、データ微分処理するのに最適なパラメータを簡便に自動で見出すことができるプラズマ処理装置、データ処理装置およびデータ処理方法を提供する。
【課題を解決するための手段】
【００２９】
本発明は、入力データと平滑化されたデータの予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いてデータを処理するデータ処理方法において、入力データと平滑化されたデータの予測値との誤差をパラメータとした、確率密度関数（正規分布型関数、ガウス分布型関数）を用いて、平滑化パラメータの応答性係数を導出することを特徴とする。
【００３０】
また、本発明は、入力データと平滑化されたデータの予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いてデータを処理するデータ処理方法において、入力データと平滑化されたデータの予測値との誤差を指数荷重移動平均し、該指数荷重移動平均の係数が、平滑化データの傾きに応じて変化する、応答型加算係数を用いることを特徴とする。
【００３１】
また、本発明は、２段指数平滑法を用いてデータを処理するデータ処理方法において、連続した複数の入力データを用いて、後退型Savitzky-Golay法でデータスムージング処理を行った結果を、２段指数平滑化法のデータの平滑化の式の入力データとして用い、連続した複数の該データの平滑化の予測値を用いて、後退型Savitzky-Golay法で１次微分処理した結果を、２段指数平滑化の平滑化データの傾きの式に、データの平滑化の予測値の傾きの入力データとして用いることを特徴とする。
【００３２】
また、本発明は、試料をプラズマ処理するプラズマ処理方法において、前記試料のプラズマ処理時の発光データを計測し、前記計測された発光データと平滑化された発光データの予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いるとともに、入力データとデータの平滑化の予測値との誤差をパラメータとした、確率密度関数（正規分布型関数、ガウス分布型関数）により、平滑化パラメータの応答性係数を導出すること、あるいは入力データとデータの平滑化の予測値との誤差を指数荷重移動平均し、該指数荷重移動平均の係数が、平滑化データの傾きの予測値に応じて変化する、応答型加算係数を用いること、あるいは連続した複数の入力データを用いて、後退型Savitzky-Golay法でデータスムージング処理を行った結果を、２段指数平滑化法のデータの平滑化の式の入力データとして用い、更に連続した複数の該データの平滑化の予測値を用いて、後退型Savitzky-Golay法で１次微分処理した結果を、２段指数平滑化の平滑化データの傾きの式に、データの平滑化の予測値の傾きの入力データとして用いることの少なくとも１つを用いた、前記２段指数平滑法を活用して求められた、１次微分処理結果、あるいは２次微分処理結果の中で、少なくとも一つの処理結果に基づいて前記試料のプラズマ処理状態あるいは前記試料のプラズマ処理状態の変化を検出し、検出結果に応じて前記試料のプラズマ処理を制御することを特徴とする。
【００３３】
また、本発明は、処理されるデータを受け取るデータ入出力装置と、前記データ入出力装置データが受け取った前記データを記憶する記憶装置と、入力データと入力データの平滑化の予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いるとともに、入力データとデータの平滑化の予測値との誤差をパラメータとした、確率密度関数（正規分布型関数、ガウス分布型関数）により、平滑化パラメータの応答性係数を導出すること、あるいは入力データとデータの平滑化の予測値との誤差を指数荷重移動平均し、該指数荷重移動平均の係数が、データの平滑化の傾きの予測値に応じて変化する、応答型加算係数を用いること、あるいは連続した複数の入力データを用いて、後退型Savitzky-Golay法でデータスムージング処理を行った結果を、２段指数平滑化法のデータの平滑化の式の入力データとして用い、更に連続した複数の該データの平滑化の予測値を用いて、後退型Savitzky-Golay法で１次微分処理した結果を、２段指数平滑化の平滑化データの傾きの式に、データの平滑化の予測値の傾きの入力データとして用いることの少なくとも１つを用いた、前記２段指数平滑法のデータ処理プログラムを記憶するデータ処理プログラム記憶装置と、前記データ処理プログラムに基づいてデータ処理を行うデータ演算処理装置とを備えたことを特徴とするデータ処理装置である。
【００３４】
また、本発明は、制御対象となる処理室と、前記処理室に関連するデータを取得する計測装置と、入力データと入力データの平滑化の予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いるとともに、入力データとデータの平滑化の予測値との誤差をパラメータとした、確率密度関数（正規分布型関数、ガウス分布型関数）により、平滑化パラメータの応答性係数を導出すること、あるいは入力データとデータの平滑化の予測値との誤差を指数荷重移動平均し、該指数荷重移動平均の係数が、データの平滑化の傾きの予測値に応じて変化する、応答型加算係数を用いること、あるいは連続した複数の入力データを用いて、後退型Savitzky-Golay法でデータスムージング処理を行った結果を、２段指数平滑化法のデータの平滑化の式の入力データとして用い、更に連続した複数の該データの平滑化の予測値を用いて、後退型Savitzky-Golay法で１次微分処理した結果を、２段指数平滑化の平滑化データの傾きの式に、データの平滑化の予測値の傾きの入力データとして用いることの少なくとも１つを用いた、前記２段指数平滑法のデータ処理プログラムに基づいてデータ処理を行うデータ演算処理装置と、前記データ演算処理装置により求められたデータ処理結果を利用して、前記処理室の状態あるいは前記処理室の状態の変化を検出し、検出結果に応じて前記処理室を制御する制御装置とを備えることを特徴とする処理装置である。
【００３５】
さらに、本発明は、試料台に載置された試料がプラズマ処理される処理室と、前記試料のプラズマ処理時の発光データを取得する計測装置と、前記計測装置により取得されたデータとデータの平滑化の予測値との誤差に応じて平滑化パラメータを変化させる応答型２段指数平滑法を用いるとともに、入力データとデータの平滑化の予測値との誤差をパラメータとした、確率密度関数（正規分布型関数、ガウス分布型関数）により、平滑化パラメータの応答性係数を導出すること、あるいは入力データと平滑化されたデータの予測値との誤差を指数荷重移動平均し、該指数荷重移動平均の係数が、平滑化データの傾きの予測値に応じて変化する、応答型加算係数を用いること、あるいは連続した複数の入力データを用いて、後退型Savitzky-Golay法でデータスムージング処理を行った結果を、２段指数平滑化法のデータの平滑化の式の入力データとして用い、更に連続した複数の該データの平滑化の予測値を用いて、後退型Savitzky-Golay法で１次微分処理した結果を、２段指数平滑化の平滑化データの傾きの式に、データの平滑化の予測値の傾きの入力データとして用いることの少なくとも１つを用いた、前記２段指数平滑法のデータ処理プログラムに基づいてデータ処理を行うデータ演算処理装置と、前記データ演算処理装置により求められたデータ処理結果を利用して、前記試料のプラズマ処理状態あるいは前記試料のプラズマ処理状態の変化を検出し、検出結果に応じて前記処理室を制御する制御装置とを備えることを特徴とする処理装置である。
【００３６】
さらに本発明は、離散化され、連続した、複数の数値データを入力データとし、データスムージング処理、あるいはデータ微分処理を行い、少なくとも１つのパラメータを有する、データ処理方法において、入力データとデータスムージング処理結果との平均二乗誤差と、データスムージング処理結果の２次微分の二乗平均値に、任意の数値の係数を乗じたものを、加算した数値を、評価関数をとし、該評価関数を最小とすることを指標し、該データ処理のパラメータを導出することを特徴とする、データ処理方法である。
【発明の効果】
【００３７】
本発明によれば、デ−タを処理するデータ処理装置を備えるプラズマ処理装置、デ−タを処理するデータ処理装置およびデータ処理方法において、データ処理結果のＳ／Ｎ比向上とデータ処理結果の応答性向上（デ−タ処理結果の遅延時間低減）を両立することができる。また、データ処理方法の最適パラメータを簡便かつ短時間に自動で見出すことができる。これにより、ユーザビリティが良い、データ処理装置、データ処理方法および処理室を制御する制御装置とを備えることを特徴とする処理装置を提供することができる。
【図面の簡単な説明】
【００３８】
【図１】実施例１に係るデ−タ処理装置１の全体構成を示す図である。
【図２】実施例１に係るデータ処理のフローを示す図である。
【図３Ａ】乱数データにおいて、平均値からの予測値のズレ量の定義を示す図である。
【図３Ｂ】図３Ａのデータにおいて、ズレ量/標準偏差σと応答係数Ftとの関係を示す図である。
【図４】データ処理の入力となる信号出力の傾きの変化例を示す図である。
【図５】実施例１に係る有磁場マイクロ波プラズマエッチング装置の縦断面図である。
【図６】エッチング終点検出のための代表的なデータ処理のフロー図である。
【図７】実施例１に係るデータ処理装置１において図２のデータ処理フローを用いて式（33）〜式（42）の応答型２段指数平滑化法(タイプIV)を用いた場合のデータ処理結果である。
【図８】式（5）〜式（14）の応答型２段指数平滑化法(タイプI)にてＮ＝0とした場合のデータ処理結果である。
【図９】式（5）〜式（14）の応答型２段指数平滑化法(タイプI)にてＮ＝１とした場合のデータ処理結果である。
【図１０】式（5）〜式（14）の応答型２段指数平滑化法(タイプI)にてＮ＝５とした場合のデータ処理結果である。
【図１１】タイプIの応答型２段指数平滑化法（Ｎ=0）の１つのパラメータを順次変化させた場合における、入力データ(Y1t)とデータスムージング処理結果(S1t)を示す（Ａ）と、２次微分処理結果（B2t）を示す(Ｂ)とからなる図である。
【図１２】図１１の場合において、平均二乗誤差であるＥ、２次微分の二乗平均であるＤおよび評価関数であるＷのパラメータ依存性を示した図である。
【図１３】図１２における、評価関数Ｗを用いて自動で導出した×印の位置のパラメータを使用した場合の、入力データ(Y1t)とデータスムージング処理結果(S1t)を示す（Ａ）と、２次微分処理結果（B2t）を示す(Ｂ)とからなる図である。
【発明を実施するための形態】
【００３９】
以下、図面を参照しながら本発明の実施例を説明する。
【実施例１】
【００４０】
まず、図１ないし図１０を用いて本発明の第1の実施例に係るデータ処理装置を説明する。ここでは、有磁場マイクロ波プラズマエッチング装置において、高精度のエッチング処理を目的としたプラズマ分光によるエッチング終点検出に本発明を適用した例について説明する。
【００４１】
図１に、第１の実施例のデ−タ処理装置１の構成図を示す。本実施例では、デ−タ処理装置１は、デ−タ入出力装置２、記憶装置であるデータ記憶装置３、データ処理プログラム記憶装置４、データ演算処理装置５で構成され、相互にデータ移動が可能なように接続されている。
【００４２】
また、必要に応じて上記に加え、データ表示装置（図示省略）が設けられる。データ処理装置１は、対象とする系６(装置、分析データ等)とデータを入出力することができる。これにより、対象とする系６を高精度に制御する。本実施例の場合、対象とする系６は、マイクロ波プラズマ処理装置である。また、データ処理装置１は、単独で使用しても良く、デ−タ解析等に利用することができる。
【００４３】
データ入出力装置２は、処理データ、データ処理プログラムのパラメータ等を入出力することができる。データ入出力装置２は、対象とする系６等より一括あるいは逐次に処理するデータを受け取る。データ入出力装置２により受け取られたデータは、RAM等のデータ記憶装置３で記憶され、RAM等のデータ処理プログラム記憶装置４に記憶されたデータ処理プログラムに従ってデータ演算処理装置５によってデータスムージング処理、データ微分処理を行われる。データ演算処理装置５によるデータの演算後、データ入出力装置２によってデータスムージング処理結果データ、データ微分処理結果データを対象とする系６等に出力して対象とする系６の制御に利用する。
【００４４】
データ処理プログラムに記憶されているデータ処理方法の全体フロー図を図２に示す。一括あるいは逐次に処理するデータを受け取り、デ−タが入力される。入力データは、時系列データY1t：t=1、2、・・・とした。次に後述する方法により初期値を導出する。次に1回目の２段指数平滑化処理を行い、１回目出力のデータの平滑化の予測値S1tと平滑化されたデータの傾きの予測値B1tを得る。
【００４５】
この場合、下記式(5)〜式(14)に示されたタイプIの応答型２段指数平滑化処理法で処理した。
【００４６】
データの平滑化    ：S1t = α1t Y1t + (1 - α1t) (S1t-1 + B1t-1)         式(5)
平滑化されたデータの傾き：B1t= γ1t (S1t - S1t-1) + (1 -γ1t) B1t-1     式(6)
平滑化係数： α1t =  (Kα-Lα) Fα+ Lα                               式(7)
応答係数：    Fαt = (｜δαt/Δαt｜+φ)N                            式(8)
相対誤差：    δαt = A1（Y1t - S1t） + (1 - A1) δαt-1              式(9)
絶対誤差：    Δαt = A1 |Y1t - S1t| + (1 - A1) Δαt-1 + φ          式(10)
平滑化係数： γ1t = (Kγ-Lγ) Fγ+ Lγ                                式(11)
応答係数：   Fγt = (｜δγt/Δγt｜+φ)N                           式(12)
相対誤差：     δγt =A2{（S1t - S1t-1）- B1t}+(1 - A2)δγt-1      式(13)
絶対誤差：     Δγt =A2|（S1t - S1t-1）- B1t|+(1 - A2) Δγt-1+φ  式(14)
ここで、入力データは、例えば、時系列データY1t：t=1、2、・・・とし、出力のデータの平滑化の予測値S1tと平滑化されたデータの傾きの予測値B1tを逐次データ処理により得ることができる。また、Kα、Lα、Kγ、Lγ、N、A1、A2、φは、任意の定数である。但し、1＞ Kα＞ Lα＞0、1＞Kγ ＞ Lγ＞0、1＞A1＞0、1＞A2＞0である。φは、絶対誤差ΔαtおよびΔγt、または、応答係数FαtおよびFγtがゼロ値となるのを回避するためにある。また、通常の演算に影響が非常に少なくなるようにφとして極微小値が選択される。
【００４７】
式(8)、式(12)においてN=0の場合、応答係数Fα=1、Fγ=1となり、平滑化係数αt＝Kα、γt＝Kγとなって一定値となる。従って、N=0の場合は、通常の２段指数平滑化法となる。N=1の場合、各平滑化係数は、各相対誤差/絶対誤差に比例して、Kα〜Lα、Kγ〜Lγの範囲を変化する応答型２段指数平滑化法となる。同様にN=5の場合、各平滑化係数は、各相対誤差/絶対誤差の5乗に比例して、Kα〜Lα、Kγ〜Lγの範囲を変化する応答型２段指数平滑化法となる。
【００４８】
応答型２段指数平滑化法は、各相対誤差/絶対誤差が小さい場合に各平滑化係数を小さくしてデータスムージング性能を向上させる。一方、各相対誤差/絶対誤差が大きい場合は、各平滑化係数を大きくしてデータ応答性能を向上させる。データスムージング性能とデータ応答性能は、トレードオフ関係にあり、両者のバランスは、上記N値で変化するため、入力データの特性に応じて、上記N値を選択する。
【００４９】
次に前述の初期値の導出方法を説明する。一般に、２段指数平滑化処理における、データの平滑化の予測値S1の初期値および平滑化されたデータの傾きの予測値B1の初期値は、例えば、次の方法で導出される。データの平滑化の予測値S1の初期値は、S11＝入力データY11(方法A1)、あるいはS11＝初期N個の入力データの平均値({Y11+Y12+・・・+Y1N}/N)等(方法A2)で導出される。
【００５０】
平滑化されたデータの傾きの予測値B1の初期値は、B11＝Y12-Y11(方法B1)、B21＝｛（Y12-Y11）+(Y14-Y13)｝/2等(方法B2)で導出される。一般に２段指数平滑化処理は、データ処理の開始直後において誤差が大きいという課題があった。その原因の１つは、上記従来の導出方法による初期値が真の初期のデータの平滑化の予測値、真の初期の平滑化されたデータの傾きの予測値との誤差が大きいことにあった。
【００５１】
ここでは、所望のN個のデータ入力開始後初期のデータY1t(t=1、2、・・・、N)を用いて最小２乗法により多項式近似式を導出する。また、初期10個の等時間の間隔の時系列データを使用した。上記で導出された多項式近似式より入力データの直前ｔ＝0の仮想データであるデータの平滑化の予測値S10、平滑化されたデータの傾きの予測値B10を導出する。
【００５２】
ここでは、多項式近似式として直線1次式を使用し、データの平滑化の予測値S10、平滑化されたデータの傾きの予測値B10は、各々、式(15)、式(16)により導出される。
データの平滑化の予測値S10 =  { 330 Y11 + 275 Y12 + 220 Y13 + 165 Y14 
+ 110 Y15 + 55 Y16 + 0 Y17 - 55 Y18 - 110 Y19 
- 165 Y110 ｝ / 825              式(15)
平滑化されたデータの傾きの予測値B10 = { -45 Y11 - 35 Y12 - 25 Y13 - 15 Y14 
- 5 Y15 + 5 Y16 + 15 Y17 + 25 Y18 + 35 Y19 
+  45 Y110 } / 825             式(16)
また、図２に記載された２回目の２段指数平滑化処理における、データの平滑化の予測値の初期値S21と平滑化されたデータの傾きの予測値の初期値B21は、各々、S21＝S11、B21＝0に設定している。上記の初期値設定方法は、データ入力開始後初期の入力データも誤差が少なく、高精度にデータスムージング処理、データ微分処理できるという効果がある。
【００５３】
ここでは、多項式近似式により入力データの直前ｔ＝0の仮想データであるデータの平滑化の予測値S10、平滑化されたデータの傾きの予測値B10を導出したが、t＝1等、任意の時間の多項式近似式による各予測値を各初期値として使用しても良い。但し、この場合、t＝0の仮想データを用いた場合よりも短時間ステップのデータスムージング処理、データ微分処理等には不利となる。
【００５４】
上述のように1回目の２段指数平滑化処理を行い、１回目出力のデータの平滑化の予測値S1tと平滑化されたデータの傾きの予測値B1tを得る。次に、1回目出力の平滑化されたデータの傾きの予測値B1tを2回目の入力データY2tとして下記式(17)及び式(18)による２回目の２段指数平滑化処理を行うことによって２回目出力のデータの平滑化の予測値S2tと平滑化されたデータの傾きの予測値B2tを得る。
【００５５】
データの平滑化    ：S2t = α2 Y2t + (1 - α2) (S2t-1 + B2t-1)         式(17)
平滑化されたデータの傾き：B2t = γ2 (S2t - S2t-1) + (1 -γ2) B2t-1    式(18)
次にデータスムージング処理結果データS1t、デ−タ１次微分処理結果データS2t、データ２次微分処理結果データB2tは、一括あるいは逐次データとして出力される。ここで、２回目の２段指数平滑化における、データの平滑化の平滑化パラメータα2、平滑化され
たデータの傾きの平滑化パラメータγ2は、予め任意の定数に設定されている。但し、0＜α２＜1、0＜γ２＜1とする。1回目出力の平滑化されたデータの傾きの予測値B1tも１次微分処理結果に相当するため、これを使用しても良いが、データ結果のばらつきが大きいため、２回目の２段指数平滑化処理により、デ−タスム−ジング処理を行っている。
【００５６】
上述のタイプIの応答型２段指数平滑化法および２段指数平滑化処理を２回実施する方法によりデータスムージング処理結果、データ１次微分処理結果、データ２次微分処理結果等における、Ｓ/Ｎ比の増加と遅延時間の短縮を両立を大きく改善することができた。しかしながら、対象とする系の変化点を１次微分データや２次微分データにより検出し、対象とする装置を更に実時間に近く、更に高精度で制御する場合には、更に時間遅延の短縮や応答性の向上、更に高S/N比の両立等のデータ処理性能を向上を図ることにより対象とする装置の制御精度を向上することができる。
【００５７】
このため、図２の１回目の２段指数平滑化処理において、本実施例のタイプIIの場合は、下記式(19)〜式(24)を用いてデータ処理を行った。
データの平滑化    ：S1t = α1t Y1t + (1 - α1t) (S1t-1 + B1t-1)           式(19)
平滑化されたデータの傾き：B1t  = α1t (S1t - S1t-1) + (1 -α1t) B1t-1     式(20)
平滑化係数： α1t =  (K - L) Ft+ L                                   式(21)
応答係数：    Ft = 1 - Exp[-δαt2/(2σt2)]                          式(22)
相対誤差：    δαt = A1（Yt - S1t） + (1 - A1) δαt-1              式(23)
予測誤差分散： σt2 = A1（Yt - S1t）2 + (1 - A1)σt-12               式(24)
ここで、K、L、A1は、任意の定数である。但し、1＞ K＞ L＞0、1＞A1＞0である。また、σtの初期値は、上述の初期値導出方法と同様に入力データにおいて、初期数点のデータを用いた多項式近似の結果と入力データとの誤差の標準偏差より算出した。また、応答係数Ftは、１から確率密度関数（正規分布型関数、ガウス分布型関数）を減算したものである。一般に、誤差は、正規分布となるケースが多い。このため、正規分布を表現する確率密度関数を利用した応答係数Ftとした。これにより誤差を最小化するのに適した平滑化係数をデータ変化に応じて設定することができる。
【００５８】
図３Ａに示すように例えば、平均０、標準偏差σ＝１の正規分布に従う乱数データにおいて、予測値が最もおこり易い値、平均値からのズレ量を定義する。図３Ｂには、前記ズレ量/標準偏差σと応答係数Ftとの関係を示す。ここで、N=0、N=1、N=5は、式(5)から式（14）に示したタイプIの応答型２段指数平滑化法における応答係数である。1- PDFは、式(19)〜式(24)に示したタイプIIの確率密度関数を用いた応答型２段指数平滑化法における応答係数である。
【００５９】
ここで、PDFは、確率密度関数（Probability Density Function）を示す。1-PDFは、Ｎ＝5と類似の応答係数の特性を示すが、1-PDFは、N=5と比較してズレ量/標準偏差σ＜約1の場合は、応答係数は大きく、ズレ量/標準偏差σ＞約1の場合は、応答係数は小さい。この応答特性によりデータスムージング性能（Ｓ/Ｎ比）とデータ応答性（遅延時間の低減）の両立できる、データスムージング処理、データ微分処理ができるという効果がある。
【００６０】
上述の１回目の２段指数平滑化処理において、本実施例のタイプIIIの場合、下記式(25)〜式(32)を用いてデータ処理を行った。
データの平滑化    ：S1t = α1t Y1t + (1 - α1t) (S1t-1 + B1t-1)           式(25)
平滑化されたデータの傾き：B1t  = α1t (S1t - S1t-1) + (1 -α1t) B1t-1   式(26)
平滑化係数： α1t =  (K - L) Ft+ L                                    式(27)
応答係数：    Ft = 1 - Exp[-δαt2/(2σt2)]                           式(28)
相対誤差：    δαt = A1t（Y1t - S1t） + (1 - A1t) δαt-1            式(29)
予測誤差分散： σt2 = A1t（Y1t - S1t）2 + (1 - A1t)σt-12             式(30)
応答型加算係数： A1t = MAX(A1、βtA1max)                            式(31)
傾き係数      ： βt = 1 - Exp[-Bt2/(2 NN C)2]                        式(32)
ここで、A1maxは、加算係数の上限値（0＜A1max＜１）であり、NNは感度係数で任意の定数である。またCは、上述の初期値導出方法と同様に入力データにおいて、初期数点のデータを用いた１次近似式の結果より算出した傾きである。タイプIIIは、タイプIIに応答型加算係数を導入したものであり、タイプIIIの式(25)〜式(30)は、タイプIIの式(19)〜式(24)とほぼ同じである。入力となる信号出力の傾きの変化の例を図４に示す。この場合には、初期値の算出期間での初期傾きＣを基準として傾きが大きくなり、再び初期と同程度の傾きに変化する。
【００６１】
本実施例では、傾きの変化に応じて加算係数を変化させる。式（29）、式（30）より、加算係数A1tは、指数荷重平均処理の係数に相当する。従って、初期傾きに応じて指数荷重平均処理の係数を変化させる。つまり、初期傾きCと比較し、平滑化データの傾きB1tが大きくなった場合には、応答型加算係数A1tを小さくし、実効的に最新のデータの係数の比率が大きくなるようにして応答性を向上させる。応答型加算係数A1tは、まず式(32)により傾き係数βtを算出する。
【００６２】
傾き係数βtは、前述の応答係数Ftと同様に確率密度関数(正規分布関数、ガウス型関数)を用いて算出している。応答型加算係数A1tは、式（31）により通常時に設定値A1を使用し、平滑化データの傾きが大きくなり、設定値A1よりもβtA1maxが大きくなった場合、この値を加算係数として使用する。但し、過度に応答性が良くなることを抑制するため、上限値A1maxを設定できる。
【００６３】
応答型加算係数A1tの導入により平滑化されたデータの傾きが初期傾きと同程度の傾きの場合、加算係数を設定値に保つ。データの状態の変化を示す平滑化されたデータの傾きが大きくなった場合には、応答型加算係数A1tを大きくすることにより最新のデータの比重が大きくなるため、データ処理の応答性が更に向上するという効果がある。
【００６４】
上述の１回目の２段指数平滑化処理において、本実施例のタイプIVの場合は、下記式(33)〜式(42)を用いてデータ処理を行った。
データの平滑化    ：S1t = α1t SGB0Dt + (1 - α1t) (S1t-1 + B1t-1)        式(33)
平滑化されたデータの傾き：B1t  = α1t SGB1Dt + (1 -α1t) B1t-1          式(34)
SGB0Dt = (83Y1t + 54Y1t-1 + 30Y1t-2 + 11Y1t-3 - 3Y1t-4 - 12Y1t-5 -16Y1t-6 
- 15Y1t-7 - 9Y1t-8 + 2Y1t-9 + 18Y1t-10 ) / 143              式(35)
SGB1Dt = (945S1t + 456S1t-1 + 67S1t-2 - 222S1t-3 - 411S1t-4 - 500S1t-5 
- 489S1t-6 - 378S1t-7 - 167S1t-8 + 144S1t-9 + 555S1t-10) / 4290  
式(36)
平滑化係数： α1t =  (K - L) Ft+ L                                    式(37)
応答係数：    Ft = 1 - Exp[-δαt2/(2σt2)]                           式(38)
相対誤差：    δαt = A1t（Yt - S1t） + (1 - A1t) δαt-1             式(39)
予測誤差分散： σt2 = A1t（Yt - S1t）2 + (1 - A1t)σt-12              式(40)
応答型加算係数： A1t = MAX(A1、βtA1max)                            式(41)
傾き係数      ： βt = 1 - Exp[-Bt2/(2 NN C)2]                        式(42)
ここで、SGB0Dtは、非特許文献３に記載の後退型Savitzky-Golay法により、連続した１１項の入力データY1t〜Y1t-10を使用すると共に最新のデータY1tの時点におけるデータスムージング処理結果である。同様にSGB1Dtは、後退型Savitzky-Golay法により、連続した１１項のデータの平滑化の予測値S1t〜S1t-10を使用すると共に最新のデータの平滑化の予測値S1tの時点における1次微分処理結果である。但し、非特許文献３には、一部誤記等があると思われ、修正し使用した。
【００６５】
式(33)、式(34)、式(37)〜式(42)は、前述のタイプIIIの式(25)〜式(32)と概ね同じである。２段指数平滑化法およびタイプIの応答型２段指数平滑化法の基本部分である式(5)、式(6)の意味は、各々、下記式(43)、式(44)である。
データの平滑化    ：S1t = α1t (入力データ)t + (1 - α1t) (S1t-1 + B1t-1)
式(43)平滑化されたデータの傾き：B1t=γ1t(データの平滑化の予測値の傾き)t+(1 -γ1t)B1t-1
式(44)
タイプIVでは、式(43)の入力データの替わりに前処理として後退型Savitzky-Golay法を用いてデータスムージング処理した結果を用いている。また、式(44)の平滑化されたデータの予測値の傾きの替わりに前処理として後退型Savitzky-Golay法を用いて1次微分処理した結果を用いている。通常、多項式適合法(Savitzky-Golay法）は、複数の連続したデータを使用し、この連続データの期間の中心点でのデータ処理結果を導出する場合が多い。ここでは、これを中心型Savitzky-Golay法とする。これに対し、最新データの時点におけるデータ処理結果を導出する場合を後退型Savitzky-Golay法とする。
【００６６】
中心型Savitzky-Golay法の場合、連続データの中心点でのデータ処理結果を導出するため、データ処理に遅延時間が生じることは避けられない。Savitzky-Golay法は、使用するデータ数が多いほどS/N比は向上するが、中心型では、遅延時間が使用するデータ数とともに増加する。後退型Savitzky-Golay法の場合、最新データの時点におけるデータ処理結果を導出するので、データ処理の遅延はほとんどない。しかし、後退型は中心型よりも、S/N比は低下する。
【００６７】
検討の結果、中心型でデータ数５ケの場合と後退型でデータ数１１ケの場合を比較した結果、ほぼ同程度のS/N比が得られると分かった。Savitzky-Golay法は、複数データを用
いて２次あるいは３次曲線に多項式近似するものであるため、データ数が増加すると、使用したデータ数の期間において３次以上の変化が生じる場合には対応できない。即ち、使用するデータ数が増加すると、高周波成分が欠落するリスクが生じる。
【００６８】
以上を鑑みてタイプIVでは、式(33)の入力データに後退型Savitzky-Golay法データ数11ケのデータスムージング処理を使用し、式(34)のデータの平滑化の予測値の傾きに後退型Savitzky-Golay法データ数11ケのデータ１次微分処理値を使用した。このため、本実施例によれば、データ遅延時間を増加させず、また、データ処理の周波数特性を劣化を少なく抑えてデータスムージング処理、データ微分処理のS/N比を向上させる効果がある。タイプIVでは、後退型Savitzky-Golay法データ数11ケを使用したが、要求されるデータ処理性能に応じて他のデータ数を使用したり、あるいは中心型と後退型の中間等の時点でデータ処理を行うSavitzky-Golay法を使用しても良い。
【００６９】
更にSavitzky-Golay法以外の前処理として入力データをデータスムージング処理した結果を式(43)の入力データとし、また、データの平滑化の予測値をデータ１次微分処理した結果を式(44)のデータの平滑化の予測値の傾きに使用しても良い。
【００７０】
図５は、本発明の第１の実施例に係る有磁場マイクロ波プラズマエッチング装置の縦断面図を示す。本実施例では、有磁場マイクロ波プラズマエッチング装置が図１の対象とする系６(装置、分析データ等)に相当する。容器７、放電管８、石英板９及び石英窓１０で区画された処理室１１の内部を排気用開閉バルブ１２を開とし、真空排気装置１３により減圧する。エッチングガスは、マスフローコントロ−ラ（図示省略）を介してガス配管１４を通り、石英板９と石英シャワープレート１５の間を通過して石英シャワープレート１５のガス孔から処理室１１に導入される。処理室１１に導入されたエッチングガスは、排気速度可変バルブ１６により処理室１１内の圧力を所望の圧力に調整される。
【００７１】
また、処理室１１は、コイル１７、１８とヨーク１９により生成される磁場領域内にある。マグネトロン２０から発振された、この場合周波数2.45GHzのマイクロ波は、アイソレータ（図示省略）、パワーモニタ（図示省略）、整合器２１を経由して矩形導波管２２
内を矩形TE10モードで伝播し、円矩形変換器２３を経由して円形導波管２４内を円形TE11モードで伝播する。その後、マイクロ波は、空洞共振器２５に導入され、石英板９、石英シャワープレート１５を透過して処理室１１内に入射される。処理室１１内には、導入される2.45GHzのマイクロ波と電子サイクロトロン共鳴を生じる磁束密度875Gaussの磁場領域が、処理室１１の中心軸およびマイクロ波の導入方向に対し垂直に、また、処理室１１の中心軸に対する断面方向に対して全面に形成されている。
【００７２】
この2.45GHzのマイクロ波と875Gaussの磁場との相互作用により主に生成されたプラズマにより試料台であるウェハ載置用電極２６に載置されたウェハ２７がエッチング処理される。また、試料であるウェハ２７のエッチング形状を制御するため、ウェハ載置用電極２６には整合器（図示省略）を介して高周波電源２８が接続され、高周波電圧を印加することが可能になっている。また、ウェハ載置用電極２６には、チラーユニット（図示省略）が接続され、ウェハ２７の温度を制御することができる。
【００７３】
処理室１１、ウェハ２７、ウェハ載置用電極２６は、各々同軸に配置されている。また、エッチングガスを導入する石英シャワープレート１５のガス孔領域、真空排気部である排気用開閉バルブ１２、排気速度可変バルブ１６、真空排気装置１３も処理室１１に対し各々同軸に配置されている。このため、ウェハ２７上でのガス流れは同軸対称である。磁場を生成するコイル１７、１８、ヨーク１９も処理室１１に対して同軸に配置されているため、処理室１１内の磁場プロファイル、磁束875Gaussの電子サイクロトロン共鳴領域は処理室１１に対し同軸に形成される。また、円形導波管２４、空洞共振器２５も処理室１１に対し同軸に配置されているため、処理室１１に導入されるマイクロ波も処理室１１に対し同軸に導入される。
【００７４】
磁場が処理室１１に対し同軸に生成され、マイクロ波も処理室１１に対し同軸に導入されるため、磁場とマイクロ波との相互作用によって形成されるプラズマは、処理室１１に対し同軸に生成される。これにより、プラズマ中の電子やイオンは、ウェハ２７に対し同軸に輸送される。また、エッチングガスの流れも処理室１１に対し同軸であるため、プラズマにより生成されたラジカルやウェハ２７のエッチングによる反応生成物もウェハ２７に対し同軸に導入、排気される。従って、エッチングレートや材料選択比やエッチング形状等のエッチングプロセス処理性能をウェハ面内で均一性良くエッチング処理することができる。
【００７５】
処理室１１で生成されたプラズマからであり処理室１１側方からの発光は、石英窓１０及び光ファイバー２９を通過して分光器３０に導入され、光強度の波長依存性の時系列データとして出力される。また、処理室１１上方からのプラズマからの発光は、石英シャワープレート１５、石英板９、空洞共振器２５、円形導波管２４、円矩形変換器２３及び光ファイバー３１を通過して分光器３２に導入され、光強度の波長依存性の時系列データとして出力される。
【００７６】
処理室１１には、エッチングガスおよびウェハ２７からのエッチング反応生成物が導入され、これらがマイクロ波と磁場の相互作用により、解離し、プラスマを生成する。このため、処理室１１で生成されたプラズマからの発光には、エッチングガスおよびエッチング反応生成物の構成される原子、分子、ラジカルおよびこれらの反応物の情報が含まれている。
【００７７】
例えば、パターンマスクの下方にpoly-Si膜とSiO2膜が配置されたSi基板の代表的なpoly-Siエッチングでは、下地SiO2と高選択比でPoly-Siエッチングすることが要求される。エッチングガスは、ハロゲン系ガスが使用され、エッチング反応生成物には、被エッチング材料であるSiとハロゲンが含まれる。エッチング反応生成物は、プラズマにより再解離されるため、プラズマからのSiに起因した波長288nmの発光の光強度を分光器３０あるいは分光器３２でモニタする。
【００７８】
この場合、poly-Si膜のエッチングが終了し、下地SiO2が出現した場合、下地SiO2のエッチレートは小さいため、Siに起因した波長288nmのプラズマ発光強度が急激に減少し、遂には一定値に近づく。このプラズマ発光の変化をモニタしてエッチング処理の終点検出を行う。
【００７９】
処理室１１側方からのプラズマの発光には、エッチングガスとエッチング反応生成物の情報が含まれるが、処理室１１上方からのプラズマの発光には、上記情報に加え、プラズマ光がウェハ２７の膜構造および段差構造で干渉を生じるため、ウェハ２７の膜構造および段差構造の情報も含んでいる。このプラズマの発光データを解析することによりエッチング中の膜厚、エッチング深さなどをモニタすることができる。本実施例では、簡単のために処理室１１側方からのプラズマの発光データをエッチング終点モニタに使用した。  エッチング終点検出のための代表的なデータ処理フローを図６に示す。入力デ−タY1tは、エッチング中のプラズマ発光強度の変化を模擬した式(45)の評価関数で作成した。
【００８０】
Y1t ＝H/[1＋exp{ -A(t - T)}]＋C t＋D＋F (R - 0.5)    式(45)
ここで、H、A、T、C、D、Fは、任意の定数、Rは0〜1の乱数である。上記評価関数を使用すれば、データスムージング処理、１次微分処理、２次微分処理の解析的な真値が既知であるため、種々のデータ処理方法において、真値との絶対誤差、データ処理に伴なう遅延時間、S/N比（シグナル/ノイズ比）等のデータ処理性能を比較、評価することができる。
【００８１】
図６に示すようにエッチング終点のための代表的なデータ処理フローは、図６(A)に示すような入力データ波形を図６(B)に示すようにデータスムージング処理を行った後、図６(C)に示すように１次微分処理、図６(D)に示すように２次微分処理を行う。ノイズが多く含まれる入力データは、データスムージング処理により変化点が明瞭となる。この変化点を１次微分処理ではピーク値の点（時間）として２次微分処理ではゼロクロスする点(時間)として検出する。ゼロクロスする点(時間)を基準にエッチング終点を判断し、エッチング装置を制御し、高精度のエッチング処理を行う。
【００８２】
変化点は、１次微分処理のピーク、２次微分処理のゼロクロスによって順次、より明確により簡便に判断できるが、シグナル強度の絶対値は、順次小さくなる。このため、S/N比の高いデータ処理が重要となる。特に、被エッチング面積が小さい低開口率のマスクパターンのエッチングの場合は、エッチング終点前後のプラズマ発光強度の変化が小さいため、更に高S/N比のデータ処理が必要となる。
【００８３】
一般にデータスムージング処理、データ微分処理では、S/N比が高くなると、遅延時間が長くなり、真値との絶対誤差が大きくなる。つまり、S/N比と遅延時間、絶対誤差は、トレードオフ関係となっており、S/N比、遅延時間、絶対誤差を同時に満足する、データスムージング処理、データ微分処理が必要である。
【００８４】
本実施例では、図２のデータスムージング処理、データ微分処理フローを用いて図６のデータ処理を行う。また、図６の入力データは、図５の有磁場マイクロ波プラズマエッチング装置におけるエッチング中のプラズマ発光をモニタする分光器３０からの出力データに相当する。図５に示すように有磁場マイクロ波エッチング装置をシステムとして制御するシステム制御装置３３（データ入出力装置、データ処理装置、データ表示装置等を含む）と本実施例のデータ処理装置１とが設けられている。データ処理装置１は、システム制御装置３３の一部として、組み込まれても良い。
【００８５】
分光器３０、分光器３２からの出力データは、データ処理装置１に伝送される。データスムージング処理結果、データ１次微分結果、データ２次微分結果は、制御装置であるシステム制御装置３３に伝送される。このデータスムージング処理結果、データ１次微分結果、データ２次微分結果に基づいてシステム制御装置３３はエッチング終点判定を行い、有磁場マイクロ波エッチング装置をシステムとして制御する。エッチング終点判定では、主にプラズマ生成を制御するため、図５において、システム制御装置３３は、マグネトロン２０、高周波電源２８への接続が図示されている。また、図５での図示は省略したが、システム制御装置３３は、システムを構成する他の機器にも接続されている。
【００８６】
図７は、本発明の実施例のデータ処理装置１において、図２のデータ処理フローを用い、式（33）〜式（42）のタイプIVの応答型２段指数平滑化法を用いた場合のデータ処理結果である。図８は、式（5）〜式（14）のタイプIの応答型２段指数平滑化法において、Ｎ＝0とした場合のデータ処理結果である。また、図９は、図８と同様でＮ＝１とした場合のデータ処理結果である。図１０は、図８と同様でＮ＝５とした場合のデータ処理結果である。
【００８７】
図７〜図１０において、（Ａ）は、入力データ（Y1t）とデータスムジング処理結果（S1t）を示し、（Ｂ）は、１次微分処理結果（B1t）を示し、（Ｃ）は、１次微分スムージング処理結果（S2t）を示し、（Ｄ）は、２次微分処理結果（B2t）を示す。また、図７〜図１０の入力データは、同一のものを使用した。
【００８８】
図７〜図１０を比較すると、まず（Ａ）のデータスムージング処理において、本発明のタイプIVの実施例である図７（Ａ）は、タイプIの、図８（Ａ）、図９（Ａ）、図１０（Ａ）よりもデータスムージング処理結果が入力データの変化に追随しており、オーバーシュート等が少なく、さらに応答性が良く、誤差が小さいことが分かる。次に（Ｃ）の１次微分スムージング処理において、本発明のタイプIVの実施例である図７（Ｃ）は、タイプIの、図８（Ｃ）、図９（Ｃ）、図１０（Ｃ）よりも1次微分のピークの時間が早く、また波形の半値幅も小さい。このため、データ処理の遅延時間が小さく、応答性が良いことが分かる。
【００８９】
（Ｄ）の２次微分処理において、本発明のタイプIVの実施例である図７（Ｄ）は、タイプIの、図８（Ｄ）、図９（Ｄ）、図１０（Ｄ）よりも、２次微分のゼロクロスの時間が
早く、また、２次微分波形の全体の幅も小さい。さらに波形もなめらかで、２次微分の信号強度も大きい。このため、データ処理の遅延時間が小さく、応答性が良く、Ｓ/Ｎ比も良いことが分かる。数値で比較すると、本発明のタイプIVの実施例である図７（Ｄ）、タイプI（Ｎ＝0）の図８（Ｄ）、タイプI（Ｎ＝1）の図９（Ｄ）、タイプI（Ｎ＝５）の図１０（Ｄ）の順で各々の２次微分波形のS/N比は、1862、68、621、957である。
【００９０】
遅延時間は、２次微分波形のゼロクロス時間において、本発明の実施例である図７（Ｄ）は、タイプI（Ｎ＝0）の図８（Ｄ）よりも3.9秒早く、タイプI（Ｎ＝1）の図９（Ｄ）よりも1.6秒早く、タイプI（Ｎ＝５）の図１０（Ｄ）よりも0.9秒早く、遅延時間が短縮されている。このため、本発明の実施例では、遅延時間とS/N比がともに、タイプIよりも性能向上しており、遅延時間短縮（応答性向上）とS/N比向上を両立させる効果がある。
【００９１】
従って、エッチング終点判定の基準となる１次微分のピークの点（時間）、２次微分のゼロクロスの点（時間）を明瞭に検出することができることが分かる。このため、本発明の実施例によれば、絶対値誤差が小さく、S/N比が高く、遅延時間が短い、データスムージング処理、データ微分処理をリアルタイムで逐次処理できるという効果がある。
【００９２】
また、（Ｂ）１次微分処理において、本発明のタイプIVの実施例である図７（Ｂ）は、タイプIの、図８（Ｂ）、図９（Ｂ）、図１０（Ｂ）よりも波形がなめらかである。これは、後退型Savitzky-Golay法による前処理の効果である。このため、図２の処理フローにおいて、２回目の入力データとなる、１回目の出力データ（S1t）、（B1t）がなめらかな波形であるため、２回目の２段指数平滑化の平滑化係数を大きくすることができる。これにより、２回目の出力である、１次微分スムージング結果（S2t）と２次微分処理結果（B2t）の遅延時間を短縮し、応答性を向上させることができる。
【００９３】
また、データスムージング処理結果およびデータ１次微分処理結果のみが必要な場合には、本実施例によれば、後退型Savitzky-Golay法による前処理の効果により１回目の２段指数平滑化処理によって１次微分のなめらかな波形が得られるため、２回目の２段指数平滑化処理を実施せずに１回目の出力データ（S1t）、（B1t）を用いても良い。この場合、データ処理プログラムが簡単となり、データ処理速度が向上するという効果がある。
【００９４】
更に式(5)〜式(14)に示すタイプIの応答型2段指数平滑化法に後退型Savitzky-Golay法を適用した場合をタイプVとし、下記式(46)〜式(57)に示す。
データの平滑化    ：S1t = α1t SGB0Dt + (1 - α1t) (S1t-1 + B1t-1)        式(46)
平滑化されたデータの傾き：B1t  = γ1t SGB1Dt + (1 -γ1t) B1t-1          式(47)
SGB0Dt = (83Y1t + 54Y1t-1 + 30Y1t-2 + 11Y1t-3 - 3Y1t-4 - 12Y1t-5 -16Y1t-6 
- 15Y1t-7 - 9Y1t-8 + 2Y1t-9 + 18Y1t-10 ) / 143              式(48)
SGB1Dt = (945S1t + 456S1t-1 + 67S1t-2 - 222S1t-3 - 411S1t-4 - 500S1t-5 
- 489S1t-6 - 378S1t-7 - 167S1t-8 + 144S1t-9 + 555S1t-10) / 4290   
式(49)
平滑化係数： α1t =  (Kα-Lα) Fα+ Lα                               式(50)
応答係数：    Fαt = (｜δαt/Δαt｜+φ)N                            式(51)
相対誤差：    δαt = A1（Y1t - S1t） + (1 - A1) δαt-1              式(52)
絶対誤差：    Δαt = A1 |Y1t - S1t| + (1 - A1) Δαt-1 + φ          式(53)
平滑化係数： γ1t = (Kγ-Lγ) Fγ+ Lγ                                式(54)
応答係数：   Fγt = (｜δγt/Δγt｜+φ)N                           式(55)
相対誤差：     δγt =A2{（S1t - S1t-1）- B1t}+(1 - A2)δγt-1      式(56)
絶対誤差：     Δγt =A2|（S1t - S1t-1）- B1t|+(1 - A2) Δγt-1+φ  式(57)
タイプVの場合は、タイプIVよりも応答性等は劣るが、データ処理プログラムが簡単となると言う効果がある。要求されるデータ処理の性能に応じて必要なレベル、複雑さの方式を選択すれば良い。
【００９５】
更に式(19)〜式(24)に示すタイプIIの応答型2段指数平滑化法に後退型Savitzky-Golay法を適用した場合をタイプVIとし、下記式(58)〜式(65)に示す。
データの平滑化    ：S1t = α1t SGB0Dt + (1 - α1t) (S1t-1 + B1t-1)        式(58)
平滑化されたデータの傾き：B1t  = α1t SGB1Dt + (1 -α1t) B1t-1          式(59)
SGB0Dt = (83Y1t + 54Y1t-1 + 30Y1t-2 + 11Y1t-3 - 3Y1t-4 - 12Y1t-5 -16Y1t-6 
- 15Y1t-7 - 9Y1t-8 + 2Y1t-9 + 18Y1t-10 ) / 143              式(60)
SGB1Dt = (945S1t + 456S1t-1 + 67S1t-2 - 222S1t-3 - 411S1t-4 - 500S1t-5 
- 489S1t-6 - 378S1t-7 - 167S1t-8 + 144S1t-9 + 555S1t-10) / 4290   
式(61)
平滑化係数： α1t =  (K - L) Ft+ L                                    式(62)
応答係数：    Ft = 1 - Exp[-δαt2/(2σt2)]                           式(63)
相対誤差：    δαt = A1（Y1t - S1t） + (1 - A1) δαt-1              式(64)
予測誤差分散： σt2 = A1（Y1t - S1t）2 + (1 - A1)σt-12               式(65)
タイプVIの場合もタイプIVよりも応答性等は劣るが、データ処理プログラムが簡単となると言う効果がある。また、タイプVにおいてN=0の場合は、平滑化係数が固定の2段指数平滑化法に後退型Savitzky-Golay法を適用した方式となり、最も簡単となる。要求されるデータ処理の性能に応じて必要なレベル、複雑さの方式を選択すれば良い。
【００９６】
タイプIIは、タイプIの応答型２段指数平滑化法に確率密度関数を用いたデータ応答型平滑化係数を導入したものである。また、タイプIIIは、タイプIIにデータ応答型加算係数を導入したものである。更にタイプIVは、タイプIIIに後退型Savitzky-Golay法の前処理を導入したものである。タイプVと同様に各々の改善要素である、確率密度関数を用いたデータ応答型平滑化係数、データ応答型加算係数を個別にタイプIの応答型２段指数平滑化法に適用しても良い。また、確率密度関数を用いたデータ応答型平滑化係数、データ応答型加算係数、後退型Savitzky-Golay法を適宜、組合わせてタイプIの応答型２段指数平滑化法に適用しても良い。前述と同様に要求されるデータ処理の性能に応じて必要なレベル、複雑さの方式を選択すれば良い。
【実施例２】
【００９７】
上述のようにデータスムージング処理、データ微分処理において、Ｓ/Ｎ比性能向上とデータ応答性能向上（遅延時間短縮）は、トレードオフ関係である。このため、上述したタイプIの応答型２段指数平滑化法、あるいは本実施例に記載したタイプII、タイプIII、タイプIV、タイプV、タイプVI等のデータスム-ジング処理方法およびデータ微分処理方法では、データ処理の対象である入力データに応じて各データ処理方法のパラメータを最適化する必要がある。
【００９８】
従来は、データ処理のパラメータを順次変更し、データ処理を行い、各データスムージング波形、データ微分波形、およびＳ/Ｎ比、遅延時間等の数値データを総合的に俯瞰し、最適なパラメータを見出していた。しかしながら、上述の方法では、最適なパラメータを見出すのに長時間を要し、また最適なパラメータを見出す時間短縮等には、データ処理の知識や経験が必要等の課題があった。
【００９９】
また、非特許文献３には、1段指数平滑化法において、1期先予測値の誤差の総和を最小とすることにより最適な平滑化パラメータを推定する方法が開示されている。しかし、この方法では、予測値の曲線のなめらかさについて配慮されていないため、変化点を検出する１次微分処理、２次微分処理において、ノイズが大きく、S/N比が良くないという課題があった。
【０１００】
このため、入力データをデータスムージング処理、データ微分処理するのに最適なパラメータを簡便かつ短時間に自動で見出すことができる方法を以下に示す。
【０１０１】
まず、下記式(66)の評価関数Ｗを使用する。
【０１０２】
評価関数Ｗ＝平均二乗誤差(Ｅ)＋係数λ×２次微分の二乗平均(Ｄ)    式（66）
ここで、平均二乗誤差であるＥは、入力データとデータスムージング処理結果の平均二乗誤差であり、入力データに対する、データスムージング処理結果の適合度を評価する。２次微分の二乗平均であるＤは、データスムージング処理波形の曲線のなめらかさを評価する。係数λは、任意の数値で、上記の入力データに対する適合度評価と曲線のなめらかさ評価の重要度の比率を調整する。
【０１０３】
入力データに対するデータスムージング処理結果の適合度が大きいほど、オーバーシュート等が無く、かつ、データ処理の応答性が良く、かつ、データ処理の遅延時間が小さい。また、データスムージング処理波形がなめらかな曲線である程、データ１次微分処理波形、データ２次微分処理波形は、なめらかな曲線となる。その結果、データ１次微分処理とデータ２次微分処理の各々のＳ／Ｎ比が向上する。曲線のなめらかさは、データスムージング処理波形の２次微分の二乗平均Ｄで評価するが、ここでは差分法で２次微分を算出した。
【０１０４】
２次微分の二乗平均Ｄは、直線に近づくほど、小さな値となる。従って、単に２次微分の二乗平均Ｄが小さいだけでなく、平均二乗誤差Ｅが小さい値であり、かつ、入力データに対するデータスムージング処理結果の適合度も良く両者が両立した場合、Ｓ/Ｎ比が良く、データ処理の応答性が良い（遅延時間が少ない）最適なデータ処理となる。本実施例で使用した評価関数Ｗを下記式（67）に示す。
【０１０５】
Ｗ＝Σ（Y1t - S1t）2 / N＋λ×Σ{（S1t+1 - 2×S1t + S1t+1）/ΔT2 }2 / N 
式(67)
ここで、Nはデータ数、ΔＴは入力データのサンプリング時間（時間間隔）である。データスムージング処理、データ微分処理の最適パラメータは、評価関数Ｗが最小値となるように、最急降下法等の勾配法等を用いて導出する。評価関数Ｗの特性を２次元グラフで示すため、簡単なタイプIの応答型２段指数平滑化法の場合を例に図１１〜図１３を用いて説明する。但し、入力データは、図７〜図１０で使用した入力データと類似しているが、違うデータを使用した。
【０１０６】
図１１は、タイプIの応答型２段指数平滑化法の１つのパラメータを順次変化させた場合における、入力データ(Y1t)及びデータスムージング処理結果(S1t)を示す（Ａ）と２次微分処理結果（B2t）を示す(Ｂ)とからなる図である。パラメータ１、２、３、４の順に数値は大きい。最も数値が大きいパラメータ４の場合、入力データ波形とデータスムジング波形はほぼ一致し、適合度が良いことがわかる。しかし、２次微分波形は、ゼロクロスの時間が早く、応答性は良いが、ノイズが多く、Ｓ/Ｎ比が良くないことが分かる。一方、数値が最も小さいパラメータ１の場合、入力波形に対してデータスムージイング波形はオーバーシュートがあり、適合度は良く無く、応答性が良くないことが分かる。また、２次微分波形は、ノイズが小さいが、信号強度も小さく、また、ゼロクロスの時間が遅く、遅延時間が大きいことが分かる。図１１の場合、（Ｂ）の２次微分処理波形からは、Ｓ/Ｎ比と遅延時間を両立しているのは、パラメータ２の場合と判断できる。
【０１０７】
図１２は、図１１の場合において、平均二乗誤差Ｅ、２次微分の二乗平均Ｄおよび評価関数Ｗのパラメータ依存性を示した図である。また、図１１に対応するパラメータの位置を矢印で図示している。また、×印は、最急降下法等の勾配法等を用いて自動的に導出した評価関数Ｗが最小値となるパラメータの位置と評価関数Ｗの数値の位置を示したものである。図１３は、図１２における、評価関数Ｗを用いて自動で導出した×印の位置のパラメータを使用した場合の、入力データ(Y1t)及びデータスムージング処理結果(S1t)を示す（Ａ）と２次微分処理結果（B2t）を示す(Ｂ)とからなる図である。
【０１０８】
図１２に示すようにパラメータの数値が増加するとともに平均二乗誤差Ｅは減少し、２次微分の二乗平均Ｄは増加することが分かる。この結果、評価関数Ｗは、最小値を有することが分かる。この結果、評価関数Ｗの最小値を求めることによりＳ/Ｎ比と応答性を両立するパラメータを見出すことができることが分かる。図１１に示す定性的なデータ処理波形の特性のパラメータ依存性を図１２において定量的に示している。
【０１０９】
図１２に示す×印のパラメータの最適値は、パラメータ２とほぼ一致する。また、図１３に示すデータスムジング処理波形及びデータ２次微分処理波形は、図１２のパラメータ２のデータスムージング処理波形とデータ２次微分処理波形と各々ほぼ一致していることが分かる。
【０１１０】
本実施例では、説明の都合上、１つのパラメータを変化させた場合について説明したが、複数のパラメータを最急降下法等の勾配法等を用いて導出することができる。この場合、局所最小値（最適値）の有無等を配慮し、パラメータの初期値および探索範囲を検討する必要がある。また、実施例では、説明の都合上、簡単なタイプIの応答型２段指数平滑化法の場合について説明したが、上述のタイプII、タイプIII、タイプIV、タイプV、タイプVIのデータ処理方式、あるいは１段指数平滑化法（指数加重移動平均：ＥＷＭＡ）、応答型1段指数平滑化法、ローパスフィルタ、カルマンフィルタ等のその他のデータスムージング処理、あるいは差分法等によるデータ微分処処理においても同様に評価関数Ｗを最小化することによってＳ/Ｎ比と応答性を両立する最適なパラメータを簡便に見出すことができる。
【０１１１】
また、本実施例では、式（66）および式（67）の２次微分の二乗平均Ｄにデータスムージング処理結果の２次微分値を使用したが、１次微分処理結果あるいは２次微分処理結果の２次微分の二乗平均を用いても良い。この場合は、各々、１次微分波形の曲線のなめらかさ、２次微分波形の曲線のなめらかさを直接、評価することになる。この場合、データスムージング処理結果と比較して１次微分処理結果、２次微分処理結果は、順に、数値の絶対値が小さくなるため、式（66）および式（67）の係数λを大きくする等の調整が必要である。
【０１１２】
非特許文献３に記載された、式(1)〜式(4)に示す応答型１段指数平滑化法では、現在の時間tでの入力データY1tと現在時間tでのデータ平滑化の予測値S1tを用いて1期先の時間t+1でのデータ平滑化の予測値St+1を導出している。一方、タイプI〜タイプVIの実施例では、現在時間tでの入力データY1tと1期前の時間t-1でのデータ平滑化の予測値S1t-1を用いて現在の時間tでのデータ平滑化の予測値S1tを導出している。前者は、「1期先予測」、後者は「現在推定」と呼ぶことができる。タイプI〜タイプVIの実施例の場合でもS1t-1→S1t、S1t→S1t+1等の変換により「1期先予測」に変更しても良い。
【０１１３】
但し、一般に「現在推定」の方が「1期先予測」よりもデータ平滑化の予測値の精度が良い。また、非特許文献３には、「1期先予測」の1段指数平滑化法において、1期先予測値の誤差の総和を最小とすることにより最適な平滑化パラメータを推定する方法が開示されている。
【０１１４】
「現在推定」の場合には、平滑化係数を1とすると、Y1t＝S1tとなり、誤差あるいは平均２乗誤差がゼロとなるため、最急降下法等を用いて最適な平滑化パラメータと推定することができない。しかし、実施例で示したように「現在推定」であっても平均二乗誤差と曲線のなめらかさを考慮した式(66)の評価関数Ｗを用いることにより最適な平滑化パラメータを導出することができる。
【０１１５】
本実施例によれば、入力データをデータスムージング処理、データ微分処理するのに最適なパラメータを短時間でデータ処理の知識や経験に頼ることなく、自動で見出すことができる効果がある。これにより、一般ユーザが使い勝手が良く、即ち、ユーザビリティが良い、データ処理装置、データ処理方法および処理室を制御する制御装置とを備える処理装置を提供することができる効果がある。
【０１１６】
以上、本発明は、特に短時間のエッチングプロセスや短時間での変化を伴うエッチングの終点検出等に有効である。半導体デバイスの高集積化、微細化に伴なって半導体エッチングでは、多層薄膜をエッチングする工程が増加しており、短時間のエッチングプロセス、短時間での変化を伴うエッチングステップでのエッチング終点検出が重要となっている。
【０１１７】
この短時間プロセス及び短時間変化プロセス対応は、本発明によるS/N比向上と応答性向上の両立した総合性能の向上により、明瞭な1次微分波形、２次微分波形を検出することができるようになる。このことにより、高精度にエッチング終点判定が可能となり、これに基づきプラズマ処理室でのプロセスを制御することによって半導体ウェハを安定した性能で高精度に微細加工できるようになった。
【０１１８】
以上、上述した実施例では、マイクロ波プラズマエッチング装置におけるエッチングの終点検出に適用して高精度にエッチングを行う場合について詳述したが、他のプラズマ生成方式（誘導結合型、平行平板型等）のエッチング装置および成膜装置、あるいは他の分野の処理装置およびその他の装置等においても装置等から得られる数値データを入力として本発明のデータ処理装置およびデータ処理方法を適用することによって装置の状態を監視し、状態の変化点を高精度に検出することができる。このため、対象とする装置を高精度に制御できるという効果がある。またその他の装置の制御においても、同様の作用効果がある。
【０１１９】
また、需要供給予測等の経済や金融分野でも本発明のデータ処理装置およびデータ処理方法を適用することにより高精度にデータを分析することができるという効果がある。
【０１２０】
さらに本発明は、逐次デ−タ処理において、デ−タスム−ジング処理、デ−タ微分処理を高S/N比でデ−タ遅延が小さく、また、データ処理開始初期の期間も信頼性が高いデ−タ処理を行うことができる。
【０１２１】
また、本発明は、データスムージング値、１次微分値、２次微分値が高いS/N比、短い遅延時間、あるいは、データ処理開始の初期も高い信頼性で逐次、リアルタイムに得られる。さらに、本発明は、このデータスムージング値、１次微分値、２次微分値を用いて対象とする系を高精度に制御することができる。
【０１２２】
なお、本発明は、上述した各実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上述した各実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることも可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【０１２３】
また、実施例１および実施例２において、「スムージング」と「平滑化」が混在しているが、実施例１および実施例２における、「スムージング」と「平滑化」は、同義語として使用している。
【符号の説明】
【０１２４】
１…デ−タ処理装置、２…データ入出力装置、３…データ記憶装置、４…データ処理プログラム記憶装置、５…データ演算処理装置、６…対象とする系、７…容器、８…放電管、９…石英板、１０…石英窓、１１…処理室、１２…排気用開閉バルブ、１３…真空排気装置、１４…ガス配管、１５…石英シャワープレート、１６…排気速度可変バルブ、１７…コイル、１８…コイル、１９…ヨーク、２０…マグネトロン、２１…整合器、２２…矩形導波管、２３…円矩形変換器、２４…円形導波管、２５…空洞共振器、２６…ウェハ載置用電極、２７…ウェハ、２８…高周波電源、２９…光ファイバー、３０…分光器、３１…光ファイバー、３２…分光器、３３…システム制御装置
【図１】
2020061569.tif 000003
【図２】
2020061569.tif 000004
【図３Ａ】
2020061569.tif 000005
【図３Ｂ】
2020061569.tif 000006
【図４】
2020061569.tif 000007
【図５】
2020061569.tif 000008
【図６】
2020061569.tif 000009
【図７】
2020061569.tif 000010
【図８】
2020061569.tif 000011
【図９】
2020061569.tif 000012
【図１０】
2020061569.tif 000013
【図１１】
2020061569.tif 000014
【図１２】
2020061569.tif 000015
【図１３】
2020061569.tif 000016