Si të aplikoni transformimin Laplace në një funksion

Autor: Ellen Moore
Data E Krijimit: 19 Janar 2021
Datën E Azhurnimit: 2 Korrik 2024
Anonim
Solving IVPs with the Laplace Transform and farewell to Freeman Dyson
Video: Solving IVPs with the Laplace Transform and farewell to Freeman Dyson

Përmbajtje

Transformimi Laplace është një transformim integral që përdoret për të zgjidhur ekuacionet diferenciale me koeficientë konstantë. Ky transformim është përdorur gjerësisht në fizikë dhe inxhinieri.

Ndërsa mund të përdorni tabelat e duhura, është e dobishme të kuptoni transformimin e Laplace në mënyrë që ta bëni vetë nëse është e nevojshme.

Informacioni paraprak

  • I është dhënë një funksion f(t){ displaystyle f (t)}përcaktuar për t0.{ displaystyle t geq 0.} Atëherë Transformimi i Laplace funksionin f(t){ displaystyle f (t)} është funksioni tjetër i secilës vlerë s{ displaystyle s}, në të cilën integrali konvergon:
    • F(s)=L{f(t)}=0f(t)estdt{ displaystyle F (s) = { mathcal {L}} {f (t) } = int _ {0} ^ { infty} f (t) e ^ {- st} mathrm {d} t}
  • Transformimi i Laplace merr një funksion nga rajoni t (shkalla kohore) në rajonin s (rajoni i transformimit), ku F(s){ stilet e ekranit F (të)} është një funksion kompleks i një ndryshoreje komplekse. Kjo ju lejon të zhvendosni funksionin në një zonë ku një zgjidhje mund të gjendet më lehtë.
  • Natyrisht, transformimi Laplace është një operator linear, kështu që nëse kemi të bëjmë me një shumë termash, secili integral mund të llogaritet veçmas.
    • 0[af(t)+bg(t)]estdt=a0f(t)estdt+b0g(t)estdt{ displaystyle int _ {0} ^ { infty} [af (t) + bg (t)] e ^ {- st} mathrm {d} t = a int _ {0} ^ { infty} f (t) e ^ {- st} mathrm {d} t + b int _ {0} ^ { infty} g (t) e ^ {- st} mathrm {d} t}
  • Mos harroni se transformimi Laplace funksionon vetëm nëse integrali konvergon. Nëse funksioni f(t){ displaystyle f (t)} ka ndërprerje, është e nevojshme të jeni të kujdesshëm dhe të vendosni saktë kufijtë e integrimit në mënyrë që të shmangni pasigurinë.

Hapa

