- หน้าแรก
- คลังความรู้
- อนุกรมฟูเรียร์
- เชิงซ้อน · สเปกตรัม
// บทที่ 14.4 · อนุกรมฟูเรียร์ · Complex Form & Spectrum
รูปเชิงซ้อน สเปกตรัมขนาด–เฟส และสมบัติของอนุกรมฟูเรียร์
สามสูตรของ a0, an, bn รวมเป็นสูตรเดียวได้ด้วยเลขชี้กำลังเชิงซ้อน รูปนี้คือภาษาของ FFT การแปลงฟูเรียร์ และการวิเคราะห์วงจรด้วยเฟเซอร์ เมื่อเข้าใจแล้วจะอ่านสเปกตรัมบนจอเครื่องวัดได้ถูก และใช้สมบัติไม่กี่ข้อหาอนุกรมของรูปใหม่ได้โดยไม่ต้องอินทิเกรต
สรุปใน 30 วินาที
- x(t)=∑n=−∞∞cnejnω0t, cn=T1∫Tx(t)e−jnω0tdt
- cn=21(an−jbn), c−n=cn∗, An=2∣cn∣, φn=∠cn
- พัลส์ดิวตี้ D: ∣cn∣=AD∣sinc(nD)∣ ศูนย์ที่ n=k/D
- เลื่อน t0: คูณ e−jnω0t0 · อนุพันธ์: คูณ jnω0 · ปริพันธ์: หาร jnω0
ระดับพื้นฐาน ม.ปลาย · ปวส. · ปี 1กลาง ป.ตรี วิศวกรรมขั้นสูง ป.โท · งานวิจัย
จากสูตรออยเลอร์สู่อนุกรมเชิงซ้อน กลาง
โคไซน์หนึ่งตัวคือเฟเซอร์สองตัวหมุนสวนกัน ตัวหนึ่งที่ความถี่ +ω อีกตัวที่ −ω ส่วนจินตภาพหักกันพอดีเหลือค่าจริง เมื่อเขียนทุกฮาร์มอนิกแบบนี้ อนุกรมทั้งก้อนกลายเป็นผลรวมของเลขชี้กำลังเพียงชนิดเดียว ทบทวนสูตรออยเลอร์ที่หน้าตรีโกณมิติ
C1 Euler's formula สูตรออยเลอร์
ejθ=cosθ+jsinθ,cosθ=2ejθ+e−jθ,sinθ=2jejθ−e−jθ
C2 Complex (exponential) Fourier series อนุกรมฟูเรียร์รูปเชิงซ้อน
x(t)=n=−∞∑∞cnejnω0t,cn=T1∫Tx(t)e−jnω0tdt
สูตรเดียวใช้ได้ทุก n รวมทั้ง n = 0 (c0 คือค่าเฉลี่ย) และ n ติดลบ ไม่ต้องแยก 1/T กับ 2/T
ความสัมพันธ์ระหว่าง cₙ, aₙ, bₙ, Aₙ และ φₙ กลาง
C3 Amplitude–phase (cosine) form รูปแอมพลิจูด–เฟส
ancosnω0t+bnsinnω0t=Ancos(nω0t+φn),An=an2+bn2,φn=atan2(−bn,an)
ใช้ atan2 (ดูเครื่องหมายทั้งสองตัว) ไม่ใช่ tan−1(b/a) ธรรมดา ซึ่งผิดในจตุภาคที่ 2 และ 3 · ตัวอย่าง bn บวก an=0 → φn=−90∘ คือ cos(⋅−90∘)=sin(⋅)
C4 cₙ from aₙ, bₙ แปลงไปมา
c0=a0,cn=2an−jbn,c−n=cn∗,∣cn∣=2An,∠cn=φn
ย้อนกลับ: an=2Recn, bn=−2Imcn
C5 Sine-referenced form รูปอ้างอิงไซน์ (นิยมในไฟฟ้ากำลัง)
Ansin(nω0t+ψn),ψn=atan2(an,bn)=φn+90∘
เมื่อแรงดันเขียนเป็น Vmsinωt มุม ψ1 ของกระแสคือมุมเฟสที่ใช้คำนวณ DPF โดยตรง
สเปกตรัมขนาดและเฟส: ด้านเดียวกับสองด้าน กลาง
ตำราคณิตศาสตร์และการแปลงฟูเรียร์วาด ∣cn∣ ทั้งความถี่บวกและลบ (สมมาตรกัน) เครื่องวัดจริงแสดงด้านเดียว และส่วนใหญ่แสดงเป็นค่า RMS หรือ dB ต้องรู้ว่าจอแสดงแบบไหนก่อนอ่านตัวเลข
C6 Reading a spectrum อ่านสเปกตรัมให้ถูกชนิด
สองด้าน: ∣cn∣ด้านเดียว (ค่ายอด): An=2∣cn∣ด้านเดียว (RMS): 2An=2∣cn∣Ln=20log101 VAn/2 dBV
FFT บนออสซิลโลสโคปแสดงด้านเดียว หน่วย dBV (RMS) เครื่องวิเคราะห์กำลังแสดงด้านเดียวเป็น A หรือ V RMS หรือ % ของฮาร์มอนิกหลัก วิธีอ่าน FFT บนสโคป
คลื่นสี่เหลี่ยม ±1: สเปกตรัมสองด้าน ∣cn∣ (บน) แบ่งพลังงานครึ่งหนึ่งไปไว้ที่ความถี่ลบ สเปกตรัมด้านเดียว An (ล่าง) คือสิ่งที่เครื่องวัดแสดง สูงเป็นสองเท่า
พัลส์และเอนเวโลป sinc: ดิวตี้กำหนดรูปสเปกตรัม กลาง
พัลส์สี่เหลี่ยมเป็นรูปคลื่นที่สำคัญที่สุดในอิเล็กทรอนิกส์กำลังและดิจิทัล ทั้ง PWM นาฬิกาลอจิก และแรงดันโหนดสวิตช์ของบั๊กคอนเวอร์เตอร์ รูปเชิงซ้อนให้คำตอบในบรรทัดเดียว
C7 Pulse train: sinc envelope พัลส์สูง A ดิวตี้ D เริ่มที่ t = 0
cn=T1∫0DTAe−jnω0tdt=ADsinc(nD)e−jπnD,sinc(u)=πusinπu,An=2AD∣sinc(nD)∣
เส้นสเปกตรัมห่างกัน f0=1/T แต่รูปเอนเวโลปขึ้นกับความกว้างพัลส์ τ=DT เท่านั้น ศูนย์แรกที่ f=1/τ พัลส์ยิ่งแคบ สเปกตรัมยิ่งกว้างและแบน
แอมพลิจูดของพัลส์ D = 0.25 (บน) และ D = 0.1 (ล่าง) เส้นประคือเอนเวโลป 2AD∣sinc(fτ)∣ แท่งที่ n = 4, 8 (บน) และ n = 10, 20 (ล่าง) หายไปเพราะตกที่ศูนย์ของ sinc
ตัวอย่าง C-1 · สเปกตรัมของ PWM
สัญญาณ PWM 0–5 V ความถี่ 20 kHz ดิวตี้ 20 % หา DC แอมพลิจูดและมุมเฟสของฮาร์มอนิก 1–4 และบอกว่าความถี่ใดไม่มีพลังงาน
c0=AD=1.00 V (ถ้าต่อ RC กรองจะได้ DC 1 V)
C7: An=2(5)(0.2)∣sinc(0.2n)∣, มุม φn=−πnD=−36∘n (ตราบที่ sinc ยังเป็นบวก)
| n | ความถี่ | sinc(0.2n) | An (V ยอด) | φn |
| 1 | 20 kHz | 0.9355 | 1.871 | −36° |
| 2 | 40 kHz | 0.7568 | 1.514 | −72° |
| 3 | 60 kHz | 0.5046 | 1.009 | −108° |
| 4 | 80 kHz | 0.2339 | 0.468 | −144° |
ศูนย์ของ sinc ที่ nD=1,2,… → n = 5, 10, 15 คือ 100, 200, 300 kHz
DC 1 V · ฮาร์มอนิก 1.87, 1.51, 1.01, 0.47 V ที่ 20–80 kHz · ไม่มีพลังงานที่ 100, 200, 300 kHz · ถ้าใช้ PWM ขับมอเตอร์หรือ LED ฮาร์มอนิกพวกนี้คือเสียงหวีดและสัญญาณรบกวน
ตรวจในแล็บ: “พัลส์ 0 ถึง A” A = 5, D = 20 % ✓
สมบัติที่ใช้แก้โจทย์โดยไม่ต้องอินทิเกรต ขั้นสูง
ถ้า x(t) มีสัมประสิทธิ์ cn แล้ว รูปที่ได้จากการบวก เลื่อน ดิฟ หรืออินทิเกรต x(t) ก็มีสัมประสิทธิ์ที่หาได้ทันที
C8 Linearity ความเป็นเชิงเส้น
αx(t)+βy(t)↔αcn+βdn
รูปซับซ้อน = ผลรวมของรูปในตาราง สัมประสิทธิ์ก็บวกกันได้ตรง ๆ
C9 Time shift เลื่อนเวลา
x(t−t0)↔cne−jnω0t0
ขนาดคงเดิม เฟสลดลง nω0t0 (ใน S6 คือข้อนี้)
C10 Differentiation อนุพันธ์
dtdx↔jnω0cn
ฮาร์มอนิกสูงถูกขยาย n เท่า ดิฟเฟอเรนชิเอเตอร์จึงขยายนอยส์ และขอบคมของสัญญาณให้สเปกตรัมกว้าง
C11 Integration (zero-mean x) ปริพันธ์
∫xdt↔jnω0cn,n=0
ต้องให้ c0=0 ก่อน ไม่อย่างนั้นอินทิกรัลโตไม่จำกัด (ไม่เป็นคาบ) ค่าเฉลี่ยของผลลัพธ์กำหนดจากเงื่อนไขเริ่มต้นแยกต่างหาก
C12 Multiplication & modulation การคูณสัญญาณ
x(t)ejmω0t↔cn−m,x(t)y(t)↔k∑ckdn−k
คูณกับไซน์ = เลื่อนสเปกตรัมไปทั้งสองข้าง (เกิดไซด์แบนด์) อธิบายฮาร์มอนิกรอบความถี่สวิตชิ่งของ PWM และวงจรชอปเปอร์ V6
ตัวอย่าง C-2 · ออปแอมป์อินทิเกรเตอร์ทำสี่เหลี่ยมให้เป็นสามเหลี่ยมป้อนสี่เหลี่ยม ±5 V 1 kHz เข้าอินทิเกรเตอร์ vo=−RC1∫vindt ที่ RC=1 ms หาแอมพลิจูดฮาร์มอนิกหลักของขาออก และเทียบกับค่ายอดของสามเหลี่ยม
สี่เหลี่ยม: A1=π4(5)=6.366 V W1
C11: แต่ละฮาร์มอนิกหารด้วย nω0RC และเลื่อนเฟส (เครื่องหมายลบกับ 1/j)
A1,out=(2π⋅1000)(0.001)6.366=6.2836.366=1.013 V
ในโดเมนเวลา: ครึ่งคาบ 0.5 ms ขาออกเปลี่ยน RC1(5)(0.5 ms)=2.5 V ยอดถึงยอด → สามเหลี่ยมยอด 1.25 V ซึ่งมีฮาร์มอนิกหลัก π28(1.25)=1.013 V ✓
ฮาร์มอนิกที่ n ของขาออกลดลงตาม 1/n2 (สี่เหลี่ยม 1/n หารด้วย n อีกครั้ง) ได้สามเหลี่ยมยอด 1.25 V ทดลองจริงได้ที่ แล็บออปแอมป์อินทิเกรเตอร์
ตัวอย่าง C-3 · ขบวนอิมพัลส์หาอนุกรมฟูเรียร์ของ δT(t)=∑kδ(t−kT)
C2 อินทิเกรตช่วง −T/2 ถึง T/2 มีอิมพัลส์ตัวเดียวที่ t = 0: cn=T1∫δ(t)e−jnω0tdt=T1
δT(t)=T1∑nejnω0t ทุกฮาร์มอนิกเท่ากันหมด สเปกตรัมแบนราบ ข้อนี้คือหัวใจของการสุ่มสัญญาณ (sampling) และที่มาของอเลียสซิง 12.13
เฟเซอร์ของแต่ละฮาร์มอนิก: สะพานสู่การวิเคราะห์วงจร กลาง
วงจร AC ที่เรียนมาใช้เฟเซอร์ที่ความถี่เดียว อนุกรมฟูเรียร์ทำให้สัญญาณที่ไม่ใช่ไซน์กลายเป็นปัญหาเฟเซอร์หลายปัญหา ความถี่ละหนึ่งปัญหา แต่ละปัญหาแก้ด้วยอิมพีแดนซ์ที่ความถี่นั้นแล้วนำมารวมกัน
C13 Harmonic phasor เฟเซอร์ของฮาร์มอนิกที่ n
xn(t)=Re{Xnejnω0t},Xn=an−jbn=2cn=An∠φn(ยอด),Xn,rms=Xn/2
แล็บในชุดนี้ใช้นิยามนี้ทุกที่ ตัวอย่างการใช้ A2: Yn=H(jnω0)Xn
แล็บ: ดูสเปกตรัมขนาดและเฟส พื้นฐาน
สลับปุ่มเหนือกราฟสเปกตรัมระหว่าง “แอมพลิจูด” “dB” และ “เฟส” เลื่อนดิวตี้ของพัลส์แล้วดูศูนย์ของ sinc ขยับ ลองพัลส์ 10 % ในโหมด dB จะเห็นเอนเวโลปชัดเจน
คำถามที่พบบ่อย
cn กับ an bn สัมพันธ์กันอย่างไร?
c0=a0, cn=(an−jbn)/2 และ c−n=cn∗ สำหรับสัญญาณจริง ขนาด ∣cn∣=An/2 แอมพลิจูดของไซน์ที่ความถี่ nf0 จึงเท่ากับ 2∣cn∣ เพราะมาจากสองเส้นคือ +nf0 และ −nf0
สเปกตรัมสองด้านมีความถี่ติดลบได้อย่างไร?
ความถี่ติดลบเป็นเรื่องทางคณิตศาสตร์ ไซน์จริงหนึ่งตัวเขียนเป็นผลรวมของเฟเซอร์หมุนทวนเข็ม ejωt กับตามเข็ม e−jωt อย่างละครึ่ง เครื่องวัดและสโคปแสดงสเปกตรัมด้านเดียว (เฉพาะความถี่บวก) ซึ่งเท่ากับสองเท่าของสองด้าน
ทำไมพัลส์ดิวตี้ 20 % ไม่มีฮาร์มอนิกที่ 5?
สเปกตรัมของพัลส์มีเอนเวโลป sinc(nD) ซึ่งเป็นศูนย์เมื่อ nD เป็นจำนวนเต็ม D = 0.2 จึงไม่มีฮาร์มอนิกที่ 5, 10, 15 ทั่วไปฮาร์มอนิกที่ k/D หายไป
สมบัติอนุพันธ์ของอนุกรมฟูเรียร์ใช้ทำอะไร?
อนุพันธ์คูณสัมประสิทธิ์ด้วย jnω0 ปริพันธ์หารด้วย jnω0 จึงหาอนุกรมของสามเหลี่ยมจากสี่เหลี่ยมได้ทันที และอธิบายว่าทำไมอินทิเกรเตอร์ทำให้ฮาร์มอนิกสูงลดลง (ขอบมน) ส่วนดิฟเฟอเรนชิเอเตอร์ขยายฮาร์มอนิกสูง (ขยายนอยส์)
มด
อาจารย์มด · ธีรวุฒิ แสวงบุญ
อาจารย์สาขาวิชาวิศวกรรมไฟฟ้า มหาวิทยาลัยราชภัฏธนบุรี สมุทรปราการ · ดูโปรไฟล์