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 2022/08/15 15:25
- 物理学 アインシュタインの質量とエネルギーの等価性(E=mc²)って間違ってますよね? 4 2023/01/14 13:29
- 物理学 大学物理に詳しい方に質問です。 ラザフォードたちが実験で知りたかったことは衝突パラメータbと原子核の 1 2023/03/16 03:39
- その他(病気・怪我・症状) 夕食後のお散歩を終えると、意識がふわふわします 3 2023/05/23 23:28
- ウォーキング・ランニング 雪国の冬の運動 2 2022/07/20 01:01
- 犬 2歳半程になる猟犬系の雄の雑種ですが、散歩の時に引っ張るのと、興味のあるもの!、他の犬や動物を見つけ 2 2022/11/11 12:45
- 物理学 【物理】地上でのヘリコプターのホバリング 12 2023/04/12 14:58
- 犬 犬に詳しい方教えてください!愛犬(1歳7ヶ月)がリードを持つ人を選びます。 パピー3ヶ月の頃から飼っ 4 2022/05/23 18:36
- 統計学 加重最小二乗法=①「変数を自然対数変換」=②「誤差項の分散の逆数を重み付け」? 8 2022/11/26 11:15
- その他(SNS・コミュニケーションサービス) SNSでのトラブル 3 2023/01/10 20:53
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
鯉はどうして水面に跳ねるので...
-
pHジャンプについて教えてくだ...
-
ある芝生の整備をするのにロボ...
-
都市銀行店頭の店員さん全てロ...
-
dynamicsとkineticsの違いって...
-
慣性がなかったら、どうなりますか
-
よく「飲み行こう」「食べ行こ...
-
自由落下する箱の中の運動につ...
-
某ファミレスのロボット
-
リン酸と水酸化ナトリウムについて
-
ヤフー知恵袋、出禁になりました。
-
AIが悪用され、ネットで出し入...
-
お礼をしない人、質問を放置す...
-
漢字を調べたいのですが
-
どうして港湾施設のドルフィン...
-
生物と無生物の最大の違いって...
-
ビー玉を41個使って遊ぶゲー...
-
ロボットの操作について
-
皆さん、おはようございます♪ ...
-
大学の学科の志望理由のところ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
鯉はどうして水面に跳ねるので...
-
pHジャンプについて教えてくだ...
-
外から帰ってきて帰宅した時に...
-
やで、って何ですか?ヤガラン...
-
steamに登録できない
-
皆さん、おはようございます♪ ...
-
dynamicsとkineticsの違いって...
-
ビー玉を41個使って遊ぶゲー...
-
ワークポイントとは何ですか(...
-
走る電車の中と外でジャンプし...
-
『ロボットのような人』とはど...
-
こんな商品あったらいいな!
-
漢字を調べたいのですが
-
走っている電車内でジャンプし...
-
よく「飲み行こう」「食べ行こ...
-
AI心療のススメ
-
某ファミレスのロボット
-
中国の人が一斉にジャンプした...
-
ヤフー知恵袋、出禁になりました。
-
慣性がなかったら、どうなりますか
おすすめ情報