第九章积分变换法
§9.1傅氏变换
9.1.1傅氏积分和傅氏积分定理
1.周期函数的傅氏级数(周期为2l)
(1)三角式:f(x)=a 0 2 +∑ n=1 ∞ (a n cosnπl x+b n sinnπl x)
a n =1l ∫ l −l f(ξ)cosnπl ξdξ
b n =1l ∫ l −l f(ξ)sinnπl ξdξ
(2)复数式⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ f(x)=∑ n=−∞ ∞ c n e iω n x c n =12l ∫ l −l f(ξ)e −iω n ξ dξ
ω n =nπl ,n=⋯,−2,−1,0,1,2,⋯,Δω n =ω n+1 −ω n =πl
2.非周期函数的傅氏积分
周期(2l)的函数−−→ l→∞ 非周期函数
Δω n =πl −−→ l→∞ 0,ω n →ω
∵f(x)=lim l→∞ ∑ n=−∞ ∞ c n e iω n x =lim l→∞ ∑ n=−∞ ∞ [12l ∫ l −l f(ξ)e −iω n ξ dξ]e iω n x
=lim Δω n →0 ∑ n=−∞ ∞ [12π ∫ ∞ −∞ f(ξ)e −iω n ξ dξ]e iω n x Δω n
f(x)=12π ∫ −∞ ∞ [∫ ∞ −∞ f(ξ)e −iωξ dξ]e iωx dω−−傅氏积分
或f(x)=12π ∫ −∞ ∞ [∫ ∞ −∞ f(x)e −iωx dx]e iωx dω
3.傅氏积分存在的条件(傅氏积分定理)
(1)满足狄氏条件;2)绝对可积;∫ −∞ ∞ |f(x)|dx<∞,
则:f(x)=⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ 12π ∫ −∞ ∞ [∫ ∞ −∞ f(ξ)e −iωξ dξ]e iωx dω,连续点处;12 [f(x 0 +0)+f(x 0 −0)],x 0 间断点。
4.三维傅氏积分
f(r ⃗ )=1(2π) 3 ∭ ∞ −∞ [∭ ∞ −∞ f(r ⃗ )e −iω ⃗ r ⃗ dr ⃗ ]e iω ⃗ r ⃗ dω ⃗
r ⃗ =e ⃗ 1 x+e ⃗ 2 y+e ⃗ 3 z
ω ⃗ =e ⃗ 1 ω 1 +e ⃗ 2 ω 2 +e ⃗ 3 ω 3
dr ⃗ =dxdydz,dω ⃗ =dω 1 dω 2 dω 3
9.1.2傅氏变换
1.(一维)傅氏变换
定义:G(ω)=∫ ∞ −∞ f(ξ)e −iωξ dξ=F[f(x)]−−f(x)的傅氏变换
则f(x)=12π ∫ ∞ −∞ G(ω)e iωx dω=F −1 [G(ω)]−−G(ω)的傅氏逆变换
2.例题
(1)f(t)={0,t<0e −βt ,t>0
1)求F[f(t)]=?2)f(t)的积分表示式;
F[f(t)]=β−iωβ 2 +ω 2
f(t)=1π ∫ ∞ 0 βcosωt+ωsinωtβ 2 +ω 2 dω
∫ ∞ 0 βcosωt+ωsinωtβ 2 +ω 2 dω=⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ 0,t<0π2 ,t=0πe −βt ,t>0
(2)F[e −ax 2 =?(a>0)
F[e −ax 2 =e −ω 2 4a πa − − √
3.三维傅氏变换
定义:G(ω ⃗ )=∭ ∞ −∞ f(r ⃗ )e −iω ⃗ r ⃗ dr ⃗ ,则f(r ⃗ )=1(2π) 3 ∭ ∞ −∞ G(ω ⃗ )e iω ⃗ r ⃗ dω ⃗
dr ⃗ =dxdydz,dω ⃗ =dω 1 dω 2 dω 3
4.傅氏变换的其它几种形式
(1)⎧ ⎩ ⎨ G(ω)=12π ∫ ∞ −∞ f(x)e −iωx dxf(x)=∫ ∞ −∞ G(ω)e iωx dω
(2)⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ G(ω)=12π − − √ ∫ ∞ −∞ f(x)e −iωx dxf(x)=12π − − √ ∫ ∞ −∞ G(ω)e iωx dω
(3)⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ G(ω)=12π − − √ ∫ ∞ −∞ f(x)e iωx dxf(x)=12π − − √ ∫ ∞ −∞ G(ω)e −iωx dω
5.傅氏变换的物理意义
G(ω)=∫ ∞ −∞ f(ξ)e −iωξ dξ
G(ω)为f(x)的频率密度函数或频谱函数,它可用来反映各种频率谐波之间振幅的相对大小,并称|G(ω)|为f(x)的频谱。因为ω是相对变化的,所以f(x)的频谱是良序普。而
f(x)=12π ∫ ∞ −∞ G(ω)e iωx dω
可解释为无穷多个振幅(复振幅)为无限小的,频率为连续的谐波的连续和。
9.1.3傅氏变换的性质
1.F[αf 1 +βf 2 ]=αF[f 1 ]+βF[f 2 ]线性
2.F[e iω 0 x f(x)]=G(ω−ω 0 )[设F[f(x)]=G(ω)]延迟
3.F[f(x±x 0 )]=e ±iωx 0 F[f(x)]位移
4.F[f (n) (x)]=(iω) n F[f(x)],f (n−1) (x)−−→ |x|→∞ 0,n=1,2,⋯微分
5.F[∫ x x 0 f(ξ)dξ]=1iω F[f(x)]积分
F[f 1 ∗f 2 ]=F[f 1 ]⋅F[f 2 ]卷积定理
F[f 1 ⋅f 2 ]=12π F[f 1 ]∗F[f 2 ]像函数卷积
其中f 1 ∗f 2 =∫ ∞ −∞ f 1 (ξ)f 2 (x−ξ)dξ→卷积
例题:
(1)F[xe −ax 2 ]=−iω2a πa − − √ e −ω 2 4a
(2)已知F[φ(x)]=G(ω),求1)F −1 [G(ω)cosaωt],2)求F −1 [G(ω)aω sinaωt].
1)F −1 [G(ω)cosaωt]=12 [φ(x+at)+φ(x−at)]
2)F −1 [G(ω)aω sinaωt]=12a ∫ x+at x−at φ(ξ)dξ
小结
定义:G(ω)=∫ ∞ −∞ f(ξ)e −iωξ dξ=F[f(x)]−f(x)的傅氏变换
则:f(x)=12π ∫ ∞ −∞ G(ω)e iωx dω=F −1 [G(ω)]−G(ω)的傅氏逆变换
显然:F −1 F[f(x)]]=f(x)
§9.2傅里叶变换法
9.2.1波动问题
⎧ ⎩ ⎨ ⎪ ⎪ u tt −a 2 u xx =0,−∞<x<∞,t>0(1)u| t=0 =φ(x)(2)u t | t=0 =ψ(x)(3)
1.对定解问题各项施行傅氏变换
记∫ ∞ −∞ u(x,t)e −iωx dx=u ~ (ω,t),∫ ∞ −∞ φ(x)e −iωx dx=φ ~ (ω)
∫ ∞ −∞ ψ(x)e −iωx dx=ψ ~ (ω),
则⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ d 2 u ~ (ω,t)dt 2 +a 2 ω 2 u ~ (ω,t)=0(4)u ~ (ω,0)=φ ~ (ω)(5)u ~ t (ω,0)=ψ ~ (ω)(6)
2.解常微分方程定解问题(4) (6)得
u ~ (ω,t)=φ ~ (ω)cosαωt+1aω ψ ~ (ω)sinaωt
3.求傅氏逆变换
u(x,t)=F −1 [u ~ (ω,t)]=F −1 [φ ~ (ω)cosaωt]+F −1 [ψ ~ (ω)aω sinaωt]
u(x,t)=12 [φ(x+at)+φ(x−at)]+12a ∫ x+at x−at ψ(α)dα
9.2.2输运问题
1.{u t −a 2 u xx =f(x,t),−∞<x<∞,t>0(7)u(x,0)=φ(x)(8)
(1)对定解问题各项施行傅氏变换
记∫ ∞ −∞ u(x,t)e −iωx dx=u ~ (ω,t),∫ ∞ −∞ f(x,t)e −iωx dx=f ~ (ω,t),
∫ ∞ −∞ φ(x)e −iωx dx=φ ~ (ω),
则⎧ ⎩ ⎨ du ~ dt +a 2 ω 2 u ~ =f ~ (ω,t)(9)u ~ (ω,0)=φ ~ (ω)(10)
(2)解常微分方程定解问题(9) (10)得
u ~ (ω,t)=e −a 2 ω 2 t φ ~ (ω)+∫ t 0 f ~ (ω,τ)e −a 2 ω 2 (t−τ) dτ
(3)求傅氏逆变换
F −1 [e −a 2 ω 2 t ]=12aπt − − √ e −x 2 4a 2 t
u(x,t)=F −1 [φ ~ (ω)e −a 2 ω 2 t ]+∫ t 0 F −1 [f ~ (ω,τ)e −a 2 ω 2 (t−τ) ]dτ
=F −1 F[φ(x)∗F −1 (e −a 2 ω 2 t )]+∫ t 0 F −1 F[f(x,τ)∗F −1 (e −a 2 ω 2 (t−τ) )]dτ
=φ(x)∗12aπt − − √ e −x 2 4a 2 t +∫ t 0 f(x,τ)∗12aπ(t−τ) − − − − − − − √ e −x 2 4a 2 (t−τ) dτ
u(x,t)=12aπt − − √ ∫ −∞ ∞ φ(ξ)e −(ξ−x) 2 4a 2 t dξ+12aπ √ ∫ 0 t 1t−τ − − − − √ ∫ −∞ ∞ f(ξ,τ)e −(x−ξ) 2 4a 2 (t−τ) dξdτ
2.{u t −a 2 Δu=0u| t=0 =φ(r ⃗ )
F −1 [e −a 2 ω 2 t ]=F −1 [e −a 2 (ω 2 1 +ω 2 2 +ω 2 3 )t ]=(12aπt − − √ ) 3 e x 2 +y 2 +z 2 −4a 2 t
记F[u(r ⃗ ,t)]=∭ ∞ −∞ u(r ⃗ ,t)e −iω ⃗ ⋅r ⃗ dr ⃗ =u ~ (ω ⃗ ,t)
F[φ(r ⃗ )]=∭ ∞ −∞ φ(r ⃗ )e −iω ⃗ ⋅r ⃗ dr ⃗ =φ ~ (ω ⃗ )
F[Δu]=∭ ∞ −∞ (∂ 2 u∂x 2 +∂ 2 u∂y 2 +∂ 2 u∂z 2 )e −iω ⃗ ⋅r ⃗ dr ⃗ =−ω 2 u ~ (ω ⃗ ,t)
⎧ ⎩ ⎨ ⎪ ⎪ du ~ (ω ⃗ ,t)dt +a 2 ω 2 u ~ (ω ⃗ ,t)=0u ~ (ω ⃗ ,0)=φ ~ (ω ⃗ )
u ~ (ω ⃗ ,t)=φ ~ (ω ⃗ )e −a 2 ω 2 t
u(r ⃗ ,t)=18a 3 (πt) 32 ∭ ∞ −∞ φ(r 1 → )e −|r ⃗ −r ⃗ 1 | 2 4a 2 t dr ⃗ 1
附:F[Δu]=∭ ∞ −∞ (∂ 2 u∂x 2 +∂ 2 u∂y 2 +∂ 2 u∂z 2 )e −iω ⃗ ⋅r ⃗ dr ⃗
=∭ ∞ −∞ (∂ 2 u∂x 2 +∂ 2 u∂y 2 +∂ 2 u∂z 2 )e −i(ω 1 x+ω 2 y+ω 3 z) dxdydz
=∬[∫ ∞ −∞ (∂ 2 u∂x 2 +∂ 2 u∂y 2 +∂ 2 u∂z 2 )e −iω 1 x dx]e −iω 2 y dye −iω 3 z dz
=∫ −∞ ∞ ∫ −∞ ∞ [(iω 1 ) 2 u ~ (ω 1 ,y,z)+∂ 2 u ~ (ω 1 ,y,z)∂y 2 +∂ 2 u ~ (ω 1 ,y,z)∂z 2 ]e −iω 2 y dye −iω 3 z dz
=∫ −∞ ∞ [−ω 2 1 u ~ (ω 1 ,ω 2 ,z)−ω 2 2 u ~ (ω 1 ,ω 2 ,z)+∂ 2 u ~ (ω 1 ,ω 2 ,z)∂z 2 ]e −iω 3 z dz
=−u ~ (ω 1 ,ω 2 ,ω 3 )[ω 2 1 +ω 2 2 +ω 2 3 ]
=−u ~ (ω 1 ,ω 2 ,ω 3 )ω 2
=−ω 2 u ~ (ω ⃗ )
9.2.3稳定问题
Δu=−1ε 0 ρ(x,y,z)
记u(x,y,z)=u(r ⃗ )
1ε 0 ρ(x,y,z)=f(r ⃗ )
∭ ∞ −∞ u(r ⃗ )e −iω ⃗ ⋅r ⃗ dr ⃗ =u ~ (ω ⃗ )
∭ ∞ −∞ f(r ⃗ )e −iω ⃗ ⋅r ⃗ dr ⃗ =f ~ (ω ⃗ )
ω 2 u ~ (ω ⃗ )=f(ω ⃗ )
u ~ (ω ⃗ )=f ~ (ω ⃗ )ω 2 =F[f(r ⃗ )]⋅F[F −1 [1ω 2 ]]=F[F −1 [1ω 2 ]∗f(r ⃗ )]
∵F[1r ]=4πω 2 ,∴F −1 [1ω 2 ]=14πr =14π|r ⃗ |
∴u(r ⃗ )=F −1 [u ~ (ω ⃗ )]
u(r ⃗ )=14π ∭ ∞ −∞ f(r ⃗ ′ )|r ⃗ −r ⃗ ′ | dr ⃗ ′
9.2.4小结
1.傅氏变换法精神:
对方程中各项施行傅氏变换,从而将偏微分方程化为常微分方程求解。
2.傅氏变换法的解题步骤:
(1)对方程中各项选择适当变量施行傅氏变换
(2)解像函数的常微分方程的定解问题
(3)求逆变换的原定问题的解
3.由傅氏变换法得到的三种典型定解问题的解:
⎧ ⎩ ⎨ ⎪ ⎪ u tt −a 2 u xx =0,−∞<x<∞,t>0u| t=0 =φ(x)u t | t=0 =ψ(x)
u(x,t)=12 [φ(x+at)+φ(x−at)]+12a ∫ x+at x−at ψ(α)dα
{u t −a 2 u xx =f(x,t)u(x,0)=φ(x)
u(x,t)=12aπt − − √ ∫ ∞ −∞ φ(ξ)e −(ξ−x) 2 4a 2 t dξ+12aπ √ ∫ t 0 1t−τ − − − − √ ∫ ∞ −∞ f(ξ,τ)e −(x−ξ) 2 4a 2 (t−τ) dξdτ
Δu=−1ε 0 ρ(x,y,z)
u(r ⃗ )=14π ∭ ∞ −∞ f(r ⃗ ′ )|r ⃗ −r ⃗ ′ | dr ⃗ ′
§9.3拉普拉斯变换法
9.3.1拉氏变换
问题引入:∫ ∞ −∞ |f(t)|dt<∞难!
若β>0;当t<0时f(t)=0,则∫ ∞ −∞ |f(t)e −βt dt<∞易,
此时:F[f(t)e −βt ]=∫ ∞ 0 f(t)e −βt e −iωt dt
而:f(t)e −βt =12π ∫ ∞ −∞ F[f(t)e −βt ]e iωt dω
1.定义:记p=β+iω,F(p)=F[f(t)e −βt ],则dp=idω
F(p)=∫ ∞ 0 f(t)e −pt dt−f(t)的拉氏变换
f(t)=12πi ∫ β+i∞ β−i∞ F(p)e pt dp−F(p)的拉氏逆变换
2.存在条件
(1)f(t)及导数除有限个第一类间断点外连续
(2)|f(t)|≤Me β 0 t (M,β 0 ≥0;β 0 是增长指数)
例:(1)L[e at ]=∫ ∞ 0 e at e −pt dt=1p−a ,Rep>Rea
(2)L[t k ]=k!p k+1 =Γ(k+1)p k+1 ,Rep>0
L[t 0 ]=∫ ∞ 0 e −pt dt=∫ ∞ 0 e 0t e −pt dt=1p ,Rep>0
L[t 1 ]=∫ ∞ 0 te −pt dt=−1p ∫ ∞ 0 tde −pt =1p 2 ,Rep>0
9.3.2拉氏变换性质
1.线性:L[αf 1 (t)+βf 2 (t)]=αL[f 1 (t)]+βL[f 2 (t)]
2.延迟:L[e p 0 t f(t)]=F(p−p 0 )记L[f(t)]=F(p)
3.位移:L[f(t−τ)]=e −pτ F[p]
4.相似:L[f(at)]=1a F(pa ),a>0
5.微分:L[f (n) (t)]=p n F(p)−p n−1 f(0)−p n−2 f ′ (0)−⋯−f (n−1) (0)
6.积分:L[∫ t 0 f(τ)dτ]=1p L[f(t)]
7.卷积:L[f 1 (t)∗f 2 (t)]=L[f 1 (t)]⋅L[f 2 (t)]
f 1 (t)∗f 2 (t)=∫ t 0 f 1 (τ)f 2 (t−τ)dτ
例:(1)L[sinkt]=L[e ikt −e −ikt 2i ]
=12i [1p−ik −1p+ik ]=kp 2 +k 2 ,Rep>0
(2)L[sin(t−2π3 )]=e −2π3 p L[sint]=e −2π3 p 1p 2 +1
(3)L[coskt]=L[1k ddt sinkt]=1k p⋅kp 2 +k 2 =pp 2 +k 2
(4)已知F(p)=p 2 (p 2 +1) 2 ,求L −1 [F(p)]=?
L −1 [F(p)]=L −1 [pp 2 +1 ⋅pp 2 +1 ]
=L −1 ⋅L[cost∗cost]=∫ t 0 cosτcos(τ−t)dτ
=∫ t 0 12 [cos(2τ−t)+cost]dτ
=14 ⋅∫ t 0 cosαdα+12 tcost
=12 (tcost+12 sint)
9.3.3原函数存在定理
拉氏反演及展开定理:
若F(p)单值,在0≤argz≤2π中,当p→∞,F(p)→0,则f(t)=∑ k res[F[p k ]e p k t ],p k →全平面奇点
例:F(p)=1(p+1)(p−3) 2 ,求f(t)=?
f(t)=res[e pt (p+1)(p−3) 2 ,−1]+res[e pt (p+1)(p−3) 2 ,3]
=e pt (p−3) 2 | p=−1 +ddp [e pt p+1 ] p=3
=e −t 16 +te 3t 4 −e 3t 16
9.3.4解数理方程
1.⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ T ′′ (t)+(nπal ) 2 T(t)=f(t)T(0)=0T ′ (0)=0
记L[T(t)]=∫ ∞ 0 T(t)e −pt dt=T ~ (p),
L[f(t)]=∫ ∞ 0 f(t)e −pt dt=f ~ (p)
则p 2 T ~ (p)−pT(0)−T ′ (0)+(nπal ) 2 T ~ (p)=f ~ (p)
∴T ~ (p)=f ~ (p)p 2 +(nπal ) 2 =L[f(t)∗lnπa sinnπal t]
T(t)=lnπa ∫ t 0 f(τ)sinnπal (t−τ)dτ
2.解混合问题
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ u tt =a 2 u xx ,0<x<∞,t>0u(0,t)=f(t),lim x→∞ u(x,t)=0(t≥0)u(x,0)=0,u t (x,0)=0
记L[u(x,t)]=u ~ (x,p),L[f(t)]=F ~ (p),则:
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ p 2 u ~ (x,p)−pu(x,0)−u t (x,0)=a 2 ∂ 2 ∂x 2 u ~ (x,p)u ~ (0,p)=F ~ (p),lim x→∞ u ~ (x,p)=0
即⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ d 2 dx 2 u ~ (x,p)−p 2 a 2 u ~ (x,p)=0u ~ (0,p)=F ~ (p)lim x→∞ u ~ (x,p)=0
u ~ =c 1 (p)e −pa x +c 2 (p)e pa x
由u ~ (0,p)=F ~ →c 1 (p)+c 2 (p)=F ~ (p)
由u ~ (∞,p)=0→c 2 (p)=0
u ~ (x,p)=F ~ (p)e −pa x →u(x,t)=L −1 [F ~ (p)⋅e −pa x ]
u(x,t)=L −1 L[f(t−xa )]=f(t−xa )
3.⎧ ⎩ ⎨ ⎪ ⎪ u tt −a 2 u xx =0,−∞<x<∞,t>0u(x,0)=φ(x)u t (x,0)=ψ(x)
(1)令F[u(x,t)]=u ~ (ω,t),F[φ(x)]=φ ~ (ω),F[ψ(x)]=ψ ~ (ω)
⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ d 2 u ~ dt 2 +a 2 ω 2 u ~ (ω,t)=0,t>0u ~ (ω,0)=φ ~ (ω),u ~ t (ω,0)=ψ ~ (ω)
(2)记L[u ~ (ω,t)]=U(ω,p)
则p 2 U(ω,p)−pu ~ (ω,0)−u ~ t (ω,0)+a 2 ω 2 U(ω,p)=0
即p 2 U(ω,p)−pφ ~ (ω)−ψ ~ (ω)+a 2 ω 2 U(ω,p)=0
(3)U(ω,p)=pφ ~ (ω)+ψ ~ (ω)p 2 +a 2 ω 2
=φ ~ (ω)pp 2 +a 2 ω 2 +ψ ~ (ω)aω aωp 2 +a 2 ω 2
(4)u ~ (ω,t)=φ ~ (ω)cosaωt+ψ(ω)aω sinaωt
u(x,t)=12 [φ(x+at)+φ(x−at)]+12a ∫ x+at x−at ψ(α)dα