• 締切済み

ランダムウォーク

酔っ払いのAさんが1直線上を前に進む（x=1)確立をpとし後退(x=-1)をqとしたときに、p+q=1として p=q=1/2、Nを無限大にしたときガウス分布となり微粒子の拡散方程式でもあることを示せと問題があります。 ここでx=mになるときの確立と、分散、期待値までは求めましたが、どうしてもガウス分布になる理由がわかりません。中心極限定理を使うと思うのですが、どう解けばいいでしょうか？ ヒントでもいいので教えてもらえないでしょうか？よろしくお願いします

みんなの回答

回答No.1

＞どうしてもガウス分布になる理由がわかりません。中心極限定理を使うと思うのですが、どう解けばいいでしょうか？ 以下のサイトをご覧下さい。尚、計算は2項分布からnを大きくしてstarlingの公式logn!=nlogn-n（ｎ大きい場合の近似式）を使います。 http://physforumsaga.hp.infoseek.co.jp/WYPsaga/Lecture/Lecture.html 講義Ver2

ahoyasu
質問者

お礼

ん～結構難しいですね。 これを見てなんとか理解してみます。ありがとうございました

  1. カテゴリ
  2. 学問・教育
  3. 自然科学
  4. 物理学

関連するQ&A

  • ランダムウォークのような問題です。教えて頂けませんでしょうか。

    ランダムウォークのような問題です。教えて頂けませんでしょうか。 1次元ランダムウォークの少し変形バージョンのようなものだと思うのですが、 以下の試行を十分大きい回数行ったあとの、粒子の確率分布について どなたか教えて頂けませんでしょうか。 (解法、物理などで似たような命題(ブラウン運動だよ、のような)、もしくは参考書など) 試行:============================================================ 初期位置x=x0にいる1次元粒子を考える。 粒子は位置xにいるとき、ガウス分布f(x)=a*exp(-(x-b)^2/c)の最大ステップ幅で、しかも φを0<φ< 2πの一様分布の乱数だとしたとき、ステップ幅 f(x)*sinφ で移動する。 ================================================================ これを十分な回数繰り返します。 x0, a, b, cの値をいろいろ変えたときの挙動が知りたいと考えています。 aが小さいときにはあまり初期位置から動きませんので、x0を中心としたunipolarな分布で 一方aが大きいときにはbipolarな分布になるのかな...とまではプログラムや直感的方法で 定性的には考えたのですが、それ以上が詰められません。 どなたか分かる方、よろしくお願い致します。

  • 拡散方程式について

    一次元の拡散方程式∂P/∂t=D∂^2P/∂x^2で初期条件がP(x,0)=δ(x)のとき、方程式の解はP(x,t)=1/√４πDtexp(-x^2/4Dt)で与えられ、これは分散が2Dtであるようなガウス分布である。「この確率分布に関する物理量Xの平均を<X>=∫∞～-∞　XP(x,t)dxとすると、<x>=0,<x^2>=2Dtとなる」ようなのですが、「」の部分が理解出来きません。どなたか教えてください。

  • 一次元のランダムウォーク

    下の問が解けなくて困っています。解説をして頂けると非常にありがたいです。お願いしますm(__)m （問）一次元のランダムウォークする粒子を考える。離散的な時間（ｔ＝0,1,2,...）を考え、時刻ｔ＝0では原点にいるとする。この粒子の各時刻での変位は確率pで１、確率p-1で-1とする。時刻tでの粒子の位置ｘ（ｔ）の期待値、分散の値を一般のpで計算しなさい。

  • 単純ランダムウォークの分散の式

    下記のサイトで8.2 <ランダムウォーク>の分散の式を解いています。 https://parco1021.hatenablog.com/entry/2020/05/24/213000 これの V(S_n) = n{ E(X_i) - E(X_i)^2 } = ​n { (p + q) - (2p - 1)^2 } = ​n { p + q - (4p^2 - 4p + 1) } = n ( p + q - 4p^2 + 4p - 1) の n{ E(X_i) - E(X_i)^2 } は、例の分散の式に従うと、正しくは n{ E(X_i^2) - E(X_i)^2 } じゃないですか？ （「2乗」を消しゴムで消した跡がうっすら見えますよね？） （上記が合っているかは分からないですけど） E(X_i)　だか　E(X_i^2) が p + q になる理由が分からないです。 問題の定義から p + q = 1 ですよね・・・いや、ただ書いてみただけです・・・。 どうか教えて下さい。お願いします。

  • 分布収束の問題

    こんにちは。元日ですが、具体的な設定における分布収束の問題についてご質問致します。 nを整数(n>4)、Xを一次元の確率変数、Qを自由度nの逆カイ二乗分布にしたがう確率変数、Cを正の定数とします。また、Xは以下のように構成されているとします。 X～nCQ 【Qにnと定数Cをかけている】 なお、QとXの期待値と分散はそれぞれ以下のようになります。 Qの平均： 1/(n-2) Qの分散： 2/{(n-2)^2×(n-4)} Xの平均： nC/(n-2)【n→∞のときCとなる】 Xの分散： 2×n^2×C^2/{(n-2)^2×(n-4)}【n→∞のとき0となる】 以上の設定において、nを大きくしたとき、√{n}×(X-C)の極限分布と漸近分散を求めたいのですが、√{n}×(X-C)が確率変数列の和という形になっていないので中心極限定理等が使えず、求め方が分からない状況となっています。 nを大きくしたとき、√{n}×(X-C)は分布収束するのでしょうか。また、分布収束する場合は、極限分布はどのように求めたらいいのでしょうか。 よろしくお願い致します。

  • 中心極限定理について

    中心極限定理で理解できないところがあります。 多分、何か勘違いをしているところがあると思うので、分かる方は教えてください。 平均値 μ*，分散 σ2* をもつ，任意の分布に従う乱数列 x1，x2， … ，xnがあるとき，その平均値 ave(xn) = (x1+X2+・・・+xn)/n の確率分布は，n が大きくなるとき，平均値 μ*，分散 σ2* / n である正規分布に収束する。 　すなわち， (ave(xn) - μ*)/(σ*/√n) は，n が大きいとき，平均値 0，分散 1 の標準正規分布に従うとみなしてよい。 はどうやって証明するのでしょうか？ また、１２個の乱数rand()を発生させた場合、分散は12*1/12=1になるようですが、これは中心極限定理から分散がσ2/nの正規分布に近づくという結果（こちらは1/12を12で割っている）に矛盾しているように思うのですが、どこが勘違いしているのでしょうか？ よろしくお願い致します。

  • 有限母集団の中心極限定理

    訳あって確率統計の勉強をしておりますが、中心極限定理について質問があります。 「母集団の平均をm、分散をvとすると、そこから抽出したn個の標本の平均の分布は、平均m、分散v/nという分布になり、標本数をn→無限大とすると、分布は母集団の分布によらず正規分布に近づく。」 とありますが、母集団が有限個（N個）の集合ならどうなるでしょうか。 その場合、標本数をnがNに等しくなった時点で平均はm、分散0、つまり標本から母集団の平均が完全に推定（決定）することになります。 ●有限母集団の場合も中心極限定理は成り立つのか？ ●成り立つならn→Nで分散が0になるという点はどう表現されるのか？ このあたりを教えてください。

  • 統計について　【ランダムウォーク？】

    教えて！gooを始めて利用させていただきます。 次の問題に関してご質問させていただきます。 問　第i期に次のような確率分布に従う株価の変化Xiを考えなさい： P(xi)=p　　for xi= a 1-p for xi=-a ただし，0<p<1, a>0。 すなわち、第i期に確率pで株価はaほど上昇し、確率1-pでaほど下落する可能性があります。 簡単化のために、1)現在（第0期）の株価S0は0です、2)各期の株価の変化Xiは独立に生じるとしましょう。 このとき、第n期の株価Snをつぎのように表すことができます： Sn=X1+X2+・・・Xn この時、確率変数Snの確率関数、積率母関数、平均、分散を求めなさい。 以上のような問題です。 一応、私は次のように考えてみました。 まず、n期あるうちのr期(r回)でXi=aであるとすると、 f(x1+x2+…xn)=nCr p^r (1-p)^(n-r)・・・(1)となる。 また、Sn=X1+X2+…+Xn=ar-a(n-r)=2ar-anであり、 rについて解くとr=(an+Sn)/2a・・・(2)となる。 よって(2)を(1)に代入することにより、 確率関数f(Sn)=nC(an+Sn)/2a×p^(an+Sn)/2a×(1-p)^(an-Sn)/2aを得られる。 平均に関しては私の導出過程が少々回りくどいようですので、 それは割愛させていただきますが、結果はan(2p-1)となりました。 ここで、本題の質問なのですが、 １．この確率関数は正しいのでしょうか？ ２．同様にこの平均は正しいのでしょうか？ ３．積率母関数および分散はどのように導出すればよいのでしょうか？ どれか1つでも構いませんので、ご回答くださいますよう何卒よろしくお願いいたします。

  • 中心極限定理　実験

    中心極限定理に、 『母集団分布の平均、分散をμ、σ^2とすると、その分布が何であっても、nが大きければ、 Sn=X1+X2+...+Xnは、N(nμ, nσ^2)に、 Xmean=(X1+X2+...+Xn)/nは、N(μ, (σ^2)/n)に、従う。』 とあります。 これを実証しようと、エクセルで乱数を作りました。 RAND関数なので、母集団の平均μ=0.5、σ^2=0.083です。 サンプル10,000個を作成しました。 その合計Sn=S10,000のデータを30件取り、平均、分散を求めました。 平均は 5001.608 となり、 中心極限定理通りn×μに近い値になりましたが、 分散は460程度となってしまいました。 定理によれば、830付近になるとのことですが、このズレはなぜ発生するのでしょうか？ よろしくお願いします。

  • ランダムウォークのマルコフ性の証明

    Y_1 , Y_2 , …, Y_n , …　を、同分布でそれぞれ独立である確率変数列であるとし、 P(Y_i = 1) = p , P(Y_i = -1) = q　(i=1,2,3,…) p + q = 1 , 0≦p,q≦1 とします。このとき、{S_n}を、 S_0 = 0 S_n = Σ[i=1～n] Y_i としたとき、{S_n}はランダムウォークと呼ばれる確率過程に成ります。 この確率過程{S_n}のマルコフ性、つまり、 P( S_n = x+1 | S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x ) = P( S_n = x+1 | S_{n-1} = x ) を示したいのですが、以下の証明法は正しいでしょうか?何だかあっさりし過ぎていて不安なのですが…。 P( S_n = x+1 | S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x ) = P( S_n = x+1 , S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x )/P( S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x ) = P( Y_n = 1 , S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x )/P( S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x ) = P( Y_n = 1 )P( S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x )/P( S_1 = s_1 , S_2 = s_2 , … , S_{n-1} = x ) = P(Y_n = 1) = p = P( S_n = x+1 | S_{n-1} = x ) 特に3個目の等号が成立するかどうかが不安です。 Y_nはS_1,S_2,…S_{n-1}と独立だから成立するのではないか、と思うのですが、独立という概念に対する理解がまだまだ甘いからなのか、何だか不安です。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する