292 บทท่ี 7 สมการอนุพนั ธแ บบเชิงเสน คณติ ศาสตรว ิศวกรรมอิเลก็ ทรอนกิ ส เอกสารอา งองิ วรี ศกั ดิ์ บุญทน. (2553). คณิตศาสตรอเิ ลก็ ทรอนิกส 2. กรงุ เทพมหานคร: สาํ นักพมิ พแหง จุฬาลงกรณม หาวิทยาลยั . ภคั คนิ ี ชติ สกลุ และคณะ. (2557). คณติ ศาสตรว ิศวกรรม Advanced Engineering Mathematics. กรุงเทพมหานคร. Erwin Kreyszig . พรชัย สารทวาหา . (2550). สมการเชิงอนุพนั ธ. กรุงเทพมหานคร : ภาควชิ าคณติ ศาสตร คณะวทิ ยาศาสตร จุฬาลงกรณมหาวิทยาลยั . นิรดั ร คาํ ประเสริ ฐ . (2553). คณิตศาสตรว ศิ วกรรมไฟฟา 4. กรุงเทพมหานคร: ศนู ยส่ือเสรมิ กรุงเทพฯ . สุวัฒน รอดผล. (2546). สมการเชงิ อนพุ นั ธสาํ หรบั วศิ วกร. กรงุ เทพมหานคร : สาํ นกั พิมพ ส.ส.ท สมาคมสง เสรมิ เทคโนโลยไี ทยญีป่ นุ . สําเรงิ ชน่ื รงั สิกลุ . (2555). สมการเชิงอนุพันธ. กรุงเทพมหานคร : สํานักพิมพแ หงจุฬาลงกรณ มหาวทิ ยาลยั . Linear Differential Equations. [ออนไลน] เขาถึงไดจาก http://www.tutorial.math.lamar.edu Classes /DE/Linear.aspx. (วันท่ีคน ขอ มูล 5 เมษายน 2556) Order and Linearity of Differential Equations. [ออนไลน] เขาถึงไดจาก http://www.analyzemath.com/calculus /Differential_Equations/order_linearity.html (วันทค่ี น ขอมลู 10 เมษายน 2556) Chapter 14 Difference Equations 1. [ออนไลน] เขาถึงไดจ าก http://www.cimt.plymouth.ac.uk /Projects/.../discrete_ch14.pdf (วนั ที่คน ขอมูล 10 เมษายน 2556) Chapter11 Differential Equation. [ออนไลน] เขา ถงึ ไดจาก http://ugrad.math.ubc.ca/ Coursedoc/ math103/keshet.notes/chapter11notes.pdf. (วนั ทค่ี น ขอมลู 1 พฤษภาคม 2556) Chapter4 Application of second Order Differential Equations in Mechanical Engineering Analysis. [ออนไลน] เขา ถึงไดจาก http://www.engr.sjsu.edu/trhsu/Chapter%204% 20second %20order%20des.pdf. (วนั ทีค่ น ขอ มูล 10 เมษายน 2556) Application second order. [ออนไลน] เขา ถงึ ไดจ าก http://www.stewartcalculus.com/data/ CALCULUS%20concepts%20and%20Contexts/upfiles/3C3-Appsof2ndOrders_stu.pdf . (วันทค่ี น ขอ มูล 10 เมษายน 2556) จุไรรตั นจ ินดา อรรคนติ ย
บทที่ 8 การแปลงลาปลาซ การแปลงลาปลาซหรือผลการแปลงลาปลาซหรือเรียกวาแคลคูลัสการดําเนินการ เปนสมการท่ีสําคัญสําหรับนักคณิตศาสตร วิศวกร นักฟสิกส ซึ่งหลักการของลาปลาซ ทรานสฟอรม เปนการแกป ญหาจากโดเมนหนง่ึ ไปเปนอีกโดนเมนหน่งึ สําหรับทางวิศวกรรมไฟฟา น้ันจะเปนการแปลงสมการจากโดเมนเวลาเปนโดเมนความถ่ีเชิงซอนหรือ โดเมน s จากนั้น จึงทําการแกสมการโดยใชหลักการทางพีชคณิต แลวจึงกลับไปแปลงจากโดเมน s กลับไปเปน โดเมนของเวลาอกี คร้ังโดยใชการแปลงลาปลาซผกผัน 1 รูปท่ี 8.1 ความสมั พนั ธข องลาปลาซทรานสฟอรมและการแปลงกลบั คา ของกระแสหรือแรงดันในวงจรไฟฟาชวงภาวะทรานเชียนตตองใชวิธีการของสมการ อนพุ ันธม าชว ยแกปญหา ถา วงจรไฟฟาประกอบดวยตัวตานทานและตวั เหน่ียวนําหรือตัวตานทาน และตวั เก็บประจุ สมการท่ีใชในการแกสมการจะอยูในรปู ของสมการอนพุ ันธอันดับหนึ่ง (First- order differential equation) ซ่ึงสามารถแกปญหาไดหลายวิธี แตถาวงจรไฟฟาประกอบดวย ตัวตานทาน ตัวเก็บประจุและตัวเหนี่ยวนํา สมการความสัมพันธที่ไดจะมีลักษณะเปนสมการ อนุพันธอันดับสอง (Second-order differential equation) ซ่ึงก็มีวิธีการแกปญหาไดอีกหลายวิธี แตกตางกันออกไป การแกสมการโดยใชวิธีการลาปลาซทรานสฟอรมถือเปนเครื่องมือ ที่มีประสิทธิภาพอยางยิ่งในการแกปญหาวงจรไฟฟาในสภาวะทรานเชียนต ซ่ึงเปนการแปลง จากฟงกชันของเวลา (Time domain) ไปเปนฟงกชันของความถ่ีเชิงซอน s (Complex frequency domain) ดงั รูปท่ี 8.1
294 บทที่ 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส รูปท่ี 8.2 การแกป ญ หาวงจรไฟฟา โดยวิธีสมการอนพุ ันธแ ละวธิ ีลาปลาซทรานสฟอรม จากรูปที่ 8.2 วิธีการของลาปลาซทรานสฟอรม จะเปลี่ยนสมการอนุพันธหรอื วงจรไฟฟา นั้นไปอยูในฟงกชันของความถี่แลวใชหลักการพีชคณิตธรรมดาจัดใหตัวแปรในรูปแบบ ทเ่ี หมาะสม แลวใชต ารางลาปลาซทรานสฟอรมเทยี บหาคา ออกมาในรปู ของฟง กชนั ของเวลาจะได คาํ ตอบท่ีตอ งการ การเทียบคากลบั มาฟงกชันของความถ่ไี ปเปน ฟงกชันของเวลานัน้ เรียกวาเปน การ ทําลาปลาซทรานสฟอรมกลับ (Inverse laplace transform) หลักของการลาปลาซทรานสฟอรม จะคลา ยกับหลักการของลอการิทึม เชน ถาตอ งการหาคาของฟง กชนั ใดใหใสลอการิทึมเขา ไปจะได คา ออกมาตามทเ่ี ปดจากตารางลอการิทึม และเมื่อตองการคําตอบใหแอนตลี้ อการิทึมโดยใชตาราง เชน เดียวกนั 8.1 ปริพันธไ มต รงแบบ (Improper integral) b บทนิยาม 8.1 ปริพันธไมตรงแบบ (Improper integral) คือ ปริพันธจํากัดเขต a f (x)dx เม่ือ ปริพัทธของ f (x) มคี า ทไ่ี มต อ เน่ืองบนชวง [a, b] ตั้งแตห นึง่ จดุ ขึ้นไป หรือ การหาลมิ ิตของปรพิ นั ธ อยางนอยคา หน่ึงเปน อนนั ต โดยแบง การหาปริพนั ธไมตรงแบบออกเปน 2 ชนดิ คอื 8.1.1 ปริพันธไมตรงแบบบนชวงปด ปริพัทธของ f (x) มีจุดท่ีไมตอเน่ืองบนชวง [a, b] ตั้งแต หนึง่ จุดข้ึนไป แบงการ พิจารณาออกเปน 3 กรณี กรณี 1 ถา f ตอเนือ่ งบนชว ง [ a, b) แตไ มตอ เนอื่ งท่ี x = b นิยามโดย จไุ รรตั นจ ินดา อรรคนติ ย
คณิตศาสตรว ิศวกรรมอิเลก็ ทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 295 b f (x)dx =l limb- l f (x)dx (8.1) a a กรณี 2 ถา f ตอเน่ืองบนชวง (a ,b] แตไ มตอเนอ่ื งท่ี x = a นยิ ามโดย b b f (x)dx =l lima+ f (x)dx (8.2) a l กรณี 3 ถา f ตอเน่ืองทท่ี กุ ๆ จุดบนชว ง [a, b] ยกเวน ท่ี x = c เมอื่ c(a, b) นิยามโดย b k b f (x)dx =k limc- f ( x)dx +l limc + f (x)dx (8.3) a a l จากสมการปริพัทธไมตรงแบบ (8.1), (8.2) และสูตร (8.3) ถาลิมิตหาคาได จะเรียกวา ปริพันธไมตรงแบบชนิดลูเขา (Convergent) และถาลิมิตหาคาไมได จะรียกวาปริพันธไมตรงแบบ ชนดิ ลอู อก (Divergent) 8.1.2 ปริพันธไมตรงแบบบนชวงเปด อนนั ต ปรพิ ทั ธของ f (x) ท่ีมีลมิ ติ ของการหาปริพันธอยาง นอยหนึง่ คา เปน อนนั ต แบง การ พจิ ารณาออกเปน 3 กรณี กรณี 1 ถา f ตอ เน่อื งบนชว ง [a, ) นยิ ามโดย l f ( x )dx =l lim f ( x )dx (8.4) a a กรณี 2 ถา f ตอเนอื่ งบนชว ง (- , b] นยิ ามโดย b b f ( x )dx =l lim f ( x )dx (8.5) - l กรณี 3 ถา f ตอเนื่องทที่ ุกๆ จดุ บนชวง (- ,) และ a (- , ) นยิ ามโดย b a f ( x )dx = f ( x )dx +a f ( x)dx (8.6) a - a l =k lim f ( x )dx +l lim f ( x)dx k a จากการปริพัทธไมตรงแบบ (8.4), (8.5) และ (8.6) จะไดวาถาลิมิตหาคาได จะเรียกวาปริพันธ ไมตรงแบบชนดิ ลเู ขา (Convergent) และถาลิมติ หาคาไมไ ดจ ะรียกวาปรพิ นั ธไ มตรงแบบชนิดลอู อก (Divergent) 8.2 คาํ จํากดั ความของลาปลาซทรานสฟอรม ถา f (t) = 0 ท่ี t < 0 ซ่งึ สามารถหาคา ไดท ี่ t > 0 แลวนัน้ จะสามารถหาคา ลาปลาซทรานสฟ อรม ของ f (t) นน้ั ไดจ าก จไุ รรัตนจ ินดา อรรคนติ ย
296 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส L[ f (t)] = 0 f (t)e-st dt (8.7) เม่ือ L เปนเครื่องหมายของลาปลาซทรานสฟอรม และ s เปนฟงกชันของความถ่ีเชิงซอน อยูในรูป s = + j และ L[f (t)] เรียกวาลาปลาซทรานสฟอรมของ f (t) และเปนลาปลาซทรานส ฟอรม ดานเดียว (Unilateral laplace transform) คือจาก 0 จากสมการท่ี (8.7) เรียกวา อนิ ทิกรัลไมตรงแบบ (Improper integral) หรือปริพันธไมตรงแบบ ซึ่งกําหนดโดย เมอื่ ทําการหาปรพิ นั ธแ ลว แทนคา ลิมติ เขา ไปจะอยใู นรูปของตวั แปร s หรอื สามารถ เขียนไดเปน L[ f (t )] = f (t )e- st dt = lim b e-st f (t )dt (8.8) b 0 0 จะพบวาผลการแปลงลาปลาซของฟงกชัน f (t) จะมีผลการแปลงหรือไมนั้นข้ึนอยูกับลิมิต ทางขวามือของสมการวาลิมิตหาคาไดหรือไม ถาลิมิตหาคาไดจะเรียกวา อินทิกรัลลูเขา และ สัญลักษณ L ซ่งึ แปลง f (t) เปน F(s) จะเรยี กวา ตัวดาํ เนนิ การแปลงลาปลาซ L[ f (t )] = F (s) (8.9) เปนการเปล่ียนแปลงฟงกชันของเวลาไปเปนฟงกชันของความถ่ีเชิงซอน จากสมการที่ (8.8) จะหาคา ไดแนน อนหรอื อยใู นชวงลูเขา (Converse) เมอื่ 0 f (t ) e- tdt < ทงั้ นเี้ นื่องจาก = f (t) e-(t+ j)t 0 f (t)e-st = f (t)e-t e jt = f (t)et 1 = f (t)et คาที่จริงบวก (Positive real) ใดๆ ของ ถา f (t) < Meat เสมอไปท่ีคา t เปนบวก และ การหาปริพันธจะไดคาท่ีแนนอนเม่ือ > น่ันคือชวงที่จะทําใหสมการลูเขาอยูท่ี < < โดยพิสูจนไดจาก จุไรรัตนจ ินดา อรรคนติ ย
คณิตศาสตรว ิศวกรรมอิเล็กทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 297 0 f (t ) e-t dt < 0 Me-t e-t dt = 0 M e(-)t d (-)t - = M e -( - ) t - 0 = M ถา > - ดังน้ันถาเลือกคา > แลวผลการหาปริพันธจะไดเปน M/( - ) และจะไดคาท่ีแนนอน เมื่อ < ดวย ดงั นนั้ คา ทเี่ หมาะสมของ คอื < < จากสมการที่ (8.9) เมอ่ื ตองการเปล่ียนกลับจาก F(s) เปน f (t) ทําไดโดยการทาํ ลาปลาซทรานส ฟอรม กลบั (Inverse laplace) เขยี นไดเ ปน f (t) = L-1[F(s)] f (t) เปนการทําลาปลาซทรานสฟอรมกลับของ F(s) และจะทําไดเมื่อ t > 0 โดยหาคาไดจาก สมการ f (t) = j12 -+jj F(s)est ds (8.10) ในการทําลาปลาซทรานสฟ อรมกลบั ที่นาํ มาใชจริงๆ จะหลีกเล่ียงการหาปริพันธตามสมการที่ (8.10) แตจะพยายามใชคาในตารางลาปลาซทรานสฟอรมเปรียบเทียบกลับใหไดคาท่ีตองการ เพียงแตใ นตารางจะตอ งมฟี งกชันที่จําเปนใหเพยี งพอเทา นนั้ ตัวอยางที่ 8.1 จงหาผลการแปลงลาปลาซของฟงกช ัน f (t) = 1 เม่อื t > 0 วธิ ีทาํ ข้นั ตอนท่ี 1 หาผลการแปลงลาปลาซจากนิยามจะได L[f (t)] = 0 f (t)e-st dt = b lim 0b e-st f (t)dt ขน้ั ตอนท่ี 2 โจทยก าํ หนดให f (t) = 1 ทาํ การแทนคาในนยิ ามจะได L[f (t)] = 0 (1)e-stdt = b lim 0b (1)e-st dt = b lim 0b e- st d (-st ) -s = b lim -es-st 0 ตอบ L[1] = 1s เมอื่ s > 0 = b lim -e-st +1 = 1s เมือ่ s > 0 จุไรรตั นจ นิ ดา อรรคนติ ย
298 บทที่ 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอเิ ล็กทรอนกิ ส 8.2 การคํานวณหาลาปลาซทรานสฟอรม ทจ่ี ริงแลว f (t) ใดๆ ท่ีหาคาไดท่ี t > 0 จะสามารถหาคา ลาปลาซทรานสฟอรม ของฟงกชันนัน้ ๆ ไดท ั้งหมด ตวั อยา งการหาคาลาปลาซทรานสฟอรมของฟง กช นั ทสี่ าํ คัญมดี งั ตอ ไปนี้ 8.2.1 ฟงกชันยูนิตสเต็ป (Unit step function) เปนฟงกชันท่ีมีลักษณะเปนเสนตรงคงท่ี ในแนวแกน y สวนคา x จะมีคาเพิม่ ขนึ้ เรอ่ื ยๆ โดยคา ในแนวแกน x จะเร่ิมท่ี 0 สว นคาคงที่ y จะมีคา เทากับ 1 โดยกําหนดให u(t) แทนฟงกชันยูนิตสเต็ปลักษณะเหมือนแหลงจายไฟกระแสตรง ในวงจรไฟฟา f (t) = u(t) ดงั รปู ที่ 8.3 u( t) f( x)=u(x) 2 1 t 2 34 -2 -1 1 -1 -2 รูปที่ 8.3 ฟง กช นั ยนู ิตสเต็ป f (t) = u(t) = 1 t 0 0 t 0 < เมื่อ f (t) = u(t) = 1 สามารถหาลาปลาซทรานฟอรม ของ u(t) ตามสมการท่ี (8.8) ไดดังนี้ L[u(t )] = 0 f (t )e - st dt = u(t )e -st dt 0 เม่อื u(t) = 1 จะได L[u(t)] = 01e-st dt แทนคา t = 0 และ t = จะได L[u(t )] = - 1s e-st = - 1s (e-s - e-s0 ) 0 e-= 0 และ e-s0 = 1 ดังนั้น L[u(t)] = - 1s (0 - 1) = - (- 1s ) = 1s จไุ รรตั นจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 299 8.2.2 ฟง กชันเอกซโพเนนเชยี ล (Exponential function) เมื่อ f (t) อยูใ นรูปของเอกซโพเนนเชยี ล (Exponential) เชน f (t) = e-at ซ่งึ ในวงจรไฟฟา ไดแกการคายประจขุ องตัวเกบ็ ประจโุ ดยสามารถหา ลาปลาซทรานสฟ อรม ของ f (t) เมอื่ a เปนคา คงทไี่ ดจ ากสมการท่ี (8.8) f(t) f( x)=2 e^- 2x t 1 1.5 3 2 1 -0.5 0.5 รปู ที่ 8.4 ฟงกชันเอกซโ พเนนเชียลของตัวเก็บประจุเมอื่ มกี ารคายประจุ ลาปลาซทรานสฟอรมของ f (t) = e-at คือ L[e-at] = 0 e-at e-st dt = 0 e-(s+a)t dt = - s+1a e -( s + a ) t 0 แทนคา t = 0 และ t = จะได = - s+1a [e-(s+a) - e-(s+a)0 ] 1 1 = - s+1a [e- - e-0 ] s+1a e e0 = - [ - ] = - s+1a [0 - 11] = - s+1a (-1) L[e-at] = s+1a 8.2.3 ฟงกชันไซน (Sine function) เม่อื f (t) อยูในรูปของตรโี กณ เชน f (t) = sin t ซ่งึ พบเสมอ ในฟง กช ันของการจา ยกาํ ลงั ไฟฟาสลบั จงหาลาปลาซทรานสฟอรมของ f (t) จุไรรตั นจ ินดา อรรคนติ ย
300 บทที่ 8 การแปลงลาปลาซ คณิตศาสตรว ิศวกรรมอิเล็กทรอนกิ ส y = sin t 3 2 1 t 246 -6 -4 -2 -1 -2 -3 รปู ท่ี 8.5 ฟง กช ันไซน ลาปลาซทรานสฟอรม ของ f (t) จะได L[sin t] = 0 sin te-st dt จาก e jt = cost + j sin t และ e-jt = cost - j sin t นํามาลบกันได e jt - e-jt = 2j sin t ดงั น้ัน sin t = 21j (e jt - e- jt ) และ L[ e jt ] = s-1j ดังนัน้ L[sin t] = 21j s-1j - s+1j = s2 +2 โดยวิธีเดยี วกนั ถา f (t) = cos t s L[cos t] = s2 +2 8.2.4 ฟงกชันของอนุพันธ ถา f (t) อยใู นรูปของสมการอนุพันธ เชน d f (t ) ซ่ึงจะพบในการ dt d หาแรงดันตกครอม L ขณะท่ีมกี ระแส i(t) ไหลผานจะได L [ dt i(t) ] จะหาลาปลาซทรานสฟอรม ของฟงกช นั แบบนไ้ี ดจาก L ddt f (t ) = 0 ddt f (t)e-stdt = 0 e-st df (t ) ใชหาปรพิ ันธบ ายพาส โดย ให u = e-st du = -se-stdt จไุ รรตั นจ ินดา อรรคนติ ย
คณิตศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 301 ให dv = df (t) f = f (t) ดงั นั้น L ddt f (t ) = e-st f (t) + s0 f (t)e-stdt 0 = -f (0) + sF(s) เมื่อ f (0) คือคาเร่ิมตน (Initial) ของฟงกชันนั้น เชน ถา f (t) = i(t) จะไดวา f (0) = i(0) ก็คือ กระแสไฟฟา ไหลในวงจรเมือ่ เวลาเปน 0 นนั่ เอง ในรปู แบบทัว่ ๆ ไปจะอยใู นรปู L dn f (t) = snF(s) - sn-1 f (0) - sn-2 f (0) - ... - s0 f (n-1) (8.11) dt n 8.2.5 ฟงกชันการหาปริพันธ ถา f (t) = f (t)dt ซ่ึงจะพบฟงกชันเชนน้ีในการหาแรงดัน ตกครอม C เม่ือกระแส i(t) ไหลผา น หาคาลาปลาซทรานสฟอรม จาก L[ f (t)dt] = 0 f (t)dt.e-st dt ใชห าปรพิ นั ธบ ายพาสโดย ให u = f (t)dt du = f (t)dt ให dv = e-stdt v = - 1s e-st ดังนน้ั L[ f (t)dt] = - 1s e-st f (t)dt + 1s 0 f (t)e-stdt 0 1s F(s) = f (0)dt + s = f -1 (0) + F(s) s s เมื่อ f -1(0) = f (0)dt ในวงจรไฟฟาคา นจ้ี ะมคี า ชดั เจน 8.2.6 ฟงกช ันอมิ พลั ส ถา f (t) เปน อิมพัลส (Impulse function) หรอื เรยี กกันวาสญั ญาณอิมพัลส เปนสัญญาณที่มีความสําคัญกับวงจรไฟฟามาก จะเขียนเปนสัญลักษณวา (t) สัญญาณพิเศษนี้ มลี ักษณะดังรูปที่ 8.6 จุไรรัตนจ ินดา อรรคนติ ย
302 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอเิ ลก็ ทรอนกิ ส รูปที่ 8.6 สญั ญาณอิมพัลส (ทีม่ า : http://cnx.org/contents/30db824d-56bc-49cf-b42d-a38ff2fffd54@3) สัญญาณอิมพัลสจะมศี ูนยก ลางอยทู ี่ t = 0 และ f (t) = 01 -2 < t < 2 t ใดๆ มคี านอยมากแตถงึ อยา งไรพ้ืนท่ีไดส ัญญาณอมิ พัลสก็จะยงั คงเทา กบั หนง่ึ เสมอไมวา จะมีคาเทาใดก็ตาม ถา มีคานอยลงเร่ือยๆ นั่นคือความกวางของสัญญาณจะลดลงความสูง ของสญั ญาณก็จะสงู ขึ้นตาม แตพ น้ื ท่ีไดส ญั ญาณพลั สก จ็ ะยงั คงเทากบั หน่ึงเสมอไป คําจํากดั ความของอมิ พัลสคือ (t) = 0 ท่ี t 0 และ (t)dt = 1 ดังนั้น f (t) (t = f (0) (t)) เสมอ ดังน้นั ในการหาลาปลาซทรานสฟอรมของ(t)จะไดโ ดย จาก L (t) = 0 e-st (t)dt แทนคา t = 0 จะได = 01 (t)dt เม่ือ a เขาสู 0 จะมีผลเทา กับ = 0 (t)dt = 1 8.2.7 ฟง กชนั แกมมา (Gamma function) ฟงกชันแกมมานั้นเปนฟงกชันทางคณิ ตศาสตรที่เปนสวนขยายของฟงกชัน แฟกทอเรียลบนจํานวนเชิงซอน โดยฟง กชนั แกมมาเปนการเตมิ เตม็ ฟงกช ันแฟกทอเรียลของคา n ท่ีไมใชจ าํ นวนเตม็ บวกหรอื ศูนย ซึง่ จํานวนเชิงซอ น z ที่เปน สวนจริงมคี าเปน บวก จุไรรตั นจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 303 y f(x)=gamma (x) 8 x 6 46 4 2 -8 -6 -4 -2 2 -2 -4 -6 รปู ที่ 8.7 ฟงกช นั แกมมา ซึ่งนยิ ามของฟง กชันแกมมาบน [0,∞) เขยี นแทนดวยสัญลักษณ มนี ิยามดงั น้ี () = t-1e-t dt , > 0 0 โดยฟง กช นั แกมมาจะลเู ขา เมื่อ - 1 > -1 หรือ > 0 สว นการคํานวณคา ของฟง กชนั แกมมา สามารถคํานวณคาไดจ าก (1) = t-1e -t dt = -e-t =1 0 0 และสาํ หรบั จํานวนเต็มทีไ่ มเ ปนลบใดๆ และ α 1 แลว ( + 1) = te-t dt 0 = -te-t + t -1e -t dt = () 0 0 ดังน้นั ถาแทน ดวยจํานวนเต็ม n ที่ไมเปน ลบใดๆ และ n 1 แลวจะได (n+1) = n(n) (8.12) ซึ่งเรียกสมการที่ 8.6 วาสูตรเวียนบังเกิด (Recurrence formula) เมื่อแทน n = 1,2 ,3…. จะได (2) = 1. (1) =1 (3) = 2. (2)=2.1=2! (4) = 3. (3)=3.2.1=3! จไุ รรัตนจ ินดา อรรคนติ ย
304 บทที่ 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส . . . สาํ หรับคา n เปนจาํ นวนเต็มท่ไี มเ ปนลบใดๆ แลว (n+1) = n! เมื่อพิจารณาผลการแปลงลาปลาซของฟงกชัน f (t) = t n จากนิยามของลาปลาซทรานส ฟอรมจะได L[t n ] = e- st t ndt 0 เมอ่ื u = st ดังนน้ั du = sdt และเม่อื t 0 , u 0 และ t ∞ , u ∞ จะได L[tn] = e -u un n dsu 0 = 1 e-u u ndu sn+1 0 = 1 (n + 1) = n! sn+1 sn+1 L[tn ] (n+1) n! ดังนั้น = sn+1 หรอื L[ tn ]= sn+1 เมอื่ n =1 ,2 ,3 ……… และ 1 = (8.13) 2 ตัวอยางที่ 8.2 จงหาผลการแปลงลาปลาซของฟง กชนั ตอไปนี้ L [ t-12 ] (n+1) sn+1 วิธีทํา ขนั้ ตอนท่ี 1 หาลาปลาซทรานสฟอรม โดยใชส ูตร L[tn ] = เมอ่ื n = M n = - 21 - ขน้ั ตอนท่ี 2 แทนคา ในสมการจะได (-12+1) s(1221+1 ) ดังน้ัน L[t-12 ] = = s12 จุไรรตั นจ ินดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 305 ขน้ั ตอนที่ 3 จากสมการท่ี (8.7) จะได 21 = แทนคาจะได = ( 12 ) = s s12 ตอบ L [ t- 12 ] = s ตัวอยางท่ี 8.3 จงหาลาปลาซทรานสฟอรม ของฟงกช นั ตอ ไปนี้ L[t 23 ] วธิ ที ํา ข้ันตอนที่ 1 หาลาปลาซทรานสฟอรมโดยใชส ตู ร L[tn ] = (n+1) sn+1 3 เมื่อ n = 2 จะได L[t 23 ] = ( 23+1) s 32+1 ขน้ั ตอนท่ี 2 จากสมการท่ี (8.6) เมื่อ (n+1) = n(n) จะได ( 23+1) 23 ( 23 L[t 32 ] = s 32+1 = s 25 ) ข้ันตอนท่ี 3 จาก 23 = 23 = 12 + 1 จะได L[t 32 ] = 23 ( 23 ) = 23 ( 21 +1) s 25 s25 ใชสมการ 8.12 จะได 23 ( 21 +1) 23 21 ) 21 L[t 32 ] = s25 = ( )( s 52 ( ) จากสมการที่ (8.13) 12 = 32 )( 12 จะได L[t 23 ] = ( s 52 ) = 3 = 3 4 4 s3 . s 52 s = 4 3 s3 3 s 4 s3 ตอบ L[t 23 ] = s จไุ รรัตนจ นิ ดา อรรคนติ ย
306 บทที่ 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอิเล็กทรอนกิ ส แลวสรุปการทําลาปลาซทรานสฟอรมของฟงกชันรูปแบบตางๆ และนําลงในตาราง ของลาปลาซทรานสฟอรม ตารางท่ี 8.1 เปนคุณสมบัติของลาปลาซทรานสฟอรมในลกั ษณะตางๆ จะไดดังน้ี ตารางที่ 8.1 ลาปลาซทรานสฟอรม อันดับท่ี F (t) ท่ี t > 0 F(s) หมายเหตุ 1 Unit step function 1 u(t) s Impulse function 1 Exponential decay 2 (t) 1 s+a function 3 e-at Ramp function n! 4 tn sn+1 Damped ramp function 5 tne-at n! 6 sin t ( s+ a )n +1 Sine function 8 cos t 9 e-at sin t Cosine function 10 e-at cos t s2 +2 Exponentially decaying s sine function 11 sin (t+) 12 cos (t+) s2 +2 Exponentially decaying cosine function 13 sinh t 14 cosh t Phase shift sine function (s+a)2 + 2 15 e(+j)t Phase shift cosine s+a function 16 e(-j)t (s+a)2 +2 Hyperbolic sine function 17 t s sin +cos s2 +2 Hyperbolic cosine function s cos-sin s2 +2 Complex exponential function Complex exponential s2 - 2 function s Unit ramp function s2 - 2 1 s-- j 1 s-+ j 1 s2 จุไรรตั นจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 307 ตารางท่ี 8.1 ลาปลาซทรานสฟ อรม (ตอ) อันดบั ที่ f (t) ที่ t > 0 F(s) หมายเหตุ 1 Complex exponential 18 e(j)t s- j function 19 e(-j)t s+1j Complex exponential 20 A A decay function s Constant function 8.3 คุณสมบัติของลาปลาซทรานสฟอรม ในการหาลาปลาซทรานสฟอรมของฟงกชันบางฟงกชนั น้ันคอนขางยุงยากในการหาปริพันธ แตสามารถคํานวณไดง า ยขน้ึ โดยอาศัยตารางท่ี 8.1 และคุณสมบตั ิตา งๆ ท่ีจะพิจารณาตอ ไปน้ีมาใช คาํ นวณ 8.3.1 คณุ สมบัติเชงิ เสน (Linearity property) ถา a1 , a2 เปน ตวั คงที่ และ L[ft (t)] = f2 (s) , L[f2 (t)] = f2 (s) ดังน้นั L[a1 f1 (t) + a2 f2 (t)] = a1 L [f1(t)] + a2 L [f2(t)] L[a1 f1 (t) + a2 f2 (t)] = a1 F1(s) + a2 F2(s) 8.3.2 คณุ สมบัติการสเกล (Scaling property) ถา L[f (t)] = F(s) (8.14) ดังน้นั L[f (at)] = 1a F as ตวั อยา งที่ 8.4 จงหา L[sin 6t] วธิ ีทํา ขน้ั ตอนที่ 1 ใชค ณุ สมบตั กิ ารแปลงลาปลาซจากตารางที่ 8.1 ลําดับที่ 6 จะได ข้นั ตอนที่ 2 1 L[sin t] = s2 +1 = F(s) จากโจทย L[sin 6t] ใชคณุ สมบตั ิการสเกลจากสมการ (8.14) จะได ข้นั ตอนท่ี 3 L[f (at)] = 1a F as เมือ่ a = 6 จะได แทน s ดว ย 6s จะไดใ นสมการจะได L[sin 3t] = 16 F 3s จุไรรตั นจ นิ ดา อรรคนติ ย
308 บทที่ 8 การแปลงลาปลาซ คณิตศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส = 16 1 (6s )2 +1 16 1 16 1 36 6(s2 +36) = 3s26 = s23+636 = +1 6 = (s2 +36) ตอบ L[sin 3t] = 6 (s2 +36) 8.3.3 คณุ สมบตั กิ ารเลือ่ นทางเวลา (Time-shifting property) ถาฟงกช ันถูกเล่ือนเวลาในลกั ษณะ ดงั รูปท่ี 8.8 y x รปู ที่ 8.8 ฟง กช นั ยนู ิตสเต็ปท่ถี ูกเลอ่ื นไปที่เวลา t0 ฟงกช ันทถี่ ูกเลอ่ื นออกไปจะเขียนไดในรูป f (t-)u(t-) เมื่อ u(t-) เปนยูนติ สเตป็ ฟง กช ัน ซ่ึงจะมคี า เปน 0 เมื่อ t < และเปน t < ดังนัน้ การหาลาปลาซทรานสฟอรม ของฟง กชันท่ีถูกเลื่อน เวลาทําไดโ ดย จาก L f (t - )u(t - ) = 0 f (t - )u(t - )e-st dt = 0 f (t - )e-st dt ถา ให t - = x = 0 f (x) e-s(+ x)dx จะเหน็ ไดวาเม่อื ฟงกชันอยใู นรปู ของ x เวลาก็ถกู เรมิ่ ที่ 0 ใหม ไมใชท่ี อีกตอ ไป L f (t - )u(t - ) = e-st 0 f (x) e-sxdx จไุ รรตั นจ ินดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 309 = e-st F(s) (8.15) ถา L[f (t)] = F(s) ดังนน้ั L[f (t - t0) = F(s) e-st0 ตัวอยางที่ 8.5 จงหาคา L[u(t-3)] วิธที าํ ขัน้ ตอนที่ 1 จากตารางคุณสมบัตกิ ารแปลงลาปลาซตารางท่ี 8.1 ขน้ั ตอนที่ 2 L[u(t)] = 1s = F(s) จากโจทย L[u(t-3)] ใชค ุณสมบตั กิ ารเลอ่ื นทางเวลา ตามสมการท่ี (8.15) เม่อื t0 = 3 จะได L[u(t - 3)] = F(s) e-3s L[u(t - 3)] = 1s e-3s ตอบ L[u(t - 3)] = 1s e-3s 8.3.4 คุณสมบตั กิ ารเลอื่ นทางความถ่ี (Frequency - shifting property) ถาฟงกชันนั้นถูกเปล่ียนความถ่ี (Frequency shift) ในลักษณะ e-at f (t) ลาปลาซทรานส ฟอรม จะหาไดจ าก L[e-at f (t)] = 0 e-at f (t)e-st dt = 0 f (t )e-(s+a)t dt = F(s + a) L[ f (t)e-at ] = F(s + a) (8.16) ตวั อยางท่ี 8.6 ถา ฟง กชัน g(t) ดงั แสดงในรปู ท่ี 8.9 จงหาลาปลาซทรานสฟอรมของ g(t) y c-12 รูปที่ 8.9 เม่ือ g(t) = e-4tu(t-3) จุไรรัตนจ นิ ดา อรรคนติ ย
310 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส วธิ ีทาํ ขัน้ ตอนท่ี 1 หาลาปลาซทรานสฟ อรมของ F(s) = L[u(t)] = 1s ข้ันตอนที่ 2 ใชหลักการเล่อื นเวลาเม่อื t = 3 จะได ขนั้ ตอนที่ 3 L[f (t - t0) = F(s) e-st0 F(s) = L[u(t - 3)] = ( 1s )e-3s จากนนั้ ใชหลกั การเปลย่ี นความถี่ เม่ือมี e- 4t จะไดว า L[e-atf (t)] = F(s + a) = e-st0 G(s) = L[e-4tu(t - 3)] = F(s + a) แทน s = s+a เมื่อ a = 4 จะได s = s + 4 ดงั น้ัน G(s) = e-3(s+4) s+14 = e-s3+s+412 ตอบ L[g(t)]= G(s) = e-s3+s+412 ตัวอยางที่ 8.7 จงหา L[sin 3te-t] วิธที ํา ขั้นตอนที่ 1 ใชตารางการเปลี่ยนคาลาปลาซทรานสฟอรม จากตารางที่ 8.1 ขัน้ ตอนที่ 2 a ขั้นตอนที่ 3 L[sin at] = s2 +a2 = F(s) จากโจทยเมือ่ a = 3 จะได L[sin 3t] = 3 = F(s) s2 +32 ใชค ณุ สมบตั กิ ารเล่อื นทางความถ่ี จากโจทยใหหา L[sin 3te-t] สาํ หรบั การเล่ือนทางความถ่ีมีคา เปน (-1) ดงั นนั้ แทน s = s - (-1) ในสมการท่ี (8.16) จะได 3 3 L[sin 3te-t] = ( s-(-1))2 +32 = (s +1)2 +9 L[sin 3te-t] = s2 + 3 9 = s2 + 3 2 s +1+ 2s+10 L[sin 3te-t] 3 = s2 + 2 s+10 ตอบ L[sin 3te-t] = s2 3 + 2 s+10 จุไรรัตนจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอเิ ล็กทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 311 8.3.5 คุณสมบัตกิ ารหาอนพุ ันธทางความถี่ (Frequency differentiation property) ในกรณีท่ีตองการหาลาปลาซทรานสฟอรมของผลคณู ของ t กบั f (t) ซงึ่ เรียกวา a ทําได ดังน้ใี หพจิ ารณาจากสมการอนุพันธของ F(s) ตอ s แต F(s) = L[ f (t)] = 0 f (t)e-stdt dF (s) df (t ) ดังนน้ั ds = 0 ds e- st dt = 0 -tf (t)e-stdt ใชก ารหาอนพุ นั ธผ ลคณู ธรรมดา = L[tf (t)] dF (s) ds เขยี นสลับใหมจ ะได L[tf (t)] = ในรูปแบบท่วั ไปจะได L[tnf (t)] = (-1)n dn F(s) dsn ถา L[f (t)] = [F(s)] ดังนนั้ L[f (t)t] = -[F(s)] และ L[f (t)(t)n] = (-1)n dn F(s) = (-1)nFn(s) dsn L[f (t)(-t)n] = dn [F(s)] (8.17) dsn ตวั อยา งที่ 8.8 จงหา L[e3tt] วธิ ที าํ ขัน้ ตอนที่ 1 หาผลการแปลงลาปลาซจากตารางที่ 8.1 ข้นั ตอนที่ 2 L[f (t)] = [e3t] = s1- 3 จากโจทย L[e3tt] คือ คา L[f (t)(-t)n] เมื่อ n = 1 ใชค ุณสมบัตจิ ากสมการท่ี (8.17) จะได L[e3tt] = (-1) dds s1-3 = - dds (s - 3)-1 1 = -(-1)(s – 3)-2 = (s-3)2 ตอบ L[e3tt]= 1 (s-3)2 8.3.6 คณุ สมบัตไิ ทมคอนโวลูช่นั (Time convolution) (8.18) L[f1(t)f2(t)] = [F1(s)F2(s)] จุไรรัตนจ ินดา อรรคนติ ย
312 บทที่ 8 การแปลงลาปลาซ คณิตศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส 8.3.7 คณุ สมบัติฟรเี ควนซีคอนโวลูชน่ั (Frequency convolution) (8.19) L[f1(t)f2(t)] = 21 [F1(s)F2(s)] (8.20) ซึ่ง F1(s)F2(s) = -+jj F1 ()F2 (s - )d (8.21) 8.3.8 คุณสมบตั กิ ารหาอนพุ ันธท างเวลา (Time - differentiation property) ถา L[f (t)] = [F(s)] L[ ddt [f (t)]] = L[f(t)] = sF(s) - f (0) dn L[ dtn [f (t)]] = L[f n (t)] = snF(s) - sn-1f (0) - sn-2f (0)… f (n-1)(0) ตัวอยางที่ 8.9 ถา f (t) = cos 2t จงหา L[f(t)] วิธีทาํ ขั้นตอนที่ 1 ใชคุณสมบัตกิ ารหาอนุพนั ธท างเวลาจะได ขัน้ ตอนที่ 2 L[ ddt [f (t)]] หา L[f(t)] = L[f(t)] = sF(s) - f (0) = L[ ddt (cos 2t) = L(-2 sin 2t) หาคาลาปลาซทรานสฟอรม ของ sin at จะได a L[sin at] = s2 +a2 = F(s) ดงั นน้ั L[f(t)] = -2L[sin 2t] = -2 2 s2 +4 L[f(t)] =- 4 s2 +4 ตอบ 4 L[f(t)] = - s2 +4 8.3.9 คณุ สมบัติการหาปริพันธท างเวลา (Time - integration property) ถา L[f (t)] = [F(s)] F(s) s ดงั นั้น L[ 0t f (t)dt ] = ตวั อยา งท่ี 8.10 จงหา L[ 0t cos 2tdt ] = s วธิ ีทาํ ข้นั ตอนท่ี 1 หาคา L[cos at] s2 +a2 จไุ รรตั นจ ินดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 313 จากโจทย L[cos 2t] = s2 s = s = F(s) +22 s2 +4 ขัน้ ตอนที่ 2 ใชค ณุ สมบตั กิ ารหาปรพิ นั ธท างเวลา จากสมการ (8.21) L[ 0t f (t)dt ] = F(s) s s 1 L[ 0t cos 2tdt ] = s(s2 +4) = s2 +4 ตอบ L[ 0t cos 2tdt ] = 1 s2 +4 8.3.10 คุณสมบตั กิ ารหาปริพนั ธท างความถี่ (Frequency- integration property) ถา L[f (t)] = [F(s)] ดังนัน้ L f (t ) = s F(s) (8.22) t ตารางท่ี 8.2 คุณสมบัติของลาปลาซทรานสฟ อรม อนั ดบั ท่ี f (t) ที่ t > 0 F(s) หมายเหตุ Linearity 1 a1f1(t)+a2f2(t) a1F1(s)+a2F2(s) Scaling Time shift 2 f (at) 1a F as Frequency shift Time differentiation 3 f (t -)u(t -) e-st F (s) Time integration 4 e-at f (t) F(s + a) Frequency 5 ddt f (t) differentiation sF(s) - f (0) Frequency Integration d2 f (t) s2 F(s) - sf (0) - f (0) dt 2 dn f (t) snF(s) - sn-1 f (0) - sn-2 f (0) - ... - f n-1 (0) dt n 6 [ ft]dt F(s) + f -1(0) s s 7 tn f (t) dnF(s) (-1)n dsn 8 f (t) s F(s)ds t จุไรรตั นจ ินดา อรรคนติ ย
314 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอิเลก็ ทรอนกิ ส ตารางท่ี 8.2 คณุ สมบตั ิของลาปลาซทรานสฟ อรม (ตอ) อนั ดบั ท่ี f (t) ท่ี t > 0 F(s) หมายเหตุ 9 f (t) = f (t + nT ) F1 ( s ) Time periodic 1-e-st 10 f (0+) lim sF(s) Initial value 11 f () Final value 12 f1(t)f2(t) s Convolution lim sF(s) s0 F1(s)F(s) ตัวอยา งท่ี 8.11 จงหาลาปลาซทรานสฟอรม ของ f (t) = sinh t วธิ ที ํา ขน้ั ตอนท่ี 1 เนอ่ื งจากในตารางการแปลงลาปลาซไมม คี าการแปลงลาปลาซของ f (t) = sinh t ดงั นน้ั จงึ เปลีย่ น sinh t = et -2e-t ข้นั ตอนท่ี 2 ทาํ การหาคา ลาปลาซทรานสฟอรม จากตาราง = L et -2e-t L[sinh t] = 12 L[et - e-t ] L[eat] = s 1- a จะได เม่ือ a = และ - จะได L[sinh t] = 12 L[et - e-t ] = 21 s-1 - s+1 12 (ss+-)-((ss+-)) = 21 ss-+)(-ss++ = ( ) = 12 2 = s2 +2 s2 + 2 ตอบ L[sinh t]= s2 + 2 จไุ รรัตนจ ินดา อรรคนติ ย
คณิตศาสตรว ศิ วกรรมอเิ ลก็ ทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 315 ตวั อยางที่ 8.12 จงหาลาปลาซทรานสฟ อรม ของ f (t) = (t) + u(t) - 2e-3t ที่ t > 0 วิธีทํา ขนั้ ตอนที่ 1 จากโจทยก ําหนดใหใ ชค ณุ สมบัตกิ ารเปน เชงิ เสนตามตารางที่ 8.2 จะได L[a1f1(t) + a2f2(t)] = a1F1(s) + a2F2(s) F(s) = L[(t)] + L[u(t)] - 2L[e-3t ] ข้นั ตอนที่ 2 หาคา ลาปลาซทรานสฟอรม จากตารางท่ี 8.1 จะได ขั้นตอนท่ี 3 L[(t)] = 1 L[u(t)] = 1s L[e-at] = s+1a แทนคา ในสมการจะได F(s) = L[(t)] + L[u(t)] - 2L[e-3t ] = 1 + 1s - 2 s+13 = s2s(+s2+s3+)3 ตอบ L[f (t)] = s2s(+s2+s3+)3 ตัวอยา งท่ี 8.13 จงหาลาปลาซทรานสฟอรม ของ f (t) = t2sin 3t วธิ ที ํา ข้ันตอนท่ี 1 จากโจทยก าํ หนดให f (t) = t2sin 3t เปนฟง กชนั ทอ่ี ยใู นรปู L[tnf (t)] ใชคณุ สมบตั ิ Frequency differentiation ตามตารางที่ 8.2 จะได L[t n f (t)] = (-1)n dn F(s) dsn ขน้ั ตอนท่ี 2 จากสูตร L[ f (t)] = L[sin at] = a ขั้นตอนท่ี 3 s2 +a2 3 จะได F(s) = L[sin 3t] = s2 +9 แทนคาใน L[t n f (t)] = (-1)n dn F(s) dsn จะได F(s) = L[t2sin 3t] เมอื่ n = 2 จะได F(s) = (-1)2 dd2s 3 s2 +9 จุไรรัตนจ นิ ดา อรรคนติ ย
316 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส = dds -6s (s2 +9)2 = 18s2 +54 (s2 +9)4 ตอบ L[f (t)] = 18s2 +54 (s2 +9)4 ตัวอยางท่ี 8.14 จงหาลาปลาซทรานสฟอรมของฟงกชนั ตอไปน้ี f (t) = t2e4t วิธที ํา ขน้ั ตอนท่ี 1 จากโจทยก าํ หนด f (t) = t2e4t เปนคณุ สมบตั กิ ารเลอื่ นทางเวลาจาก ตารางที่ 8.2 จะได L[f (t)eat] = F(s - a) ข้นั ตอนท่ี 2 หาคา L[f (t)] = L[t n] = n! จากตารางที่ 8.1 ขนั้ ตอนที่ 3 sn+1 2! 2 ดังนน้ั L[f (t)] = L[t 2] = s2+1 = s3 = F(s) เมอ่ื n = 2 แทนคา จากตารางท่ี 8.2 L[f (t)eat] = F(s - a) เมอื่ a = 4 และ s = s - 4 L[f (t)] = L[t 2e4t] = F(s - 4) L[f (t)] = L[t2e4t] = 2 (s-4)3 2 ตอบ L[t 2e4t] = (s - 4)3 ตัวอยา งท่ี 8.15 จงหาลาปลาซทรานสฟอรมของฟง กชันตอ ไปนี้ f (t) = 3e7t + 6t 3- 4 sin 3t + 3cos 3t วิธีทาํ ข้นั ตอนท่ี 1 จากโจทยให L[f (t)] = L[3e7t + 6t 3 - 4 sin 3t + 3cos 3t] ใชคุณสมบตั กิ ารเปนเชิงเสน ตามตารางท่ี 8.2 จะได L[f (t)] = L[3e7t + 6t 3 - 4 sin 3t + 3cos 3t] = L[3e7t] + L[6t3] - L[4 sin 3t] + L[3cos 3t] ขน้ั ตอนที่ 2 ใชค ุณสมบัตกิ ารแปลงลาปลาซจากตารางท่ี 8.1 จะได L[eat] = s 1- a ดงั นน้ั L[3e7t] = s3- 7 จุไรรตั นจ นิ ดา อรรคนติ ย
คณิตศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 317 L[tn]= n! ดังนัน้ L[6t 3] =6 3! sn+1 s3+1 L[sin at] = a ดังน้นั L[4 sin 3t] = 4 3 s2 +a2 s2 +32 L[cos at] = s ดังน้นั L[3 cos 3t] = 3 s2 s s2 +a2 +32 ขั้นตอนที่ 3 แทนคา ลงในสมการจะได L[f (t)] = L[3e7t] + L[6t3] – L[4 sin 3t] + L[3 cos 3t] = s3-7 +6 3! -4 s2 3 +3 s2 s s3+1 +32 +32 = s3-7 + 36 - 12 + 3s s4 s2 +32 s2 +32 ตอบ 36 12 3s L[f (t)] = s3-7 + s4 - s2 +32 + s2 +32 ตัวอยา งท่ี 8.16 จงหาลาปลาซทรานสฟอรม ของสัญญาณตามรปู ท่ี 8.10 4 3 2 1 1234 5678 -1 รปู ท่ี 8.10 รปู คล่ืนสญั ญาณแบบสามเหลย่ี ม วิธที ํา ขั้นตอนที่ 1 วาดรูปสญั ญาณที่เกิดเปน รูปสามเหลีย่ ม ประกอบดว ยผลรวมของสาม ฟง กช ันดว ยกนั คอื f (t) = f1(t) + f2(t) +f3(t) 1.1 สมการเสน ตรงตัดแกน x และ แกน y ท่ีจุด (0,0) ซง่ึ เปน ฟง กช นั แรม มีสมการเสน ตรงคอื y = mx+c เมอ่ื m คอื ความชันของกราฟจะได y = f (t) = 45 t + 0 f1(t) = 45 t จไุ รรตั นจ นิ ดา อรรคนติ ย
318 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอเิ ล็กทรอนกิ ส f 1( t ) 5 4 3 2 1 t 12345 1.2 สมการเสนตรงท่ีมีจุดเริ่มตน คือ (4 , 0) คาความชันของกราฟมีคาเปน - 45 t ซงึ่ กราฟมีการเลอ่ื นทางเวลาไป 2 วินาทีจะได (t - 4) สมการมคี าเปน f2(t) = - 45 t (t-4) f2 (t) x y f2(t) = - 45 t (t - 4) 1.3 กราฟรูปที่ 3 เปนเสนตรง ขนานแกน y เมื่อ y = -5 จุดเร่ิมตนอยูท่ี (0,4) เปน ฟงกชันยูนิตสเต็ป ที่คาเวลาเลื่อนไปเริ่มตนท่ี t = 4 ดังนั้น f3(t) มีสมการเปน f3(t) = -5u(t - 4) f3 (t ) y 2 x -2 2 4 6 8 -2 -4 -6 ข้นั ตอนที่ 2 ทําการรวมสัญญาณจะไดเ ปน ข้ันตอนท่ี 3 f (t) = 45 t - 45 t(t - 4) - 5u(t - 4) ทําการแปลงลาปลาซทรานสฟอรม โดยใชต ารางคูแปลงลาปลาซ L[f1(t)] = L[t] = 1! = 1 s2 s2 จไุ รรัตนจ ินดา อรรคนติ ย
คณิตศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 319 ใชคณุ สมบตั กิ ารเลื่อนทางเวลา L[f (t - t0)] = F(s) e-st0 L[f2(t)] = L[t(t - 4)] = 1 e-4s s2 = L[u(t - 4)] = 1 e-4s L[f3(t)] s e-4 s -5 e-s4s จะได F(s) = 45 1 - 45 s2 s2 5 = 4s2 (1- e-4s- se-4s) ตอบ F(s) = 5 (1- e-4s- se-4s) 4s2 ตัวอยางที่ 8.17 จากวงจรสญั ญาณรูปฟน เลอ่ื ย ใหห าลาปลาซทรานสฟอรม ของวงจรในรูปดงั กลาว รูปที่ 8.11 รูปคลืน่ ของฟน เล่ือย y x รูปท่ี 8.12 รปู คล่ืน จุไรรัตนจ นิ ดา อรรคนติ ย
320 บทที่ 8 การแปลงลาปลาซ คณิตศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส วธิ ีทาํ ขนั้ ตอนท่ี 1 หาสมการของสญั ญาณรูปฟนเลือ่ ยเพยี งรูปเดียวซงึ่ จะประกอบดวย สัญญาณ 3 สญั ญาณรวมกนั ไดแ ก สญั ญาณแรมปส องสัญญาณและ ยนู ิตสเตป็ 1 สัญญาณ ดงั น้ี f (t) = f1(t) + f2(t) + f3(t) 1.1 สมการเสนตรงตัดแกน x และ แกน y ท่ีจุด (0,0) ซ่ึงเปนฟงกชันแรม มีสมการเสนตรง คือ y = mx + c เม่ือ m คือ ความชันของกราฟ และ C จุดตดั แกน y จะได y = f1(t) = At 1.2 สมการเสนตรงทม่ี จี ดุ เริม่ ตน คอื (T,0) คาความชันของกราฟมคี า เปน - TA t ซึง่ กราฟมกี ารเลื่อนทางเวลาไป T วินาที จะได (t - T) สมการมีคา เปน f2(t) = -t TA t(t-T) f2 (t) yx f2(t) = -t TA t(t - T) 1.3 กราฟรูปท่ี 3 เปนเสนตรง ขนานแกน y เมื่อ y = -A จุดเร่ิมตนอยูท่ี (0,4) เปนฟงกชันยูนิตสเต็ปท่ีคาเวลาเลื่อนไปเร่ิมตนที่ t = T ดังน้ัน f3(t) มีสมการเปน f3(t) = -Au(t-T) จไุ รรัตนจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอเิ ล็กทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 321 f3 (t) ขน้ั ตอนท่ี 2 ทําการรวมสัญญาณจะไดเ ปน ข้ันตอนท่ี 3 f (t) = TA t - TA t(t - T) -Au(t -T) ขั้นตอนท่ี 4 ทาํ ลาปลาซทรานสฟอรม โดยใชต ารางคูแปลงลาปลาซไดเ ปน e-2 s -A e-s2s F(s)= TA 1 - TA s2 s2 A = Ts2 (1- e-TS-Tse-TS) จากรปู ที่ 8.10 สัญญาณน้ีจะซ้าํ เร่ือยไปในเวลา T จากหลกั การของ ฟงกช ันเลอื่ นเวลาจะได f (t) = f1(t) + f1(t -T)t(t-T) + f1(t-2T)t(t-2T) +… แปลงลาปลาซจะได L[f (t)] = L[f1(t)] + L[f1(t - T)t(t-T)] + L[f1(t - 2T)t(t-2T)} +… และจาก L[f (t)] = F(s) L[f (t-T)t(t -T)] = F(s)e-Ts ดังน้นั จากสมการจะไดเปน F(s) = F1(s)+F1(s)e-TS+F1(s)e-2Ts+F1(s)e-3Ts+… และจาก Infinite series คอื 11-x = 1 +x + x2 + x3+ … ดงั น้ัน 1 = 1+e- TS+e-2TS+e-3TS +…+e-nTS 1-e-Ts จากสมการ F(s) = F1(s)+F1(s)e-TS+ F1(s)e-2Ts+F1(s)e-3Ts+… จะได F(s) = F1 ( s ) 1-e-Ts = A (1 - e -Ts - Tse-Ts ) (1-e-Ts )Ts2 จไุ รรตั นจ ินดา อรรคนติ ย
322 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอิเล็กทรอนกิ ส = A (1 - Tse-Ts ) Ts2 1-e-Ts = A 1 - e-Ts ) Ts2 s(1-e-Ts ตอบ F(s) = A 1 - e-Ts Ts2 s(1-e-Ts ) ตัวอยา งที่ 8.17 จงหาลาปลาซทรานสฟอรม ของรปู คลน่ื แรงดันขาออกของวงจรเรยี งกระแสแบบ ครงึ่ คลื่นเมอื่ แรงดนั อนิ พตุ มสี ัญญาณรูปไซน รปู ท่ี 8.13 วงจรเรยี งกระแสแบบครึ่งคล่นื และสัญญาณเอาตพ ตุ (ที่มา : www.CircuitToday.com) รูปที่ 8.14 สญั ญาณพลั สร ปู ไซน วิธีทาํ ขัน้ ตอนที่ 1 วาดรปู คลืน่ ไซนที่มจี ดุ เริ่มตนของสญั ญาณที่ตาํ แหนง (0,0) จะได สมการของฟงกชนั f1(t) = A sin t จไุ รรัตนจ ินดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอเิ ลก็ ทรอนกิ ส บทท่ี 8 การแปลงลาปลาซ 323 f1(1) y f1(t) = Asin t 6 4 x 2 -2 -4 -6 ขั้นตอนที่ 2 วาดรูปคลื่นไซนท่ีมจี ดุ เรมิ่ ตน ของสญั ญาณทต่ี ําแหนง (,0) จะไดสมการ ของฟง กชนั f2(t) = A sin (t -) f2 (t) f2 (t) = A sin (t - ) ข้นั ตอนท่ี 3 นํารปู สญั ญาณทัง้ สองมาบวกกนั จะได f (t) = f1(t) + f2(t) = A sin t + A sin(t - ) L[f (t)] = L[A sin t +A sin(t - )] = L[A sin t] + L[A sin(t - )] ใชคุณสมบัตกิ ารเปลี่ยนกลับของฟงกชนั ไซนแ ละการเล่อื นทางความถี่ จะได a L[sin at]= s2 +a2 เมอ่ื a = 1 จะได L[sin t] = 1 s2 + 1 1 e -Ts และ L[sin(t -T)] = s2 +t2 L[sin(t -)] = 1 e -s s2 +1 A A .e-s F(s) = s2 +1 + s2 +1 จไุ รรตั นจ นิ ดา อรรคนติ ย
324 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ศิ วกรรมอเิ ลก็ ทรอนกิ ส = A (1 + e-s ) s2 +1 =A1+e-s s2 +1 ตอบ 1+e-s F(s) = A s2 +1 ตัวอยางท่ี 8.18 จงหาลาปลาซทรานสฟอรมของสญั ญาณตามรูปท่ี 8.15 0.5 1 1.5 2 2.5 3 3.5 รูปท่ี 8.15 รปู คล่นื สัญญาณแบบสเ่ี หลี่ยม วิธที าํ ขัน้ ตอนท่ี 1 วาดรปู สญั ญาณท่ีเกิดเปน รูปสเี่ หล่ียม ประกอบดว ยผลรวมของสอง ฟงกช ันดว ยกนั คอื f (t) = f1(t) + f2(t) 1.1 สมการเสนตรงขนานแกน x จุดเริ่มตน ที่ (0,5) หรอื ฟงกชนั ยนู ติ สเตป็ y = f1(t) = 5u(t) y 8 6 4 2 x -1 12 34 567 1.2 สมการเสน ตรงขนานแกน x จดุ เรม่ิ ตน ท่ี (4,-5) หรือฟง กช ันยูนติ สเตป็ ทมี่ คี า เวลาเล่ือนไป 4 วนิ าที ความสงู ของคลื่น -5 y = f2(t) = - 5u(t - 4) จุไรรัตนจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอเิ ลก็ ทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 325 y x -2 2 4 6 8 10 12 14 -2 -4 -6 ขัน้ ตอนที่ 2 ทาํ การรวมสญั ญาณจะไดเ ปน f (t) = 5u(t) - 5u(t - 4) ข้ันตอนท่ี 3 ทําการแปลงลาปลาซทรานสฟอรม โดยใชตารางคแู ปลงลาปลาซ L[f1(t)] = L[u(t)] = 1s ใชคุณสมบตั กิ ารเล่อื นทางเวลา L[f (t - t0) = F(s) e-st0 L[f2(t)] = L[u(t - 4)] = 1s e-4s จะได F(s) = 5s - 5s e-4s = 5s (1- e-4s) ตอบ F(s) = 5s (1- e-4s) 8.4 สรุป การแปลงลาปลาซหรือผลการแปลงลาปลาซ เปนสมการที่สําคัญในการแกปญหา ทางวิศวกรรมไฟฟา ซึ่งเปนสมการอนพุ ันธอนั ดับหนึง่ และอันดับท่ีสูงกวา 1 การแกส มการอนุพนั ธ สามารถทําไดโดยวิธีแบบตางๆ ซ่ึงมีความยุงยากในการแกสมการ โดยการแปลงลาปลาซจะเปน การแปลงสมการจากโดเมนเวลาเปนโดเมนความถี่เชิงซอนหรือโดเมน S จากนั้นจึงทําการ แกสมการโดยใชหลักการทางพีชคณิต ซ่ึงสามารถทําการแปลงลาปลาซได 2 วิธี คอื โดยใชตาราง การแปลงลาปลาซโดยตรง และโดยออมซ่ึงจะตองทําการจัดรูปสมการใหอยูในรูปแบบ ตามคุณสมบัติของการแปลงลาปลาซกอน เชน คุณสมบัติการเปนเชิงเสน คุณสมบัติการเล่ือน ทางเวลา คุณสมบัติการเลอื่ นทางความถี่ คุณสมบัติการสเกล เปนตน แลวจึงทําการแปลงลาปลาซ ตามตารางคุณสมบัติ จุไรรตั นจ นิ ดา อรรคนติ ย
326 บทที่ 8 การแปลงลาปลาซ คณิตศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส แบบฝกหัดทายบทท่ี 8 จงหาลาปลาซทรานฟอรมของฟงกชันตอไปนี้ F(s) = 3s + s+24 1. f (t) = 3u(t) + 2e-4tu(t) ตอบ 2. f (t) = 3e5t+ 8 sinh 3t + 6 ตอบ F(s) = s3-5 + 24 + 6s s2 -9 3. f (t) = 7t – 5t4 +4 te8t 7 4! 4 ตอบ F(s) = s2 - 5 s5 + ( s-8)2 4. f (t) = 3t2 (t) + 2e-t ตอบ F(s) = s3 + 2 s+2 + s+24 5. f (t) = 2cos 3t ตอบ s3 ( s +1) 6. f (t) = 7sin 6t - 9cos 4t +10u(t) ตอบ 2s 7. f (t) = 5sin 2t ตอบ F(s) = s2 +9 8. f (t) = 2sin 2t - cos 2t ตอบ F(s) = 42 - 9s + 1s0 s2 +36 s2 +16 10 F(s) = s2 + 4 F(s) = 4-s s2 +4 s2 9. f (t) = (t2 + 2)2 ตอบ F(s) = 4 +8s+24 10. f (t) = 2cos2 2t ตอบ s5 11. f (t) = e-tcos 2t ตอบ F(s) = 1s s 12. f (t) = t3e-3t ตอบ + s2 +16 ตอบ 13. f (t) = 5(t) + 7u(t) - 9e-7t ตอบ F(s) = s+1 14. f (t) = (t - 4)2 (s2 +2s+5) 6 F(s) = ( s+3)4 F(s) = 5 + 7s - s+97 2! e-4 s F(s) = s3 15. f (t) = t2 sin 2t ตอบ F(s) = 12 s2 -16 16. f (t) = 3(t) - 5u(t) - 3e-4t ตอบ (s2 +4)3 F(s) = 3- 5s - s+34 จไุ รรตั นจ นิ ดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอิเลก็ ทรอนกิ ส บทที่ 8 การแปลงลาปลาซ 327 17. จงหาลาปลาซทรานสฟ อรม ของรปู คลืน่ สญั ญาณตอไปนี้ 18. จงหาลาปลาซทรานสฟ อรม ของรปู คลืน่ สญั ญาณตอไปนี้ f1 (t) y 19. จงหาลาปลาซทรานสฟ อรมของรูปคลนื่ สญั ญาณตอ ไปนี้ f (t) 7 05 t 20. จงหาลาปลาซทรานสฟอรม ของรูปคล่ืนสญั ญาณตอไปน้ี จุไรรตั นจ ินดา อรรคนติ ย
328 บทท่ี 8 การแปลงลาปลาซ คณติ ศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส เอกสารอางอิง มงคล ทองสงคราม. (2542). การแปลงลาปลาซและการแปลงฟูเรยี ร. กรงุ เทพฯ : รามาการพิมพ. พมิ พค ร้ังที่ 2 ศรบี ตุ ร แววเจรญิ , ชนศกั ดิ์ บา ยเทีย่ ง. (2542). สมการเชิงอนุพนั ธ 1 Differential equations 1: คณติ ศาสตรส าํ หรับวศิ วกรรมและวทิ ยาศาสตร. กรงุ เทพฯ: สาํ นักพิมพว งตะวัน. ศรบี ตุ ร แววเจรญิ , ชนศักด์ิ บา ยเท่ยี ง. (2543). สมการเชิงอนุพันธ 2 และการแปลงลาปลาซ Differential equations 2 & laplace transform : คณติ ศาสตรส ําหรบั วิศวกรรมและ วทิ ยาศาสตร. กรงุ เทพฯ : สาํ นกั พมิ พว งตะวัน. นริ ันดร คาํ ประเสริฐ. (2538). คณติ ศาสตรวศิ วกรรมไฟฟา 4 : สมการเชงิ อนุพนั ธส ามญั สมการเชงิ อนุพนั ธยอ ย . กรุงเทพฯ : ศนู ยส ่ือเสรมิ กรงุ เทพฯ. Goyal Gupta . (2011). Laplace Transform and Fourier Transform .Pragati Prakashan. Joel L. Schiff . (1999). The Laplace Transform: Theory and Applications.Springer . จนิ ดา ไชยชว ย. Engineering Mathematics2. เขา ถงึ ไดจาก http://www.kmitl.ac.th/~kcchinda/ PDFs/inverseLP29_43.pdf .(วนั ทค่ี น ขอ มูล 5 มีนาคม 2556) การแปลงลาปลาซ. [ออนไลน] . เขา ถึงไดจ าก http://www.crc.ac.th/online/22201/312ch09_2550 _3rd.pdf. (วนั ทค่ี น ขอ มลู 6 พฤษภาคม 2556). การแปลงลาปลาซ. [ออนไลน] . เขา ถึงไดจาก http://ssru.ac.th/teacher/ubol/file.php/1/ Math_for_science /30_7.pdf. (วนั ท่ีคนขอ มลู 3 มนี าคม 2556). การแปลงลาปลาซ. [ออนไลน] . เขาถึงไดจาก http://www.vcharkarn.com/uploads/252/252771.pdf. (วันทคี่ นขอ มลู 8 กรกฎาคม 2556). ลาปลาซทรานสฟอรม. [ออนไลน] เขาถึงไดจ าก http://montri.rmutl.ac.th/old/ee/04210204mat /Presentation8.pdf. (วนั ที่คน ขอมูล7 มถิ นุ ายน 2556). Chater7 Laplace Transform . [ออนไลน] เขา ถึงไดจ าก http://www.math.utah.edu/~gustafso /laplace Transform .pdf. (วนั ทค่ี น ขอมูล12 พฤษภาคม 2556). S.boyd. Laplace Transform. [ออนไลน] เขาถึงไดจาก https://web.stanford.edu/~boyd/ ee102/laplace.pdf. (วนั ทคี่ น ขอมูล 5 เมษายน 2556) จุไรรตั นจ นิ ดา อรรคนติ ย
บทท่ี 9 ลาปลาซทรานสฟอรมกลับ 9.1 ความหมายของลาปลาซทรานสฟ อรม กลบั ลาปลาซทรานสฟอรมกลับหรือการแปลงลาปลาซกลบั เปนการหาฟงกช ัน f (t) เดมิ จากฟงกชัน F(s) L[f (t)] = F(s) (9.1) L[f (t)] = 0 f (t )e-st dt = F(s) (9.2) L-1[F(s)]= f (t) (9.3) ซึง่ L-1 เรียกวา ตัวดาํ เนินการลาปลาซทรานสฟอรมกลบั ( Inverse laplace transform) ซงึ่ f (t) อาจหาไดจากสมการ L-1[F(s)]= f (t) = 21j -+jj F(s)est ds (9.4) ตัวอยางท่ี 9.1 จงหาผลการแปลงกลบั ของลาปลาซทรานสฟอรมตอ ไปน้ี F(s) = s+14 วธิ ที าํ ขนั้ ตอนที่ 1 ใชส มการท่ี (9.4) ในการแกสมการจะได ขัน้ ตอนที่ 2 L-1[F(s)]= f (t) = 21j -+jj F(s)est ds แทน F(s) = s+14 ในสมการจะได f (t) = 21j -+jj s+14 estds = e-4t ตอบ L-1[F(s)] = e-4t
330 บทท่ี 9 ลาปลาซทรานสฟอรม กลบั คณิตศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส ตารางท่ี 9.1 ลาปลาซทราสฟอรมกลับ L-1[F(s)] = f (t) 1 ลาํ ดบั ที่ L[f (t)] = F(s) t 1 1s tnn! 2 1 eat 3 s2 4 1 sin at 5 s n +1 cos at 6 s1-a sinh at 7 a cosh at 8 e(a+j)t 9 s2 +a2 e(a-j)t 10 s 11 tn 12 s2 +a2 teat 13 a t neat 14 e-atsin 0t s2 -a2 15 s e-atcos 0t 16 s2 -a2 (t) s-a1- j s-a+1 j n! sn+1 1 (s-a)2 n! (s-a)n+1 0 (s+a)2 +02 ) s+a (s+a)2 +02 ) 1 จุไรรัตนจ ินดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเลก็ ทรอนกิ ส บทที่ 9 ลาปลาซทรานสฟอรมกลบั 331 9.2 คณุ สมบตั ขิ องลาปลาซทรานสฟอรมกลบั 9.2.1 คณุ สมบัตกิ ารเปนเชิงเสน (Linearity property) L-1[aF1(s) + bF2(s)] = aL-1[F1(s)] + bL-1[F2(s)] (9.5) L-1[aF1(s)+bF2(s)] = af1(t) + bf2(t) ตวั อยา งที่ 9.1 จงหาการแปลงกลับของฟง กชนั ตอไปน้ี L-1 s5-3 - 2s + 4 s2 +36 s2 +16 วธิ ีทาํ ข้นั ตอนที่ 1 ใชค ณุ สมบัตกิ ารเปน เชงิ เสน จากสมการท่ี (9.5) จะได L-1[aF1(s)+ bF2(s)] = af1(t)+bf2(t) L-1 s5-3 - 2s + 4 = 5L-1 s1-3 - 2L-1 s2 s s2 +36 s2 +16 +36 + L-1 s2 4 +16 ขนั้ ตอนที่ 2 ใชต ารางคณุ สมบัตกิ ารแปลงกลับลาปลาซจากตารางท่ี (9.1) จะไดพ จน แรก เมอ่ื 5L-1 s1-3 ใชตารางการแปลงกลบั ลําดับท่ี 4 L-1 s1-a = eat เม่ือ a = 3 จะได 5L-1 s1-3 = 5e3t พจนที่สอง เมอื่ 2L-1 s2 s ใชต ารางการแปลงกลับลําดบั ท่ี 6 +36 L-1 s2 s = cos at เมือ่ a = 6 จะได +a2 2 L-1 s2 s = 2cos 6t +36 พจนทีส่ าม เมอื่ L-1 s2 4 ใชต ารางการแปลงกลบั ลําดบั ที่ 5 +16 L-1 s2 s = sin at เม่อื a = 4 จะได +a2 L-1 s2 4 = sin 4t +16 จุไรรตั นจ นิ ดา อรรคนติ ย
332 บทที่ 9 ลาปลาซทรานสฟ อรม กลบั คณติ ศาสตรว ิศวกรรมอิเล็กทรอนกิ ส ตอบ L-1 s5-3 - 2s + 4 = 5e3t+2cos 6t + sin 4t s2 +36 s2 +16 9.2.2 คณุ สมบตั ิการเลอ่ื นทางความถี่ (Frequency - shifting property) (9.5) L-1[F(s - a)] = f (t)eat ตัวอยา งที่ 9.2 จงหา L-1 s2 1 -8s+20 วิธีทํา ขนั้ ตอนที่ 1 ทาํ สมการใหเ หมือนกับฟงกช นั F(s) ในตารางที่ (9.1) โดยใชสมการ กําลัง 2 สมบูรณได (s - a)2 = s2- 2as +a2 ข้ันตอนท่ี 2 จากโจทย s2- 8s+20 = (s2-2(4)s + 42) + 4 เมื่อ a = 4 จะได s2- 8s+20 = (s - 4)2+ 4 L-1 s2 - 1 = L-1 ( 1 + 4 8s+20 s-4)2 ข้ันตอนท่ี 3 เมอ่ื ทําเปน กาํ ลังสองสมบูรณเรยี บรอ ยแลว จดั สมการใหอยใู นรูป L-1 s2 s = sin at +a2 นําคา 22 คณู เขาไปในสมการจะได L-1 ( s - 1 + 4 = 22 L-1 ( 1 + 4 4)2 s-4)2 = 12 L-1 1 + 22 (s-4)2 คา a = 2 และจาก s = s - 4 ใชคณุ สมบตั กิ ารเล่ือนทางความถจ่ี ะได L-1 ( 1 +4 = 12 e4tsin 2t s-4)2 ตอบ L-1 s2 1 = 12 e4t sin 2t -8s+20 จุไรรตั นจ ินดา อรรคนติ ย
คณิตศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทท่ี 9 ลาปลาซทรานสฟ อรมกลับ 333 9.2.3 คุณสมบตั ิการเล่ือนทางเวลา (Time - shifting property) (9.6) L-1[F(s)e-as] = f (t - a) (9.7) ตวั อยา งที่ 9.3 จงหา L-1 e-( / 3) s s2 +1 วิธที าํ ข้ันตอนท่ี 1 จากสตู ร L-1[F(s)] = L-1 1 = sin t = f (t) ข้ันตอนที่ 2 s2 +1 ใชค ณุ สมบัตกิ ารเล่ือนทางเวลา L-1[F(s)e-as] = f (t - a) เมอ่ื - a = 3 แทนคา ในสมการจะได 1 -( / 3) s L-1 s2 +1 e = sin t - 3 1 -( / 3) s ตอบ s2 +1 L-1 e = sin t - 3 9.2.4 คุณสมบัตกิ ารสเกล (Scaling property) [F(as)] = a1 f at ตัวอยา งที่ 9.4 จงหา L-1 3s (3s)2 +36 s วธิ ที าํ ข้ันตอนท่ี 1 หา L-1[F(s)] เม่อื F(s) = s2 +62 จากตารางที่ 9.1 ลาํ ดับที่ 6 L-1 s2 s a2 = cos at + เมื่อ a=6 จะได L-1 s2 s = cos 6t = f (t) +62 ขั้นตอนท่ี 2 ใชคณุ สมบตั กิ ารสเกล L-1[F(as)] = a1 f at เม่ือ f (t) = cos 6t และ a = 3 จะได จุไรรตั นจ ินดา อรรคนติ ย
334 บทที่ 9 ลาปลาซทรานสฟอรม กลับ คณติ ศาสตรว ศิ วกรรมอเิ ล็กทรอนกิ ส L-1 (3s 3s = 13 cos 63t = 13 cos 2t )2 +36 ตอบ L-1 (3s 3s = 13 cos 2t )2 +36 9.2.5 ลาปลาซทรานสฟอรมกลบั อนุพันธ (Differential inverse laplace transform) ถา L-1[F(s)] = f (t) ดังน้ัน L-1[Fn(s)] = L-1 dn [ F ( s)] = (-1)nt nf (t) (9.8) dsn ตัวอยางท่ี 9.5 จงหา L-1 1 s(s2 +4) วธิ ที ํา ข้นั ตอนที่ 1 L-1[F(s)] = L-1 1 = L-1 s2 1 s2 +4 +22 นาํ คา 2 คูณท้ังเศษและสว น จะได L-1 s2 1 = 12 L-1 s2 2 +22 +22 f (t) = 12 sin 2t ขนั้ ตอนท่ี 2 ทําการหาอินเวอรสลาปลาซจากสตู ร L-1 1s [F ( s )] = t f ( t ) dt 0 L-1 1s 1 = 0t 21 sin 2tdt s2 +4 - 21 21 cos (2t) t เมอื่ 0 u = 2t จะได = แทนคา t = 0 และ t จะได = - 41 (cos 2t - cos 0) = - 41 (cos 2t - 1) L-1 1s 1 = 41 (1 - cos 2t) s2 +4 ตอบ L-1 1s 1 = 14 (1 - cos 2t) s2 +4 จุไรรัตนจ ินดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทที่ 9 ลาปลาซทรานสฟ อรม กลับ 335 ตารางท่ี 9.2 คุณสมบตั ขิ องลาปลาซผกพนั ลาํ ดบั ท่ี คุณสมบตั ิ F(s) L-1[F(s)] = f (t) aF1(s)+ bF2(s) af1(t)+ bf2(t) 1 เชงิ เสน F(s- a) f (t)eat 2 การเลอื่ นทางความถ่ี F(s)e-at f (t - a) 3 การเลอื่ นทางเวลา F(as) a1 f at 4 การสเกล (-1)ntnf (t) เมื่อ n =1 ,2 ,3… 5 อนุพนั ธ F n(s) ตัวอยา งที่ 9.6 จงหาคา ลาปาซทรานสฟอรมกลับของฟงกชนั ตอ ไปนี้ [F(s)] = 2 s6-3 - 3+4s 9s2 -16 วธิ ที าํ ข้นั ตอนที่ 1 หาคาลาปาซทรานสฟอรม กลับของแตละฟง กช ันจะได 6 เม่อื [F1(s)] = 2s-3 จะได L-1[F1(s)] = L-1 2s6-3 = L-1 22 s-323 = L-1 s -323 = 3e 23 t และ L-1[F2(s)] = L-1 - 3+4s 9s2-16 ทาํ การแยกเศษสว นจะได L-1[F2(s)] = L-1 - 3 - 4s 9s2 -16 9s2-16 = L-1 - 19 s2 3 - 19 4s -196 s2 -196 ขัน้ ตอนท่ี 2 ทาํ การแปลงกลับโดยใชตาราง จากตารางที่ 9.1 จะได L-1[F(s)] = L-1 a = sinh at s2- a2 L-1 - 19 s2 3 คา a2 = 196 -196 จไุ รรตั นจ นิ ดา อรรคนติ ย
336 บทที่ 9 ลาปลาซทรานสฟ อรมกลบั คณิตศาสตรว ิศวกรรมอเิ ล็กทรอนกิ ส ดังน้ัน a = 43 นาํ คา 43 คณู เขาไปในสมการทั้งเศษและสวนจะได L-1 - 91 s2 3 = L-1 - 43439 3 - 196 s2 -196 = L-1 - 4339 43 s2 -196 = L-1 - 4339 43 s2 -196 = L-1 - 41 43 s2 -196 = - 41 sinh 43 t สว นการหาคา L-1 - 91 4s คา a2 = 196 s2 -196 ดงั นนั้ a = 43 จากตารางท่ี 9.10 จะได L-1[F(s)] = L-1 s = cosh at จัดรูปสมการ s2- a2 L-1 - 19 4s = L-1 - 49 s2 s s2 -196 -196 = - 49 cosh 43 t ตอบ L-1 2 s6-3 - 3+4s = 3e 23t - 14 sinh 43 t - 49 cosh 43 t 9s2-16 ตวั อยางท่ี 9.7 จงหาคาลาปลาซทรานสฟอรม กลับของฟงกชนั ตอไปนี้ [F(s)] = s2 6s-4 -4s+20 วิธที ํา ขน้ั ตอนที่ 1 ใชก ารแกสมการแบบกาํ ลงั สองสมบรู ณ a2 - 2ab + b2 = (a - b)2 จากเศษของสมการ คือ s2- 4s+20 จะได a = 1 และ b = 2 ดังนน้ั s2- 4s + 20 = (s2 - 2 2s + 22) + 16 จุไรรตั นจ ินดา อรรคนติ ย
คณิตศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส บทที่ 9 ลาปลาซทรานสฟ อรมกลบั 337 = (s - 2)2+ 16 L-1[F(s)] = L-1 6s-4 s2-4s+20 = L-1 6s-4 (s-2)2 +16 ขั้นตอนที่ 2 เมื่อทาํ เปนกาํ ลังสองสมบรู ณเ รียบรอยแลว จัดสมการใหอ ยูในรูป ตอบ L-1 a = sin at s2 +a2 L-1 ( 6s-4 = L-1 6(s-2)+8 s-2)2 +16 (s-2)2 +16 ทําการแยกเศษสว นยอยเปน 2 พจน จะได = L-1 ( 6(s-2) + L-1 ( 8 +16 s-2)2 +16 s-2)2 = 6 L-1 ( s-2 4 2 + 2 L-1 ( 4 + 42 s-2)2 + s-2)2 จากตารางท่ี 9.1 จะได L-1 a = sin at และ L-1 s = cos at s2 +a2 s2 +a2 และใชค ณุ สมบัตกิ ารเล่ือนทางความถจี่ ะได L-1[F(s)] = f (t)eat เมื่อ a = 2 จะได = 6 L-1 ( s-2 4 2 + 2 L-1 ( 4 +4 2 s-2)2 + s-2)2 = 6e2tcos 4t + 2e2tsin 4t = 2e2t(3cos 4t + sin 4t) L-1 6s-4 = 2e2t(3cos 4t + sin 4t) s2 -4 s+ 20 ตัวอยางท่ี 9.8 L-1[F(s)] = L-1 s2 5s+30 +10s+25 วธิ ีทาํ ขนั้ ตอนที่ 1 หาสมการกาํ ลังสองสมบรู ณจาก s2+10s+25 = s2+2(5)s+52= (s+5)2 จุไรรัตนจ นิ ดา อรรคนติ ย
338 บทที่ 9 ลาปลาซทรานสฟอรม กลบั คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส ขัน้ ตอนท่ี 2 แทนคา ในสมการจะได L-1 s2 5s+30 = L-1 5s+30 +10s+25 (s+5)2 = L-1 5( s+5)+5 ( s+5)2 ขัน้ ตอนท่ี 3 ทําการแยกเศษสวนจะได L-1 5( s+5)+5 = L-1 5(s+5) + 5 ( s+5)2 (s+5)2 (s+5)2 = L-1 (s+55) + ( 5 s+5)2 = 5e-5t- 5te-5t = 5e-5t(1- t) ตอบ L-1 5s+30 = 5e-5t(1- t) s2 +10 s+25 ตัวอยา งท่ี 9.9 จงหาคา ลาปลาซทรานสฟอรม กลับ L-1[F(s)] = L-1 4s+3 s2 -2 s-5 วธิ ีทํา ขัน้ ตอนที่ 1 ใชส มการกําลังสองสมบรู ณ จาก s2- 2s - 5 = s2-2(1)s+12 = (s - 1)2- 6 ขน้ั ตอนท่ี 2 ทําการหาคา การแปลงกลบั จะได L-1[F(s)]= L-1 4s+3 s2 -2 s-5 = L-1 ( 4s+3 s-1)2 -6 ขน้ั ตอนท่ี 3 จัดสมการตวั เศษใหเ หมอื นกับตวั สว น จะได = L-1 4(s-1)+7 (s-1)2 - 6 ข้นั ตอนที่ 4 ทาํ การแยกพจนหาร 2 พจน จะได = L-1 ( 4(s-1) + 7 -6 s-1)2 -6 (s-1)2 จไุ รรตั นจ ินดา อรรคนติ ย
คณติ ศาสตรว ิศวกรรมอิเลก็ ทรอนกิ ส บทที่ 9 ลาปลาซทรานสฟ อรม กลบั 339 = L-1 ( 4(s-1) +7 s-1)2- 62 ( s-1)2 - 62 = 4 L-1 ( (s-1) 62 +7 s-1)2 - ( s-1)2 - 62 ขนั้ ตอนท่ี 5 จากตารางท่ี 9.10 จะได ตอบ L-1[F(s)]= L-1 s = cosh at s2-a2 L-1[F(s)]= L-1 a = sinh at s2-a2 จดั รปู สมการ ใชค ณุ สมบัตกิ ารเล่ือนทางความถ่ีจะได L-1[F(s - a)] = f (t)eat เมอ่ื a = 1 จะได = 4etcosh ( 6t ) + 7etsinh ( 6t ) = et[4cosh( 6t ) + 7sinh ( 6t )] L-1 4s+3 = et[4cosh( 6t ) + 7sinh( 6t )] s2 -2 s-5 9.3 การแยกเศษสว นยอ ย (Partial – fraction expansion) ในการหาคาลาปลาซทรานสฟอรมกลับสามารถทําไดโดยใชตารางที่ 9.1 และตารางท่ี 9.2 แตหากคาฟงกชัน F(s) ไมอยูในรูปแบบท่ีเทียบกลับไดทันที จะตองทําเปล่ียนเปลี่ยนรูป F(s) ใหสามารถอยใู นรปู แบบทสี่ ามารถใชตารางได ซ่ึงวิธีที่นํามาใช คือ การแยกเศษสวนยอย พิจารณา ฟงกชัน F(s) ในรูป F(s) = N (s) = bmsm +bm-1sm-1+...+b1s+b0 (9.9) D(s) ansn +an-1sn-1+...+a1s+a0 เมื่อ N(s) คือ ตัวตั้ง (Numerator) และ D(s) คือ ตัวหาร (Denominator) ซ่ึงอยูในรูปแบบของ สมการแบบโพลิโนเมียลของฟงกช นั s โดยใชตัวหาร คอื D(s) นาํ มาแยกเปนเศษสว นยอย ซ่งึ n > m โดยพิจารณา D(s) ไดเ ปน 2 กรณี คือ 9.3.1 เมื่อ D(s) มีรากปกติ (Simple roots) เมื่อ D(s) สามารถแยกตัวประกอบเปนหลายวงเล็บ ซ่งึ อยูใ นรปู ของ s ยกกําลังหนึง่ จไุ รรตั นจ ินดา อรรคนติ ย
340 บทที่ 9 ลาปลาซทรานสฟอรม กลบั คณิตศาสตรว ิศวกรรมอเิ ลก็ ทรอนกิ ส F(s) = N (s) = N (s) (9.10) D(s) (s - a1)(s - a2 )...(s - an ) เมื่อ a1 ,a2 ,….,an เปนรากของ D(s) และไมมีรากตัวใดที่ซํ้ากัน คือ a1 a2 …. an ซึ่งสามารถแยกเปนเศษสว นยอยไดดังน้ี s A-1a1 s A- 2a2 F(s) = N (s) = + + ... + s A- nan (9.11) D(s) ซึ่ง A1 , A2 ,……,An เปนคาคงท่ี สามารถหาคาไดโดย 2 วิธี คือ การเทียบสัมประสิทธ์ิ และการใชส ตู ร ซง่ึ มีสูตรการหาคา An ดงั น้ี An = F (s)(s - an ) s=an หรอื An = sA-a11 + sA-a22 + ... + sA-ann N (s) An = dds [D(s)] s = an , n = 1,2,3…. An = N(s) [D(s)] s = an ตวั อยา งที่ 9.10 จงหาลาปลาซทรานสฟอรมกลับของ L-1 3s-5 s2 -2 s-3 3s-5 วิธีทํา ขัน้ ตอนท่ี 1 เมอื่ F(s) = s2 -2 s-3 ทาํ การแยกตัวประกอบของสวน คอื s2-2s-3 = (s-3)(s+1) มีรากของสมการ 2 ตัว คอื (s- 3) = 0 จะได a1 = 3 (s+1) = 0 จะได a2 = -1 ข้ันตอนที่ 2 ทําการแยกตวั แปรจะได = (s-33)s(-s5+1) 3s-5 s2 -2 s-3 3s-5 sA-3 + sB+1 (เมอื่ A และ B เปนคา คงท)ี่ s2 -2 s-3 = (9.12) จุไรรตั นจ ินดา อรรคนติ ย
คณติ ศาสตรว ศิ วกรรมอิเล็กทรอนกิ ส บทที่ 9 ลาปลาซทรานสฟ อรม กลบั 341 ขน้ั ตอนท่ี 3 3.1 ทาํ การแกส มการโดยใชวิธีเทยี บสมั ประสทิ ธ์ิ จะได นาํ (s2- 2s - 3) หรือ (s-3)(s+1) คณู สมการท่ี (9.12) ทงั้ ซายมือและขวามอื ดงั น้ัน (เมอ่ื A และ B เปนคา คงที่) 3s-5 s2 -2 s-3 s2-2s-3 = sA-3 (s - 3)(s + 1) + sB+1 (s - 3)(s + 1) จะได 3s - 5 = A(s + 1) + B(s - 3) 3s - 5 = As + A + Bs - 3B 3s - 5 = (A + B)s + (A - 3B) A + B =3 จะได A = 3 - B (9.13) A - 3B = -5 (9.14) แทนคา A ในสมการที่ (9.14) จะได 3 - B - 3B = -5 3 - 4B = - 5 -4B = -8 ดงั น้นั B = 2 และ A =1 3s-5 ดังนนั้ s2 -2 s-3 = s1-3 - s+21 L-1 3s-5 = L-1 s1-3 - 2L-1 s+11 s2 -2 s-3 แทนคาในตารางที่ 9.1 จะได L-1 ( 3s-5 = e3t - 2e-t s-3)(s+1) 3.2 ทําการหาคา โดยการใชสตู รจะได A = F(s)(s - 3) s=3 A = (s-33)s(-s5+1) (s - 3) s=3 A = (3ss+-15) s=3 3(3)-5 A = 3+1 =1 B = F (s)(s + 1) s=-1 B = (s-33)s(-s5+1) (s + 1) s=-1 B = (3ss--35) s=-1 จไุ รรตั นจ ินดา อรรคนติ ย
Search
Read the Text Version
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
- 42
- 43
- 44
- 45
- 46
- 47
- 48
- 49
- 50
- 51
- 52
- 53
- 54
- 55
- 56
- 57
- 58
- 59
- 60
- 61
- 62
- 63
- 64
- 65
- 66
- 67
- 68
- 69
- 70
- 71
- 72
- 73
- 74
- 75
- 76
- 77
- 78
- 79
- 80
- 81
- 82
- 83
- 84
- 85
- 86
- 87
- 88
- 89
- 90
- 91
- 92
- 93
- 94
- 95
- 96
- 97
- 98
- 99
- 100
- 101
- 102
- 103
- 104
- 105
- 106
- 107
- 108
- 109
- 110
- 111
- 112
- 113
- 114
- 115
- 116
- 117
- 118
- 119
- 120
- 121
- 122
- 123
- 124
- 125
- 126
- 127
- 128
- 129
- 130
- 131
- 132
- 133
- 134
- 135
- 136
- 137
- 138
- 139
- 140
- 141
- 142
- 143
- 144
- 145
- 146
- 147
- 148
- 149
- 150
- 151
- 152
- 153
- 154
- 155
- 156
- 157
- 158
- 159
- 160
- 161
- 162
- 163
- 164
- 165
- 166
- 167
- 168
- 169
- 170
- 171
- 172
- 173
- 174
- 175
- 176
- 177
- 178
- 179
- 180
- 181
- 182
- 183
- 184
- 185
- 186
- 187
- 188
- 189
- 190
- 191
- 192
- 193
- 194
- 195
- 196
- 197
- 198
- 199
- 200
- 201
- 202
- 203
- 204
- 205
- 206
- 207
- 208
- 209
- 210
- 211
- 212
- 213
- 214
- 215
- 216
- 217
- 218
- 219
- 220
- 221
- 222
- 223
- 224
- 225
- 226
- 227
- 228
- 229
- 230
- 231
- 232
- 233
- 234
- 235
- 236
- 237
- 238
- 239
- 240
- 241
- 242
- 243
- 244
- 245
- 246
- 247
- 248
- 249
- 250
- 251
- 252
- 253
- 254
- 255
- 256
- 257
- 258
- 259
- 260
- 261
- 262
- 263
- 264
- 265
- 266
- 267
- 268
- 269
- 270
- 271
- 272
- 273
- 274
- 275
- 276
- 277
- 278
- 279
- 280
- 281
- 282
- 283
- 284
- 285
- 286
- 287
- 288
- 289
- 290
- 291
- 292
- 293
- 294
- 295
- 296
- 297
- 298
- 299
- 300
- 301
- 302
- 303
- 304
- 305
- 306
- 307
- 308
- 309
- 310
- 311
- 312
- 313
- 314
- 315
- 316
- 317
- 318
- 319
- 320
- 321
- 322
- 323
- 324
- 325
- 326
- 327
- 328
- 329
- 330
- 331
- 332
- 333
- 334
- 335
- 336
- 337
- 338
- 339
- 340
- 341
- 342
- 343
- 344
- 345
- 346
- 347
- 348
- 349
- 350
- 351
- 352
- 353
- 354
- 355
- 356
- 357
- 358
- 359
- 360
- 361
- 362
- 363
- 364
- 365
- 366
- 367
- 368
- 369
- 370
- 371
- 372
- 373
- 374
- 375
- 376
- 377
- 378
- 379
- 380
- 381
- 382
- 383
- 384
- 385
- 386
- 387
- 388
- 389
- 390
- 391
- 392
- 393
- 394
- 395
- 396
- 397
- 398
- 399
- 400
- 401
- 402
- 403
- 404
- 405
- 406
- 407
- 408
- 409
- 410
- 411
- 412
- 413
- 414
- 415
- 416
- 417
- 418
- 419
- 420
- 421
- 422
- 423
- 424
- 425
- 426
- 427
- 428
- 429
- 430
- 431
- 432
- 433
- 434
- 435
- 436
- 437
- 438
- 439
- 440
- 441