เหตุใด DFT ถือว่าสัญญาณที่แปลงเป็นระยะ?


10

ในหนังสือประมวลผลสัญญาณหลายเล่มอ้างว่า DFT ถือว่าสัญญาณที่แปลงเป็นระยะ (และนี่คือเหตุผลว่าทำไมการรั่วไหลของสเปกตรัมเช่นอาจเกิดขึ้นได้)

ทีนี้ถ้าคุณดูนิยามของ DFT ก็ไม่มีข้อสันนิษฐานแบบนั้น อย่างไรก็ตามในบทความ Wikipediaเกี่ยวกับการแปลงฟูริเยร์แบบไม่ต่อเนื่อง (DTFT) มีการระบุไว้ว่า

เมื่อมีการป้อนข้อมูลตามลำดับ x[n] คือ N-periodic, Eq.2 สามารถลดลงเป็นการแปลงฟูริเยร์โดยสิ้นเชิง (DFT)

  • ดังนั้นสมมติฐานนี้เกิดขึ้นจาก DTFT หรือไม่
  • ที่จริงแล้วเมื่อคำนวณ DFT ฉันจริง ๆ แล้วทำการคำนวณ DTFT ด้วยการสันนิษฐานว่าสัญญาณเป็นระยะ?

เนื่องจาก DFT X [k] ของ x [n] เป็นช่วงแรกของ Discrete Fourier Series (DFS) ของสัญญาณคาบ xp [n] ซึ่งมีคาบแรกเป็น x [n]
Fat32

1
ดูเหมือนว่าฉันจะต้องเขียนคำตอบที่ไม่เห็นด้วยกับเรื่องนี้ DFT ถือว่าสัญญาณที่แปลงเป็นคาบเพราะมันเป็นชุดของฟังก์ชั่นพื้นฐานที่เหมาะสมกับสัญญาณที่ถูกแปลงซึ่งทั้งหมดนั้นเป็นคาบ
robert bristow-johnson

1
DFT เป็นเพียงการแสดงออกอย่างง่ายของ DFS ดังนั้นจึงมีสมมติฐานเป็นระยะ ๆ
lxg

คำตอบ:


12

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

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

(1)X[k]=n=0N1x[n]ej2πnk/N,k=0,1,,N1

จาก (1) เป็นที่ชัดเจนว่ามีเพียงตัวอย่างของ x[n] ในช่วงเวลา [0,N1]ได้รับการพิจารณาจึงไม่มีการสันนิษฐานเป็นระยะ ในทางกลับกันค่าสัมประสิทธิ์X[k] สามารถตีความได้ว่าเป็นค่าสัมประสิทธิ์ฟูริเยร์ของความต่อเนื่องเป็นระยะของสัญญาณ x[n]. สามารถเห็นได้จากการแปลงผกผัน

(2)x[n]=k=0N1X[k]ej2πnk/N

ซึ่งคำนวณ x[n] อย่างถูกต้องในช่วงเวลา [0,N1]แต่มันยังคำนวณความต่อเนื่องเป็นระยะนอกช่วงเวลานี้เนื่องจากด้านขวามือของ (2) เป็นคาบที่มีคาบ N. คุณสมบัตินี้มีอยู่ในคำจำกัดความของ DFT แต่ไม่จำเป็นต้องรบกวนเราเพราะโดยปกติเราสนใจเฉพาะช่วงเวลาเท่านั้น[0,N1].

พิจารณา DTFT ของ x[n]

(3)X(ω)=n=x[n]ejnω

เราสามารถดูได้โดยเปรียบเทียบ (3) กับ (1) ว่าถ้า x[n]เป็นลำดับที่แน่นอนในช่วงเวลา[0,N1]ค่าสัมประสิทธิ์ DFT X[k] เป็นตัวอย่างของ DTFT X(ω):

(4)X[k]=X(2πk/N)

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

เป็นคำพูดสุดท้ายโปรดทราบว่า DTFT ของความต่อเนื่องเป็นระยะ x~[n] ของลำดับที่แน่นอน x[n] สามารถแสดงในรูปของสัมประสิทธิ์ DFT ของ x[n]:

(5)x~[n]=k=x[nkN]
(6)X~(ω)=2πNk=X[k]δ(ω2πk/N)

แก้ไข:ความจริงที่ว่าx~[n] และ X~(ω)ที่ระบุข้างต้นเป็นคู่แปลง DTFT สามารถแสดงได้ดังนี้ สิ่งแรกที่ทราบว่า DTFT ของหวีแรงกระตุ้นเวลาไม่ต่อเนื่องเป็นหวี Dirac:

(7)k=δ[nkN]2πNk=δ(ω2πk/N)

เป็นลำดับ x~[n] สามารถเขียนได้ตามความเชื่อของ x[n] ด้วยหวีแรงกระตุ้น:

