説明01 - 版の履歴

2022年1月4日 (火) 01:33にKamimuraによる

2022-01-04T01:33:55Z

← 古い版		2022年1月4日 (火) 01:33時点における版
27行目:		27行目:
	: 大きければ、ローレンツ関数を1つ増やす。この時、1つ目の初期値は今回の最小二乗法で求まった値を用いる。		: 大きければ、ローレンツ関数を1つ増やす。この時、1つ目の初期値は今回の最小二乗法で求まった値を用いる。
	: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left \| r \right \vert_{calc}^2</math>を<math>\left \| r \right \vert_{exp}^2</math>から引いた値計算する。この差の関数のうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。絶対値が最も大きな点が極大値の場合と極小値の場合の2パターンがある。模式的に示すと以下のようになる。		: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left \| r \right \vert_{calc}^2</math>を<math>\left \| r \right \vert_{exp}^2</math>から引いた値計算する。この差の関数のうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。絶対値が最も大きな点が極大値の場合と極小値の場合の2パターンがある。模式的に示すと以下のようになる。
	[[File:theory.png\|~~200px~~]]		[[File:theory.png\|500px]]
	: <math>A</math>の大きさは1つ目で求まったAの半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。今回はAの符号も検討する必要がある。これは、差の値が入力ファイルのように反射率が下に凸とは限らないからだ。よってνに対応する値が最小値なら正、最大値なら負というように与えている。最小値の場合、正として与えるのは1つ目の時と同じで反射率が最小の時吸収が最大になると考えるからである。逆の場合は、物理的な意味は考えづらいが数学的に符号を反対にすることが妥当だと考えたからである。		: <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}^2</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>に対応する値が最小値なら正、最大値なら負というように与えている。

2022年1月4日 (火) 01:32にKamimuraによる

2022-01-04T01:32:18Z

← 古い版		2022年1月4日 (火) 01:32時点における版
26行目:		26行目:
	: これが入力の⑩で与えられた値より小さければ結果を返す。		: これが入力の⑩で与えられた値より小さければ結果を返す。
	: 大きければ、ローレンツ関数を1つ増やす。この時、1つ目の初期値は今回の最小二乗法で求まった値を用いる。		: 大きければ、ローレンツ関数を1つ増やす。この時、1つ目の初期値は今回の最小二乗法で求まった値を用いる。
	: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left \| r \right \vert_{calc}^2</math>を<math>\left \| r \right \vert_{exp}^2</math>~~から引いた値計算する。その様子を模式的に図にすると以下のようになる。~~		: 2つ目の初期値の作り方を説明する。これは、1つ目のローレンツ関数を用いて計算される <math>\left \| r \right \vert_{calc}^2</math>を<math>\left \| r \right \vert_{exp}^2</math>から引いた値計算する。この差の関数のうち絶対値が最も大きい値に対応する波数を2つ目の<math>\nu_l</math>の初期値として与える。絶対値が最も大きな点が極大値の場合と極小値の場合の2パターンがある。模式的に示すと以下のようになる。
	~~[[File:差の関数.png\|200px]]~~		[[File:theory.png\|200px]]
	~~: このうち絶対値が最も大きい値に対応する波数を2つ目の~~<math>\nu_l</math>~~の初期値として与える。~~
	~~: つまり~~
	[[File:~~nuの与え方~~.png\|200px]]
	: <math>A</math>の大きさは1つ目で求まったAの半分としている。これは屈折率の分散の概形が一つ目で与えられていると考え、足りない部分を補うという考えから、一つ目のAより小さいと考えたためである。今回はAの符号も検討する必要がある。これは、差の値が入力ファイルのように反射率が下に凸とは限らないからだ。よってνに対応する値が最小値なら正、最大値なら負というように与えている。最小値の場合、正として与えるのは1つ目の時と同じで反射率が最小の時吸収が最大になると考えるからである。逆の場合は、物理的な意味は考えづらいが数学的に符号を反対にすることが妥当だと考えたからである。		: <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}^2</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>に対応する値が最小値なら正、最大値なら負というように与えている。

2021年12月1日 (水) 03:09にKamimuraによる

2021-12-01T03:09:09Z

