第二十七讲积分变换 介绍求解偏微分方程定解问题的另一种做法,即将积分变换应用于在求解偏微分方 程定解问题. 常用的积分变换有 Laplace变换和 Fourier变换两种 §27.1应用 Laplace变换求解偏微分方程定解问题 Laplace变换可用于求解含时间的偏微分方程定解问题.变换后,自变量的个数比原来减少 一个.例如,原来是x和t两个自变量的偏微分方程定解问题,变换后就只需求解常微分方程(自 变量为x)的定解问题.一般说来,后者总比较容易求解.这样求得的是原始的定解问题的解的象 函数,还必须反演,才能得到原始问题的解 例271求无界杆的热传导问题 -6x2=f(x,,-∞<r<∞,t>0 ∞<x<∞ 解 解在这种无界区间的定解问題中,往往并不明确列出边界条件.实际上,无界区间,只是 一个物理上的抽象,它只是表明在所考察的限度(时间,精度,…)內,两端的影响可以忽略.因 此,如果要完整地列出定解问题的话,则还应当有边界条件 作 Laplace变换.令 u(a, t)=U(a, p) 利用初始条件,有 oU(E, P) 把变换后的象函数只看成是x的函数,p是参数,所以 a2u d2U d r2 微商运算就是一元函数的微商.再进一步令 f(r, t): F(a, p) 这样,在经过 Laplace变换后,定解问题就变成 pU(E, P)-K d2U(x,p)∠F(x,p) 0的条件下可以得到 U(r, p) F(a, p)ex
Wu Chong-shi ✁✂✄☎ ✆✝✞✟ ✠✡☛☞✌✍✎✏✑✒☞✓✔✕✖✗✘✙✚✛✜✢✣✎✤✥✦✧★✩☛☞✌✍✎✏ ✑✒☞✓✔✪ ✫✧✕✣✎✤✥✬ Laplace ✤✥✭ Fourier ✤✥✮✘✪ §27.1 ✯✰ Laplace ✱✲✳✴✵✶✷✸✹✺✴✻✼ Laplace ✤✥✽✧★☛☞✾✿❀✕✌✍✎✏✑✒☞✓✔✪✤✥❁✛ ❂✤❃✕❄❅ ❆❇❈❉❊ ✗❄✪❋●✛❇❈❍ x ✭ t ✮❄ ❂✤❃✕✌✍✎✏✑✒☞✓✔✛✤✥❁■❏❑☛☞✫✍✎✏✑ ( ❂ ✤❃▲ x) ✕✒☞✓✔✪✗▼◆❈✛❁❖P ❆◗❘❙☛☞✪❚❯☛❱✕❍❇❲✕✒☞✓✔✕☞✕❳ ❨❅✛❩❬❭❪❫✛❴❵❱❛❇❲✓✔✕☞✪ ❜ 27.1 ☛❝❞❡✕❢❣❤✓✔ ∂u ∂t − κ ∂ 2u ∂x2 = f(x, t), − ∞ < x < ∞, t > 0; u t=0 = 0, − ∞ < x < ∞ ✕☞✪ ✐ ❥❦❧♠♥ ♦♣qrs t✉ ✈✛✇✇①② ③④⑤ ⑥⑦♥⑧⑨✪⑩❶❷✛♠♥ ♦♣✛❸❹ ❺❻❼❽❷q❾❿✛➀ ❸❹➁ ③❥➂➃➄q➅➆ (➇ ♣✛➈➆✛ · · ·) ➉ ✛➊➋q➌➍➎ ➏➐➑✪➒ ➓✛➔→➣ ↔↕➙⑤ ⑥rs t✉q➛✛➜➝➞ ➟➠⑦♥⑧⑨ u x→±∞ → 0. ➡ Laplace ✤✥✪➢ u(x, t) ; U(x, p) = Z ∞ 0 u(x, t)e−ptdt, ➤✧➥❲➦➧✛✬ ∂u ∂t ; pU(x, p). ➨✤✥❁✕❳❨❅❏➩➫❍ x ✕❨❅✛ p ❍➭❅✛➯➲ ∂ 2u ∂x2 ; d 2U(x, p) dx 2 , ✍➳➵➸■❍✗➺❨❅✕✍➳✪➻➼✗➽➢ f(x, t) ; F(x, p), ❚❯✛✩➾➚ Laplace ✤✥❁✛✒☞✓✔■✤➫ pU(x, p) − κ d 2U(x, p) dx 2 = F(x, p). ✩ U x→±∞ → 0 ✕➦➧➪✽➲❱❛ U(x, p) = 1 2 1 √κp Z ∞ −∞ F(x 0 , p) exp − r p κ |x − x 0 | dx 0
§27.1应用 Laplace变换求解偏微分方程定解问题 再根据 Laplace变换的反演公式 以及卷积定理,就能够最后得到 u(a, t) dr'/exp (-r)Id 用 Laplace变换求解偏微分方程定解间题,除了可以减少自变量的数目以外,某些已 知函数的象函数(例如方程的非齐次项,它的形式可能很复杂)甚至都不必具体求出 在求反演时只需应用卷积定理即可 例27.2用 Laplace变换求解无界弦的波动问题 2 0, φ(x) 解设在 Laplace变换之 (ar, t)=U(r, p). 于是,原来的定解问题就化为 P-U(r, p)-a d-U,2=po(r)+v(r) 可以求得此方程的解(实际上还考虑了U(x,p)在x→士∞的行为 U(a, p) po(r)+v()exp P 因为 6(t-a) 7(t-a) 所以 o(x)6( v(x)( p(a)(at-r-a'Ddx+ N v(a)n(at-lr-r'Ddr 注意到 6(at-|z-x|) 0,1x-x1≠at 0,|x-xl> at 1, r-r'<at 就可以求出 (x,t) 2a v(a')dr
Wu Chong-shi §27.1 ➶➹ Laplace ➘➴➷➬➮➱✃❐❒❮➬❰Ï Ð 2 Ñ ➻ÒÓ Laplace ✤✥✕❪❫ÔÕ 1 √ p e −α √p : 1 √ πt exp − α 2 4t ➲Ö×✣✒Ø✛■❵ÙÚ❁❱❛ u(x, t) = 1 2 √ κπ Z ∞ −∞ dx 0 Z t 0 exp − (x − x 0 ) 2 4κ(t − τ) f(x 0 , τ) √ t − τ dτ. ✧ Laplace ✤✥☛☞✌✍✎✏✑✒☞✓✔✛ÛÜ✽➲❉❊ ❂✤❃✕❅ Ý➲Þ✛ßà á â❨❅✕❳❨❅ (❋●✏✑✕ãäåæ✛ç✕èÕ✽❵éêë) ìíîï❬ðñ☛ò✛ ✩☛❪❫✿❏❑✦✧×✣✒Ø✜✽✪ ❜ 27.2 ✧ Laplace ✤✥☛☞❝❞ó✕ôõ✓✔ ∂ 2u ∂t2 − a 2 ∂ 2u ∂x2 = 0, − ∞ < x < ∞, t > 0; u t=0 = φ(x), ∂u ∂t t=0 = ψ(x), − ∞ < x < ∞. ✐ ö✩ Laplace ✤✥÷➪✛ u(x, t) ; U(x, p), ★❍✛❇❈✕✒☞✓✔■ø▲ p 2U(x, p) − a 2 d 2U(x, p) dx 2 = pφ(x) + ψ(x). ✽➲☛❱ù✏✑✕☞ (úûü❩ýþÜ U(x, p) ✩ x → ±∞ ✕ÿ▲) U(x, p) = 1 2ap Z ∞ −∞ h pφ(x 0 ) + ψ(x 0 ) i exp n − p a x − x 0 o dx 0 = 1 2a Z ∞ −∞ h φ(x 0 ) + ψ(x 0 ) p i expn − p a x − x 0 o dx 0 . ▲ e −αp : δ(t − α), 1 p e −αp : η(t − α), ➯➲ u(x, t) = 1 2a Z ∞ −∞ φ(x 0 ) δ t − |x − x 0 | a dx 0 + 1 2a Z ∞ −∞ ψ(x 0 ) η t − |x − x 0 | a dx 0 = 1 2 Z ∞ −∞ φ(x 0 ) δ(at − |x − x 0 |) dx 0 + 1 2a Z ∞ −∞ ψ(x 0 ) η(at − |x − x 0 |) dx 0 , ✁✂❛ δ(at − |x − x 0 |) = ( 0, |x − x 0 | 6= at, ∞, |x − x 0 | = at, η(at − |x − x 0 |) = ( 0, |x − x 0 | > at, 1, |x − x 0 | < at, ■✽➲☛ò u(x, t) = 1 2 Z ∞ −∞ φ(x 0 )δ t − |x − x 0 | a dx 0 + 1 2a Z x+at x−at ψ(x 0 )dx 0
第二十七讲积分变换 第3页 -at)+o(=+at)+ v(adr x一at 用 Laplace变换求解偏微分方程定解问题还有一个优点,这就是不必将非齐次的边 界条件齐次化,因为这时原有的偏微分方程定解问题的非齐次边界条件将转化为常微 分方程的非齐次边界条件,这并不会带来原则的困难
Wu Chong-shi ✄☎✆✝✞ ✟ ✃➘➴ Ð 3 Ñ = 1 2 h φ(x − at) + φ(x + at) i + 1 2a Z x+at x−at ψ(x 0 )dx 0 . ✧ Laplace ✤✥☛☞✌✍✎✏✑✒☞✓✔❩✬✗❄✠✡✛❚■❍ï ❬✢ãäå✕☛ ❞➦➧äåø✛▲❚✿❇✬✕✌✍✎✏✑✒☞✓✔✕ãäå☛❞➦➧✢☞ø▲✫✍ ✎✏✑✕ãäå☛❞➦➧✛❚✌ ï✍✎❈❇✏✕✑✒✪
§272 Fourier变换 第4页 8272 Fourier变换 Fourier变换可对空间变量进行.根据空间变量的变化区间,可以选用 · curler变换对于无界区间(-∞,∞)上的函数f(x),如果在任意有限区间上只有有限个极 大极小和有限个第一类间断点,且 积分/f(x)dx绝对收敛,则它的 Fourier变换存在 F(k)=U()]=6/f()- dr, 而逆变换(反演)是 f(x)=[F(k)≡ k)eirik 这里的 Fourier变换和逆变换的形式可能和读者熟悉的形式略有不同.形式更加对 称,更多地为物理学家所采用 ·正弦变换和余弦变换如果f(x)是定义在半无界区间,∞)上,则可根据x=0端边界条件 的不同类型,选用正弦变换 F()=v#/ f(a) sin kidz, f()=V元F(6)skdk 或余弦变换 F(k) f(ar)cos kadz f(x)= F()cos cdk ·有限正弦、余弦变换如果f(x)是定义在有界区间上,则应采用有限正弦或余弦变换 本节介绍无界空间上的 Fourier变换 用 Fourier变换来求解上一节的例271和例272
Wu Chong-shi §27.2 Fourier ➘➴ Ð 4 Ñ §27.2 Fourier ✱✲ Fourier ✤✥✽✓✔❀✤❃➼ÿ✪ÒÓ✔❀✤❃✕✤ø✕❀✛✽➲✖✧ • Fourier ✤✥ ✓★ ✗✘✙✚(−∞, ∞) ü ✕❨❅ f(x) ✛●✛✩✜✂✬✢✕❀ ü ❏✬✬✢❄✣ ✤✣✥✭✬✢❄✦✗✧❀★✡✛✩ ✣✎ Z ∞ −∞ f(x)dx ✪ ✓✫✬✛✏ç✕ Fourier ✤✥✭✩✛ F(k) = ✮ [f(x)] ≡ 1 √ 2π Z ∞ −∞ f(x)e−ikxdx, ✯✰✤✥ (❪❫) ❍ f(x) = ✮ −1 [F(k)] ≡ 1 √ 2π Z ∞ −∞ F(k)eikxdk. ❚✱✕ Fourier ✤✥✭✰✤✥✕èÕ✽❵✭✲❖✳✴✕èÕ✵✬ ï✶ ✪èÕ✷✸✓ ✹✛✷✺✻▲✼Ø✽✾➯✿✧✪ • ❀ ó✤✥✭❁ó✤✥ ●✛ f(x) ❍✒❂✩ ❃✗✘✙✚[0, ∞) ü ✛✏✽ÒÓ x = 0 ❄ ☛❞➦➧ ✕ ï✶ ✧❅✛✖✧ ❀ ó✤✥ F(k) = r 2 π Z ∞ 0 f(x) sin kxdx, f(x) = r 2 π Z ∞ 0 F(k) sin kxdk, ❆❁ó✤✥ F(k) = r 2 π Z ∞ 0 f(x) cos kxdx, f(x) = r 2 π Z ∞ 0 F(k) cos kxdk.. • ✬✢ ❀ ó❇❁ó✤✥ ●✛ f(x) ❍✒❂✩✬❞✕❀ ü ✛✏✦✿✧✬✢ ❀ ó❆❁ó✤✥✪ ❈❉✠✡❝❞✔❀ ü ✕ Fourier ✤✥✪ ✧ Fourier ✤✥❈☛☞ü ✗❉✕❋ 27.1 ✭❋ 27.2 ✪
第二十七讲积分变换 第5页 ★对于例27.1,即无界杆的热传导问题 0-b2f(x,t),-∞<x<∞,t>0 0<x<0. 假设u(x,t)的 Fourier变换存在 、-山,并设F(=层0a 这样,在作 Fourier变换后,定解问题就变为 dU(h, t dt+Ek U(k, t)=F(k, U(k, t) 用常数变易法求解这个一阶常微分方程的初值问题,就得到 U(k, t)=ekAt/F(k, r)erk4Tdr 再求反演 u()=/(.dk 家/F(k(=ka 利用 t cos 2xt dt=-VTe 可以算出 :k(t-r)eikrd -kk(t-T) cos krak exp 一 2 一c T) 再利用 f(a, t)= F(k, t)edk, 根据 Fourier变换的卷积公式 f1(xf2(x)= 27 fi1()f2(x-5)d 就能最后得到 u(a, t) { f(,r) (x-) V2J∞√2k(t-r) k(t- T) 和上一节中得到的解式完全一样 从解法上看, Fourier变换的反演问题似乎要比 Laplace变换简单一些,往往不需要 用留数定理来计算反演中出现的定积分.就本例而言,两种方法都要用到卷积公式
Wu Chong-shi ✄☎✆✝✞ ✟ ✃➘➴ Ð 5 Ñ F ✓★❋ 27.1 ✛✜❝❞❡✕❢❣❤✓✔ ∂u ∂t − κ ∂ 2u ∂x2 = f(x, t), − ∞ < x < ∞, t > 0; u t=0 = 0, − ∞ < x < ∞. ❊ö u(x, t) ✕ Fourier ✤✥✭✩✛ U(k, t) = 1 √ 2π Z ∞ −∞ u(x, t)e−ikxdx, ✌ö F(k, t) = 1 √ 2π Z ∞ −∞ f(x, t)e−ikxdx, ❚❯✛✩➡ Fourier ✤✥❁✛✒☞✓✔■✤▲ dU(k, t) dt + κk2U(k, t) = F(k, t), U(k, t) t=0 = 0, ✧✫❅✤❙✚☛☞❚❄✗❋✫✍✎✏✑✕➥●✓✔✛■❱❛ U(k, t) = e−κk2 t Z t 0 F(k, τ)eκk2 τ dτ. ➻☛❪❫✛ u(x, t) = 1 √ 2π Z ∞ −∞ U(k, t)eikxdk = Z t 0 1 √ 2π Z ∞ −∞ F(k, τ)e−κk2 (t−τ) e ikxdk dτ. ➤✧ Z ∞ 0 e −t 2 cos 2xt dt = 1 2 √ πe −x 2 , ✽➲➸ò 1 √ 2π Z ∞ −∞ e −κk2 (t−τ) e ikxdk = 1 √ 2π Z ∞ −∞ e −κk2 (t−τ) cos kxdk = 1 p 2κ(t − τ) exp − x 2 4κ(t − τ) , ➻➤✧ f(x, t) = 1 √ 2π Z ∞ −∞ F(k, t)eikxdk, ÒÓ Fourier ✤✥✕×✣ÔÕ✛ ✮ [f1(x)] ✮ [f2(x)] = ✮ 1 √ 2π Z ∞ −∞ f1(ξ)f2(x − ξ)dξ , ■❵Ú❁❱❛ u(x, t) = Z t 0 ( 1 √ 2π Z ∞ −∞ f(ξ, τ) p 2κ(t − τ) exp − (x − ξ) 2 4κ(t − τ) dξ ) dτ = 1 2 √ κπ Z t 0 Z ∞ −∞ f(ξ, τ) exp − (x − ξ) 2 4κ(t − τ) dξ dτ √ t − τ . ✭ ü ✗❉ ❍❱❛✕☞Õ■❏✗❯✪ ❑☞✚ü ➩✛ Fourier ✤✥✕❪❫✓✔▲▼◆ ❆ Laplace ✤✥❖P✗à✛◗◗ï ❑◆ ✧ ❘❅✒Ø❈❙➸❪❫ ❍ò❚✕✒✣✎✪■❈❋✯❯✛✮✘✏✚î ◆✧❛×✣ÔÕ✪