(8)x~[n]=x[n]k=δ[nkN]

เนื่องจากการแปลงสอดคล้องกับการคูณในโดเมน DTFT, DTFT X~(ω) ของ x~[n] ได้รับจากการคูณของ X(ω) ด้วยหวี Dirac:

(9)X~(ω)=X(ω)2πNk=δ(ω2πk/N)=2πNk=X(2πk/N)δ(ω2πk/N)

รวม (9) กับ (4) สร้างผลลัพธ์ (6).


ลูกศรลงคำตอบนี้ด้วยเหตุผลเดียวกันฉันมีคำตอบล่าสุดของ @ hotpaw2 ในคำสั่งนี้: "จาก (1) เป็นที่ชัดเจนว่ามีเพียงตัวอย่างของx[n] ในช่วงเวลา [0,N1]ได้รับการพิจารณาดังนั้นจึงไม่มีการสันนิษฐานระยะเวลา "ข้อสรุปไม่ได้ติดตามจากสถานที่ตั้ง
robert bristow-johnson

4
@ robertbristow-johnson: ใช่ ให้ฉันNตัวอย่างต่อเนื่องและฉันให้ DFT กับคุณ ฉันไม่จำเป็นต้องคิดอะไรเกี่ยวกับสัญญาณนอกช่วง[0,N1]ไม่แม้แต่มีอยู่จริง นี่เป็นสิ่งเดียวที่ฉันอ้างในประโยคนั้นและมันเป็นความจริงอย่างชัดเจน สำหรับการคำนวณ DFT ฉันไม่จำเป็นต้องรู้อะไรนอกจากค่าในช่วงเวลา[0,N1]. ไม่แน่ใจว่าคุณจะเข้าใจผิดหรืออ่านคำสั่งของฉันผิดได้อย่างไร หากเป็นปัญหาการกำหนดสูตรฉันยินดีที่จะชี้แจงประโยคของฉัน แต่เนื้อหาที่ชาญฉลาดมันเป็นเรื่องเล็กน้อย
Matt L.

4
อ่านคำตอบอื่น ๆ ด้านล่างและคำตอบของฉันที่หัวข้ออื่น มันไม่เกี่ยวกับสิ่งที่คุณคิดx[n] ด้านนอกของ 0nN1. มันเกี่ยวกับการแปลง "สมมติ" (ถ้าเราอนุญาตให้มนุษย์เปลี่ยนรูปเล็กน้อย) เกี่ยวกับx[n] ด้านนอกของ 0nN1. เราสามารถทราบได้ว่าการแปลงรูปแบบใดเมื่อเราเรียกใช้การดำเนินการในโดเมนหนึ่งที่เปลี่ยนโดเมนอื่นด้วยจำนวนเต็ม
robert bristow-johnson

@MattL (9) ควรอ่าน
=2πNk=X[k]δ(ω2πk/N)
แทน
=2πNk=X(2πk/N)δ(ω2πk/N)
jomegaA

@jomegaA: ไม่มีทั้งสองกรณี ตามที่ระบุไว้ในประโยคสุดท้ายของคำตอบของฉันผลลัพธ์สุดท้าย (6) สรุปได้จากการรวม (9) กับ (4) ดังนั้นแน่นอนX[k]=X(2πk/N)แต่ใน (9) มันมาจาก DTFT X(ω). และเกี่ยวกับปัจจัยการปรับ2π/Nแน่นอนมันต้องอยู่ที่นั่น อย่าสับสนนิพจน์ที่ใช้ω และ fพวกเขามีปัจจัยที่แตกต่างกัน
Matt L.

8

มันมาจากคำจำกัดความของสัญญาณโดเมนเวลา:

x[n]=k=0N1X[k]e2πinkN

คุณสามารถดูได้จากนิยามที่ว่า x[n]=x[n+N].
ในทางกลับกัน DFT สร้างตัวอย่าง N ของสัญญาณได้อย่างสมบูรณ์แบบ
ดังนั้นคุณสามารถสรุปได้ว่ามันจะมีความต่อเนื่องเป็นระยะ

อีกมุมมองหนึ่งน่าจะมองว่า DFT เป็นอนุกรมฟูเรียร์แบบไม่ต่อเนื่องแบบ จำกัด (อันที่จริงแล้วดูที่Discrete Fourier Series - DFS ) ซึ่งแน่นอนว่าสัญญาณเป็นระยะ (สรุปรวมของสัญญาณที่มีระยะเวลาT เป็นสัญญาณที่มีจุด T)


2
ฉันไม่เห็นว่ามันมาจากคำจำกัดความ
user10839