Pjesa 1 nga 3: Bazat

  1. 1 Zëvendësoni funksionin në formulën e transformimit Laplace. Teorikisht, transformimi Laplace i një funksioni është shumë i lehtë për t'u llogaritur. Si shembull, merrni parasysh funksionin f(t)=eat{ displaystyle f (t) = e ^ {at}}, ku a{ displaystyle a} është një konstante komplekse me Re(s)Re(a).{ displaystyle operatorname {Re} (s) operatorname {Re} (a).}
    • L{eat}=0eatestdt{ displaystyle { mathcal {L}} {e ^ {at} } = int _ {0} ^ { infty} e ^ {at} e ^ {- st} mathrm {d} t}
  2. 2 Vlerësoni integralin duke përdorur metodat në dispozicion. Në shembullin tonë, vlerësimi është shumë i thjeshtë dhe mund të arrini me llogaritjet e thjeshta. Në raste më komplekse, mund të nevojiten metoda më komplekse, për shembull, integrimi sipas pjesëve ose diferencimi nën shenjën integrale. Kusht kufizues Re(s)Re(a){ displaystyle operatorname {Re} (s) operatorname {Re} (a)} do të thotë që integrali konvergon, domethënë vlera e tij tenton në 0 si t.{ displaystyle t to infty.}
    • L{eat}=0e(as)tdt=e(as)tas|0=1sa{ displaystyle { filloj {në linjë} { mathcal {L}} {e ^ {at} } & = int _ {0} ^ { infty} e ^ {(si) t} mathrm {d } t & = { frac {e ^ {(si) t}} {as}} Bigg _ {0} ^ { infty} & = { frac {1} {sa}} fund {lidhur}}}}
    • Vini re se kjo na jep dy lloje të transformimit Laplace, me sinus dhe kosinus, pasi sipas formulës së Euler euneat{ displaystyle e ^ {iat}}... Në këtë rast, në emërues marrim sunea,{ displaystyle s-ia,} dhe mbetet vetëm për të përcaktuar pjesët reale dhe imagjinare. Ju gjithashtu mund të vlerësoni rezultatin drejtpërdrejt, por kjo do të zgjasë pak më shumë.
      • L{cosat}=Re(1sunea)=ss2+a2{ displaystyle { mathcal {L}} { cos at } = emri i operës {Re} majtas ({ frac {1} {s-ia}} djathtas) = ​​{ frac {s} {s ^ {2} + a ^ {2}}}}
      • L{mëkatat}=Une jam(1sunea)=as2+a2{ displaystyle { mathcal {L}} { sin at } = emri i operës {Im} majtas ({ frac {1} {s-ia}} djathtas) = ​​{ frac {a} {s ^ {2} + a ^ {2}}}}
  3. 3 Konsideroni transformimin Laplace të një funksioni fuqie. Së pari, ju duhet të përcaktoni transformimin e funksionit të fuqisë, pasi vetia e linearitetit ju lejon të gjeni transformimin për nga te gjitha polinome. Një funksion i formës tn,{ displaystyle t ^ {n},} ku n{ displaystyle n} - çdo numër i plotë pozitiv. Mund të integrohet pjesë për pjesë për të përcaktuar një rregull rekursive.
    • L{tn}=0tnestdt=nsL{tn1}{ displaystyle { mathcal {L}} {t ^ {n} } = int _ {0} ^ { infty} t ^ {n} e ^ {- st} mathrm {d} t = { frac {n} {s}} { mathcal {L}} {t ^ {n-1} }}
    • Ky rezultat shprehet në mënyrë të nënkuptuar, por nëse zëvendësoni disa vlera n,{ displaystyle n,} mund të krijoni një model të caktuar (përpiquni ta bëni vetë), i cili ju lejon të merrni rezultatin e mëposhtëm:
      • L{tn}=n!sn+1{ displaystyle { mathcal {L}} {t ^ {n} } = { frac {n!} {s ^ {n + 1}}}}
    • Ju gjithashtu mund të përcaktoni transformimin Laplace të fuqive fraksionale duke përdorur funksionin gama. Për shembull, në këtë mënyrë mund të gjeni transformimin e një funksioni si p.sh f(t)=t.{ displaystyle f (t) = { sqrt {t}}.}
      • L{tn}=Γ(n+1)sn+1{ displaystyle { mathcal {L}} {t ^ {n} } = { frac { Gama (n + 1)} {s ^ {n + 1}}}}
      • L{t1/2}=Γ(3/2)s3/2=π2ss{ displaystyle { mathcal {L}} {t ^ {1/2} } = { frac { Gama (3/2)} {s ^ {3/2}}} = { frac { sqrt { pi}} {2s { sqrt {s}}}}}
    • Edhe pse funksionet me fuqi thyesore duhet të kenë shkurtime (mbani mend, çdo numër kompleks z{ displaystyle z} dhe α{ displaystyle alpha} mund të shkruhet si zα{ displaystyle z ^ { alpha}}, sepse eαRegjistrohuz{ displaystyle e ^ { alpha emri i operës {Log} z}}), ato gjithmonë mund të përcaktohen në mënyrë të tillë që prerjet të qëndrojnë në gjysmën e majtë, dhe kështu të shmangin problemet me analitizmin.

