ทดสอบเพื่อแยกความแตกต่างเป็นระยะจากข้อมูลเกือบเป็นระยะ


9

สมมติว่าฉันมีบางฟังก์ชั่นที่ไม่รู้จักโดเมนซึ่งผมทราบว่าจะปฏิบัติตามเงื่อนไขที่เหมาะสมบางอย่างเช่นความต่อเนื่อง ฉันรู้ค่าที่แน่นอนของf (เพราะข้อมูลมาจากการจำลอง) ที่จุดสุ่มตัวอย่างบางระยะเท่ากันt_i = t_0 + iΔtด้วยi∈ \ {1, …, n \}ซึ่งฉันคิดว่าเพียงพอที่จะจับภาพทั้งหมด ด้านที่เกี่ยวข้องของfเช่นฉันสามารถสรุปได้ว่ามี extremum ท้องถิ่นมากที่สุดหนึ่งแห่งfในระหว่างจุดสุ่มสองจุด ฉันกำลังมองหาการทดสอบที่บอกฉันว่าข้อมูลของฉันสอดคล้องกับfเป็นระยะอย่างแน่นอนหรือไม่นั่นคือ∃τ: f (t + τ) = f (t) \, ∀ \, tffti=t0+iΔti{1,,n}fffτ:f(t+τ)=f(t)tด้วยความยาวของช่วงเวลาที่สามารถสะท้อนกลับได้เช่นΔt<τ<n·Δt (แต่เป็นไปได้ว่าฉันสามารถสร้างข้อ จำกัด ที่แข็งแกร่งกว่านี้ได้ถ้าจำเป็น)

จากมุมมองอื่นผมมีข้อมูลx0,,xnและกำลังมองหาการทดสอบที่ตอบคำถามไม่ว่าจะเป็นฟังก์ชั่นเป็นระยะf (เงื่อนไขดังกล่าวตอบสนอง) อยู่เช่นว่าf(ti)=xii∀ฉัน

จุดสำคัญคืออย่างน้อยfใกล้เคียงกับช่วงเวลาอย่างมาก (อาจเป็นเช่นf(t):=sin(g(t)·t)หรือf(t):=g(t)·sin(t)ด้วยg(t)g(t0)/Δt ) เท่าที่การเปลี่ยนแปลงจุดข้อมูลหนึ่งด้วยจำนวนเล็กน้อยอาจพอเพียงเพื่อให้ข้อมูลสอดคล้องกับfเป็นระยะอย่างแน่นอน ดังนั้นเครื่องมือมาตรฐานสำหรับการวิเคราะห์ความถี่เช่นการแปลงฟูริเยร์หรือการวิเคราะห์การข้ามศูนย์จะไม่ช่วยอะไรมาก

โปรดทราบว่าการทดสอบที่ฉันกำลังมองหาอาจไม่น่าจะเป็นไปได้

ฉันมีความคิดบางอย่างเกี่ยวกับวิธีการออกแบบการทดสอบด้วยตัวเอง แต่ต้องการหลีกเลี่ยงการสร้างใหม่ของวงล้อ ดังนั้นฉันกำลังมองหาการทดสอบที่มีอยู่


4
ระบุว่าคุณมีข้อมูลคุณสามารถอธิบายความหมายของการทดสอบที่ไม่เป็น "สถิติ" ได้หรือไม่ ตอนนี้คุณมีการทดสอบแบบไหน
whuber

1
โดยวิธีการที่คุณอาจต้องการเริ่มต้นที่นี่ในกรณีที่คุณกำลังมองหาการทดสอบทางสถิติของช่วงเวลา
tchakravarty

2
มีการกำหนดจุดสุ่มตัวอย่างอย่างไร เนื่องจากคุณไม่ทราบแน่ชัดว่าคืออะไรดังนั้นถ้ามีคนอื่นตัวอย่างพวกเขาจะไม่ใช้ "เวลา" ที่แตกต่างกันดังนั้นจึงได้ค่าต่างกันหรือไม่ นั่นคือความแปรปรวน อนึ่งไม่มีสิ่งดังกล่าวเป็นที่แน่นอนข้อมูลจนกว่าคุณจะได้ดำเนินการออกกำลังกายทางคณิตศาสตร์ทฤษฎีดังนั้นมันจะเป็นความคิดที่ดีที่จะอธิบายวิธีการที่คุณได้พบค่าของฉfff
whuber

