「ComplexRI」の版間の差分

2021年11月29日 (月) 00:54時点における版

概要

本ComplexRIは全反射実験から得られた反射率あるいは吸光度のデータをもとにして媒質の複素屈折率の分散を出力するWebアプリケーションである。 1⃣実験の結果、 2⃣ファイル形式に関する入力、 3⃣実験の条件、 4⃣解析に関する条件という四つの情報を入力として与えると、解析が行える。

チュートリアル

チュートリアルでは4種類全10個の入力について、実際に解析を行い、その結果をもとに説明する。

目次
チュートリアル01：ファイルを入れて解析する(入力1⃣を与える練習)
チュートリアル02：解析に関する条件を指定する(入力4⃣を与える練習)
チュートリアル03：ファイルの形式を変える(入力2⃣を与える練習)
チュートリアル04：実験条件を正しく与える(入力3⃣を与える練習)

チュートリアル01：ファイルを入れて解析してみよう➡チュートリアル01の解説へ

ここでは1⃣実験の結果を入力として与えて解析結果を出力させてみよう。

実験結果のファイルはComplexRIのMANUALページのFile Listからダウンロードできる。

今回はDownload the File01のボタンを押してダウンロードしてこれを使用してほしい。

ダウンロードが出来たらLSRページに行きファイルを入れてみよう。

入力は①から⑩の10個だ。今回は②Fileにダウンロードしたファイルを与えるだけでFよい。

それではexecute LSRボタンをクリックして解析してみよう。

チュートリアル01の解説

解析結果は以下の画像のようなったはずである。

[図1-1]

ここまでComplexRIの基本的な流れを体験してもらった。

ここで入れたファイルについて簡単に説明する。

この様になっているはずである。

一列目には波数、二列目に一列目の波数の光が入射した時の反射率が(％)単位が並んだ構成になっている。

この様に解析に必要な実験結果は波数とそれに対応する反射率の二種類だけであり、これだけの実験結果から解析が行えるということである。

結果については次のチュートリアルで説明する。

チュートリアル02：解析に関する条件を指定しよう➡チュートリアル02の解説へ

ここでは4⃣解析に関する条件を入力として与えてみよう。

まず、先程と同様にLSRページに行き②Fileにファイルを与える。これは前回と同じものを使用してほしい。

先程はこれで実行していたが、以下のように④と⑩に追加の設定をしてみよう。

それではexecute LSRボタンをクリックして解析してみよう。

チュートリアル02の解説

解析結果は以下のようなったはずである。

[図2-1]

④と⑩の入力の意味を説明する前に、結果の見方について簡単に説明する。

(1)数字、(2)表、(3)グラフが並んでいる。

(2)表の数値はフィッティング関数のパラメータである。このフィッティング関数を足し合わせることで、左側の(3)グラフを描いている。また、このグラフの各点の値から反射率を計算してプロットしたのが右側の(3)グラフのcalculateである。実際の回帰ではこの表の数値をいろいろ試して、右側の(3)グラフのcalculateとexperimentと比較するということを行っている。ここで、experimentは入力ファイルとして与えられた反射率の値である。(1)数字は解析が終了した時点で、この比較の差がどれくらい小さくなったかを示している。具体的なフィッティング関数の形や方法については内部処理の説明01：フィッティング方法を参考にしてほしい。

これを踏まえて、④と⑩の入力の意味を説明していこう。

説明は今回の結果をチュートリアル01の時と比較しながら行っていく。

図の比較は以下のようになる。

左がチュートリアル01の結果。右がチュートリアル02➡チュートリアル01の結果を直接見る。の結果。

④の入力をしたことは、グラフの横軸の幅を設定することなっている。確かに、[入力4の図]と同じ範囲になっている。デフォルトの値で行われたチュートリアル01の時と比べると変化がわかりやすいと思う。

このグラフの幅の設定は、言い換えると、解析する波数の範囲設定したことになっている。これは直前に説明した解析の原理から考えればわかる。(3)グラフのcalculateとexperimentを比較するのであった。このときグラフの横幅を変えれば、比べる範囲つまり解析の範囲を変えることに相当するということだ。

次に⑩について説明する。今回は[入力１０の図]のように0.01した。解析は(1)数字がこの値より小さくなったら終了するようになっている。つまり解析の正確さもしくは厳しさを設定することになっている。これは小さくしすぎないように注意してほしい。なぜなら、比較の差分がこれよりも小さくならず、解析がとても長くなることがあるからだ。目安としてはデフォルトの0.02から0.01がよい。これは適宜試してみてほしい。

