การกำหนด "ความขาว" ของเสียงรบกวน


14

คนเราคิดว่าเสียง "สีขาว" เป็นอย่างไร? มีมาตรการทางสถิติใด ๆ หรือมาตรการอื่น ๆ (เช่น FFT) ที่สามารถวัดจำนวนว่าเสียงใดเป็นตัวอย่างที่ใกล้เคียงกับเสียงสีขาว?


3
คุณสนใจคำแนะนำเกี่ยวกับวิธีเปรียบเทียบแหล่งสัญญาณ / สัญญาณรบกวนที่แตกต่างกันหรือคุณกำลังมองหาตัวชี้วัด "มาตรฐานอุตสาหกรรม" ที่ใช้กับจำนวน "สี" ในแหล่งเสียง ฉันไม่ทราบถึงตัวชี้วัดทั่วไปที่ใช้ แต่คุณสามารถเปรียบเทียบปริมาณการระบายสีได้โดยดูที่การกระจายพลังงานเสียงรบกวนใน FFT หรือ PSD (ประจบประแจง = ขาว) หรือคุณสามารถเปรียบเทียบ fucitons autocorrelation (แคบลง = ประจบ)
2718

ถ้าฉันเข้าใจคุณถูกต้องคุณกำลังมองหาเครื่องคำนวณกล่องดำอัตโนมัติที่มี 'ความขาว' ใช่ไหม?
Spacey

1
+1 สำหรับการคำนวณPower Spectral Densityของต้นทาง สำหรับเร็กคอร์ดนั้นฉันต้องการที่จะเพิ่มสัญญาณรบกวนสีขาวที่ไม่สามารถสุ่มตัวอย่างในทางปฏิบัติได้เนื่องจาก PSD นั้นแบนใน-∞ <f <∞
Serge

@ Mohammad - ไม่จำเป็นต้องมีกล่องดำให้คำนวณ ฉันแค่อยากรู้ว่ามันมีตัวประมาณค่าทางคณิตศาสตร์ของความขาวหรือไม่
Kitchi

1
@ BruceZenone - สำหรับตัวอย่างจริงของข้อมูลดังที่ Serge ชี้ให้เห็นว่า PSD จะไม่แบนอย่างสมบูรณ์ใช่ไหม? แต่ฉันก็ยังคาดเดาได้ว่าเสียงประจบมันยิ่งเข้ามาใกล้กับเสียงสีขาว "จริง"
Kitchi

คำตอบ:


12

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

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

สิ่งนี้อาจถูกนำไปใช้ในสลิปดังนี้

การใช้งานฟังก์ชั่นนี้กับลำดับสัญญาณรบกวนสองชุด: เสียงรบกวนสีขาวและเสียงรบกวนสีขาวที่ถูกกรองเบา ๆ จากนั้นจึงแสดงผลการลงจุดต่อไปนี้ สคริปต์สำหรับการสร้างการรับรู้ลำดับเสียงแต่ละครั้งจะสิ้นสุดลง

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

ค่าเฉลี่ยของสถิติสำหรับเสียงสีขาวคือ 9.79; ค่าเฉลี่ยของสถิติสำหรับเสียงที่ถูกกรองคือ 343.3

เมื่อมองดูตารางไคสแควร์เพื่ออิสรภาพ 10 องศาเราจะได้รับ:

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

และเราเห็นว่าไม่มีระดับนัยสำคัญที่ 9.79 (ในตาราง) ว่าเสียงสีขาวไม่ใช่สีขาว เรายังเห็นว่าค่าของ 343.3 นั้นมีแนวโน้มที่จะไม่ใช่แบบสีขาว (เปรียบเทียบกับค่า 23.2093 ในคอลัมน์นัยสำคัญ)พี=0.01


function R = whiteness_test(x,m)
    N = length(x);
    XC = xcorr(x);
    len = length(XC);
    lags = len/2+1 + [1:m];
    R = N*sum(XC(lags).^2)/XC(len/2+1).^2;
 endfunction

X = rand(1,1000,'normal');
Y = filter(1,[1 -0.5],X)
R = [R; whiteness_test(X,10)];
R2 = [R2; whiteness_test(Y,10)];