Pjesa 2 nga 3: Vetitë e transformimit të Laplasit

  1. 1 Le të gjejmë transformimin Laplace të funksionit të shumëzuar me eat{ displaystyle e ^ {at}}. Rezultatet e marra në pjesën e mëparshme na lejuan të zbulojmë disa veti interesante të transformimit të Laplace. Transformimi Laplace i funksioneve të tilla si kosinusi, sinusi dhe funksioni eksponencial duket të jetë më i thjeshtë sesa transformimi i funksionit të fuqisë. Shumëzimi me eat{ displaystyle e ^ {at}} në rajonin t korrespondon me zhvendosje në rajonin s:
    • L{eatf(t)}=0f(t)e(sa)tdt=F(sa){ displaystyle { mathcal {L}} {e ^ {at} f (t) } = int _ {0} ^ { infty} f (t) e ^ {- (sa) t} mathrm {d} t = F (sa)}
    • Kjo pronë ju lejon menjëherë të gjeni transformimin e funksioneve të tilla si f(t)=e3tmëkat2t{ displaystyle f (t) = e ^ {3t} sin 2t}, pa pasur nevojë të llogarisim integralin:
      • L{e3tmëkat2t}=2(s3)2+4{ displaystyle { mathcal {L}} {e ^ {3t} sin 2t } = { frac {2} {(s-3) ^ {2} +4}}}
  2. 2 Le të gjejmë transformimin Laplace të funksionit të shumëzuar me tn{ displaystyle t ^ {n}}. Së pari, merrni parasysh shumëzimin me t{ displaystyle t}... Sipas përkufizimit, mund të dalloni një funksion nën një integral dhe të merrni një rezultat çuditërisht të thjeshtë:
    • L{tf(t)}=0tf(t)estdt=0f(t)sestdt=dds0f(t)estdt=dFds{ displaystyle { filloj {në linjë} { mathcal {L}} {tf (t) } & = int _ {0} ^ { infty} tf (t) e ^ {- st} mathrm { d} t & = - int _ {0} ^ { infty} f (t) { frac { partial} { partial s}} e ^ { - st} mathrm {d} t & = - { frac { mathrm {d}} { mathrm {d} s}} int _ {0} ^ { infty} f (t) e ^ { - st} mathrm {d} t & = - { frac { mathrm {d} F} { mathrm {d} s}} në fund {në linjë}}}
    • Duke përsëritur këtë operacion, marrim rezultatin përfundimtar:
      • L{tnf(t)}=(1)ndnFdsn{ displaystyle { mathcal {L}} {t ^ {n} f (t) } = (- 1) ^ {n} { frac { mathrm {d} ^ {n} F} { mathrm {d} s ^ {n}}}}
    • Megjithëse rirregullimi i operatorëve të integrimit dhe diferencimit kërkon një justifikim shtesë, ne nuk do ta paraqesim këtu, por vetëm vëmë re se ky operacion është i saktë nëse rezultati përfundimtar ka kuptim. Ju gjithashtu mund të merrni parasysh faktin se variablat s{ displaystyle s} dhe t{ displaystyle t} nuk varen nga njeri -tjetri.
    • Duke përdorur këtë rregull, është e lehtë të gjesh transformimin e funksioneve të tilla si p.sh t2cos2t{ displaystyle t ^ {2} cos 2t}, pa ri-integrim sipas pjesëve:
      • L{t2cos2t}=d2ds2ss2+4=2s324s(s2+4)3{ displaystyle { mathcal {L}} {t ^ {2} cos 2t } = { frac { mathrm {d} ^ {2}} { mathrm {d} s ^ {2}}} { frac {s} {s ^ {2} +4}} = { frac {2s ^ {3} -24s} {(s ^ {2} +4) ^ {3}}}}
  3. 3 Gjeni transformimin Laplace të funksionit f(at){ displaystyle f (at)}. Kjo mund të bëhet lehtë duke zëvendësuar ndryshoren me u duke përdorur përkufizimin e transformimit:
    • L{f(at)}=0f(at)estdt,  u=at=1a0f(u)esu/adu=1aF(sa){ displaystyle { begin {përputhje} { mathcal {L}} {f (at) } & = int _ {0} ^ { infty} f (at) e ^ {- st} mathrm { d} t, u = at & = { frac {1} {a}} int _ {0} ^ { infty} f (u) e ^ {- su / a} mathrm {d } u & = { frac {1} {a}} F majtas ({ frac {s} {a}} djathtas) në fund {në linjë}}}
    • Më sipër, gjetëm transformimin e funksioneve Laplace mëkatat{ displaystyle sin at} dhe cosat{ displaystyle cos at} drejtpërdrejt nga funksioni eksponencial. Duke përdorur këtë pronë, mund të merrni të njëjtin rezultat nëse gjeni pjesët e vërteta dhe imagjinare L{eunet}=1sune{ displaystyle { mathcal {L}} {e ^ {it} } = { frac {1} {s-i}}}.
  4. 4 Gjeni transformimin Laplace të derivatit f(t){ displaystyle f ^ { prime} (t)}. Ndryshe nga shembujt e mëparshëm, në këtë rast duhet të integroni pjesë për pjesë:
    • L{f(t)}=0f(t)estdt,  u=est, dv=f(t)dt=f(t)est|0+s0f(t)estdt=sF(s)f(0){ displaystyle { begin {përputhje} { mathcal {L}} {f ^ { prime} (t) } & = int _ {0} ^ { infty} f ^ { prime} (t ) e ^ {- st} mathrm {d} t, u = e ^ {- st}, mathrm {d} v = f ^ { prime} (t) mathrm {d} t & = f (t) e ^ {- st} Big _ {0} ^ { infty} + s int _ {0} ^ { infty} f (t) e ^ {- st} mathrm {d } t & = sF (s) -f (0) fund {rreshtuar}}}
    • Meqenëse derivati ​​i dytë ndodh në shumë probleme fizike, ne gjejmë transformimin e Laplace edhe për të:
      • L{f(t)}=s2F(s)sf(0)f(0){ displaystyle { mathcal {L}} {f ^ { prime prime} (t) } = s ^ {2} F (s) -sf (0) -f ^ { prime} (0) }
    • Në rastin e përgjithshëm, transformimi Laplace i derivatit të rendit të nëntë përcaktohet si më poshtë (kjo lejon zgjidhjen e ekuacioneve diferenciale duke përdorur transformimin Laplace):
      • L{f(n)(t)}=snF(s)k=0n1snk1f(k)(0){ displaystyle { mathcal {L}} {f ^ {(n)} (t) } = s ^ {n} F (s) - sum _ {k = 0} ^ {n -1} s ^ {nk-1} f ^ {(k)} (0)}