2
เนื่องจาก @whuber และอะมีบากำลังขับรถไปที่คำถามนี้จะยังคงตอบยากจนกว่าจะได้คำจำกัดความที่น่าพอใจเกี่ยวกับการทดสอบเป็นระยะและ / หรือการทดสอบ เมื่อกำหนดจุดสุ่มโดยไม่มีข้อผิดพลาดจะมีฟังก์ชั่นเป็นระยะอย่างต่อเนื่องมากมาย (ใช้การกำหนดตามตัวอักษร) ซึ่งจะพอดีกับจุด มันเป็นแบบฝึกหัดง่ายๆในการแก้ไข แต่เห็นได้ชัดว่านี่ไม่ใช่คำตอบสำหรับคำถามของคุณนอกจากความจริงที่ว่าชุดการทำนายแบบสุ่มจำนวนตัวจะพอดีกับจุดอย่างสมบูรณ์แบบผ่านการถดถอยเชิงเส้น ดังนั้นเรารอด้วยลมหายใจซึ้งน้อยลงสำหรับการชี้แจงของคุณ nnn
พระคาร์ดินัล

1
สำหรับการใด ๆที่ไม่ได้เป็นหลายเหตุผลของ , ข้อมูลที่คุณมีสามารถเสมอถูกมองว่าเป็นตัวอย่างของการทำงานเป็นระยะ ๆ อย่างต่อเนื่องเป็นระยะเวลาเพราะคุณไม่มีข้อสังเกตว่าผลคูณหนึ่งของออกจากกัน สิ่งนี้นำไปสู่การสังเกตของ @ cardinal ซึ่งเป็นจำนวนที่สังเกตได้ว่าข้อสรุปนี้เล็กน้อยที่จะมีประโยชน์ แต่คุณก็ยังไม่ได้กำหนดหลักเกณฑ์ใด ๆ τΔtττ
whuber

คำตอบ:


0

ที่ผมกล่าวว่าผมมีความคิดวิธีการทำเช่นนี้ซึ่งผมรู้ได้จากการกลั่นและเขียนกระดาษเกี่ยวกับซึ่งเป็นที่เผยแพร่ขณะนี้: Chaos 25, 113,106 (2015) - preprint บน arXiv

เกณฑ์ที่ตรวจสอบเกือบจะเหมือนกับร่างในคำถาม: จากข้อมูลเก็บตัวอย่าง ณ จุดเวลาการทดสอบจะตัดสินว่ามีฟังก์ชันและ aนั้น:x1,,xnt0,t0+Δt,,t0+nΔtf:[t0,t0+Δt]τ[2Δt,(n1)Δt]

  • f(t0+iΔt)=xii{1,,n}
  • f(t+τ)=f(t)t[t0,t0+Δtτ]
  • fไม่มี extrema ในท้องที่มากไปกว่าซีเควนซ์โดยมีข้อยกเว้นที่เป็นไปได้ของ extremum สูงสุดหนึ่งอันใกล้กับจุดเริ่มต้นและจุดสิ้นสุดของแต่ละอันxf

การทดสอบสามารถแก้ไขได้สำหรับข้อผิดพลาดเล็ก ๆ เช่นข้อผิดพลาดเชิงตัวเลขของวิธีการจำลอง

ฉันหวังว่ากระดาษของฉันยังตอบว่าทำไมฉันจึงสนใจในการทดสอบดังกล่าว


-1

แปลงข้อมูลเป็นโดเมนความถี่โดยใช้การแปลงฟูริเยร์แบบไม่ต่อเนื่อง (DFT) หากข้อมูลเป็นระยะอย่างสมบูรณ์จะมีช่องความถี่หนึ่งช่องที่มีค่าสูงและถังขยะอื่นจะเป็นศูนย์ (หรือใกล้กับศูนย์ให้ดูที่การรั่วไหลของสเปกตรัม)

หมายเหตุว่ามติความถี่จะได้รับจาก{จำนวนตัวอย่าง}} ดังนั้นนี่จึงเป็นการจำกัดความแม่นยำในการตรวจจับsampling frequencyNumber of samples


1
ดังที่ฉันได้กล่าวไปแล้วในคำถามการแปลงฟูริเยร์ (อย่างน้อยที่สุดทั้งหมด) นั้นไม่แม่นยำพอที่จะตรวจจับความแตกต่างที่ฉันสนใจและจะตรวจจับความแตกต่างระหว่างและ(x) นอกจากนี้สิ่งที่คุณอ้างว่ามีไว้สำหรับข้อมูลไซน์เท่านั้น สำหรับข้อมูลอื่น ๆ ซับฮาร์โมนิกจะปรากฏขึ้น sin(x)(1+εx)·sin(x)
Wrzlprmft

-2

หากคุณรู้สัญญาณตามระยะเวลาจริงให้คำนวณ

difference=|theoretical datameasured data|

แล้วสรุปองค์ประกอบของ{} หากอยู่เหนือขีด จำกัด (พิจารณาข้อผิดพลาดจากการคำนวณเลขทศนิยม) ข้อมูลจะไม่เป็นงวดdifference


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