チュートリアル03：ファイルの形式を変えてみよう➡チュートリアル03の解説へ

ここでは2⃣ファイル形式に関する条件を入力として与えてみよう。

チュートリアル01の時と同様にMANUALページのFile Listからファイルをダウンロードしてほしい。

今回はFile02を使う。

ダウンロード出来たら、LSRページに行き②Fileにファイルを与える。

今回は③，⑤，⑥の入力を以下のように変更してみよう。

それではexecute LSRボタンをクリックして解析してみよう。

チュートリアル03の解説

解析結果は以下のようなったはずである。

じつはこれはチュートリアル01と同じ結果になっている。

左がチュートリアル01の結果。右がチュートリアル03の結果。➡チュートリアル01の結果を直接見る。

これは全く同じになっている。

これは同じ実験結果を入力したことを意味している。しかし、チュートリアル01と02で使用した入力ファイルと今回使用した入力ファイルを比較すると、以下のように一見同じには見えない。これが同じであることを説明する。（一部分のみ示す。）

左がチュートリアル01の入力ファイル。右がチュートリアル02の入力ファイル。

異なる点は２つある。

1つは行が3行になっていることである。これは波数、反射率、吸光度の順に並んでいる。ここで新しく出てきた吸光度というのは反射率と等価な情報を表している。なぜなら、

A=-\log _{10}R

(証明)という関係式が成り立っており、これによって吸光度を反射率に変換することができるからだ。

もう1つは波数と反射率の並びが昇順から降順になっていることだ。ここで、吸光度はこの様に並び替えた後で反射率から計算したので波数に対応している。

このように形は異なるが等価な実験結果を与えた。だから結果が同じだったのだ。

③，⑤，⑥の入力はこのファイルの形の違いをアプリ側に考慮させるために入力した。

⑤では読み込む列の番号を与えたことになっている。今回は入力として1行目の波数と3行目の吸光度を使うのでこのように設定した。

吸光度を使うときは、忘れずに③をabsorptanceに設定してほしい。すると、解析時に反射率に変換して解析する。

この設定はデフォルトではreflectanceになっているのでこれまでは変更する必要がなかった。

最後に⑥はデータの波数のの並び方昇順か降順かを入力している。今回は降順なのでこのように設定した。デフォルトは昇順なので今までは設定する必要がなかった。

チュートリアル04：実験条件を変えてみよう➡チュートリアル04の解説へ

ここでは3⃣実験の条件を入力として与えてみよう。

LSRページに行き②Fileにファイルを与える。これはチュートリアル01と02で使ったものと同じものを使ってほしい。

今回は⑦，⑧，⑨の入力を以下のように変更してみよう。

それではexecute LSRボタンをクリックして解析してみよう。

チュートリアル04の解説

解析結果は以下のようなったはずである。

チュートリアル01とグラフを比較してみると以下のようになる。

左がチュートリアル04の結果。右がチュートリアル01の結果。➡チュートリアル01の結果を直接見る。

同じファイル、つまり同じ実験結果を与えているが、実験条件が違えば、出力も異なることがわかる。正しい分散を求めるには自身の行った実験の条件を正しく与える必要があるということだ。

各入力を説明しよう。

⑦では全反射実験のときに使った基質の種類を指定している。

⑧では全反射実験のときの入射光の角度を指定している。

⑨では分散のバックグラウンドを与えている。これは測定する基質の赤外領域での屈折率を与える必要がある。

これは可視光領域の光に対する屈折率からコーシの方程式のパラメータを求め、これを外挿することによって求めることができる。

詳しい求め方は内部処理の説明02：屈折率のバックグラウンドを参考にしてほしい。

チュートリアルtopへ

マニュアル

目次
入力01：タイトル
入力02：ファイル
入力03：反射率か吸光度か？
入力04：波数範囲
入力05：データの列の指定
入力06：昇順か降順か？
入力07：実験で用いた基質の種類
入力08：入射角
入力09：吸収がないときの屈折率
入力10：残差の指定
入力11：フィッティング関数のパラメータの設定の有無

チュートリアル01：タイトル

チュートリアル02：ファイル

チュートリアル03：反射率か吸光度か？

チュートリアル04：波数範囲

チュートリアル05：データの列の指定

チュートリアル06：昇順か降順か？

チュートリアル07：実験で用いた基質の種類

チュートリアル08：入射角

