- หน้าแรก
- คลังความรู้
- การแปลงลาปลาซ
- โพล-ซีโร่และ ROC
// บทที่ 8.2 · การแปลงลาปลาซ · Poles, Zeros & ROC
แผนภาพโพล-ซีโร่และบริเวณลู่เข้า (ROC) ของการแปลงลาปลาซ
X(s) บอกว่าน้ำหนักของ est แต่ละตัวเป็นเท่าไร แต่บอกไม่ครบ ต้องมี ROC บอกว่าอินทิกรัลลู่เข้าที่ไหน หน้านี้พาหาโพล ซีโร่ และ ROC ทีละขั้นจากนิยาม สรุปสมบัติ 8 ข้อที่ใช้หา ROC โดยไม่ต้องอินทิเกรต และอ่านเสถียรภาพกับรูปคลื่นจากตำแหน่งโพลได้ในพริบตา สูตรใช้รหัส R
สรุปใน 30 วินาที
- โพล = รากของตัวส่วน (× บนระนาบ) ซีโร่ = รากของตัวเศษ (○) อันดับ = จำนวนซ้ำ
- ROC เป็นแถบแนวตั้งเสมอ ไม่มีโพลข้างใน ขอบเขตผ่านโพล
- ด้านขวา → ขวาของโพลขวาสุด · ด้านซ้าย → ซ้ายของโพลซ้ายสุด · สองด้าน → แถบระหว่างโพล
- ROC ครอบแกน jω ⇔ เสถียร · เหตุผล + เสถียร ⇔ โพลทุกตัวอยู่ครึ่งซ้าย
ระดับพื้นฐาน ปี 2กลาง ป.ตรี วิศวกรรมขั้นสูง สัญญาณและระบบ · ป.โท
โพลและซีโร่คืออะไร พื้นฐาน
ฟังก์ชันตรรกยะใด ๆ แยกตัวประกอบได้เป็นผลคูณของ (s−zi) บนเศษและ (s−pk) บนส่วน รู้ตำแหน่งราก + ค่าคงที่หนึ่งตัวก็รู้ฟังก์ชันทั้งหมด แผนภาพโพล-ซีโร่จึงเป็น “ลายมือชื่อ” ของระบบ
R1 Pole–zero–gain form รูปโพล-ซีโร่-เกน
X(s)=D(s)N(s)=K(s−p1)(s−p2)⋯(s−pn)(s−z1)(s−z2)⋯(s−zm)
ซีโร่ zi: X(zi)=0 · โพล pk: ∣X(s)∣→∞ · K = อัตราส่วนสัมประสิทธิ์นำหน้า · ซีโร่ที่ s=0 หรือโพลที่ s=0 มักพิมพ์ว่า s เฉย ๆ
R2 Proper / improper เศษส่วนแท้และไม่แท้
m<n: แท้อย่างเข้ม,m=n: แท้,m>n: ไม่แท้
ระบบกายภาพเกือบทั้งหมดเป็นแท้ m≤n ส่วนไม่แท้มีอิมพัลส์และอนุพันธ์ของอิมพัลส์ในคำตอบ
R3 Conjugate symmetry สมมาตรสังยุค
สัมประสิทธิ์จริง⇒p เป็นโพล⟺pˉ เป็นโพล
แผนภาพโพล-ซีโร่ของระบบจริงจึงสมมาตรกับแกนจริงเสมอ
ตัวอย่าง R-1 · วาดแผนภาพโพล-ซีโร่H(s)=s3+5s2+12s+82s+6
เศษ: 2s+6=2(s+3) ซีโร่ s=−3 และ K=2
ส่วน: ลองรากตรรกยะ s=−1: −1+5−12+8=0 ✓ หารสังเคราะห์ได้ s2+4s+8
s2+4s+8=0⇒s=2−4±16−32=−2±2j
- H(s)=(s+1)(s+2−2j)(s+2+2j)2(s+3)
ซีโร่ ○ ที่ −3 · โพล × ที่ −1 และ −2±2j ทุกโพลอยู่ครึ่งซ้าย ระบบเหตุผลนี้เสถียร
ROC คืออะไร: หาจากนิยามทีละขั้น พื้นฐาน
อินทิกรัล ∫x(t)e−stdt ลู่เข้าหรือไม่ขึ้นกับ ∣x(t)e−st∣=∣x(t)∣e−σt เท่านั้น ส่วน ω ไม่มีผลเลย ROC จึงขึ้นกับ σ=Res อย่างเดียว และเป็นแถบแนวตั้งเสมอ
R4 Region of convergence บริเวณลู่เข้า
ROC={s:∫−∞∞∣x(t)∣e−σtdt<∞},σ=Res
R5 Right-sided exponential เอกซ์โพเนนเชียลด้านขวา
e−atu(t) ⟷ s+a1,Res>−a
R6 Left-sided exponential เอกซ์โพเนนเชียลด้านซ้าย
−e−atu(−t) ⟷ s+a1,Res<−a
ตัวอย่าง R-2 · สัญญาณด้านซ้ายหาการแปลงสองด้านและ ROC ของ x(t)=−e−atu(−t)
สัญญาณมีค่าเฉพาะ t<0
X(s)=∫−∞0−e−ate−stdt=−∫−∞0e−(s+a)tdt
เมื่อ t→−∞ พจน์ e−(s+a)t=e−(σ+a)te−jωt ต้องหายไป ซึ่งเกิดเมื่อ σ+a<0
X(s)=−[s+a−e−(s+a)t]−∞0=−(s+a−1−0)=s+a1
X(s)=s+a1 เหมือนของ e−atu(t) ทุกตัวอักษร แต่ ROC เป็น Res<−a (ครึ่งซ้าย) แทน Res>−a
โพลเดียวกันที่ s=−a ROC ทางขวาให้สัญญาณด้านขวา e−atu(t) ROC ทางซ้ายให้สัญญาณด้านซ้าย −e−atu(−t) — X(s) ตัวเดียวไม่พอ ต้องบอก ROC เสมอ
สมบัติของ ROC ทั้ง 8 ข้อ กลาง
แปดข้อนี้ทำให้หา ROC ได้โดยไม่ต้องอินทิเกรตเลย แค่รู้ตำแหน่งโพลและรู้ว่าสัญญาณอยู่ด้านไหนของเวลา (เรียงตาม Oppenheim & Willsky)
R7 Property 1 เป็นแถบขนานกับแกน jω
ROC={s:σ1<Res<σ2}
ขอบอาจเป็น −∞ หรือ +∞ ก็ได้ เพราะการลู่เข้าขึ้นกับ σ เท่านั้น
R8 Property 2 ไม่มีโพลอยู่ใน ROC
p เป็นโพล⇒p∈/ROC
ที่โพล X(s) เป็นอนันต์ อินทิกรัลจึงลู่เข้าไม่ได้
R9 Property 3 สัญญาณช่วงจำกัด → ทั้งระนาบ
x(t)=0 นอก [T1,T2] และอินทิเกรตได้สัมบูรณ์⇒ROC=C
R10 Property 4 ด้านขวา → ครึ่งระนาบขวา
x(t)=0, t<T1: σ0∈ROC⇒{Res>σ0}⊂ROC
R11 Property 5 ด้านซ้าย → ครึ่งระนาบซ้าย
x(t)=0, t>T2: σ0∈ROC⇒{Res<σ0}⊂ROC
R12 Property 6 สองด้าน → แถบ (หรือว่างเปล่า)
ROC=ROCขวา∩ROCซ้าย
ถ้าส่วนขวาต้องการ σ>σa แต่ส่วนซ้ายต้องการ σ<σb ที่ σb≤σa ROC ว่าง: ไม่มีการแปลงลาปลาซ
R13 Property 7 ตรรกยะ + ด้านขวา
ROC={Res>Repขวาสุด}
การแปลงด้านเดียวทั้งหมดอยู่ในกรณีนี้
R14 Property 8 ตรรกยะ + ด้านซ้าย
ROC={Res<Repซ้ายสุด}
X(s)=(s+1)(s+2)1 มี ROC ได้สามแบบ: ขวาของ −1 (ด้านขวา) แถบระหว่าง −2 กับ −1 (สองด้าน) และซ้ายของ −2 (ด้านซ้าย) ไม่มีแบบอื่นเพราะ ROC ห้ามมีโพลข้างใน
ตัวอย่าง R-3 · X(s) เดียว ROC สามแบบ สัญญาณสามตัวหา x(t) ของ X(s)=(s+1)(s+2)1 สำหรับทั้งสาม ROC
เศษส่วนย่อย (ไม่ขึ้นกับ ROC): X(s)=s+11−s+21
กฎ: โพลที่อยู่ซ้ายของ ROC ให้พจน์ด้านขวา eptu(t) R5 · โพลที่อยู่ขวาของ ROC ให้พจน์ด้านซ้าย −eptu(−t) R6
(ก) Res>−1: ทั้งสองโพลอยู่ซ้าย
x(t)=(e−t−e−2t)u(t)
(ข) −2<Res<−1: โพล −1 อยู่ขวา โพล −2 อยู่ซ้าย
x(t)=−e−tu(−t)−e−2tu(t)
(ค) Res<−2: ทั้งสองโพลอยู่ขวา
x(t)=(−e−t+e−2t)u(−t)
สามสัญญาณต่างกันโดยสิ้นเชิง มีเพียง (ก) ที่เป็นเหตุผล และมีเพียง (ก) ที่เสถียร (ROC ครอบ σ=0) · ลองเลือก ROC ทั้งสามแบบในแล็บด้านล่างด้วยปุ่ม “สัญญาณสองด้าน”
สัญญาณสองด้านและสัญญาณช่วงจำกัด กลาง
ตัวอย่าง R-4 · ผลรวมด้านขวา + ด้านซ้ายx(t)=e−2tu(t)+e3tu(−t)
ส่วนขวา: e−2tu(t)→s+21, Res>−2
ส่วนซ้าย: e3tu(−t)=−(−e3tu(−t))→−s−31, Res<3 R6 ด้วย a=−3
รวมและหาส่วนร่วมของ ROC R12
X(s)=s+21−s−31=(s+2)(s−3)−5,−2<Res<3
ROC เป็นแถบ −2<σ<3 ซึ่งครอบแกน jω จึงมีการแปลงฟูเรียร์ X(jω)=(jω+2)(jω−3)−5 แม้จะมีโพลครึ่งขวา (เพราะส่วนที่โตขึ้นอยู่ฝั่ง t<0 ซึ่งลดลงเมื่อย้อนเวลา)
ตัวอย่าง R-5 · เอกซ์โพเนนเชียลสองด้านสมมาตรx(t)=e−b∣t∣, b>0
- x(t)=e−btu(t)+ebtu(−t) → s+b1−s−b1
- X(s)=s2−b2(s−b)−(s+b)=s2−b2−2b,−b<Res<b
โพลที่ ±b ROC เป็นแถบระหว่างโพลทั้งสอง ถ้า b≤0 ส่วนขวาต้องการ σ>b ส่วนซ้ายต้องการ σ<−b ไม่มีส่วนร่วม จึงไม่มีการแปลงลาปลาซ
ตัวอย่าง R-6 · ช่วงจำกัด: ภาวะเอกฐานที่ถอดได้พัลส์สี่เหลี่ยม x(t)=u(t)−u(t−T)
- X(s)=∫0Te−stdt=s1−e−sT
ดูเหมือนมีโพลที่ s=0 แต่ lims→0s1−e−sT=T (โลปิตาล) จึงเป็นภาวะเอกฐานที่ถอดได้ ไม่ใช่โพล Z19
ROC คือทั้งระนาบ s ตามสมบัติ 3 R9 ซีโร่อยู่ที่ s=j2πk/T, k=0 บนแกน jω ซึ่งคือศูนย์ของฟังก์ชัน sinc ในสเปกตรัมนั่นเอง
ตารางคู่การแปลงสองด้านพร้อม ROC กลาง
| สัญญาณ x(t) | X(s) | ROC |
| δ(t) | 1 | ทุก s |
| δ(t−T) | e−sT | ทุก s |
| u(t) | s1 | Res>0 |
| −u(−t) | s1 | Res<0 |
| (n−1)!tn−1u(t) | sn1 | Res>0 |
| −(n−1)!tn−1u(−t) | sn1 | Res<0 |
| e−atu(t) | s+a1 | Res>−a |
| −e−atu(−t) | s+a1 | Res<−a |
| (n−1)!tn−1e−atu(t) | (s+a)n1 | Res>−a |
| −(n−1)!tn−1e−atu(−t) | (s+a)n1 | Res<−a |
| cos(ω0t)u(t) | s2+ω02s | Res>0 |
| sin(ω0t)u(t) | s2+ω02ω0 | Res>0 |
| e−atcos(ω0t)u(t) | (s+a)2+ω02s+a | Res>−a |
| e−atsin(ω0t)u(t) | (s+a)2+ω02ω0 | Res>−a |
| e−a∣t∣, a>0 | s2−a2−2a | −a<Res<a |
| u(t)−u(t−T) | s1−e−sT | ทุก s |
ทุกคู่ที่มี −(…)u(−t) คือฝาแฝดด้านซ้ายของคู่ด้านขวาบรรทัดบน: X(s) เหมือนกัน ROC อยู่คนละฝั่งของโพล
ตำแหน่งโพลบอกรูปคลื่น พื้นฐาน
แต่ละโพล p ของระบบเหตุผลให้พจน์ ept ในผลตอบสนองอิมพัลส์ (คู่เชิงซ้อนรวมกันเป็น eσtcos(ωt+ϕ)) การอ่านแผนภาพโพลจึงคือการอ่านรูปคลื่นโดยไม่ต้องแปลงกลับ
ตำแหน่งโพล (ซ้าย) กับผลตอบสนองอิมพัลส์ (ขวา) ระยะจากแกน jω (|σ|) บอกความเร็วในการลดลง ระยะจากแกนจริง (ω) บอกความถี่ของการแกว่ง โพลที่อยู่ใกล้แกน jω ที่สุดเรียกว่าโพลเด่น (dominant pole) เพราะลดลงช้าที่สุดจึงกำหนดรูปร่างหลัก
R15 Reading a pole pair อ่านคู่โพล −σ±jωd
τ=σ1,ωn=∣p∣=σ2+ωd2,ζ=cosθ=ωnσ
θ คือมุมระหว่างเส้นจากจุดกำเนิดไปยังโพลกับแกนจริงลบ เส้นมุมคงที่ = ζ คงที่ (%OS คงที่) วงกลมรัศมีคงที่ = ωn คงที่
R16 What zeros do ซีโร่ทำอะไร
เรซิดิวที่โพล pk=K∏j=k(pk−pj)∏i(pk−zi)
ซีโร่ไม่สร้างพจน์ใหม่ แต่เปลี่ยน “น้ำหนัก” ของแต่ละโพล ซีโร่ใกล้โพลไหนจะลดพจน์นั้นลง ซีโร่ครึ่งขวาทำให้ผลตอบสนองขั้นบันไดวิ่งสวนทางก่อน (undershoot)
ตัวอย่าง R-7 · อ่านแผนภาพโพลระบบมีโพลที่ −2±4j และ −10 ประมาณรูปผลตอบสนองโดยไม่แปลงกลับ
โพลเด่นคือ −2±4j (ใกล้แกน jω กว่า) โพล −10 ลดลงเร็วกว่า 5 เท่า (τ=0.1 s เทียบกับ 0.5 s) แทบไม่มีผลหลังผ่านไปไม่กี่ร้อยมิลลิวินาที
จาก R15: ωn=4+16=20≈4.47 rad/s, ζ=2/4.47≈0.447
แกว่งที่ ωd=4 rad/s (คาบ ≈1.57 s) พุ่งเกิน %OS=100e−ζπ/1−ζ2≈100e−π/2≈20.8% เข้าที่ใน ts≈4/σ=2 s
ประมาณได้เป็นระบบอันดับสองที่ ζ≈0.45, ωn≈4.5 rad/s พุ่งเกินราว 21% เข้าที่ราว 2 วินาที ดูวิธีเลือกตำแหน่งโพลให้ได้สเปกที่ต้องการในผลตอบสนองทางเวลา
ROC กับเสถียรภาพ ความเป็นเหตุเป็นผล และการแปลงฟูเรียร์ ขั้นสูง
R17 Causality ความเป็นเหตุเป็นผล
h(t)=0, t<0⟺ROC ของ H(s) เป็นครึ่งระนาบขวา
สำหรับ H(s) ตรรกยะ: ROC อยู่ขวาของโพลขวาสุด
R18 BIBO stability เสถียรภาพ
∫−∞∞∣h(t)∣dt<∞⟺jω-axis⊂ROC
R19 Causal and stable ทั้งเหตุผลและเสถียร
Repk<0 ทุกโพล
กฎ “โพลครึ่งซ้าย = เสถียร” ในวิชาควบคุมเป็นจริงเพราะระบบจริงเป็นเหตุผลเสมอ
R20 Fourier transform การแปลงฟูเรียร์จากลาปลาซ
X(jω)=X(s)s=jωเมื่อ ROC ครอบแกน jω
ตัวอย่าง R-8 · ระบบเดียว สามสมมติฐานH(s)=(s+1)(s−2)s−1 จงหา ROC ถ้า (ก) ระบบเหตุผล (ข) ระบบเสถียร (ค) ทั้งสองอย่างได้ไหม
โพล −1 และ 2 ROC ที่เป็นไปได้: σ<−1, −1<σ<2, σ>2
(ก) เหตุผล → ขวาของโพลขวาสุด: Res>2 (ไม่ครอบแกน jω → ไม่เสถียร)
(ข) เสถียร → ต้องครอบ σ=0: −1<Res<2 (สองด้าน → ไม่เหตุผล)
(ค) ต้องเป็นครึ่งระนาบขวาและครอบแกน jω พร้อมกัน แต่โพล s=2 ขวางอยู่
เป็นไปไม่ได้ที่จะทั้งเหตุผลและเสถียร เพราะมีโพลครึ่งขวา ส่วนซีโร่ที่ s=1 (ครึ่งขวา) ไม่มีผลต่อเสถียรภาพ แต่ทำให้เป็นระบบ non-minimum phase
ROC เปลี่ยนอย่างไรเมื่อใช้คุณสมบัติ ขั้นสูง
| การดำเนินการ | X(s) ใหม่ | ROC ใหม่ (R = ROC เดิม) |
| ax1+bx2 | aX1+bX2 | อย่างน้อย R1∩R2 (กว้างขึ้นได้ถ้าโพลตัดกับซีโร่) |
| x(t−T) | e−sTX(s) | R เท่าเดิม |
| es0tx(t) | X(s−s0) | R+Res0 (เลื่อนทั้งแถบ) |
| x(at) | ∣a∣1X(s/a) | aR (ถ้า a<0 แถบกลับข้าง) |
| x(−t) | X(−s) | −R (สะท้อนข้ามแกน jω) |
| dx/dt | sX(s) | อย่างน้อย R |
| ∫−∞txdτ | X(s)/s | อย่างน้อย R∩{Res>0} |
| x1∗x2 | X1X2 | อย่างน้อย R1∩R2 |
รายละเอียดและพิสูจน์ของทุกคุณสมบัติอยู่ในหน้า 8.3
แล็บโพล-ซีโร่-ROC: ลากแล้วดูสัญญาณเปลี่ยน
แล็บนี้ใช้การแปลงสองด้าน เลือก ROC ได้ทุกแถบที่เป็นไปได้ แล้วดูว่าสัญญาณ x(t) เปลี่ยนเป็นด้านขวา ด้านซ้าย หรือสองด้าน พร้อมป้ายบอกความเป็นเหตุเป็นผลและเสถียรภาพ ถ้า ROC ครอบแกน jω จะแสดง ∣X(jω)∣ ให้ด้วย
ลองสิ่งเหล่านี้: (1) ลากโพลจริงข้ามแกน jω แล้วดูป้ายเสถียรภาพ (2) กด “สัญญาณสองด้าน” แล้วเลือก ROC แต่ละแถบ (3) เพิ่มคู่โพลเชิงซ้อนแล้วลากเข้าใกล้แกน jω ดูยอดของ |X(jω)| สูงขึ้นและการแกว่งช้าลงที่จะหาย (4) วางซีโร่ทับบนโพลเกือบพอดี ดูพจน์นั้นหายไปจากสัญญาณ
ดูภูมิประเทศ |X(s)| พร้อมกัน
คำถามที่พบบ่อย
ROC ของการแปลงลาปลาซหาอย่างไรเร็วที่สุด?
ถ้า X(s) เป็นฟังก์ชันตรรกยะ ให้หาโพลก่อน ROC ไม่มีโพลอยู่ข้างในและขอบเขตเป็นเส้นแนวตั้งที่ผ่านโพล สัญญาณด้านขวา (เช่นคูณ u(t)) ROC อยู่ขวาของโพลขวาสุด สัญญาณด้านซ้าย ROC อยู่ซ้ายของโพลซ้ายสุด สัญญาณสองด้าน ROC เป็นแถบระหว่างโพลสองตัว
ทำไมต้องบอก ROC ด้วย ในเมื่อ X(s) ก็รู้แล้ว?
เพราะในการแปลงสองด้าน X(s) เดียวกันมาจากสัญญาณได้หลายตัว เช่น e−atu(t) และ −e−atu(−t) มีการแปลงเท่ากันคือ 1/(s+a) ต่างกันเพียง ROC ซึ่งเป็น Res>−a กับ Res<−a ถ้าไม่บอก ROC การแปลงกลับจะกำกวม
ROC บอกเสถียรภาพได้อย่างไร?
ระบบเสถียรแบบ BIBO ก็ต่อเมื่อ ROC ของ H(s) ครอบแกน jω ถ้าระบบเป็นเหตุเป็นผลด้วย ROC เป็นครึ่งระนาบทางขวาของโพลขวาสุด เงื่อนไขทั้งสองรวมกันจึงเป็นกฎที่คุ้นเคย: ระบบเหตุผลเสถียรเมื่อโพลทุกตัวอยู่ครึ่งซ้ายของระนาบ s
ในวิชาวงจรใช้การแปลงด้านเดียว ต้องสนใจ ROC ไหม?
ไม่ค่อยต้อง เพราะการแปลงด้านเดียวมี ROC เป็นครึ่งระนาบขวาของโพลขวาสุดเสมอ จึงไม่กำกวม แต่ยังต้องใช้แนวคิดนี้ในการตรวจว่าใช้ทฤษฎีบทค่าสุดท้ายได้ไหม และในวิชาสัญญาณและระบบซึ่งใช้การแปลงสองด้าน
มด
อาจารย์มด · ธีรวุฒิ แสวงบุญ
อาจารย์สาขาวิชาวิศวกรรมไฟฟ้า มหาวิทยาลัยราชภัฏธนบุรี สมุทรปราการ · ดูโปรไฟล์