- ベストアンサー
Monte Carlo法の誤差
Monte Carlo法の計算誤差に関して詳しい本かwebサイトを 教えてください。 知りたい例の一つとして、空間に256個の球がちらばっている状態 にてz軸下から上方向に 数 nrayの光線を入射した時にその面積 計算の誤差(およそ?%のオーダー)を知りたいと思います。 上記は一例ですが、入射光線の数と誤差の関係が他の例でも 知ることができましたらと思います。
- 数学・算数
- 回答数1
- ありがとう数1
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
- ベストアンサー
一般にMonte Carlo法の精度は、 （１）モデリングの精度（どのくらいはしょらないでモデル化したか） （２）領域を分ける細かさ （３）試行回数（精度は一般に√Nに比例して向上） （４）現実の確率過程と発生させる乱数の適合度 （一様乱数/正規乱数/Poisson乱数などのうちどれを使うべきか ＆ 不規則性・偏りなどに関しどのくらい"良い"乱数か） などによって決まります。 なので、理論一発で「Monte Carlo法で光線何本通したら誤差何％」なんて 簡単には出ませんよ。
関連するQ&A
- 誤差と最小二乗法
単振り子の振れの大きさを小さくし、 振り子の長さLと周期Tを5通りほど計測しました。 そこでT=2π√(L/g)⇔g=4π^2L/T^2から 重力加速度を5通り求めました。 結果、どれも相対誤差0.3％以内で、 どれもマイナスの誤差になりました。 次にT=2π√(L/g)⇔T=(2π/√g)*√Lとし、 √Lをx軸に、Tをy軸にとり、(2π/√g)を傾きとするグラフを描きました。 ここで、最小二乗法よりグラフの傾きを求め、 そこから重力加速度gを求めたところ、 相対誤差が3％になり、さらにそれはプラスの誤差でした。 そこで質問です。 最小二乗法を用いるのは今回が初めてなんですが、 このように測定値を単体で考えたときと、 測定値をまとめて最小二乗法で考えたときとで、 こんなにも誤差が変わるのはよくあることですか？ それとも計算に誤りがあるのでしょうか。 よろしくお願いします。
- 締切済み
- 物理学
- 複屈折材の異常光線について質問です
複屈折材の異常光線について質問です 書きのように複屈折材の光学軸がφ傾いた素子に面内方向の直線偏光を入れると 異常屈折して異常光線になりますが、 添付画像の赤の光線のように素子に垂直に入射する場合は 書籍などにも載っており、光線の軌跡を計算できるのですが、 素子に角度θをもって入射する場合の軌跡 （図中では緑の光線）を計算したのですが、 考え方からしてわかりません。 この場合もスネルの法則には従わないのでしょうか？ どのようにこの問題を考えればよいのか ご教授願えませんでしょうか？
- ベストアンサー
- 物理学
- 磁界中での運動において
xyz空間において+z方向に磁界が発生している（ここでは磁束密度の強さをＢとしてください）状況で、 xz平面でz軸からある角（ここではθ）をなす方向に速さ（v)で電子を入射した場合、螺旋運動をしますよね。 そこまでは良いのですが、この状況で電子はz軸上（例えば座標が(0,0,z)）にどのような条件で出現するのでしょうか？ 私は、-z側から見て円運動を始め、その円運動を続けるという風に考えたのですが・・・。
- 締切済み
- 物理学
- 誤差を最小化する相似変換行列の求め方
三次元の変換行列について質問です。 空間上のある3点について、座標変換前と変換後の座標が与えられているとします。 その座標を元に、行った変換の表す行列を知りたいと考えております。 ただし変換は並進と回転のみとし、剪断変形等は行わないという条件です。 （行列は同次変換の形などで得られれば良いです） ただし、変換前と変換後では三点の相対的な位置関係は完璧には一致しておりません。 そのため得られた変換を元の3点に行った際に、変換後の3点の座標（既知）との 誤差を最小化したいという問題です。 変換が相似変換でなければ疑似逆行列により計算が可能と思いますが、 並進・回転のみなのでどうしようか聞きたいと思い質問をしました。 パラメータは「x, y, z 軸方向の並進移動量 ＋ x, y, z 軸周りの回転量」の6つなので、 プログラムを組んで6パラメータを少しずつ変化させ、 誤差の二乗和が最小となる箇所を探すという方法は思いついております。 ですが、その他に行列計算により解析的に答えを得る方法はないでしょうか。 また、変換前後の3点を元に誤差を最小化する変換行列を得ることが可能であるなら、 座標の情報が3点でなくそれ以上の場合でも計算は可能でしょうか。 よろしくお願いいたします。
- 締切済み
- 数学・算数
- 誤差を最小化する相似変換行列の求め方
三次元の変換行列について質問です。 空間上のある3点について、座標変換前と変換後の座標が与えられているとします。 その座標を元に、行った変換の表す行列を知りたいと考えております。 ただし変換は並進と回転のみとし、剪断変形等は行わないという条件です。 （行列は同次変換の形などで得られれば良いです） ただし、変換前と変換後では三点の相対的な位置関係は完璧には一致しておりません。 そのため得られた変換を元の3点に行った際に、変換後の3点の座標（既知）との 誤差を最小化したいという問題です。 変換が相似変換でなければ疑似逆行列により計算が可能と思いますが、 並進・回転のみなのでどうしようか聞きたいと思い質問をしました。 パラメータは「x, y, z 軸方向の並進移動量 ＋ x, y, z 軸周りの回転量」の6つなので、 プログラムを組んで6パラメータを少しずつ変化させ、 誤差の二乗和が最小となる箇所を探すという方法は思いついております。 ですが、その他に行列計算により解析的に答えを得る方法はないでしょうか。 また、変換前後の3点を元に誤差を最小化する変換行列を得ることが可能であるなら、 座標の情報が3点でなくそれ以上の場合でも計算は可能でしょうか。 質問は数学のようですが、この手の計算はCGや画像処理等で多く行われていると 感じたことから、情報工学カテゴリーで質問をさせて頂きました。 よろしくお願いいたします。
- 締切済み
- 情報工学
- 空間上の複数の点から球の中心を求めるときの誤差
空間上の複数の点から球の中心を求めるとき、ｒが既知のとき、 (x-a)^2+(y-b)^2+(z-c)^2=r^2 から変数が３つあるので、３点以上あれば、a,b,cが求まるのかなと思うのですが、計測誤差はどう考えればよいのでしょうか？ 感覚的に、x,y,zの計測点が近ければ球の中心からの立体角が小さいので誤差が乗りやすく、しかも推定した中心位置がオフセットしやすく、離れていれば誤差が小さい気がするのですが、この式を見ていてもピンと来ません。計測誤差は同じなのでしょうか？ その場合はｒを大きくしても同じなのでしょうか？ 最終的な目的は１点の計測精度はあがらないので、多点計測で計測精度を良くしたいのですが、理論式で傾向をつかみたいと思ったのですが・・・。 どなたかわかる方が居ましたらご教授願えないでしょうか？ よろしくお願いいたします。
- 締切済み
- 物理学
- 重みつき最小二乗法の誤差
まず、条件を書かせていただきます。 ln(DZ)=ax+B a=-1/T x:条件(x=1、x=5) Z:xを決めると決まる値。（例：x=1ならZ=20、x=5ならZ=68） B:定数 D：条件xの時に行なった測定の測定値 (x=1のときに行なった測定値をD1、x=5のときに行なった測定値をD5) ln(DZ)=yとすれば、y=ax+Bとなります。 一回の測定すると、x=1とx=5に対するD1とD5を測定できます。 一回の測定操作をηとします。 上記の条件で、下記に示す2種類の方法でTの平均値とTの誤差を求めます。 (1) ηを行なうことにより、aとTを求めることができます。 ηを10回行なうことにより、10個のTを求めるます。 10個のTからTの平均値とTの誤差を求めます。 (2) ηを10回行なうことにより、10個のD1と10個のD5を求めるます。 10個のD1と10個のD5を使い、重みつき最小二乗法によって Tの平均値とTの誤差を求めます。 どちらの方法でも平均値と誤差は出ますが、 どちらの方がより正確にTの平均値と誤差を求めることができるのでしょうか。 根拠を示した上で教えてください。 じっくり考えてみましたが、よくわかりません。 文章がおかしいために 、たぶん質問の内容自体が理解できない可能性が高いので分からないことがあれば教えてください。 もし、参考になるホームページや本などがあればご紹介ください。 よろしくお願いいたします。
- ベストアンサー
- 科学
- 逆空間のメッシュについて
逆空間のメッシュについて 数値計算において実空間と同様、逆空間での演算を考えています。 その際逆空間も実空間と同様にメッシュに切って離散的に扱うのですが、 その際のメッシュの切り方を教えてもらいたいです。 1次元では 実空間を r(i) = Δr*i とグリッドを切ったとき 逆空間は k(I) = π*I/(N*Δr) としました。(i,Iは自然数、Δrはグリッド幅、Nはグリッド点の数) すると3次元では 実空間は r(i,j,k) = i*Δx*a + j*Δy*b + k*Δz*c (i,j,kは自然数、Δx,Δy,Δzは各軸方向におけるグリッド幅、 a,b,cはそれぞれ3つの方向の単位ベクトル) で 逆空間は k(I,J,K)=π*I*A/(Nx*Δx) + π*J*B/(Ny*Δy) + π*K*C/(Nz*Δz) (I,J,Kは自然数、Δx,Δy,Δzは各軸方向におけるグリッド幅、 A,B,C 、Nx,Ny,Nzはそれぞれ3つの方向の単位ベクトル、とグリッド点の数) としていいのでしょうか？ 教えてください よろしくお願いします。
- ベストアンサー
- 物理学
- 光学 レトロリフレクタ内の光路長
立方体の一部を切り取ったタイプのレトロリフレクタ（コーナーキューブプリズム）は入射光線を元の方向に反射させるプリズムですが、リフレクタに入射してから光がリフレクタから出てくるまでの光路長（プリズム内部での光路長）はどのように計算したらいいのでしょうか？ また、どの方向から入射してもプリズム内部での光路長は変わりませんよね？ 側面から見れば、単純な幾何学で計算できそうですが、反射の様子を見るとそこまで単純ではなさそうですが。一様自分なりに計算してみたところ、開口部の一辺と同じ長さになりました。分かる方、お願いします。
- 締切済み
- 物理学
お礼
回答ありがとうございました。 精度について参考になりました。
補足
回答ありがとうございます。上記のことなど記載している 本を購入したいのですが、お薦めなどありますでしょうか?