Bağımsız bir Poisson sürecinin diğerini geçme olasılığı


9

Bu soruyu daha önce diğer yığın değişimleri üzerinde başka bir şekilde sordum, biraz repost için çok üzgünüm.

Profesöruma ve birkaç doktora öğrencisine kesin bir cevap vermeden sordum. Önce sorunu, sonra potansiyel çözümümü ve çözümümle ilgili sorunu, metin duvarı için çok üzgünüm.

Sorun:

İki bağımsız Poisson süreci varsayalım Mve , ve ile aynı aralıkta, . Zamanın herhangi bir noktasında, sonsuzluğa eğilim gösterdikçe, işleminin toplam çıktısının, artı işleminin , yani toplam çıktısından daha büyük olma olasılığı nedir . Bir örnekle göstermek için , ortalama ve köprülerinin , ortalama ve araçlarının aralık başına sırasıyla ve köprüsü üzerinden veRλRλMλR>λMMRDP(M>R+D)RMλRλMRMλR>λM. otomobilleri zaten köprüsü üzerinde sürdü, zaman içinde herhangi bir noktada toplamda daha fazla otomobilin köprüsü üzerinde daha fazla sürme olasılığı nedir .DRMR

Bu sorunu çözme yolum:

İlk önce iki Poisson süreci tanımlarız:

M(I)Poisson(μMI)R(I)Poisson(μRI)

Bir sonraki adım tanımlayan bir fonksiyon bulmaktır aralıkları, belirli bir süre geçtikten sonra . Bu, tüm negatif olmayan değerleri için nin çıktısına bağlı olması durumunda olur . Toplam çıkış ise, göstermek için olduğu daha sonra toplam çıkış daha büyük olması gerekmektedir . Aşağıda gösterildiği gibi.P(M>R+D)IM(I)>k+DR(I)=kkRXMX+D

P(M(I))>R(I)+D)=k=0n[P(M(I)>k+DR(I)=k)]

n

Bağımsızlık nedeniyle, bu, iki elementin ürünü olarak yeniden yazılabilir, burada birinci eleman Poisson dağılımının 1-CDF'si ve ikinci eleman Poisson pmf'dir:

P(M(I)>R(I)+D)=k=0n[P(M(I)>k+D)1Poisson CDFP(R(I)=k)Poisson pmf]

n

Bir örnek oluşturmak için , ve olduğunu varsayalım , aşağıda bu fonksiyonun üzerinden grafiği var :D=6λR=0.6λM=0.4I

resim açıklamasını buraya girin

Bir sonraki adım, bunun herhangi bir zamanda gerçekleşme olasılığını bulmak, buna . Benim düşüncem, bunun 1 eksi asla üzerinde olma olasılığını bulmaya eşdeğer olduğudur . Yani sonsuzluğa , koşullu olanı da önceki tüm değerleri için de geçerlidir .QMR+DNP(R(N)+DM(N))N

P(R(I)+DM(I)) aynıdır, g (I) fonksiyonu olarak tanımlayalım:1P(M(I)>R(I)+D)

g(I)=1P(M(I)>R(I)+D)

Olarak sonsuza eğilimi, bu da fonksiyon üzerinde geometrik integrali olarak yeniden yazılabilir .Ng(I)

Q=1exp(0Nln(g(I))dI)

Q=1exp(0Nln(1P(M(I)>R(I)+D))dI)

N

Yukarıda fonksiyonuna sahibiz .P(M(I)>R(I)+D)

Q=1exp(0Nln(1k=0n[P(M(I)>k+D)1Poisson CDFP(R(I)=k)Poisson pmf])dI)

N

n

Şimdi bana bunu bana son değerini vermelidir verilen herhangi için, , ve . Bununla birlikte, bir sorun var, lambdasları istediğimiz gibi yeniden yazabilmeliyiz çünkü önemli olan tek şey birbirleriyle orantılı olması. Daha önce , ve ile örnek oluşturmak için, aralıkları bölündüğü sürece bu etkili bir şekilde , ve . 10. Yani her 10 dakikada bir 10 araba, her dakikada 1 araba ile aynıdır. Ancak bunu yapmak farklı bir sonuç doğurur. ,QDλRλMD=6λR=0.6λM=0.4D=6λR=0.06λM=0.04D=6λR=0.6 ve , verir ve , ve , verir . Anında gerçekleşme, ve iki sonucun grafiklerini karşılaştırırsak bunun nedeni oldukça basittir, aşağıdaki grafik , ve .λM=0.4Q0.5856116D=6λR=0.06λM=0.04Q0.99985071(10.5856116)10=0.9998507D=6λR=0.06λM=0.04

resim açıklamasını buraya girin