ฉันไม่ได้เป็นสถิติขนาดใหญ่ ... แต่ฉันมีความกังวลเกี่ยวกับความถูกต้องทั่วไปของการทดสอบ metioned ข้างต้นสำหรับกระบวนการเสียงสีขาวที่ไม่ใช่เกาส์: เท่าที่ฉันเข้าใจเสียงสีขาวเท่านั้นหมายความว่าไม่มีความสัมพันธ์ในเวลาและ autocorrelation คือ แรงกระตุ้นที่ 0 lag สีขาวไม่ได้แปลว่าแอมพลิจูดโดยปกติจะมีการกระจายซึ่งเป็นสิ่งที่การทดสอบถือว่า ... ดังนั้นเท่าที่ฉันเข้าใจว่าการทดสอบนั้นใช้ได้กับเสียงเกาส์เซียนสีขาว (เพราะผลรวมของการแจกแจงแบบเกาส์ Squared คือ Chi-squared) สำหรับเสียงสีขาวทั่วไป? ฉันพูดถูกหรือมีอะไรเหรอ
เฟเบียน

@Fabian: ใช่และไม่ใช่ คุณถูกต้องในการทดสอบถือว่าค่าความสัมพันธ์อัตโนมัติเป็นแบบเกาส์เซียน ถ้าเสียงรบกวนดั้งเดิมนั้นเกี่ยวกับการกระจายตัวใด ๆ ทฤษฎีบทขีด จำกัด กลางก็หมายความว่าการกระจายตัวของการประมาณค่าสัญญาณอัตโนมัติจะเป็นแบบเกาส์เซียน มีบางกรณีที่พยาธิวิทยาที่ค่าสัมประสิทธิ์ออโตคอร์เรชั่นจะไม่เป็นแบบเกาส์เซียน แต่โดยทั่วไปแล้วจะมีอยู่ไม่มากและอยู่ไกล (และอาจเป็นการวิเคราะห์แบบออโตคอร์เรชั่นไม่ใช่สิ่งที่ดีที่สุดในกรณีเหล่านั้น)
Peter K.

@ PeterK จะไม่มีการทดสอบที่ยากขึ้นเพื่อกำหนดความเรียบของ PSD ใช่ไหม ด้วยวิธีนี้ไม่มีการตั้งสมมติฐานและการกระจายตัวของตัวอย่างเสียงไม่เกี่ยวข้อง
Envidia

@Envidia: ทั้งสองเทียบเท่ากันใช่มั้ย PSD เป็นเพียง DFT ของลำดับความสัมพันธ์อัตโนมัติ
Peter K.

@PeterK ในตัวอย่างของคุณใช่พวกมันเทียบเท่ากันเป็นหลัก อย่างไรก็ตามโพรซีเดอร์จะถือว่า iid โดยทั่วไปตัวอย่างสามารถกระจายในรูปแบบใดก็ได้ ฉันเข้าใจว่าทฤษฎีบทขีด จำกัด กลางเข้ามามีบทบาทและใช้ได้ดังนั้นเหตุใดฉันจึงใช้คำว่า "ยากขึ้น" บางทีคำที่ดีกว่าน่าจะเป็น "ทั่วไป"
Envidia

5

ฉันจะใช้คุณสมบัติความสัมพันธ์อัตโนมัติของสัญญาณหรือความเรียบของ PSD เพื่อพิจารณาสิ่งนี้ ความสัมพันธ์เชิงทฤษฎีของสัญญาณรบกวนสีขาวในทางทฤษฎีนั้นเป็นแรงกระตุ้นที่เกิดความล่าช้า 0 นอกจากนี้ PSD ของการแปลงฟูริเยร์ของฟังก์ชั่น autocorrelation

สิ่งเหล่านี้ควรให้ความคิดที่ดีเกี่ยวกับความขาวของเสียงของคุณ


3

ความขาวนั้นเทียบเท่ากับความเป็นอิสระ

คุณสามารถดู diehard https://en.m.wikipedia.org/wiki/Diehard_tests

เล่มที่ 2 ของอัลกอริธึม Seminumerical ของ Knuth มีส่วนเกี่ยวกับเครื่องกำเนิดหมายเลขแบบสุ่มและการทดสอบ

ปัญหาของการทดสอบโดยใช้ DFT คือว่ามีการรั่วไหลของสเปกตรัมเล็กน้อยเทคนิคแนะนำความสัมพันธ์บางอย่างซึ่งถ้าคุณทำให้การแปลงของคุณ "ยาว" มักจะถูกทอดทิ้ง

มีการทดสอบบิตสตรีมแบบสุ่มเช่นกันที่ NIST


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