1
@ user10839: แค่ประเมิน x[n+N] และคุณจะเห็นว่ามันเท่ากับ x[n]. ตามที่ระบุไว้ในคำตอบ DFT เป็นเพียงชุดฟูริเยร์ของสัญญาณโดเมนเวลา ความยาวที่แน่นอนของสัญญาณโดเมนเวลาถือเป็นช่วงเวลาพื้นฐาน
Matt L.

@ user10839 เพียงแค่เสียบเข้ากับสมการ เลขชี้กำลังสามารถกำหนดได้ด้วยฟังก์ชัน Cosine และ Sine ซึ่งสามารถดูได้ในช่วงเวลาnkN.
Royi

1
DFT ไม่ใช่ DFS นี่เป็นเรื่องอวดดี แต่ DFT ให้ค่าสัมประสิทธิ์อนุกรมฟูริเยร์แก่คุณ สิ่งสำคัญคือให้สังเกตว่า DFT นั้นเหมือนกับการแปลงเชิงเส้นอื่น ๆ มันคือการคูณเมทริกซ์ เมทริกซ์เป็นแบบออโธเทนไนซ์ซึ่งทำให้ดี นอกจากนี้ยังสามารถแสดงให้เห็นว่าค่าสัมประสิทธิ์การส่งออกเท่ากันของการขยายตัวอนุกรมฟูริเยร์ที่สอดคล้องกันของข้อมูล แต่การแปลงฟูริเยร์ไม่ได้เป็นอนุกรมฟูริเยร์ (ประเภทไม่ตรงกัน: p)
Thang

@thang ฉันไม่รู้ว่าคุณหมายถึงอะไร DFT คือ DFS พวกเขาก็เหมือน ๆ กัน. มันง่ายที่จะเห็นว่า ให้ความสนใจนี่คือ Discrete Fourier Series ไม่ใช่ซีรี่ส์ Fourier (พร้อมอินทิกรัล) ดูที่นี่en.wikipedia.org/wiki/Discrete_Fourier_seriesและดูว่าเป็น DFT
Royi

5

มันเป็นสมมติฐานที่ไม่จำเป็น (และมักจะเป็นเท็จ) DFT เป็นเพียงการแปลงพื้นฐานของเวกเตอร์ จำกัด

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

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

ช่วงเวลาอาจทำให้คณิตศาสตร์ที่เกี่ยวข้องกับ DFT ไปยัง DTFT เวไนยมากขึ้น แต่ความสัมพันธ์ใด ๆ กับ DTFT อาจหรืออาจไม่จำเป็นเมื่อใช้ FFT จริงสำหรับการประมวลผลสัญญาณ (ขึ้นอยู่กับคุณสมบัติการแปลงฟูริเยร์ที่จำเป็นสำหรับการวิเคราะห์วิธีการประมวลผลเพิ่มเติม)


ลูกศรลงด้วยเหตุผลเดียวกันกับที่ลูกศรชี้ลงคำตอบล่าสุดของคุณเกี่ยวกับเรื่องนี้
robert bristow-johnson

5

ตกลงคำตอบของฉันจะค่อนข้างแตกต่างจากคำตอบอื่น ๆ คำตอบของฉันยอมรับหลักฐานของคำถามแทนที่จะปฏิเสธหลักฐานของคำถาม

เหตุผลที่ DFT "ถือว่า" สัญญาณอินพุต (สัญญาณที่จะถูกแปลงสิ่งที่ฉันถือว่า OP หมายถึงโดย "transformed สัญญาณ") เป็นระยะคือเพราะ DFT เหมาะกับการรวบรวมฟังก์ชั่นพื้นฐานของสัญญาณอินพุตซึ่งทั้งหมดนั้น เป็นระยะ ๆ

พิจารณาฟังก์ชั่นพื้นฐานที่แตกต่าง:

gk(u)uk0k<N

และได้รับ N ตัวอย่างอินพุต:

x[n]0n<N

เราสามารถใส่ผลรวมเชิงเส้นของฟังก์ชันพื้นฐานเหล่านี้ gk(n) ไปยังลำดับการป้อนข้อมูล

x[n]=k=0N1X[k]gk(n)=k=0N1X[k]nk

ด้วยการเลือกสัมประสิทธิ์อย่างรอบคอบ X[k]. กำลังคำนวณทั้งหมดX[k] ต้องแก้ N สมการเชิงเส้นด้วย Nราชวงศ์ คุณสามารถใช้การกำจัดแบบเกาส์เพื่อทำมันได้

กับ N ค่าที่ถูกต้องสำหรับ X[k] สำหรับ 0kN1เราสามารถตรวจสอบให้แน่ใจว่าผลรวมของการใช้พลังงานเหล่านี้ (ซึ่งก็คือ (N1)พหุนามลำดับที่สาม) จะประเมินตรงกับ x[n] แต่ละ n ดังนั้น 0nN1.

