운동량을 가진 무작위 도보


18

다음 조건에서 0부터 시작하는 정수 랜덤 워크를 고려하십시오.

  • 첫 번째 단계는 동일한 확률로 플러스 또는 마이너스 1입니다.

  • 미래의 모든 단계는 60 %가 이전 단계와 같은 방향 일 가능성이 40 %는 반대 방향 일 가능성이 높습니다

이로 인해 어떤 종류의 분포가 나옵니까?

나는 운동량이 아닌 무작위 걷기가 정규 분포를 산출한다는 것을 알고 있습니다. 운동량이 분산을 변경하거나 분포의 특성을 완전히 변경합니까?

나는 일반적인 대답을 찾고 있으므로 60 %와 40 % 이상이면 실제로 p1-p를 의미합니다.


실제로 @Dilip에는 순서 쌍 및 , 로 상태가 색인화 된 Markov 체인이 필요합니다 . 전환은(i,i+1)(i,i1)iZp ( i , i + 1 ) ( i + 1 , i ) ( , - 1 ) ( -(i,i+1)(i+1,i+1)(i,i1)(i1,i) 확률과 및 확률이 및 입니다 . p(i,i+1)(i+1,i)1 - p(i,i1)(i1,i2)1p
whuber

단계 크기는 에 Markov 체인을 형성 하며 고정 분포에서 시작하기 위해 발생합니다 (?!). {1,+1}
추기경

당신은에 대한 제한 (한계) 유통 원하는 있습니까 도보의 단계는? X n{ 1 , + 1 }Sn=i=1nXnXn{1,+1}
추기경

또 다른 접근법은 기하학적 랜덤 변수의 교호 합을 살펴본 다음 일부 마틴 게일 이론을 적용하는 것입니다. 문제는 일종의 중지 시간을 정의해야한다는 것인데, 이는 까다로울 수 있습니다.
shabbychef

답변:


8

결론으로 바로 넘어 가기 위해 "모멘텀"은 정규 분포가 랜덤 워크 분포의 점근 근사라는 사실을 변경하지는 않지만 분산은 4np(1p) 에서 np/(1p) . 이것은이 특별한 경우에 비교적 기본적인 고려 사항에 의해 도출 될 수있다. Markov 체인의 유한 상태 공간에 대해 아래의 인수를 CLT로 일반화하는 것은 끔찍한 일이 아니지만 가장 큰 문제는 실제로 분산 계산입니다. 특정 문제에 대한이 아래의 주장이 독자에게 그것이 올바른 분산이라고 확신시킬 수 있기를 바랍니다.

추기경이 주석에서 제공하는 통찰력을 사용하여, 랜덤 보행은

Sn=k=1nXk
로 주어지며 여기서 Xk{1,1} 이고 Xk 는 천이 확률 행렬을 가진 Markov 사슬을 형성합니다 n ∞ 일
(p1p1pp).
때의 점근선 고려 사항의 경우 X 1 의 초기 분포 는 아무런 역할을하지 않으므로 수정하겠습니다.nX1다음 인수를 위해 X 1 = 1 이고 0 < p < 1로 가정합니다. 매끄러운 기술은 Markov 체인을 독립적 인 사이클로 분해하는 것입니다. 하자 σ 1 나타내고 처음 시간 한 후, 즉 1 행의 마르코프 체인 복귀 경우 X 2 = 1 다음 σ 1 = 2 , 그리고 만약 X 2 = X 3 = - 1 X 4 = 1 다음 σ 1 = 4X1=10<p<1σ1X2=1σ1=2X2=X3=1X4=1σ1=4. 일반적으로, i 번째 리턴 시간을 1로 나타내고 τ i = σ i - σ i - 1 이 리턴 간 시간을 나타냅니다 ( σ 0 = 1 ). 이러한 정의를 통해σiiτi=σiσi1σ0=1
  • 함께 다음 S σ N = X 1 + N Σ= 1 U I .나는=케이=σ나는1+1σ나는엑스케이
    에스σ=엑스1+나는=1나는.
  • 이후 값이 얻어 - 1 위한 K = σ I - 1 + 1 , ... , σ I을 - 1X는 σ = 1 이 보유하고 있음 U = 2 - τ .엑스케이1케이=σ나는1+1,,σ나는1엑스σ나는=1
    나는=2τ나는.
  • Markov 체인의 반환 시간 는 공식적으로 강력한 Markov 속성으로 인해 iid이며이 경우 평균 E ( τ i ) = 2 이고 분산 V ( τ i ) = 2 p입니다.τ나는이자형(τ나는)=2 . 아래 평균과 분산을 계산하는 방법이 표시됩니다.V(τ나는)=21
  • iid 변수에 대한 일반적인 CLT는
    SσnasympN(0,2np1p).
  • 마지막으로 주목해야 할 것은 약간의 믿음의 도약이 필요하기 때문입니다. 세부 사항을 생략했기 때문에 이므로 S n asymp N ( 0 , n pσn=1+i=1nτi2n
    SnasympN(0,np1p).

모멘트의 계산하도록 주의 수를 그 P ( τ 1 = 1 ) = P 및 대한 m 2 , P ( τ 1 = m ) = ( 1 - P ) (2) (P)의 m - 2 . 그리고 기하학적 분포를위한 모멘트를 계산할 때 사용 된 것과 유사한 기술이 적용될 수있다. 또는 X 가 성공 확률 1 - p 이고 Z =τ1P(τ1=1)=pm2P(τ1=m)=(1p)2pm2X1p 다음 1 + X ( 1 - Z ) τ 1 과동일한 분포를 가지며,이 후자의 표현에 대한 평균과 분산을 쉽게 계산할 수 있습니다.Z=1(τ1=1)1+X(1Z)τ1


+1 좋아요 나는 1 / √에 대한 점근 분포를 작성했을 것입니다., CLT는 일반적인 방법으로 적용된다는 것을 명확하게 표시합니다. 그러나 그것은 맛의 문제 일뿐입니다. 1/nSn
mpiktas

2

Van Belle의 'Rule of Thumb'8.7 (그의 저서 제 2 판 )은 혁신이 자기 상관 관계 가질 때 평균의 표준 오차에 대한 근사값을 포함합니다 . 이 사용하여 번역 ρ = 2 - 1 제공 의 진정한 표준 오류  ˉ X√를ρρ=21 여기서nˉxn단계후 랜덤 보행의 위치이고s는 표본 표준 편차입니다 (n무조건n,임)

의 진정한 표준 오차 엑스¯1에스,
엑스¯에스 . 결론은 대략적인 근사치로n ˉ x 의 표준 편차가 약 이어야한다고 예상합니다.1엑스¯2엑스¯ ./(1)

편집 : 나는 자기 상관이 잘못 되었습니다 (또는 는 다르게 해석되어야합니다). 이제 일관성이 있습니다.


p=0

ρ=2p1,12p
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.