如何研究變步長隨機遊走?

考慮R上的隨機遊走,每一步有1/2的概率選擇往左走,有1/2的概率往右走.但是第n步走的距離為1/n.當n-&>∞時,終點的概率分布是否存在?如果存在如何求出?

感覺一個調和級數一半隨機取負之後還是會發散的,終點的分布應該不存在...

但是跑模擬跑出來是這個樣子的,跑了1000步,100W次試驗,好像形成了一個類似鍾型的曲線,是因為發散比較慢沒有表現出來嗎?還是確實能收斂?

data=Inner[Times,RandomChoice[{0.5,0.5}-&>{-1,1},{100000,100}],1./Range[100]];
Histogram[data,40,"Probability",ChartStyle-&>60]
FindDistribution[data,TargetFunctions-&>{NormalDistribution}]


極限分布的存在性可以用下面的定理來證明(該定理是證明Kolmogorovs three-series theorem的中間步驟).

顯然現在條件(1)是滿足的,所以隨機級數almost surely收斂.

至於極限分布如何求,也許可以考慮用特徵函數,但是我覺得應該沒有解析解.


我是題主...回來填坑....http://www.stat.ualberta.ca/people/schmu/preprints/rhs.pdf

這個問題屬於一類更加龐大的問題,隨機級數:[sumlimits_{n = 0}^infty {{s_n}{a_n}{x^n}} ]

其中s_n 是個隨機變數,這類級數有一套自己的審斂法則.

有空寫個專欄...

===========================

原題用生成函數描述就是求下列級數的係數展開

[prodlimits_{n > 0} {left( {frac{1}{2}{z^{frac{1}{n}}} + frac{1}{2}{z^{ - frac{1}{n}}}}<br />
ight)} ] 展不開的.....

從特徵函數來說就是求:[f(t) = prodlimits_{n = 1}^infty {cos } left( {frac{t}{n}} 
ight)] 的逆傅里葉變換...

但還是沒法求啊,於是就要用下面的目測大法了:

[prodlimits_{n = 1}^infty {cos } left( {frac{t}{n}} 
ight) sim {e^{ - frac{{{t^2}}}{6}}}{	ext{sinc}}(2t)]

我也不知道為什麼可以這樣....我也很絕望啊......

然後做逆傅里葉變換...

[P(x) = frac{1}{2}left( {frac{1}{2}sqrt {frac{pi }{2}} {	ext{erf}}left( {sqrt {frac{3}{2}} x + sqrt 6 } 
ight) - frac{1}{2}sqrt {frac{pi }{2}} {	ext{erf}}left( {sqrt {frac{3}{2}} (x - 2)} 
ight)} 
ight)]

不知道為什麼歸一化常數總是要調整一下...

=================================

連續型好像簡單一點.

[{x_i} sim Uleft[ { - frac{1}{n},frac{1}{n}} 
ight]]

其特徵函數為: [{f_n} = frac{n}{t}sin frac{t}{n}]

但是說的我好像會求 [prodlimits_{n = 0}^infty {frac{n}{t}sin frac{t}{n}} ] 一樣...

於是我就描點作圖想找個函數擬合...

然後群里大佬和我說這個圖像目測長的很像[f(t) = {e^{ - frac{{{t^2}}}{{{pi ^2}}}}}{	ext{sinc}}(t)]

神?目測..........

Ok,那就變換回概率分布函數...求完傅里葉後就是

[P(x) = frac{1}{2}sqrt {frac{pi }{2}} left( {{	ext{erf}}left( {frac{pi }{2}(x + 1)} 
ight) - {	ext{erf}}left( {frac{pi }{2}(x - 1)} 
ight)} 
ight)]

這個擬合程度Perfect吶...看來我只要躺屍就可以了.........


把特徵函數寫開之後歸結為要求prod_{n=1}^infty cos(t/n) 。試著找解析解未果之後,Google "Harmonic Series With Random Signs" 有不少結果,不過還沒翻到有解析解。

Byron Schmuland. Random Harmonic Series. The American Mathematical Monthly, Vol. 110, No. 5 (May, 2003), pp. 407-416:

https://www.jstor.org/stable/pdf/3647827.pdf

Kent E. Morrison. Cosine products, Fourier transforms, and random sums: https://arxiv.org/pdf/math/0411380.pdf


當且僅當向左向右的概率一樣的時候不會發散。


把你這個問題簡化成連續型想了一下,感覺n步之後的分布大概是一個標準差為1/sqrt{n} 的高斯分布,無窮步後的話就分布不存在了。


推薦閱讀:

一個醉漢,從原點開始每秒50%的概率向左走一米,50%概率向右走一米,問一個小時後,醉漢大概率離原點多遠?
圖靈斑圖:生命圖案的奧秘
可以解釋一下密碼學中什麼叫雙線性配對嗎?
西爾維婭《美麗心靈》中,納什提出的一個關於π的小數部分形成的子數列的問題?
如果博弈雙方中的一方的收益變成之前的10%,納什均衡會不會改變?

TAG:數學 | 隨機過程 |