pdf của một sản phẩm của hai biến ngẫu nhiên thống nhất độc lập


12

Đặt ~ U ( 0 , 2 )Y ~ U ( - 10 , 10 ) là hai biến ngẫu nhiên độc lập với các phân phối đã cho. Phân phối của V = X Y là gì?XU(0,2)YU(10,10)V=XY

Tôi đã thử tích chập, biết rằng

h(v)=y=y=+1yfY(y)fX(vy)dy

Chúng ta cũng biết rằng , fY(y)=120

h(v)=1

h(v)=120y=10y=101y12dy
h(v)=140y=10y=101ydy

Một cái gì đó cho tôi biết, có một cái gì đó kỳ lạ ở đây vì nó không liên tục ở mức 0. Xin hãy giúp đỡ.


1
Nếu đây là một câu hỏi bài tập về nhà, bạn có thể vui lòng thêm thẻ tự học không? Cảm ơn bạn!
Andy

đây có thể không phải là một RV thống nhất?
Yair Daon

Nó không giống như đồng phục. có lẽ một cái gì đó với log? Nhưng tôi không biết làm thế nào để viết nó ra vì số 0 nằm giữa các giới hạn và hàm không được xác định bằng 0.
CGO

Câu trả lời:


14

Một câu trả lời tốt, nghiêm ngặt, thanh lịch đã được đăng. Mục đích của việc này là để lấy được kết quả tương tự trong một cách mà có thể là một chút tiết lộ nhiều về cấu trúc cơ bản của . Nó cho thấy tại sao hàm mật độ xác suất (pdf) phải là số ít ở 0 .XY0


Nhiều thứ có thể được thực hiện bằng cách tập trung vào các hình thức phân phối thành phần :

  • là hai lần mộtbiến ngẫu nhiên U ( 0 , 1 ) . U ( 0 , 1 ) là một dạng tiêu chuẩn, "đẹp" đặc trưng của tất cả các phân phối đồng đều.XU(0,1)U(0,1)

  • là mười lần mộtbiến ngẫu nhiên U ( 0 , 1 ) .|Y|U(0,1)

  • Các dấu hiệu của tuân theo phân phối Rademacher: nó bằng - 1 hoặc 1 , mỗi với xác suất 1 / 2 .Y111/2

(Bước cuối cùng này chuyển đổi một biến không âm thành phân phối đối xứng quanh , cả hai đuôi của chúng trông giống như phân phối ban đầu.)0

Do đó (a) đối xứng khoảng 0 và (b) giá trị tuyệt đối của nó là 2 × 10 = 20 lần tích của hai biến ngẫu nhiên U ( 0 , 1 ) độc lập .XY02×10=20U(0,1)

