No.1ベストアンサー
- 回答日時:
blue_monkeyと言います。
すでに解決済みだと思いますが、アドバイスさせていただきます。参考にしていただければ幸いです。
【参考資料】
統計物理学 岩波講座
キーポイント確率・統計 岩波書店
「拡散方程式」、「酔歩」のキーワードで質問検索を行えば、参考になる回答がいくつかあります。
上記以外にも、統計力学、ブラウン運動、確率の参考書を探せば記述があると思います。
下記に蛇足の記述をさせていただきます。読み捨ててください。
【蛇足】
位置は
x=nx*a
y=ny*a
z=nz*a
時間は
t=N*τ
またτ毎にx,y,z方向に対して、独立に、確率0.5で、+a,-aだけ移動すると考えます。
1次元
(1)……<a*nx>=0
Δx=a*nx-<a*nx>
(2)……<Δx*Δx>=N*a^2
Nステップで位置nxにある確率をP(t,x)=P(τ*N,a*nx)とすると
(3)……P(τ*(N+1),a*nx)=P(τ*N,a*(nx+1))+P(τ*N,a*(nx-1))
を満たすことが確認でき、
右辺についてτについて1次まで、
左辺についてaについて2次まで
展開すると下記の拡散方程式が求まります。
(4)……(∂/∂t)P(t,x)=0.5*a^2/τ*(∂/∂x)(∂/∂x)P(t,x)
(4)式より拡散係数は
(5)……D=0.5*a^2/τ
(2)式の<Δx*Δx>=N*a^2=t/τ*a^2を(5)に代入すると
D=0.5*<Δx*Δx>/t
2次元の拡散方程式は、ステップNで、位置nx,nyにある確率P(τ*N,nx*a,ny*a)が
(6)……P(τ*(N+1),a*nx)
=P(τ*N,a*(nx+1),a*(ny-1))+P(τ*N,a*(nx-1),a*(ny+1))
+P(τ*N,a*(nx+1),a*(ny+1))+P(τ*N,a*(nx-1),a*(ny-1))
を満たすことが確認でき、1次元の時と同様に、τについて1次、aについて2次まで展開すると、(7)式の拡散方程式が得られます。
(7)……(∂/∂t)P(t,x)
=0.5*a^2/τ*[(∂/∂x)(∂/∂x)+(∂/∂y)(∂/∂y)]P(t,x)
拡散係数は1次元のときと同じD=0.5*a^2/τとなります。3次元についても同様な計算をおこない拡散方程式を導出すれば、拡散係数はD=0.5*a^2/τになることが予想されます(Blue_monkeyは1次元,2次元までしか導出の確認を行っていません)。
以上の拡散方程式から導出される拡散係数Dは
1次元,2次元,3次元でも全て同じく
D=0.5*a^2/τ
で記述されます。
1次元の場合の距離の分散は
(8)……<Δr*Δr>
=<Δx*Δx>
=N*a^2
=a^2*t/τ
2次元の場合の距離の分散は
(9)……<Δr*Δr>
=<Δx*Δx+Δy*Δy>
=2*N*a^2
=2*a^2*t/τ
3次元の場合の距離の分散は
(10)……<Δr*Δr>
=<Δx*Δx+Δy*Δy+Δz*Δz>
=3*N*a^2
=3*a^2*t/τ
(8),(9),(10)の結果を用いると、拡散係数Dは
1次元
(11)……D=(1/2)*<Δr*Δr>/t
2次元
(12)……D=(1/4)*<Δr*Δr>/t
3次元
(13)……D=(1/6)*<Δr*Δr>/t
質問中にあります拡散係数について考えますと
D=0.5*a^2/τ=(1/6)*(3*a^2/τ)
τ秒毎に±a移動するということから、1/τは、chibittzさんのf:拡散ジャンプ頻度に相当します。
aが移動するステップより、chibittzさんのr:ジャンプ距離に相当します。
残りのファクタ3がつじつま合わせ的にはchibittzさんのF:相関係数に対応することになりますが、1より大きいこと、また導出の上でも相関係数の意味を持たないように思われます。
chibittzさんの記載された関係式は、現在のところ、得られていません。
【その他】
2次元の酔歩問題で、
1)(r,φ)の極座標で考えたときの取り扱い。
2)1ステップで、x,y方向を独立にしているが、確率を0.25として、(0,a),(0-a),(a,0),(-a,0)としたときの取り扱い。
がどうなるのか、暇なときに考えようと思っています。
誤記、計算間違い、ウソがありましたらゴメンナサイ。
この回答へのお礼
お礼日時:2002/01/14 18:12
ご丁寧に教えていただき、ありがとうございます。
図書を探す際に「ブラウン運動」のほうで探していたためよい参考書に
めぐり合えませんでした。
そうですね、拡散や確率のほうから探せばよかったんですね。
どうもありがとうございました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
鯉はどうして水面に跳ねるので...
-
pHジャンプについて教えてくだ...
-
よく「飲み行こう」「食べ行こ...
-
ビー玉を41個使って遊ぶゲー...
-
魚は、空を飛びたいんでしょうか?
-
ペットロボットに愛着を感じますか
-
【物理】地上でのヘリコプター...
-
ヴィーナスのえくぼとは?
-
とびうおはなぜ飛ぶのか?
-
人造人間を英語で言うと・・・...
-
月の重力
-
これではいけない。農業は集団...
-
バスの中でジャンプしても空間...
-
全知全能の神になりたいんです...
-
慣性がなかったら、どうなりますか
-
dynamicsとkineticsの違いって...
-
自分のブログやHPがGoogleのロ...
-
夢のロボット・こんなロボット...
-
やで、って何ですか?ヤガラン...
-
地面にドリルで穴を掘りたい 地...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
鯉はどうして水面に跳ねるので...
-
pHジャンプについて教えてくだ...
-
「人間味」って何ですか?
-
ビー玉を41個使って遊ぶゲー...
-
よく「飲み行こう」「食べ行こ...
-
dynamicsとkineticsの違いって...
-
ワークポイントとは何ですか(...
-
『ロボットのような人』とはど...
-
中国の人が一斉にジャンプした...
-
こんな商品あったらいいな!
-
漢字を調べたいのですが
-
もしAIロボットが自分の考えで...
-
ロボットの操作について
-
話したことないけど、なんか雰...
-
走っている電車内でジャンプし...
-
AIは、いらないと思います。大...
-
慣性がなかったら、どうなりますか
-
走る電車の中と外でジャンプし...
-
コントローラビリティとはなん...
-
頭おかしい人
おすすめ情報