← 古い版		2021年12月1日 (水) 03:09時点における版
19行目:		19行目:
	: <math>\gamma</math>は15とする。これは参考文献[2]<ref name = "ref2" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。		: <math>\gamma</math>は15とする。これは参考文献[2]<ref name = "ref2" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。
	: <math>A</math>は7.5とする。これは参考文献[2]<ref name = "ref2" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。		: <math>A</math>は7.5とする。これは参考文献[2]<ref name = "ref2" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。
	: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。屈折率の虚部が吸収を表すことは[[~~詳細#~~詳細03\|詳細03]]で説明している。		: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。屈折率の虚部が吸収を表すことは[[詳細03\|詳細03]]で説明している。
	1つずつ与える。		1つずつ与える。
	: 最小二乗法を行う際、パラメータの範囲はそれぞれ<math>A</math>は-200から200。<math>\nu_l</math>入力④で与えられた範囲。<math>\gamma</math>は0から500で行っている。		: 最小二乗法を行う際、パラメータの範囲はそれぞれ<math>A</math>は-200から200。<math>\nu_l</math>入力④で与えられた範囲。<math>\gamma</math>は0から500で行っている。

2021年12月1日 (水) 03:07にKamimuraによる

2021-12-01T03:07:39Z

← 古い版		2021年12月1日 (水) 03:07時点における版
19行目:		19行目:
	: <math>\gamma</math>は15とする。これは参考文献[2]<ref name = "ref2" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。		: <math>\gamma</math>は15とする。これは参考文献[2]<ref name = "ref2" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。
	: <math>A</math>は7.5とする。これは参考文献[2]<ref name = "ref2" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。		: <math>A</math>は7.5とする。これは参考文献[2]<ref name = "ref2" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。
	: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>~~であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。~~		: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。屈折率の虚部が吸収を表すことは[[詳細#詳細03\|詳細03]]で説明している。
	1つずつ与える。		1つずつ与える。
	: 最小二乗法を行う際、パラメータの範囲はそれぞれ<math>A</math>は-200から200。<math>\nu_l</math>入力④で与えられた範囲。<math>\gamma</math>は0から500で行っている。		: 最小二乗法を行う際、パラメータの範囲はそれぞれ<math>A</math>は-200から200。<math>\nu_l</math>入力④で与えられた範囲。<math>\gamma</math>は0から500で行っている。

2021年12月1日 (水) 03:01にKamimuraによる

2021-12-01T03:01:51Z

← 古い版		2021年12月1日 (水) 03:01時点における版
1行目:		1行目:
	<span id="説明01" style="font-size: 150%~~; color:Aqua~~;">01:フィッティング方法</span>		<span id="説明01" style="font-size: 150%;">01:フィッティング方法</span>
	: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[詳細01\|詳細01]]</u>参照。		: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[詳細01\|詳細01]]</u>参照。
	: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。		: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。

2021年11月29日 (月) 08:30にKamimuraによる

2021-11-29T08:30:51Z

← 古い版		2021年11月29日 (月) 08:30時点における版
1行目:		1行目:
	<span id="説明01" style="font-size: 150%; color:Aqua;">01:フィッティング方法</span>		<span id="説明01" style="font-size: 150%; color:Aqua;">01:フィッティング方法</span>
	: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[#詳細01\|詳細01]]</u>参照。		: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[詳細01\|詳細01]]</u>参照。
	: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。		: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。
	: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "ref1" />参照		: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "ref1" />参照
	: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。		: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。
	: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>		: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>
	: のようにp波とs波の平均値で与えられる。詳細は<u>[[#詳細02\|詳細02]]</u>		: のようにp波とs波の平均値で与えられる。詳細は<u>[[詳細02\|詳細02]]</u>
	: また、各<math>r</math>はフレネルの公式		: また、各<math>r</math>はフレネルの公式
	: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>		: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>

2021年11月29日 (月) 08:26にKamimuraによる

2021-11-29T08:26:36Z

← 古い版		2021年11月29日 (月) 08:26時点における版
40行目:		40行目:
	<ref name = "ref1">遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー</ref>		<ref name = "ref1">遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー</ref>
	<ref name = "ref2">論文</ref>		<ref name = "ref2">論文</ref>
	<references />		</references>

2021年11月29日 (月) 08:23にKamimuraによる

2021-11-29T08:23:40Z

← 古い版		2021年11月29日 (月) 08:23時点における版
2行目:		2行目:
	: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[#詳細01\|詳細01]]</u>参照。		: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[#詳細01\|詳細01]]</u>参照。
	: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。		: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。
	: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "~~ref2~~" />参照		: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "ref1" />参照
	: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。		: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。
	: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>		: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>
8行目:		8行目:
	: また、各<math>r</math>はフレネルの公式		: また、各<math>r</math>はフレネルの公式
	: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>		: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>
	によって求まる。公式の詳細は<ref name = "~~ref2~~" />を参照。		によって求まる。公式の詳細は<ref name = "ref1" />を参照。
	: ここで、<math>\sin{\theta}_t</math>はスネルの法則<math>n_1\sin{\theta}_i = n_2\sin{\theta}_t</math>によって計算される。公式の詳細は<ref name = "~~ref2~~" />を参照。		: ここで、<math>\sin{\theta}_t</math>はスネルの法則<math>n_1\sin{\theta}_i = n_2\sin{\theta}_t</math>によって計算される。公式の詳細は<ref name = "ref1" />を参照。
	: また、<math>{\theta}_i, n_1</math>は入力として与えられる定数である。		: また、<math>{\theta}_i, n_1</math>は入力として与えられる定数である。
	: したがって、<math>\left \| r \right \vert_{calc}^2</math>は<math>n_2</math>の関数になっている。		: したがって、<math>\left \| r \right \vert_{calc}^2</math>は<math>n_2</math>の関数になっている。
17行目:		17行目:
	: <math>n_j^0, A, \nu_l, \gamma</math>について初期をそれぞれ与える。		: <math>n_j^0, A, \nu_l, \gamma</math>について初期をそれぞれ与える。
	: <math>n_j^0</math>は吸収がない場合の屈折率の値である。コーシーの方程式<math>n = A + \frac{B}{{\lambda}^2} + \frac{C}{{\lambda}^4}</math>によって近似的に求めることができるので、物質及び解析範囲に合わせてパラメータと<math>\lambda</math>を用いて適切な値を各自計算して、入力⑨として与えてほしい。		: <math>n_j^0</math>は吸収がない場合の屈折率の値である。コーシーの方程式<math>n = A + \frac{B}{{\lambda}^2} + \frac{C}{{\lambda}^4}</math>によって近似的に求めることができるので、物質及び解析範囲に合わせてパラメータと<math>\lambda</math>を用いて適切な値を各自計算して、入力⑨として与えてほしい。
	: <math>\gamma</math>は15とする。これは参考文献[1]<ref name = "~~ref1~~" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。		: <math>\gamma</math>は15とする。これは参考文献[2]<ref name = "ref2" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。
	: <math>A</math>は7.5とする。これは参考文献[1]<ref name = "~~ref1~~" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。		: <math>A</math>は7.5とする。これは参考文献[2]<ref name = "ref2" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。
	: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。		: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。
	1つずつ与える。		1つずつ与える。
38行目:		38行目:
	==参考文献==		==参考文献==
	<references>		<references>
	<ref name = "ref1">論文</ref>		<ref name = "ref1">遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー</ref>
	<ref name = "ref2">~~遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー~~</ref>		<ref name = "ref2">論文</ref>
	<references />		<references />

2021年11月29日 (月) 08:20にKamimuraによる

2021-11-29T08:20:50Z

← 古い版		2021年11月29日 (月) 08:20時点における版
35行目:		35行目:
	: ここまで、入力の⑪がNOの場合の説明をしてきた。次に、入力の⑪がYESの場合の説明をしよう。		: ここまで、入力の⑪がNOの場合の説明をしてきた。次に、入力の⑪がYESの場合の説明をしよう。
	: これは、NOの場合とほとんど変わらない。変わるのは、最初の初期値が入力値を採用するということだけである。そのほかの、残差を求めて、それが十分小さくなるまでローレンツ関数を増やしながら、最小二乗法を繰り返すのは同じである。		: これは、NOの場合とほとんど変わらない。変わるのは、最初の初期値が入力値を採用するということだけである。そのほかの、残差を求めて、それが十分小さくなるまでローレンツ関数を増やしながら、最小二乗法を繰り返すのは同じである。

			==参考文献==
	<references>		<references>
	<ref name = "ref1">論文</ref>		<ref name = "ref1">論文</ref>
	<ref name = "ref2">遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー</ref>		<ref name = "ref2">遠藤雅守：電磁場の物理学ーその発生・伝搬・吸収・増幅・共振を電磁気学で理解するー</ref>
	<references />		<references />

2021年11月29日 (月) 08:17にKamimuraによる

2021-11-29T08:17:50Z

← 古い版		2021年11月29日 (月) 08:17時点における版
2行目:		2行目:
	: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[#詳細01\|詳細01]]</u>参照。		: 屈折率の分散をパラメータを使って表したい。ローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>の重ね合わせによって<math>n_j=n_j^0+\sum_{l=1}^{lmax} \frac{A}{\nu_l-\nu-i\gamma}</math>のように表せると考えている。理由は<u>[[#詳細01\|詳細01]]</u>参照。
	: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。		: 光の強度反射率について最小二乗法を行ってパラメータを決める。つまり、<math>[\left \| r \right \vert_{calc}^2 - \left \| r \right \vert_{exp}^2]^2</math>という式について最小二乗法を行っている。
	: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "ref2">参照		: ここで<math>r</math>は両方とも電場の反射率である。光の強度反射率が電場の反射率の2乗になる理由は<ref name = "ref2" />参照
	: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。		: <math>\left \| r \right \vert_{exp}</math>については入力ファイルの値(単位が%)を100で割って用いる。
	: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>		: <math>\left \| r \right \vert_{calc}^2</math>は<math>\frac{\left \| r^s \right \vert_{calc}^2 + \left \| r^p \right \vert_{calc}^2}{2}</math>
8行目:		8行目:
	: また、各<math>r</math>はフレネルの公式		: また、各<math>r</math>はフレネルの公式
	: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>		: <math>r_s=\frac{n_1\cos{\theta}_i -n_2\cos{\theta}_t}{n_1\cos{\theta}_i +n_2\cos{\theta}_t}</math>, <math>r_p=\frac{n_2\cos{\theta}_i -n_1\cos{\theta}_t}{n_2\cos{\theta}_i +n_1\cos{\theta}_t}</math>
	によって求まる。公式の詳細は<ref name = "ref2">を参照。		によって求まる。公式の詳細は<ref name = "ref2" />を参照。
	: ここで、<math>\sin{\theta}_t</math>はスネルの法則<math>n_1\sin{\theta}_i = n_2\sin{\theta}_t</math>によって計算される。公式の詳細は<ref name = "ref2">を参照。		: ここで、<math>\sin{\theta}_t</math>はスネルの法則<math>n_1\sin{\theta}_i = n_2\sin{\theta}_t</math>によって計算される。公式の詳細は<ref name = "ref2" />を参照。
	: また、<math>{\theta}_i, n_1</math>は入力として与えられる定数である。		: また、<math>{\theta}_i, n_1</math>は入力として与えられる定数である。
	: したがって、<math>\left \| r \right \vert_{calc}^2</math>は<math>n_2</math>の関数になっている。		: したがって、<math>\left \| r \right \vert_{calc}^2</math>は<math>n_2</math>の関数になっている。
17行目:		17行目:
	: <math>n_j^0, A, \nu_l, \gamma</math>について初期をそれぞれ与える。		: <math>n_j^0, A, \nu_l, \gamma</math>について初期をそれぞれ与える。
	: <math>n_j^0</math>は吸収がない場合の屈折率の値である。コーシーの方程式<math>n = A + \frac{B}{{\lambda}^2} + \frac{C}{{\lambda}^4}</math>によって近似的に求めることができるので、物質及び解析範囲に合わせてパラメータと<math>\lambda</math>を用いて適切な値を各自計算して、入力⑨として与えてほしい。		: <math>n_j^0</math>は吸収がない場合の屈折率の値である。コーシーの方程式<math>n = A + \frac{B}{{\lambda}^2} + \frac{C}{{\lambda}^4}</math>によって近似的に求めることができるので、物質及び解析範囲に合わせてパラメータと<math>\lambda</math>を用いて適切な値を各自計算して、入力⑨として与えてほしい。
	: <math>\gamma</math>は15とする。これは参考文献[1]<ref name = "ref1">から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。		: <math>\gamma</math>は15とする。これは参考文献[1]<ref name = "ref1" />から、<math>\gamma</math>は15付近であるという傾向がわかっているからだ。
	: <math>A</math>は7.5とする。これは参考文献[1]<ref name = "ref1">から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。		: <math>A</math>は7.5とする。これは参考文献[1]<ref name = "ref1" />から、<math>A</math>と<math>\gamma</math>の間に<math>\frac{A}{\gamma} = {\kappa}_{max}</math>という関係があり、さらに<math>{\kappa}_{max}</math>が0.5付近の値を多くとるということがわかっているからだ。
	: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。		: 最後に<math>\nu_l</math>についてこれは入力④で与えられた解析範囲において入力⑤で読み込みを指定した列の反射率が最小になるところとしている。これはローレンツ関数<math>\frac{A}{\nu_l-\nu-i\gamma}</math>を実部と虚部に分けたとき、<math>\frac{A\gamma}{{\nu_l-\nu}^2 + {\gamma}^2} + i\frac{{\nu_l-\nu}^2}{{\nu_l-\nu}^2 + {\gamma}^2}</math>となり、虚部が吸収を表すことが関係している。虚部において最大になるところこそが<math>\nu_l</math>であり、この時吸収は最大である。これは反射が最小に対応する。ゆえにこのように与えている。
	1つずつ与える。		1つずつ与える。