Sản phẩm thường được đơn giản hóa bằng cách lấy logarit. Thật vậy, ai cũng biết rằng nhật ký âm của biến có phân phối mũ (vì đây là cách đơn giản nhất để tạo ra các biến thiên hàm mũ ngẫu nhiên), từ đó nhật ký âm của sản phẩm của hai trong số đó có phân phối tổng của hai số mũ. Các Exponential là một Γ ( 1 , 1 ) phân phối. Các bản phân phối gamma có cùng tham số tỷ lệ rất dễ thêm: bạn chỉ cần thêm tham số hình dạng của chúng. Một Γ ( 1 , 1 ) cộng với một Γ ( 1U(0,1)Γ(1,1)Γ(1,1) do đó cóphân phối Γ ( 2 , 1 ) . hậu quả làΓ(1,1)Γ(2,1)

Biến ngẫu nhiên là phiên bản symmetrized của 20 lần so với hàm mũ của tiêu cực của một Γ ( 2 , 1 ) biến.XY20Γ(2,1)

Nhân vật

Việc xây dựng PDF của từ đó của một U ( 0 , 1 ) phân phối được thể hiện từ trái sang phải, tiến hành trên đồng phục, đến mũ, với Γ ( 2 , 1 ) , với hàm mũ của tiêu cực của nó , với điều tương tự được nhân rộng bằng 20 , và cuối cùng là phiên bản đối xứng của điều đó. PDF của nó là vô hạn ở 0 , xác nhận sự gián đoạn ở đó.XYU(0,1)Γ(2,1)200

Chúng tôi có thể hài lòng để dừng lại ở đây. Chẳng hạn, việc mô tả đặc điểm này cho chúng ta một cách để tạo ra các nhận thức trực tiếp về , như trong biểu thức này :XYR

n <- 1; 20 * exp(-rgamma(n, 2, scale=1)) * ifelse(runif(n) < 1/2, -1, 1)

Phân tích này cũng cho thấy lý do tại sao pdf tăng lên . 0 Kỳ dị mà lần đầu tiên xuất hiện khi chúng ta coi là mũ của (các tiêu cực của) một phân phối, tương ứng với nhân một U ( 0 , 1 ) variate bằng nhau. Giá trị trong (nói) ε của 0 nảy sinh bằng nhiều cách, bao gồm (nhưng không giới hạn) khi (a) một trong những yếu tố nhỏ hơn ε hoặc (b) tất cả những yếu tố nhỏ hơn Γ(2,1)U(0,1)ε0ε . Đó là căn bậc hai là vô cùng lớn hơnεbản thân khiεgần0. Điều này buộc rất nhiều xác suất, với số tiền lớn hơnεεε0 , được vắt vào một khoảng thời gian dàiε. Để điều này có thể xảy ra, mật độ của sản phẩm phải trở nên lớn tùy ý ở mức0. Các thao tác tiếp theo - thay đổi kích thước theo hệ số20và đối xứng - rõ ràng sẽ không loại bỏ được điểm kỳ dị đó.εε020

XYΓ(2,1)

f(t)dt=tetdt, 0<t<.

t=log(z)dt=d(log(z))=dz/z0<z<1tz

f(t)dt=(log(z)e(log(z))(dz/z))=log(z)dz, 0<z<1.

20

log(z/20)d(z/20)=120log(z/20)dz, 0<z<20.

z|z|20202(20,0)(0,20)

fXY(z)dz=12120log(|z|/20), 20<z<20;fXY(z)dz=0 otherwise.

Cảm ơn bạn đã cố gắng làm cho nó trở nên "dễ tiếp cận hơn. Tôi vẫn thấy điều này hơi phản cảm nên tôi mới thực hiện điều này (tương tự như" mô phỏng "của Xi'an): Tạo plot( density( outer(seq(-10,10,length=10),seq(0,2,length=10), "*") ) )độ dài lên đến 100 để tránh một số vật phẩm cho mật độ phân phối giới hạn.
DWin

9

XXU(0,2)

fX(x)=12I(0,2)(x)
h(v)=12010101|y|12I(0,2)(v/y)dy
h(v)=1401001|y|I0v/y2dy+1400101|y|I0v/y2dy=1401001|y|I0v/2y10dy+1400101|y|I0v/2y10dy=140I20v010v/21|y|dy+140I20v0v/2101|y|dy=140I20v0log{20/|v|}+140I0v20log{20/|v|}=log{20/|v|}40I20v20
nhập mô tả hình ảnh ở đây

thu được như

   hist(runif(10^6,0,2)*runif(10^6,10,10),prob=TRUE,
   nclass=789,border=FALSE,col="wheat",xlab="",main="")
   curve(log(20/abs(x))/40,add=TRUE,col="sienna2",lwd=2,n=10^4)

Xin chào, cảm ơn. Tôi muốn hỏi tại sao giới hạn thay đổi từ -10 thành 10 thành -10 thành v / 2?
cgo

Có nên có một tiêu cực ở đâu đó? Cảm ơn
cgo

1
20<v<20log(20/|v|)20/|v|>1log(|v|/20)
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.