AB2.8: Laplacetransformation av derivator och integraler. Differentialekvationer Laplacetransformen som an analytisk funktion SATS 1 Om Laplaceintegralen F (s) = L (f ) = Z ∞ e−st f (t)dt 0 ¨r konvergent f¨or s = s0 , s˚ a a konvergerar den f¨or alla s i halvplanet Re s >Re s0 (och f¨or alla reela s > s0 ) och F (s) a¨r i detta halvplan en regulj¨ ar analytisk funktion, dvs en funktion som har derivator av godtyckligt h¨og ordning. Derivatorna f˚ as genom derivering under integraltecknet och allm¨ ant Z ∞ d e−st (−t)f (t)dt F (s) = ds 0 (1) Z ∞ dn F (s) = e−st (−t)n f (t)dt n ds 0 (2) Vi kan multipliciera b˚ ada leden i (2) med (−1)n och f˚ ar d˚ a L (tn f (t)) = (−1)n dn F (s) (n = 1, 2, . . .). dsn (3) F¨or den speciella funktionen f (t) = 1 ¨ar F (s) = L (f ) = Z ∞ e−st dt = s−1 0 vilket ger L (tn ) = (−1)n dn −1 s = (−1)n (−1)(−2) . . . (−n)s−n−1 = n!s−n−1 dsn och allts˚ a L (tn ) = n! sn+1 (n = 0, 1, 2, . . .). (n = 0, 1, 2, . . .) (4) (5) T ex 1 2 6 , L (t2 ) = 3 , L (t3 ) = 4 . (6) 2 s s s En funktion f (t) ¨ar av exponentiell ordning k om det existerar positiva konstanter k och M sadana att f¨or n˚ agot T olikheten L (t) = |f (t)| ≤ M ekt , t≥T (7) g¨aller. Exponentialfunktionen ect ¨ar en funktion av exponentiell ordning k = c (eftersom ect ≤ ct e , t ≥ 0). sin t och cos t ¨ar funktioner av exponentiell ordning k = 0 (eftersom t ex | sin t| ≤ e0·t = 1). SATS 2 L˚ at f (t) vara en styckvis kontinuerlig funktion i intervallet t ≥ 0 (dvs med m¨ojligt undantag f¨or isolerade punkter.) Antag att f (t) ¨ ar en funktion av exponentiell ordning k, s˚ a existerar dess Laplacetransformen for alla komplexa s sadana att Re s > k (och f¨ or alla reela s > k). Bevis.R Betrakta reela s. f (t) har en Laplacetransform eftersom enligt (7) motsvarande integralen 0∞ e−st f (t)dt ¨ar konvergent: Z ∞ 0 |e −st f (t)|dt ≤ Z ∞ 0 −st e |f (t)|dt ≤ Z ∞ e −st kt M e dt = M 0 Z ∞ e(k−s)t dt = 0 M s−k (s > k). SATS 5.2.1 L˚ at f (t) vara en styckvis kontinuerlig funktion i intervallet t ≥ 0 och en funktion av exponentiell ordning k (dvs satisfierar villkoret |f (t)| ≤ M ekt (t ≥ T ), (8) d¨ ar k och M ¨ar vissa positiva konstanter,) och derivatan f 0 (t) en styckvis kontinuerlig funktion i intervallet t ≥ 0. D˚ a existerar Laplacetransformen av derivatan f 0 (t) om s > k och L (f 0 ) = sL (f ) − f (0), s > k. (9) Bevis. Antag att derivatan f 0 (t) ¨ar en kontinuerlig funktion i intervallet t ≥ 0. D˚ a , enligt definitionen av Laplacetransform, kan man v¨alja s s˚ a stor att integralen Z ∞ e−st f 0 (t)dt 0 ¨ar konvergent och lim e−st f (t) = 0. t→∞ L¨agg m¨arke till att f (t) har en Laplacetransform (se SATS 2). Genom att anv¨anda definitionen av Laplacetransform och partiell integration, f˚ ar vi 0 L (f ) = Z ∞ 0 e −st 0 −st f (t)dt == e ∞ f (t) 0 +s Z ∞ 0 e−st f (t)dt = 0 − f (0) + sL (f ), s > k. Best¨am Laplacetransformen av derivatan av andra ordningen L (f 00 ) = sL (f 0 ) − f 0 (0) = s[sL (f ) − f (0)] − f 0 (0) = s2 L (f ) − sf (0) − f 0 (0). (10) Vidare, L (f 000 ) = s3 L (f ) − s2 f (0) − sf 0 (0) − f 00 (0). Genom induktion (och upprepad anv¨andning av (9)), f˚ ar vi Laplacetransformer av deriva(n) torna f av en godtycklig ordning n = 1, 2, . . . L (f (n) ) = sn L (f ) − sn−1 f (0) − sn−2 f 0 (0) − . . . − sf (n−2) (0) − f (n−1) (0). (11) EXEMPEL 1 Laplacetransformen av f (t) = t2 L˚ at f (t) = t2 , t ≥ 0. Best¨am F (s). L¨ osning. Vi har f 0 (t) = 2t, f 00 (t) = 2, d˚ a f (0) = 0, f 0 (0) = 0, f 00 (0) = 2. och Laplacetransformen ¨ar 2 = s2 L (f ), s Betrakta den sista likheten i (12)) som en ekvation och best¨am Laplacetransformen L L (f 00 ) = L (2) = L (t2 ) = 2 . s3 EXEMPEL 2 Laplacetransformen av cos ωt Best¨am Laplacetransformen av f (t) = cos ωt. L¨ osning. Vi har f 00 (t) = −ω 2 cos ωt = −ω 2 f (t), och f (0) = 1, f 0 (0) = 0. (10) ger L (f 00 ) = −ω 2 L (f ) = s2 L (f ) − sf (0) − f 0 (0) = s2 L (f ) − s. Betrakta den sista likheten som en ekvation och best¨am Laplacetransformen L L (f ) = s2 s . + ω2 EXEMPEL 2’ Laplacetransformen av sin ωt Best¨am Laplacetransformen av f (t) = sin ωt. L¨ osning. Vi har f 0 (t) = ω cos ωt, f 00 (t) = −ω 2 sin ωt = −ω 2 f (t). Ber¨akna derivator i punkten t = 0 f (0) = 0, f 0 (0) = ω. D˚ a ger (10) L (f 00 ) = −ω 2 L (f ) = s2 L (f ) − sf (0) − f 0 (0) = s2 L (f ) − ω. Betrakta den sista likheten som en ekvation och best¨am Laplacetransformen L L (f ) = s2 ω . + ω2 (12) EXEMPEL 3 Best¨am Laplacetransformen av f (t) = sin2 t. L¨ osning. Vi har f 0 (t) = 2 sin t cos t = sin 2t, f 00 (t) = 2 cos 2t. Ber¨akna derivator i punkten t = 0 f 0 (0) = 0, f (0) = 0, f 00 (0) = 2. D˚ a ger (10) 2s = s2 L (f ) − sf (0) − f 0 (0) = s2 L (f ). 2 s +4 L (f 00 ) = 2L (cos 2t) = Betrakta likheten 2s = s2 L (f ) +4 som en ekvation och best¨am Laplacetransformen L s2 L (sin2 t) = 2 . + 4) s(s2 EXEMPEL 4 Best¨am Laplacetransformen av f (t) = t sin ωt. L¨ osning. Vi har f (0) = 0 och f 0 (t) = sin ωt + ωt cos ωt, f 0 (0) = 0. Ber¨akna derivatan av andra ordningen f 00 (t) = 2ω cos ωt − ω 2 t sin ωt = 2ω cos ωt − ω 2 f (t), Enligt (10) L (f 00 ) = 2ωL (cos ωt) − ω 2 L (f ) = s2 L (f ) − sf (0) − f 0 (0) = s2 L (f ). Betrakta likheten s − ω 2 L (f ) = s2 L (f ) 2 +ω som en ekvation och best¨am Laplacetransformen L 2ω s2 L (f ) = L (t sin ωt) = (s2 2ωs . + ω 2 )2 Differentialekvationer. Begynnelsev¨ ardesproblem Grundl¨ aggande begrepp Den enklaste formen f¨or en differentialekvation av f¨orsta ordningen ¨ar y 0 = h(x). (13) En s˚ adan ekvation kan l¨osas direkt. Om H(x) = Z h(x)dx [H 0 (x) = h(x)] ¨ar en primitiv funktion till h(x) s˚ a ¨ar ju y(x) = H(x) + C den alm¨anna l¨osningen till (13). Konstanten C best¨ams av n˚ agot begynnelsevillkor. En linj¨ar differentialekvation av f¨orsta ordningen ¨ar D(y) ≡ y 0 + g(x)y = h(x). (14) H¨ar ¨ar g och h givna kontinuerliga funktioner i ett ¨oppet intervall p˚ a reela axeln x. D kallas en linj¨ ar differentialoperator (av f¨orsta ordningen) eftersom D(y1 + y2 ) = (y1 + y2 )0 + g(x)(y1 + y2 ) = y10 + y20 + g(x)y1 + g(x)y2 = D(y1 ) + D(y2 ), D(αy) = (αy)0 + g(x)(αy) = αy 0 + αg(x)y = αD(y) och D(αy1 + βy2 ) = αL(y1 ) + βL(y2 ). (15) Om y1 och y2 l¨oser de tv˚ a ekvationerna D(y) = h1 (x) respektive D(y) = h2 (x) s˚ a l¨oser y1 + y2 ekvationen D(y) = h1 (x) + h2 (x) och αy1 l¨oser ekvationen D(y) = αh1 (x). Detta kallas superpositionprincipen. En differentialekvation av f¨orsta ordningen s¨ages vara separabel eller ha separabla variabler om den kan skrivas p˚ a formen g(y)y 0 = h(x). En s˚ adan ekvation kan l¨osas direkt. Betrakta t ex begynnelsev¨ ardesproblemet ( √ dy = −k y, y = y(t), t > 0, dt y(0) = h, (16) Differentialekvationen kan skrivas som 1 dy = −k √ y dt (y > 0). (17) Den har separabla variabler. Vi har d √ 1 (2 y) = √ : dy y d q 1 dy (2 y(t)) = √ , dt y dt (18) vilket betyder att differentialekvationen kan skrivas som d √ (2 y) = −k dt Detta ¨ar ekvivalent med √ 2 y= (19) Z (−k)dt = −kt + C. √ Begynnelsev¨ardet y(0) = h ger att C = 2 h och allts˚ a √ y= √ −k t+ h 2 (20) I detta fall kan vi enkla l¨osa y genom kvadrering: y= !2 √ −k t+ h 2 . (21) En linj¨ar differentialekvation av andra ordningen ¨ar M (y) ≡ y 00 + a(x)y 0 + b(x)y = h(x). (22) H¨ar ¨ar a, b och h givna kontinuerliga funktioner. M kallas en linj¨ar differentialoperator (av andra ordningen) eftersom M satisfierar (15). Ekvationen y 00 + a(x)y 0 + b(x)y = 0 (23) kallas den till (22) h¨orande homogena ekvationen. (22) kallas inhomogena ekvationen L˚ at yp vara en given l¨osning till (22). D˚ a ¨ar funktionen y l¨osning till (22) om och endast om y ¨ar av formen y = yh + yp , d¨ar funktionen yh ¨ar en l¨osning till motsvarande homogena ekvationen (23). L¨osningen yp kallas en partikul¨ ar l¨osning. Betrakta ekvationen y 00 + ay 0 + by = 0 (24) med konstanta (komplexa) koefficienterna a och b. L¨osningen till homogena ekvationen (24) ¨ar av formen y = C1 er1 x + C2 er2 x , r1 6= r2 , (25) eller y = (C1 + C2 x)erx , r1 = r2 = r, (26) d¨ar r1 och r2 ¨ar nollst¨allena till motsvarande karakteristiska polynomet r2 + ar + b. (27) y 00 − 4y 0 + 3y = 0 (28) r2 − 4r + 3 (29) Exempel Ekvationen har karakteristiska polynomet med nollst¨allena r1 = 1 och r2 = 3. Den fullst¨andiga l¨osningen till homogena ekvationen (28) ¨ar y = C1 ex + C2 e3x . (30) y 00 + y = 0 (31) r2 + 1 (32) Ekvationen har karakteristiska polynomet med komplexa nollst¨allena r1 = i och r2 = −i. Den fullst¨andiga l¨osningen till (31) ¨ar y = C1 eix + C2 e−ix = C1 cos x + iC1 sin x + C2 cos x − iC2 sin x = C˜1 cos x + C˜2 sin x. (33) Randv¨ardesproblemet f¨or linj¨ara differentialekvationen av andra ordningen y 00 + ay 0 + by = r(t) skrivas som y 00 + ay 0 + by = r(t), y(0) = K0 , y 0 (0) = K1 (a, b = const). (34) d¨ar a, b, K0 och K1 ¨ar givna tal och r(t) ¨ar en given kontinuerlig funktion (indata). L¨os (34) med h¨ajlp av Laplaces metod. Steg 1. S¨att Y = L (y), R = L (r) Laplacetransformationen av b˚ ada leden i (34) ger L (y 00 ) + aL (y 0 ) + bL (y) = R, som skrivas [s2 Y − sy(0) − y 0 (0)] + a[sY − y(0)] + bY = R. (35) (35) kallas sidoekvationen (subsidiary equation) och reduceras till (s2 + as + b)Y = (s + a)y(0) + y 0 (0) + R(s). Steg 2. Definiera transferfunktionen Q = Q(s) = s2 1 + as + b (36) och l¨os sidoekvationen Y (s) = [(s + a)y(0) + y 0 (0)]Q(s) + R(s)Q(s). (37) Om y(0) = y 0 (0) = 0, d˚ a f˚ ar vi Y = RQ och Q ¨ar kvot av tv˚ a polynom (en rationell funktion) Q= Y L(output) = R L(input) Steg 3. Best¨am y(t) = L−1 (Y ) med h¨ajlp av inverstransformationen, partialbr˚ aksuppdelning i (37) och handp˚ al¨aggningsmetoden. EXEMPEL 5 L¨os begynnelsev¨ardesproblemet y 00 − y = t, y(0) = 1, y 0 (0) = 1. L¨ osning. Enligt (34) ¨ar indata a = 0, b = −1, r(t) = t, K0 = K1 = 1. Vi l¨oser begynnelsev¨ardesproblemet (38) med h¨ajlp av Laplaces metod. Steg 1. S¨att 1 Y = L (y), R = L (r) = L (t) = 2 . s 00 Laplacetransformationen av b˚ ada leden i ekvationen y − y = t ger L (y 00 ) − L (y) = R; den sidoekvationen [s2 Y − sy(0) − y 0 (0)] − Y = eller (s2 − 1)Y = s + 1 + 1 . s2 Steg 2. Definiera transferfunktionen Q = Q(s) = s2 1 −1 1 , s2 (38) och l¨os sidoekvationen 1 s+1 1 1 1 1 Q(s) = 2 + 2 2 = + 2 − 2. 2 s s − 1 s (s − 1) s−1 s −1 s Y (s) = [s + 1]Q(s) + Steg 3. Best¨am l¨osningen till begynnelsev¨ardesproblemet (38) −1 y(t) = L −1 (Y ) = L 1 1 1 + L−1 2 − L−1 2 = et + sinh t − t. s−1 s −1 s Laplacetransformation av integral SATS 5.2.2 L˚ at F (s) vara Laplacetransformen av f (t). Om f (t) ¨ ar en styckvis kontinuerlig funktion i intervallet t ≥ 0 och satisfierar villkoret |f (t)| ≤ M ekt (t ≥ 0), d¨ ar k och M ¨ar positiva tal, d˚ a L Z t f (τ )dτ 0 eller Z t 1 = F (s) (s > 0, s > k), s 1 F (s) . s 0 Bevis. Om f (t) ¨ar en styckvis kontinuerlig funktion s˚ adan att f (τ )dτ = L−1 |f (t)| ≤ M ekt (t ≥ 0), (39) d¨ar k > 0 och M > 0, d˚ a ¨ar integralen g(t) = Z t f (τ )dτ 0 en kontinuerlig funktion som satisfierar |g(t)| ≤ Z t 0 |f (τ )|dτ ≤ M Z t ekτ dτ = 0 M kt M kt (e − 1) ≤ e k k (k > 0), dvs (37) g¨aller f¨or g(t). Vi har g 0 (t) = f (t); d˚ a ¨ar g 0 (t) en styckvis kontinuerlig funktion i intervallet t ≥ 0. Vi f˚ ar Laplacetransformen L [f (t)] = L [g 0 (t)] = sL (g) − g(0) (s > k) och F (s) = L (f ) = sL (g) eftersom i (40), g(0) = 0. Laplacetransformen av integralen kan d¨arf¨or skrivas L (g) = L Z t 0 f (τ )dτ 1 = F (s). s (40) EXEMPEL 7 L˚ at L (f ) = s(s2 1 . + ω2) Best¨am f (t). L¨ osning. Inverstransformationen ger L och L −1 " −1 1 1 = sin ωt 2 2 s +ω ω # 1 1Zt 1 Zt = sin ωτ dτ = (1 − cos ωt). s(s2 + ω 2 ) ω 0 ω2 0 PROBLEM 5.2.1 L¨os begynnelsev¨ardesproblemet y 0 + 3y = 10 sin t, y(0) = 0. (41) L¨ osning. Vi ska l¨osa begynnelsev¨ardesproblemet (41) med h¨ajlp av Laplaces metod. Steg 1. Skriv 10 Y = L (y), R = L (10 sin t) = 2 s +1 0 och Laplacetransformera y + 3y = 10 sin t L (y 0 ) + 3L (y) = 10 s2 + 1 Skriv motsvarande sidoekvationen [sY − y(0)] + 3Y = 10 . +1 s2 Begynnelsevillkoret y(0) = 0 i (41) ger (s + 3)Y = 10 . +1 s2 Steg 2. L¨os sidoekvationen med h¨ajlp av transferfunktionen Q = Q(s) = 1 . s+3 Vi f˚ ar 10 10 1 A Bs + C Q(s) = 2 = + 2 . +1 s +1s+3 s+3 s +1 Anv¨and handp˚ al¨aggningsmetoden f¨or att best¨amma talen A, B och C Y (s) = s2 10 = A(s2 + 1) + (Bs + C)s + 3, A = 1, B = −1, C = 3. Steg 3. Best¨am l¨osningen till begynnelsev¨ardesproblemet (41) y(t) = L−1 (Y ) = L−1 s 1 1 − L−1 2 + 3L−1 2 = e−3t − cos t + 3 sin t. s+3 s +1 s +1 PROBLEM 5.2.5 L¨os begynnelsev¨ardesproblemet y 00 + ay 0 − 2a2 y = 0, y(0) = 6, y 0 (0) = 0 (a = const). (42) L¨ osning. Vi ska l¨osa begynnelsev¨ardesproblemet (42) med h¨ajlp av Laplaces metod. I (42), ¨ar indata b = −2a2 , r(t) = 0, K0 = 6, K1 = 0. Steg 1. S¨att Y = L (y), R = L (r) = 0 och Laplacetransformera y 00 + ay 0 − 2a2 y = 0 L (y 00 ) + aL (y 0 ) − 2a2 L (y) = 0. Skriv motsvarande sidoekvationen [s2 Y − sy(0) − y 0 (0)] + a[sY − y(0)] − 2a2 Y = 0. Begynnelsevillkoren y(0) = 6 och y 0 (0) = 0 i (42) ger (s2 + as − 2a2 )Y = 6(s + a). Steg 2. L¨os sidoekvationen med h¨ajlp av transferfunktionen Q = Q(s) = 1 . s2 + as − 2a2 Vi f˚ ar s+a + as − 2a2 Steg 3. Best¨am l¨osningen till begynnelsev¨ardesproblemet (41) Y (s) = 6(s + a)Q(s) = 6 s2 s a +6 2 = 2 + as − 2a s + as − 2a2 s a 6 +6 = 2 2 2 2 (s + a/2) − (a/2) − 2a (s + a/2) − (a/2)2 − 2a2 s a 6 +6 = 2 2 2 (s + a/2) − (3a/2) (s + a/2) − (3a/2)2 y(t) = L−1 (Y ) = 6 s2 6e−(a/2)t cosh (3a/2)t + 4e−(a/2)t sinh (3a/2)t = 2e−2at + 4eat . PROBLEM 5.2.11 Best¨am Laplacetransformen F (s) = L (cos2 t). L¨ osning. S¨att f = cos2 t = 1 − sin2 t. Vidare, f (0) = 1, f 0 = −2 cos t sin t = −2 sin 2t. Laplacetransformera − sin 2t L (− sin 2t) = − s2 2 = sL (f ) − 1. +4 (43) Betrakta den sista likheten i (43)) som en ekvation och best¨am Laplacetransformen L av cos2 t 1 2 1 L (f ) = 1− 2 = s s +4 s ! s2 + 2 . s2 + 4 PROBLEM 5.2.12a Best¨am Laplacetransformen av f (t) = t cos ωt. L¨ osning. Eftersom f (0) = 0, f˚ ar vi f 0 (t) = cos ωt − ωt sin ωt, f 0 (0) = 1; f 00 (t) = −2ω sin ωt − ω 2 t cos ωt = −2ω sin ωt − ω 2 f (t), och L (f 00 ) = −2ωL (sin ωt) − ω 2 L (f ) = s2 L (f ) − sf (0) − f 0 (0) = s2 L (f ) − 1. Betrakta likheten 2ω 2 − ω 2 L (f ) = s2 L (f ) − 1 s2 + ω 2 som en ekvation och best¨am Laplacetransformen L av f (t) − L (f ) = L (t cos ωt) = s2 − ω 2 . (s2 + ω 2 )2 PROBLEM 5.2.12b Visa att −1 L ! 1 2 (s + ω 2 )2 = 1 (sin ωt − ωt cos ωt). 2ω 3 L¨ osning. Vi har ω(s2 − ω 2 ) , (s2 + ω 2 )2 ω L (sin ωt) = 2 , s + ω2 −L (ωt cos ωt) = − och ω(s2 − ω 2 ) 2ω 3 ω − = . s2 + ω 2 (s2 + ω 2 )2 (s2 + ω 2 )2 L (sin ωt − ωt cos ωt) = PROBLEM 5.2.13 Best¨am originalfunktionen f (t) till Laplacetransformen F (s) = L (f ) = s2 1 . + 4s L¨ osning. Vi har 1 1 1 = L (e−4t ). ss+4 s F (s) = Laplacetransformen av integral ¨ar L Z t e 0 och f (t) = −4τ Z t 0 dτ 1 = L (e−4t ), s 1 e−4τ dτ = (1 − e−4t ). 4
© Copyright 2024