การกระจายความน่าจะเป็นสำหรับคลื่นไซน์ที่มีเสียงดัง


12

ฉันต้องการวิเคราะห์การกระจายความน่าจะเป็นของจุดสุ่มตัวอย่างจากฟังก์ชันการแกว่งเมื่อมีข้อผิดพลาดในการวัด ฉันได้คำนวณการแจกแจงความน่าจะเป็นสำหรับส่วน "ไม่มีเสียง" แล้ว (ฉันจะใส่ท้ายนี้) แต่ฉันไม่สามารถหาวิธีรวม "เสียง" ได้

การประมาณเชิงตัวเลข

เพื่อให้ชัดเจนยิ่งขึ้นลองจินตนาการว่ามีฟังก์ชั่นซึ่งคุณสุ่มเลือกคะแนนจากในรอบเดียว หากคุณได้รับคะแนนในฮิสโตแกรมคุณจะได้รับบางสิ่งที่เกี่ยวข้องกับการแจกแจงy(x)=sin(x)

ไม่มีเสียงดังรบกวน

ตัวอย่างเช่นนี่คือและฮิสโตแกรมที่เกี่ยวข้องsin(x)

ป้อนคำอธิบายรูปภาพที่นี่

พร้อมเสียงดัง

ตอนนี้หากมีข้อผิดพลาดในการวัดบางอย่างมันจะเปลี่ยนรูปร่างของฮิสโตแกรม ตัวอย่างเช่น

ป้อนคำอธิบายรูปภาพที่นี่

การคำนวณเชิงวิเคราะห์

ดังนั้นหวังว่าฉันจะทำให้คุณมั่นใจว่ามีความแตกต่างระหว่างสองอย่างนี้ตอนนี้ฉันจะเขียนวิธีคำนวณกรณี "ไม่มีเสียง":

ไม่มีเสียงดังรบกวน

y(x)=sin(x)

จากนั้นถ้าเวลาที่เราสุ่มตัวอย่างกระจายอย่างสม่ำเสมอการแจกแจงความน่าจะเป็นสำหรับyจะต้องเป็นไปตาม:

P(y)dy=dx2π

ตั้งแต่นั้นมา

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

และอื่น ๆ

P(y)=12π1y2

ซึ่งการปรับมาตรฐานให้เหมาะสมนั้นเหมาะกับฮิสโตแกรมที่สร้างขึ้นในกรณี "ไม่มีเสียงรบกวน"

พร้อมเสียงดัง

ดังนั้นคำถามของฉันคือฉันจะวิเคราะห์เสียงรบกวนในการกระจายได้อย่างไร ฉันคิดว่ามันเหมือนกับการรวมการแจกแจงในวิธีที่ฉลาดหรือรวมถึงเสียงรบกวนในคำจำกัดความของแต่ฉันไม่มีความคิดและวิธีที่จะก้าวไปข้างหน้าดังนั้นคำแนะนำ / เคล็ดลับหรือแม้แต่การอ่านที่แนะนำจะมาก ชื่นชมy(x)

คำตอบ:


10

ขึ้นอยู่กับว่าโครงสร้างเสียงนั้นเป็นอย่างไร

สมมติว่าผมเคยเข้าใจสถานการณ์ของคุณได้อย่างถูกต้องถ้าเสียงเป็นสารเติมแต่งที่เป็นอิสระและกระจายเหมือนกันคุณก็จะใช้เวลาบิดของความหนาแน่นของเสียงที่มีความหนาแน่นของYY

ถ้าสุ่มชุดตลอดวัฏจักรกระบวนการไม่มีเสียงของคุณบนคือซึ่งเลวลงโดยมีค่าเฉลี่ยและความแปรปรวน 0 การกระจายตัวเล็กน้อยของคือการรวมกันของการแจกแจงแบบถดถอยเหล่านั้น ; ดูเหมือนว่าคุณได้ทำการแจกจ่ายออกอย่างถูกต้องแล้ว ขอเรียกว่าความหนาแน่นกรัมXixYi|Xi=xisin(xi)Yg

ตัวอย่างเช่นหากเสียงของคุณคือซึ่งจะกล่าวถึงจากนั้นคือความหนาแน่นของผลรวมของเสียงด้วยการรวมกันของตัวแปรที่ไม่มีเสียงϵiN(0,σ2)f(ϵ)=12πσexp(ϵ22σ2)fg

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

ป้อนคำอธิบายรูปภาพที่นี่

(การบิดนี้ทำด้วยตัวเลขฉันไม่รู้ว่าอินทิกรัลที่อินทิกรัลอยู่ในตัวอย่างนี้เพราะฉันไม่ได้ลอง)


สิ่งมหัศจรรย์ฉันหายไปจากความคิด "convolution" เนื่องจากตัวเลขของคุณแสดงว่านี่เป็นจุด เพิ่งลองรวมเข้าด้วยกัน ขอบคุณ
Greg

2
คุณอาจพบว่ามันยาก แต่มักจะไม่ยากที่จะประมาณผลลัพธ์ที่เป็นตัวเลข
Glen_b -Reinstate Monica

0

ฉันคิดว่านิพจน์ที่ได้รับสำหรับ P (x) นั้นปิดด้วยสองปัจจัย เวลาตัวอย่างที่กระจายแบบสม่ำเสมอนั้นเทียบเท่ากับการกระจายเฟสอย่างสม่ำเสมอในช่วง -pi, pi ฟังก์ชันตรีโกณมิติกระจายความน่าจะเป็นในช่วง y {-1,1} การรวม P (y) ในช่วงเวลานี้จะต้อง = 1 ไม่ใช่ 2 ตามที่ได้รับโดยใช้การรวมของคุณด้านบน ฉันคิดว่า P (y) = 1 / (pi Sqrt (1-y ^ 2))


อาจเป็นไปได้ว่าเป็นเหตุผลที่ฉันกล่าวว่า "ด้วยการปรับสภาพที่เหมาะสม" เนื่องจากฉันขี้เกียจเกินกว่าที่จะคิดในเวลานั้น ขอบคุณ
เกร็ก
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.