Phân phối xác suất cho một sóng hình sin ồn ào


12

Tôi đang tìm cách phân tích xác suất phân phối xác suất của các điểm lấy mẫu từ hàm dao động khi có một số lỗi đo lường. Tôi đã tính phân phối xác suất cho phần "không có tiếng ồn" (tôi sẽ đặt phần này ở cuối), nhưng tôi không thể tìm ra cách bao gồm "tiếng ồn".

Ước tính số

Để rõ ràng hơn, hãy tưởng tượng có một số hàm mà bạn chọn ngẫu nhiên các điểm trong một chu kỳ; nếu bạn bin các điểm trong biểu đồ, bạn sẽ nhận được một cái gì đó liên quan đến phân phối.y(x)=sin(x)

Không có tiếng ồn

Ví dụ ở đây là và biểu đồ tương ứngsin(x)

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

Với tiếng ồn

Bây giờ nếu có một số lỗi đo lường thì nó sẽ thay đổi hình dạng của biểu đồ (và do đó tôi nghĩ rằng phân phối cơ bản). Ví dụ

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

Tính toán phân tích

Vì vậy, hy vọng tôi đã thuyết phục bạn rằng có một số khác biệt giữa hai điều này, bây giờ tôi sẽ viết ra cách tôi tính toán trường hợp "không có tiếng ồn":

Không có tiếng ồn

y(x)=sin(x)

Sau đó, nếu thời gian mà chúng tôi lấy mẫu được phân phối đồng đều thì phân phối xác suất cho phải thỏa mãn:y

P(y)dy=dx2π

sau đó kể từ

dxdy=ddy(arcsin(y))=11y2

và vì thế

P(y)=12π1y2

mà với chuẩn hóa phù hợp phù hợp với biểu đồ được tạo ra trong trường hợp "không nhiễu".

Với tiếng ồn

Vì vậy, câu hỏi của tôi là: làm thế nào tôi có thể phân tích tiếng ồn trong phân phối? Tôi nghĩ rằng đó là một cái gì đó giống như kết hợp các bản phân phối một cách thông minh, hoặc bao gồm cả tiếng ồn trong định nghĩa của , nhưng tôi không có ý tưởng và cách để chuyển tiếp nên bất kỳ gợi ý / mẹo hoặc thậm chí đọc khuyến nghị sẽ nhiều đánh giá cao.y(x)

Câu trả lời:


10

Nó phụ thuộc vào cách quá trình tiếng ồn được cấu trúc.

Giả sử rằng tôi đã hiểu tình hình của bạn một cách chính xác, nếu tiếng ồn là phụ gia, độc lập và phân phối giống nhau, bạn sẽ chỉ cần đi chập của mật độ tiếng ồn với mật độ của .Y

Nếu là thống nhất ngẫu nhiên trong một chu kỳ, quá trình ồn ào của bạn có điều kiện trên là , là suy biến, với trung bình và phương sai 0. Phân phối biên của là một hỗn hợp đồng nhất của các phân phối suy biến đó ; có vẻ như bạn đã làm việc phân phối chính xác; Hãy gọi mật độ đó là .XixYi|Xi=xisin(xi)Yg

Ví dụ: nếu tiếng ồn của bạn là , nghĩa là , sau đó là mật độ của tổng nhiễu với hỗn hợp đồng nhất của các biến không nhiễu.ϵiN(0,σ2)f(ϵ)=12πσexp(ϵ22σ2)fg

fY+ϵ(z)=(fg)(z)=fY(y)fϵ(zy)dy=fY(zw)fϵ(w)dw

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

(phép tích chập này được thực hiện bằng số; tôi không biết tính tích phân đó có thể biến đổi như thế nào trong ví dụ này, vì tôi đã không thử nó.)


Công cụ tuyệt vời, tôi đã bỏ lỡ ý tưởng "tích chập", vì số của bạn cho thấy điều này là đúng chỗ. Chỉ cần cố gắng tích hợp. Cảm ơn
Greg

2
Bạn có thể thấy nó khó hiểu, nhưng thường không khó để ước tính kết quả bằng số.
Glen_b -Reinstate Monica

0

Tôi nghĩ rằng biểu thức dẫn xuất cho P (x) bị tắt bởi hệ số hai. Thời gian mẫu phân bố đồng đều tương đương với pha phân phối đồng đều trong khoảng -pi, pi. Hàm lượng giác phân phối xác suất trong khoảng y {-1,1}. Tích hợp P (y) trong khoảng này phải = 1, không phải 2 như thu được bằng cách sử dụng tích phân của bạn ở trên. Tôi nghĩ P (y) = 1 / (pi Sqrt (1-y ^ 2))


Cũng có thể, đó là lý do tại sao tôi tuyên bố "với sự bình thường hóa phù hợp" vì tôi quá lười biếng để nghĩ về điều đó vào thời điểm đó. Cảm ơn.
Greg
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.