ตอนนี้จะเกิดอะไรขึ้นถ้าคุณใช้การรวมนั้นเพื่อก้าวข้ามช่วงเวลาของ 0nN1? คุณสามารถประเมินมันเพื่ออะไรก็ได้ n. คุณจะสังเกตเห็นว่าพฤติกรรมของฟังก์ชั่นนั้นจะเป็นของ(N1)พหุนามลำดับที่สามเพราะนั่นคือสิ่งที่มันเป็น สำหรับn ใหญ่เพียงพอเฉพาะกำลังสูงสุดที่มีค่าสัมประสิทธิ์ไม่เป็นศูนย์เท่านั้นที่จะกำหนดแนวโน้มสำหรับการประมาณค่า x[n].

ดังนั้นตอนนี้ด้วย DFT เราจึงปรับฟังก์ชั่นพื้นฐานให้เข้ากับลำดับการป้อนข้อมูลของเรา:

gk(u)1Ne+j2πku/N0k<N

x[n]=k=0N1X[k]gk(n)=1Nk=0N1X[k]e+j2πnk/N

และค่าสัมประสิทธิ์ X[k]สามารถแก้ไขได้และเป็น:

X[k]=n=0N1x[n] ej2πnk/N

ตำแหน่งของที่ 1Nเป็นเรื่องของการประชุม ฉันกำลังวางไว้ในที่ที่วรรณกรรมส่วนใหญ่ใส่1Nปัจจัย. มันอาจถูกลบออกจากx[n] สมการและใส่ใน X[k]สมการแทน หรือ "ครึ่ง" ของมัน (1N) สามารถวางด้วยสมการทั้งสอง มันเป็นเพียงเรื่องของการประชุม

แต่ที่นี่เรามีชุดฟังก์ชั่นพื้นฐานที่มีคาบเป็นระยะ N ไปที่เดิม x[n]. ดังนั้นแม้ว่าx[n] มาจากลำดับที่นานกว่านั้นไม่ได้เป็นระยะ DFT กำลังพิจารณาว่า x[n]คือผลรวมของฟังก์ชันพื้นฐานแต่ละอย่างที่มีคาบเป็นระยะN. หากคุณรวมฟังก์ชันตามช่วงเวลาจำนวนมากฟังก์ชันทั้งหมดที่มีช่วงเวลาเดียวกันผลรวมต้องเป็นงวดด้วยช่วงเวลาเดียวกัน


สำหรับน้อยทะเลาะเพิ่มเติมที่ฉันโต้แย้งความคิดที่ว่าผิวเผินไม่จำเป็นต้องขยายระยะข้อมูลที่ส่งผ่านไปมันโปรดดูคำตอบก่อนหน้านี้จากฉัน ฉันจะไม่ทำซ้ำที่นี่
robert bristow-johnson

1

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


DTFT ต่อเนื่องหรือไม่ มาทำไม
jojek

2
ผลลัพธ์ของ DTFT นั้นต่อเนื่อง (เป็นความถี่)
Deve

แน่นอน - คุณควรระบุให้ชัดเจนเพื่อหลีกเลี่ยงความเข้าใจผิดและให้สมการที่เพียงพอ
jojek

@jojek Thats จริงผมยังคิดว่าคำตอบนี้อาจจะดีขึ้นโดยสมการบาง
Deve

1
ฉันจะเพิ่มรายละเอียดเพิ่มเติมในไม่ช้า
ผู้เรียน

0

เฉพาะ DFT เท่านั้นที่ใช้งานได้จริงในโลกดิจิทัลที่ไม่ต่อเนื่องเนื่องจากมีการสันนิษฐานเป็นระยะ ๆ ในทั้งสองโดเมน (ถ้าคุณเรียกมันว่าอย่างนั้น) เพราะสัญญาณที่ไม่เป็นระยะในโดเมนหนึ่งทำให้เกิดสัญญาณต่อเนื่องในอีกโดเมนหนึ่งและคุณสามารถเก็บสัญญาณไม่ต่อเนื่องในหน่วยความจำดิจิตอล ดังนั้นคุณต้องสมมติว่าสัญญาณเป็นระยะ ๆ ในทั้งสองโดเมนเพื่อให้แยกออกจากทั้งสองโดเมน

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


0

เนื่องจากสัญญาณเป็นระยะสัญญาณเปลี่ยนเวลาจึงไม่เปลี่ยนขนาดที่แน่นอนของโดเมนความถี่

X[k]=k=0N1x[n]e2πinkN

e2πiDkNX[k]=k=0N1x[nD]e2πinkNe2πiDkN

ยังไงก็ไม่มีอะไรหยุดคุณจากการรับ FFT ของสัญญาณที่ไม่ได้เป็นระยะ ๆ แต่ก็มีการใช้งานได้จริงเพียงเล็กน้อยหากไม่มีการแปลงใด ๆ ทำงาน

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