Görüldüğü gibi, olasılık değişmemektedir, ancak şimdi aynı olasılığa ulaşmak için on kat daha fazla zaman almaktadır. Olarak fonksiyonunun aralığına bağlıdır, bu doğal olarak bir etkisi söz konusudur. Bu açıkça bir şeyin yanlış olduğu anlamına gelir, çünkü sonuç başlangıç ​​lambdama bağlı olmamalıdır, özellikle aralık ve kadar doğru olan başlangıç ​​lambda olmadığı için ve veya ve vb. buna göre ölçeklendirilir. Kolayca olasılığını ölçeklendirilebilir iken Bu nedenle, örneğin, giden ve için veQ0.040.060.40.611.50.40.60.040.06 , olasılığı 10 faktörle ölçeklemekle aynıdır. Bu açıkça aynı sonucu verir, ancak bu lambdaların hepsi eşit derecede geçerli başlangıç ​​noktaları olduğundan, bu doğru değildir.

Bu etkiyi göstermek için bir fonksiyonu olarak grafikledim ; burada , lambdaların ölçeklendirme faktörüdür, başlangıç ve . Çıktı aşağıdaki grafikte görülebilir:QttλM=0.4λR=λM1.5

resim açıklamasını buraya girin

Burada takıldığım yer, bana göre yaklaşım iyi ve doğru görünüyor, ancak sonuç açıkça yanlış. İlk düşüncem, bir yerlerde temel bir yeniden ölçeğin eksik olduğudur, ancak hayatım boyunca nerede olduğunu anlayamıyorum.

Okuduğunuz için teşekkür ederiz, her türlü yardım büyük beğeni topluyor.

Ayrıca, R kodumu isteyen varsa lütfen bize bildirin, ben de yükleyeceğim.


MathJax kodunuzu oldukça kapsamlı bir şekilde temizledim. Bir göz atarsanız, standart ve doğru kullanım hakkında birkaç şey görürsünüz. (Daha fazla iş yapılabilir; belki daha sonra.)
Michael Hardy

Müthiş! Çok teşekkür ederim, bundan habersizdim, izlemem gereken özel bir rehber var mı?
hayır nein

Yaptıklarınıza uygun olarak bazı ekstra şeyler düzenledim.
hayır nein

@nonein Düzenleme yardımında ufacık bir bit var, ancak bunun ötesinde matematik var.SE'nin MathJax temel öğretici ve hızlı referansı . LaTeX'te matematik yazmayla ilgili kılavuzlar (google için kolay olan), oradaki hızlı referansta yer almayan bir şey bulmaya çalışıyorsanız (şu anda MathJax'ın alt kümesinin oldukça kapsamlı bir kapsamına sahip olsa da) genellikle yardımcı olur.
Glen_b

Yanıtlar:


3

Süreçlerin kolektif zamanları olsun T=(0=t0<t1<t2<). Bunlar bağımsız Poisson süreçleri olduğu için, bu zamanların her birinde neredeyse kesinlikle bunlardan biri gözlenmektedir. İçini>0, tanımlamak

B(i)={+1if R(ti)=11if M(ti)=1

ve biriktirir B(i) sürece W: yani, W(0)=0 ve W(i+1)=W(i)+B(i) hepsi için i>0. W(i) kaç kez daha sayar R daha ortaya çıktı M hemen sonra ti.

Şekil: simülasyon

Bu şekil, R (kırmızı) ve M(orta mavi renkte) üstte "halı parselleri" olarak. Noktalar,(ti,W(i)). Her kırmızı nokta fazlalıktaki bir artışı temsil ederR(ti)M(ti) her mavi nokta fazlalıkta bir azalma gösterir.

İçin b=0,1,2,, İzin Vermek Eb şansın en az biri Wi küçük veya eşit b ve bırak f(b) olasılığı olabilir.

Soru soruyor f(D+1).

İzin Vermek λ=λR+λM. Bu, birleştirilmiş süreçlerin oranıdır. W Binom rastgele bir yürüyüş, çünkü

Pr(B(i)=1)=λRλ and Pr(B(i)=1)=λMλ.

Böylece,

Cevap, bu Binom rasgele yürüyüş emici bir engelle karşılaşma şansına eşittirWD1.

Bu şansı bulmanın en temel yolu,

f(0)=1

çünküve tüm iki olası sonraki adımı tekrar tekrarW(0)=0;b>0,±1

f(b)=λRλf(b+1)+λMλf(b1).

Varsayarak eşsiz çözümler olduğuλRλM,b0

f(b)=(λMλR)b,

bunu yukarıdaki tanımlayıcı denklemlere takarak kontrol edebilirsiniz. Böylece,

Cevap

Pr(ED+1)=f(D+1)=(λMλR)D+1.

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.