Pjesa 3 nga 3: Gjetja e transformimit të Laplace sipas zgjerimit të serisë

  1. 1 Le të gjejmë transformimin Laplace për një funksion periodik. Funksioni periodik plotëson kushtin f(t)=f(t+nT),{ displaystyle f (t) = f (t + nT),} ku T{ displaystyle T} është periudha e funksionit, dhe n{ displaystyle n} është një numër i plotë pozitiv. Funksionet periodike përdoren gjerësisht në shumë aplikime, përfshirë përpunimin e sinjalit dhe inxhinierinë elektrike. Duke përdorur transformime të thjeshta, marrim rezultatin e mëposhtëm:
    • L{f(t)}=0f(t)estdt=n=0nT(n+1)Tf(t)estdt=n=00Tf(t+nT)es(t+nT)dt=n=0esnT0Tf(t)estdt=11esT0Tf(t)estdt{ displaystyle { filloj {në linjë} { mathcal {L}} {f (t) } & = int _ {0} ^ { infty} f (t) e ^ {- st} mathrm { d} t & = sum _ {n = 0} ^ { infty} int _ {nT} ^ {(n + 1) T} f (t) e ^ {- st} matematikë {d} t & = sum _ {n = 0} ^ { infty} int _ {0} ^ {T} f (t + nT) e ^ {- s (t + nT)} mathrm {d} t & = sum _ {n = 0} ^ { infty} e ^ {- snT} int _ {0} ^ {T} f (t) e ^ {- st} mathrm {d} t & = { frac {1} {1-e ^ {- sT}}} int _ {0} ^ {T} f (t) e ^ {- st} mathrm {d} t end { rreshtuar}}}
    • Siç mund ta shihni, në rastin e një funksioni periodik, mjafton të kryeni transformimin Laplace për një periudhë.
  2. 2 Kryeni transformimin Laplace për logaritmin natyror. Në këtë rast, integrali nuk mund të shprehet në formën e funksioneve elementare. Përdorimi i funksionit gama dhe zgjerimi i serisë së tij ju lejon të vlerësoni logaritmin natyror dhe shkallët e tij. Prania e konstantes Euler-Mascheroni γ{ displaystyle gama} tregon se për të vlerësuar këtë integral, është e nevojshme të përdoret një zgjerim serik.
    • L{lnt}=γ+lnss{ displaystyle { mathcal {L}} { ln t } = - { frac { gamma + ln s} {s}}}
  3. 3 Konsideroni transformimin Laplace të funksionit të pa normalizuar të sin. Funksioni i sinqerte(t)=mëkattt{ displaystyle operatornname {sinc} (t) = { frac { sin t} {t}}} përdoret gjerësisht për përpunimin e sinjalit, në ekuacionet diferenciale është ekuivalente me funksionin sferik të Beselit të llojit të parë dhe rendit zero j0(x).{ displaystyle j_ {0} (x).} Transformimi Laplace i këtij funksioni gjithashtu nuk mund të llogaritet me metoda standarde. Në këtë rast, bëhet transformimi i anëtarëve individualë të serisë, të cilët janë funksione të fuqisë, kështu që transformimet e tyre konvergojnë domosdoshmërisht në një interval të caktuar.
    • Së pari, ne shkruajmë zgjerimin e funksionit në një seri Taylor:
      • mëkattt=n=0(1)nt2n(2n+1)!{ displaystyle { frac { sin t} {t}} = sum _ {n = 0} ^ { infty} { frac {(-1) ^ {n} t ^ {2n}} {(2n +1)!}}}
    • Tani ne përdorim transformimin tashmë të njohur të funksionit të fuqisë Laplace. Factorials anulohen, dhe si rezultat marrim zgjerimin Taylor për arctangent, domethënë, një seri alternative që i ngjan serisë Taylor për sinusin, por pa faktoriale:
      • L{mëkattt}=n=0(1)n(2n)!(2n+1)!1s2n+1=n=0(1)n2n+11s2n+1=nxi11s{ displaystyle { fillo {në linjë} { mathcal {L}} majtas {{ frac { sin t} {t}} djathtas } & = sum _ {n = 0} ^ { infty } { frac {(-1) ^ {n} (2n)!} {(2n + 1)!}} { frac {1} {s ^ {2n + 1}}} & = sum _ {n = 0} ^ { infty} { frac {(-1) ^ {n}} {2n + 1}} { frac {1} {s ^ {2n + 1}}} & = tan ^ {- 1} { frac {1} {s}} në fund {në linjë}}}