Xác suất của một quá trình Poisson độc lập vượt qua một quá trình khác


9

Tôi đã hỏi câu hỏi này trước đây trong một thời trang khác trên các stackexchanges khác, rất xin lỗi vì đã hơi repost.

Tôi đã hỏi giáo sư của tôi và một vài sinh viên tiến sĩ về, mà không có câu trả lời dứt khoát. Trước tiên tôi sẽ nêu vấn đề, sau đó là giải pháp tiềm năng của tôi và vấn đề với giải pháp của tôi, rất xin lỗi vì bức tường văn bản.

Vấn đề:

Giả sử hai độc lập Poisson xử lý R , với λ Rλ M cho cùng một khoảng thời, tùy thuộc vào λ R > λ M . Xác suất mà tại bất kỳ thời điểm nào, khi thời gian có xu hướng vô cùng, thì sản lượng tổng hợp của quá trình M lớn hơn sản lượng tổng hợp của quá trình R cộng với D , tức là P ( M > R + D ) . Để minh họa bằng một ví dụ, giả sử hai cầu RM , trung bình λMRλRλMλR>λMMRDP(M>R+D)RM λ M ô tô lái xe qua cầu R M tương ứng cho mỗi khoảng thời gian, và λ R > λ M . D xe đã lái xe qua cầu R , xác suất mà tại bất kỳ điểm nào trong thời gian hơn xe trong tổng số đã điều khiển qua cầu là những gì M hơn R .λRλMRMλR>λMDRMR

Cách giải quyết vấn đề này của tôi:

Đầu tiên, chúng tôi xác định hai quá trình Poisson:

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

Bước tiếp theo là tìm một hàm mô tả sau một số khoảng I nhất định . Điều này sẽ xảy ra trong trường hợp M ( I ) > k + D có điều kiện trên đầu ra của R ( I ) = k , cho tất cả các giá trị không âm của k . Để minh họa, nếu sản lượng tổng hợp của RX sau đó sản lượng tổng hợp của M cần phải được lớn hơn X + D . Như hình dưới đây.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

Do tính độc lập, điều này có thể được viết lại dưới dạng sản phẩm của hai phần tử, trong đó phần tử đầu tiên là 1-CDF của bản phân phối Poisson và phần tử thứ hai là phần tử Poisson:

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

n

Để tạo một ví dụ, giả sử , λ R = 0,6λ M = 0,4 , dưới đây là đồ thị của hàm trên tôi :D=6λR=0.6λM=0.4I

nhập mô tả hình ảnh ở đây

Bước tiếp theo là tìm xác suất xảy ra điều này tại bất kỳ thời điểm nào, cho phép cuộc gọi mà . Suy nghĩ của tôi là điều này tương đương với việc tìm kiếm 1 trừ khả năng M không bao giờ là trên R + D . Tức là chúng ta hãy N tiếp cận vô cực gì là P ( R ( N ) + D M ( N ) ) có điều kiện về vấn đề này cũng là đúng cho tất cả các giá trị trước đó của N .QMR+DNP(R(N)+DM(N))N

giống với 1 - P ( M ( I ) > R ( I ) + D ) , hãy xác định đó là hàm g (I):P(R(I)+DM(I))1P(M(I)>R(I)+D)

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

Khi có xu hướng vô cùng, điều này cũng có thể được viết lại dưới dạng tích phân hình học trên hàm g ( I ) .Ng(I)

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

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

N

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

QDλRλMD=6λR=0.6λM=0.4D=6λR=0.06λM=0.04D=6λR=0.6λM=0.4Q0.5856116D=6λR=0.06λM=0.04Q0.99985071(10.5856116)10=0.9998507D=6λR=0.06λM=0.04

nhập mô tả hình ảnh ở đây

Q0.040.060.40.611.50.40.60.040.06

QttλM=0.4λR=λM1.5

nhập mô tả hình ảnh ở đây

Đây là nơi tôi đang bị mắc kẹt, với tôi cách tiếp cận có vẻ tốt và đúng, nhưng kết quả rõ ràng là sai. Suy nghĩ ban đầu của tôi là tôi đang thiếu một quy mô lại cơ bản ở đâu đó, nhưng tôi không thể tìm ra cuộc sống của mình ở đâu.

Cảm ơn bạn đã đọc, bất kỳ và tất cả sự giúp đỡ được đánh giá rất cao.

Ngoài ra, nếu bất cứ ai muốn mã R của tôi xin vui lòng cho tôi biết và tôi sẽ tải nó lên.


Tôi đã thực hiện một số dọn dẹp khá rộng rãi mã MathJax của bạn. Nếu bạn xem bạn sẽ thấy một vài điều về cách sử dụng tiêu chuẩn và phù hợp. (Nhiều việc có thể được thực hiện; có thể sau này.)
Michael Hardy

Tuyệt vời! Cảm ơn bạn rất nhiều, tôi đã không biết điều đó, có hướng dẫn cụ thể nào tôi nên làm theo không?
không có nein

Tôi đã chỉnh sửa một số điều bổ sung phù hợp với những gì bạn đã làm.
không có nein

@nonein Có một chút thiếu sót trong trợ giúp chỉnh sửa, nhưng ngoài đó còn có toán học. Hướng dẫn cơ bản về MathJax & tham khảo nhanh . Các hướng dẫn về cách viết toán học trong LaTeX (rất dễ google) thường giúp ích nếu bạn đang cố gắng tìm thứ gì đó không có trong tài liệu tham khảo nhanh ở đó (mặc dù hiện tại nó đã có một phạm vi bao quát khá đầy đủ về tập hợp con của MathJax).
Glen_b -Reinstate Monica

Câu trả lời:


3

T=(0=t0<t1<t2<).i>0,

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

B(i)W:W(0)=0W(i+1)=W(i)+B(i)i>0. W(i)RMti.

Hình: mô phỏng

RM(ti,W(i))R(ti)M(ti)

b=0,1,2,,EbWibf(b)

f(D+1).

λ=λR+λM.W

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

Như vậy

WD1.

Cách cơ bản nhất để tìm cơ hội này quan sát rằng

f(0)=1

W(0)=0;b>0,±1

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

λRλM,b0

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

như bạn có thể kiểm tra bằng cách cắm cái này vào các phương trình xác định đã nói ở trên. Như vậy

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

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.