チュートリアル09：吸収がないときの屈折率

チュートリアル10：残差の指定

チュートリアル11：フィッティング関数のパラメータの設定の有無

マニュアルtopへ

内部処理

目次
説明01：フィッティング方法
説明02：屈折率のバックグラウンド
説明03
説明04

説明01:フィッティング方法

屈折率の分散をパラメータを使って表したい。ローレンツ関数

{\frac {A}{\nu _{l}-\nu -i\gamma }}

の重ね合わせによって

n_{j}=n_{j}^{0}+\sum _{l=1}^{lmax}{\frac {A}{\nu _{l}-\nu -i\gamma }}

のように表せると考えている。理由は詳細01参照。

光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、

[\left|r\right\vert _{calc}^{2}-\left|r\right\vert _{exp}^{2}]^{2}

という式について最小二乗法を行っている。

ここでrは両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は引用エラー: <ref> タグに対応する </ref> タグが不足しています

引用エラー: 無効な <ref> タグです。名前に単なる整数は使用できません。説明的なものにしてください

ComplexRIページtopへ]

@@ 202行目: / 202行目: @@
 : これが入力の⑩で与えられた値より小さければ結果を返す。
 : 大きければ、ローレンツ関数を1つ増やす。この時、1つ目の初期値は今回の最小二乗法で求まった値を用いる。
-: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left | r \right \vert_{calc}^2</math>を<math>\left | r \right \vert_{exp}</math>から引いた値計算する。このうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。<math>A</math>の大きさは1つ目で求まったAの半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。今回はAの符号も検討する必要がある。これは、差の値が入力ファイルのように反射率が下に凸とは限らないからだ。よってνに対応する値が最小値なら正、最大値なら負というように与えている。最小値の場合、正として与えるのは1つ目の時と同じで反射率が最小の時吸収が最大になると考えるからである。逆の場合は、物理的な意味は考えづらいが数学的に符号を反対にすることが妥当だと考えたからである。
+: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left | r \right \vert_{calc}^2</math>を<math>\left | r \right \vert_{exp}^2</math>から引いた値計算する。このうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。<math>A</math>の大きさは1つ目で求まったAの半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。今回はAの符号も検討する必要がある。これは、差の値が入力ファイルのように反射率が下に凸とは限らないからだ。よってνに対応する値が最小値なら正、最大値なら負というように与えている。最小値の場合、正として与えるのは1つ目の時と同じで反射率が最小の時吸収が最大になると考えるからである。逆の場合は、物理的な意味は考えづらいが数学的に符号を反対にすることが妥当だと考えたからである。
-: 2回目以降も同じように残差が入力値より小さいかを判定して、小さければ結果を返し、大きければ1つ増やして、初期値を与えて、最小二乗法行う。これを結果が出るまで繰り返す。2回目以降の残差が大きかった場合の初期値の与え方は1回目と同じである。求まったパ:ラメータから計算される <math>\left | r \right \vert_{calc}^2</math>を<math>\left | r \right \vert_{exp}</math>から引いた値計算する。このうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。<math>A</math>の大きさは1つ目で求まった<math>A</math>の半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。<math>A</math>の符号は<math>\nu_l</math>に対応する値が最小値なら正、最大値なら負というように与えている。
+: 2回目以降も同じように残差が入力値より小さいかを判定して、小さければ結果を返し、大きければ1つ増やして、初期値を与えて、最小二乗法行う。これを結果が出るまで繰り返す。2回目以降の残差が大きかった場合の初期値の与え方は1回目と同じである。求まったパ:ラメータから計算される <math>\left | r \right \vert_{calc}^2</math>を<math>\left | r \right \vert_{exp}^2</math>から引いた値計算する。このうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。<math>A</math>の大きさは1つ目で求まった<math>A</math>の半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。<math>A</math>の符号は<math>\nu_l</math>に対応する値が最小値なら正、最大値なら負というように与えている。
 : ここまで、入力の⑪がNOの場合の説明をしてきた。次に、入力の⑪がYESの場合の説明をしよう。
 : これは、NOの場合とほとんど変わらない。変わるのは、最初の初期値が入力値を採用するということだけである。そのほかの、残差を求めて、それが十分小さくなるまでローレンツ関数を増やしながら、最小二乗法を繰り返すのは同じである。

「ComplexRI」の版間の差分

2021年11月29日 (月) 00:54時点における版

目次

概要

チュートリアル

マニュアル

内部処理

案内メニュー

検索