数学物理方程第四章 积分变换法
第四章积分变换法详解
第二十五页,共47页。
L(eat ) 1 , pa
p L(cos at ) p2 a2
4.3 拉普拉斯变换的概念和性质
2)线性性质
L f g L f L g
3) 微分性质
若 F ( p) L[ f (t)], 则
L[ f 't ] pF p f 0 ,
L[ f ''t ] p2F p p f 0 f '0, L[ f n t ] pnF p pn1 f 0 pn2 f '0
u x, y FouxrierU , y
由傅立叶变换的线性性质
u y
x,
y
Fouxrier
U
,
y
y
d dy
U
,
y
是参数
同理,
第九页,共47页。
2u y2
x
,
y
Fouxrier
d2 dy 2
U
,
y
4.2 傅立叶变换的应用
第十页,共47页。
4.2 傅立叶变换的应用
例 用积分变换法解方程:
1 t fˆ (, ) sina(t ) d
a 0
由初始条件 U (,t) ()cosat () sinat a
第十八页,共47页。
1 t fˆ (, ) sina(t ) d
a 0
4.2 傅立叶变换的应用
注意到 ()cosat 1 [()eiat ()eiat ]
2
取傅立叶逆变换,得
u 2u
t
x 2
,
t 0, x R .
u x, 0 f x
解:由自变量的取值范围 ,对 x 进行傅立叶变换,设
数学物理方程课件 积分变换法
设F[ f1(x)] F1(), F[ f2 (x)] F2 (),
则F[ f1(x) f2 (x)] F1() F2 ()
(5)
其中,为常数,逆变换也成立,即
F-1[ F1() F2 ()] f1(x) f2 (x)
(6)
试证明Fourier正弦变换和Fourier余弦变换的公式分别为
Fs1[Fs ()]
f (x)
2
0 fs (x) sin xdx
Fc1[Fc ()]
f
(x)
2
0 fc (x) cos xdx
§4.1.1 Fourier变换法
证明:F () F[ f (x)] f (x)eixdx
i
2
0
Fs
(
)
ei
x
d
(欧拉公式)
即Fourier正弦变换的公式为
f (x) 2
0 Fs () cos xd
§4.1.1 Fourier变换法
例9:证明
x 0 1 x2
sin xdx
2
e
(
0)。
证明:本题直接积分不易计算,考虑到fs
1 l
l l
f (x) cos n
l
xdx, n 0,1, 2,...
bn
1 l
l l
f (x) sin n
l
xdx, n 1, 2,...
§4.1.1 Fourier变换法
二、Fourier变换
设f (x)在(-, )上满足
i)逐段光滑(可导);
高等应用数学方法
高等应用数学方法
1. 积分变换法:积分变换法是一种用于解决复杂微分方程的数学方法,它通过将原始微分方程转化为一系列积分方程来求解。
2. 广义矩阵反转法:广义矩阵反转法是一种用于求解线性方程组的数学方法,它利用矩阵的反转来求解线性方程组。
3. 广义逆矩阵法:广义逆矩阵法是一种用于求解线性方程组的数学方法,它利用矩阵的逆来求解线性方程组。
4. 拉格朗日乘子法:拉格朗日乘子法是一种用于求解非线性方程组的数学方法,它利用拉格朗日乘子来求解非线性方程组。
5. 拉格朗日方程法:拉格朗日方程法是一种用于求解最优化问题的数学方法,它利用拉格朗日方程来求解最优化问题。
6. 高斯消元法:高斯消元法是一种用于求解线性方程组的数学方法,它利用高斯消元法来求解线性方程组。
7. 广义逆矩阵法:广义逆矩阵法是一种用于求解复杂线性方程组的数学方法,它利用矩阵的逆来求解复杂线性方程组。
第四章积分变换法
即:由三角函数组成的函项级数成为三角级数。
三角函数系的正交性
(1)三角函数系
1,cos x,sin x,cos 2x,sin 2x, cos kx,sin kx,
( 2)正交 :
任意两个不同函数在[ , ]上的积分等于零。即
i)
cos kxdx 0,
sin kxdx 0,
16
ii)
sin kx cos nxdx 0.
3
特别是对于无界或半无界的定解问题,用积分变换来 求解,最合适不过了。(注明:无界或半无界的定解问题 也可以用第三章方法求解)
4
所谓积分变换,就是把某函数类A中的任意一个函数 f (t)
,经过某种可逆的积分方法(即为通过含参变量 的积分)
b
F( ) f (t)K(t, ) d t
a
变为另一函数类 B中的函数 F ( ), 这里 K (t, ) 是一个确
u(t) 4 (sin t 1 sin 3t 1 sin5t 1 sin7t )
3
Байду номын сангаас
5
7
( t , t 0)
由以上可以看到:一个比较复杂的周期函数可以看 作是许多不同频率的简谐函数的叠加
14
2 三角级数 三角函数系的正交性
三角级数
引例中的简谐振动函数
f (t ) A0 Ak sin(k t k )
傅立叶的两个最主要的贡献:
• “周期信号都可表示为谐波关 系的正弦信号的加权和”—— 傅里叶的第一个主要论点
• “非周期信号都可用正弦信号 的加权积分表示” ——傅里叶的第二个主要论点
10
(一) 周期函数的傅里叶展开 1.傅里叶级数的引进
在物理学中,我们已经知道最简单的波是谐波(正弦
数学物理方法3-4积分变换法
§3.4.1
第三章 偏微分方程的定解问题 第四节 积分变换法
直线上的初值问题
例3.4.1求解热传导 问题
dU(, t) 2 2 a U(, t), t 0 解:利用傅立 dt 叶变换的性质 U(, 0) (), t a22 a22t C () U(, t) e C F(, ) e d
思考 利用积分变换方法求解问题的好处是什么?
第三章 偏微分方程的定解问题 第四节 积分变换法
傅立叶变换的定义
U ( , t ) u ( x, t )e
j x
1 dx , u ( x , t ) 2
U ( , t )e j x d
傅立叶变换的性质 微分性 位移性 f ( n ) (x) ( j ) n F ( )
e
d d
1 2a
t
( )e
2 x
4 a 2t
d
第三章 偏微分方程的定解问题 第四节 积分变换法
§3.4.2
半无界直线上的问题
半无界区域上的热传导(扩散)问题 2 u 2 u 0 x , t 0 t a x 2 0, 例3.4.4 求解 t 0 u (0, t ) u0 , u ( x, 0) 0, 0 x 做代换 u ( x, t ) v( x, t ) u0 转化为直线上热传导方程 2 v v 2 对称延拓法(奇延拓) a , 0 x , t 0 2 x t x0 u0 , v(0, t ) 0, t0 ( x) u0 , x0 v( x, 0) u0 , 0 x 考虑与无界区域上 波传播问题的差别
数学物理方程行波法与积分变换
常见数学物理方程
波动方程
描述波动现象的数学模型,如声波、光波和水波 等。
热传导方程
描述热量传递过程的数学模型,如温度场的变化 和热传导等。
弹性力学方程
描述弹性物体变形的数学模型,如物体的应力和 应变等。
数学物理方程的解法
行波法
通过将方程转化为行波方程,利用行波的特性求解原 方程。
分离变量法
将多变量问题转化为单变量问题,通过求解单变量方 程得到原问题的解。
拉普拉斯变换
01
拉普拉斯变换的定 义
将一个时域函数转换为复平面上 的函数。
02
拉普拉斯变换的性 质
线性、时移、复频移、微分、积 分等。
03
拉普拉斯变换的应 用
控制系统分析、电路分析等领域。
积分变换的性质和应用
积分变换的性质
线性性质、时移性质、频移性质、微 分性质等。
积分变换的应用
求解偏微分方程、求解常微分方程、 求解积分方程等。
应用
一维波动方程的行波法广泛应用于求解一维波动问题,如弦振动、 波动传播等。
高维波动方程的行波法
方法
转化
应用
对于高维波动方程,行波法同样适用。 设解为多个行波的叠加形式,利用波 的传播性质和叠加原理,将高维波动 方程转化为多个一维或低维的常微分 方程或代数方程。
通过行波变换,将高维波动方程分解 为多个一维或低维的方程,简化求解 过程。
。
03
对于某些问题,可能需要复杂的积分变换和逆变换计
算。
行波法与积分变换的联系
行波法和积分变换都是求解数学物理方程的方法,它们之间存在一定的联 系。
在某些情况下,行波法可以通过适当的变量替换转化为积分变换的形式。
《数理方程》积分变换法解析
x2
x2
1 p2
dU dx
2x p
x2 p3
.
而 u |x1 cos y
变为
U
x,
p
|x1
1
p p2
,
解常微分方程得
U x, p
1 3 p3
x3
1 p
x2
p 1 p2
1 3 p3
1 p
.
取拉普拉斯逆变换,得
L(t n )
n! pn1 , n 0,1,
u
|x
0
f
t.
思考:需要对哪一个自变量进行哪一种积分变换?
对 t 进行拉普拉斯变换,设
u x,t U x, p, f t F p
于是方程变为
a2
d 2U x,
dx 2
p
pU
x,
p,
U x, p |x0 F p
这是二阶常微分方程的边值问题,它的通解为
根据傅里叶变换的微分性质,
方程转化为
dU ,
t
2U , t
dt
U , t |t0 F
于是 U ,t F e2t .
为了求出原方程的解,下面对 U ,t 关于 进行
傅立叶逆变换.
U ,t F e2t .
再由边值条件 U x, p |x0 F p 可知,C = F(p).
U
x,
p
F
pe
p a
x
.
为求出 u(x,t), 需要对 U(x,p) 进行拉普拉斯 逆变换。
积分变换法
dU (,t) a22U (,t) G(,t),
dt
它满足初值条件
U (, t) |t0 ().
(39) (40)
为了求解常微分方程初值问题(39)(40),记
19
例1 求解下列问题的解 ut a 2uxx f (x,t) ( x , t 0), (37)
u |t0 (x).
(38)
10
例3 求fˆ() e2t 的傅里叶逆变换,其中t 0.
解 由定义知
f (x) 1 fˆ ()eixd 1 e2t eix d
2
2
1 e2t (cosx i sin x)d,
2
1 e2t cos xd,
0
对 f (x) 求导,并利用一次分部积分得
df (x) x f (x) 0. dx 2t
( ) L1
s
1
2a 2
L1
s
1
2 a 2
G
(, s)
L[eat ] 1 sa
()ea22t
G(,t) ea22t
()ea22t t G(, )ea22 (t ) d . 0
(42)
为了求出问题(37)(38)的解,还需要对U (,t)
取傅氏逆变换。
22
例1 求解下列问题的解 ut a 2uxx f (x,t) ( x , t 0), (37)
t t0 t t0
证明 由拉氏变换的定义知
L[ f (t t0 )u(t t0 )]
0
f
(t
t0 )u(t
t0 )est dt
t0
f
(t
t0
)e st
dt
令 y t t0 , 则上式变为
“数学物理方程”课程教学大纲
“数学物理方程”课程教学大纲英文名称:Mathematical and physical equation课程编号:math2029学时:32 学分: 2适用对象:全校二年级本科生先修课程:高等数学,线性代数,复变函数,积分变换。
使用教材及参考书:申建中刘峰编,《数学物理方程》,2010年,西安交通大学出版社一、课程性质、目的和任务“数学物理方程”是高等学校工科本科有关专业的一门基础课。
本课程旨在使学生初步掌握数学物理方程的基本理论和基本方法,为学习有关后继课程和进一步扩大数学知识面而奠定必要的基础。
本课程的内容包括:弦振动方程、热传导方程和拉普拉斯方程等定解问题的提出,达朗贝尔法、分离变量法、贝塞尔函数与勒让德多项式的基本性质和应用。
二、教学基本要求本课程的内容按教学要求的不同,分为两个层次。
文中用黑体字排印的,属较高要求,必须使学生深入理解,牢固掌握,熟练应用。
其中,概念、理论用“理解”一词表述,方法、运算用“掌握”一词表述。
非黑体字排印的,也是教学中必不可少的,只是在要求是低于前者。
其中,概念、理论用“了解”一词表述,方法、运算用“会”或“了解”表述。
也是教学中必不可少的,属基本要求。
第一章数学建模与基本原理介绍了解三个典型方程(弦振动方程、热传导方程和拉普拉斯方程)的建立,了解定解条件的物理意义及三种定解问题(初值问题、边值问题和混合问题)的提法,了解偏微分方程的一些基本概念(解、阶、维数、线性与非线性、齐次与非齐交),理解线性问题的叠加定理。
第二章分离变量法掌握有界弦自由振动问题和有限长杆上热传导问题的分离变量解法,掌握圆域内拉普拉斯方程的狄利克雷问题的分离变量解法,会用固有函数法解非齐次方程的定解问题,会用辅助函数和叠加原理处理非齐次边值问题。
第三章贝塞尔函数了解贝塞尔(Bessel)方程的幂级数解法,掌握整数阶贝塞尔函数的一些性质(递推公式、零点、正交性),了解傅里叶——贝塞尔展开式,会用贝塞尔函数解有关的定解问题。
数学物理方程-第四章积分变换法
第四章 积分变换法积分变换法是求解偏微分方程的一种基本方法. 不仅如此,在自然科学和工程技术的许多领域也有着广泛应用. 本章介绍Fourier 变换在求解偏微分方程定解问题中的应用. 主要以一维热传导方程,一维波动方程及平面上的Laplace 方程为主. 对于高维情形,由于计算过程要复杂一些,故只做简单介绍,也不做过多要求.§4⋅1 热传导方程Cauchy 问题4.1.1 一维热传导方程Cauchy 问题 考虑如下问题2(,), , 0 (1. 1)(,0)(), (1. 2)t xx u a u f x t x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩ 下面利用Fourier 变换求解该定解问题.设0>β为常数,函数2x e β-的Fourier 变换为()224()x F e e ωββπωβ--=(1.3)为书写方便起见,引入记号ˆ()(())(),f F f x ωω=,如果f 为二元函数),(t x f , ),))(,((),(ˆt t x f F t fωω=表示对),(t x f 中的空间变量x 作Fourier 变换的像函数,此时t 作为参数对待.对(1.1)—(1.2)关于空间变量x 作Fourier 变换得22ˆ(,)ˆˆ(,)(,), 0ˆˆ(,0)() .du t a u t f t t dtuωωωωωϕω⎧+=>⎪⎨⎪=⎩ 上面是一阶线性常微分方程的初值问题,解之可得2222()0ˆˆˆ(,)()(,)t a t a t ut e f e d ωωτωϕωωττ---=+⎰ (1.4) 利用(1.3)得),)( 21(22224t eta F e ta x t a ωπω--=),)()(21()(4)(2222τωτπττω--=----t e t a F e t a x t a记2241Γ(,)()2x a tx t eu t a πt-=(1.5)其中)(t u 为单位阶跃函数. 则有22ˆ((,))()(,)a t e F x t t ωωω-=Γ=Γ22()ˆ((,))()(,)a t e F x t t ωττωωτ--=Γ-=Γ-利用上面结果将(1.4)改写为ˆˆˆˆˆ(,)()(,)(,)(,)tu t t f t d ωϕωωωτωττ=Γ+Γ-⎰ (1.6) 对(1.6)两边取Fourier 逆变换,并利用Fourier 变换卷积公式 ))(()))((ˆ)(ˆ(21211x f f x f fF *=-ωω 便得0(,)()Γ(,(,)*Γ(,)tu x t x x t)f x x t d ϕτττ=*+-⎰()(,)(,)(,)t x t d d f x t d ϕξξξτξτξτξ∞∞-∞-∞=Γ-+Γ--⎰⎰⎰2222()()4()401()(,)2 2( )x x t a t a td ed f ed at at ξξττϕξξξτξππτ----∞∞--∞-∞=+-⎰⎰⎰(1.7)(1.7)即为定解问题(1.1)—(1.2)的解.在),(t x u 的表达式(1.7)中,函数(;)x t Γ起着一个基本作用. 如果令0≡f ,)()(x x δϕ=,则有(,)()(,)(;).u x t x x t x t δ=*Γ=Γ因此,(;)x t Γ是如下问题的解20, , 0 (1. 8)(,0)(), . (1. 9)t xx u a u x t u x x x δ⎧-=-∞<<∞>⎨=-∞<<∞⎩而(,)x t ξΓ-和(,)x t ξτΓ--分别是下面两问题的解20, , 0 (1. 10)(,0)(), . (1. 11)t xx u a u x t u x x x δξ⎧-=-∞<<∞>⎨=--∞<<∞⎩ 2()(),, 0 (1. 12)(,0)0, . (1. 13)t xx u a u x t x t u x x δξδτ⎧-=---∞<<∞>⎨=-∞<<∞⎩ 由于知道了(;)x t Γ,就可直接写出(1.1)—(1.2)的解(1.7)式. 类似于求解线性方程组0=Ax ,其中A 为n m ⨯矩阵. 如果知道该齐次方程组的一个基解组,则方程的任一解可由基解组的线性组合表出. 因此,),(t x Γ的作用就相当于向量空间中的基,故称),(t x Γ为定解问题(1.1)—(1.2)的基本解(fundamental solution).基本解是线性微分方程的一个很重要的概念,不仅可以表示Cauchy 问题的解,也可用来构造Green 函数表示边值问题的解.基本解有明确的物理解释. 若在初始时刻0t =时在0x =处置放一单位点热源,则此单位点热源在x 轴上产生的温度分布便是(,)x t Γ. 类似地,若在初始时刻0t =时在x ξ=处置放一单位点热源,则此点热源在x 轴上产生的温度分布为(,)x t ξΓ-. 而将初始时刻0t =变为t τ=时,其温度分布就是(,)x t ξτΓ--.注1 在(1.1)—(1.2)解的表达式(1.7)中,如果将其中的第一项和第二项分别记为1(,)u x t 和2(,)u x t ,则1(,)u x t 是相应于(,)0f x t =时齐次方程的解,而2(,)u x t 是相应于0)(=x ϕ时非齐次方程的解.若记1(,)()*Γ(,)(,)u x t x x t M x t ϕϕ==,则由齐次化原理可知20(,)(,)tf u x t M x t d τττ=-⎰.另外,和1(,)u x t 表达式中的卷积形式类似,2(,)u x t 也可表示成某种卷积形式,请同学们试给出这一表示形式. 例1.1 求解如下定解问题20, , 0 (1.14)(,0)(), . (1.15)t xx x u a u bu cu x t u x x x ϕ⎧---=-∞<<∞>⎨=-∞<<∞⎩ 其中,,a b c 均为常数.解 对(1.14)-(1.15)关于x 作Fourier 变换得22ˆ(,)ˆˆˆ(,)(,)(,), 0ˆˆ(,0)()dut a u t bi u t cu t t dtu ωωωωωωωϕω⎧=-++>⎪⎨⎪=⎩解之可得22() ˆˆ(,)().a bi c tut e ωωωϕω---= (1.16)为了求函数22()a bi c teωω---的Fourier 逆变换,利用配方法将其改写为222222224()()42.b a c bi t a t a bi c taaeeeωωω-------=由于222241()(),2x a t a tF eea tωωπ--=利用Fourier 变换的位移性质得000ˆ(())()()()() ,i x F f x e F f f ωωωωωω=-=- 取022biaω=得222222()4221()().2x bibi ixa t a ta aF ee ea tωωπ---=故有2222224()42((,))()b a c bi t a t aaF g x t eeωω----=22(),a bi c teωω---=其中22222244421(,)2b a c x bx taa ta g x t eeea tπ----=22()4.2x bt cta tee a tπ+-=记22()41Γ(,)()2x bt c ta tex t e u t a πt+-=其中)(t u 为单位阶跃函数. 1(;)x t Γ即为定解问题(1.14)—(1.15)的基本解.将(1.16)改写为1ˆˆˆ(,)()(,) ,u t t ωϕωω=Γ.,求Fourier 逆变换得1(,)()Γ(,)u x t x x t ϕ=*1()(,)x t d ϕξξξ∞-∞=Γ-⎰ 22()4() .2x bt cta te ed a tξϕξξπ-+-∞-∞=⎰如果将(1.15)中的齐次方程改为非齐次方程 ,考虑如下定解问题2(,),, 0 (,0)0, . t xx x u a u bu cu f x t x t u x x ⎧=+++-∞<<∞>⎨=-∞<<∞⎩请同学们写出该定解问题的解.例1.2 求解如下定解问题20, , 0(,0)(), .t xx u a u x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩ 其中0, ()0, .A x x x x x ϕ>⎧=⎨<⎩解 由(1.7)可得该问题的解为22220()()441(,)(),2 2 x x a ta tx A u x t ed ed at at ξξϕξξξππ----∞∞-∞==⎰⎰对积分作变量代换 2x a tξα-=得 02222202020(,) [][]2x x a t x x a t x x a t Au x t e d Aed e d Ae d αααααπααππαπ---∞----∞--==+=+⎰⎰⎰⎰引入下面函数22()xx e d ααπ-Φ=⎰(1.17)该函数称为误差函数. 利用误差函数可得(,)()222x x A A u x t a t-=+Φ. 4.1.2* 二维热传导方程Cauchy 问题为加深对线性微分方程基本解的进一步理解,下面再求解二维热传导方程Cauchy 问题222()(,,), (,)R , 0 (1.18)(,,0)(,), (,)R . (1.19)t xx yy u a u u f x y t x y t u x y x y x y ϕ⎧-+=∈>⎪⎨=∈⎪⎩ 为求解(1.19)—(1.20),先求二维热传导方程的基本解,即如下定解问题的解222()0, (,)R , 0 (1.20)(,,0)()(), (,)R . (1.21)t xx yy u a u u x y t u x y x y x y δδ⎧-+=∈>⎪⎨=∈⎪⎩引入二元函数的Fourier 变换12()12()(,)(,)i x y F f f x y e dxdy ωωωω∞∞-+-∞-∞=⎰⎰和一元函数Fourier 变换的性质相对应,二元函数的Fourier 变换也有类似性质.对(1.20)-(1.21)关于空间变量作Fourier 变换得22ˆ(,)ˆ(,)0, 0ˆ(,0) 1.dut a u t t dtuωωωω⎧+=>⎪⎨⎪=⎩其中2221212(,) , ωωωωωω==+. 解之可得22222212ˆ(,)a ta t a t ut e e e ωωωω---==.故有2212222211222222222()1121221244421(,,)()(,)(2)11=2211=221=.(2)a t i x y a t ix a t iy xy a t a tx y a tu x y t F f eed d eed e e d eea ta t e a t ωωωωωωωωωωωπωωπππππ∞∞-+--∞-∞∞∞---∞-∞--+-==⎰⎰⎰⎰即(1.18)-(1.19)的基本解为222421Γ(,,)().(2)x y a tx y t eu t a πt +-=与(1.7)相对应,(1.20)—(1.21)的解为(,,)(,)*Γ(,,)(,,)*Γ(,,)tu x y t x y x y t f x y x y t d ϕτττ=+-⎰(,)(,,)x y t d d ϕξηξηξη∞∞-∞-∞=Γ--+⎰⎰(,,)(,,).t d f x y t d d τξητξητξη∞∞-∞-∞Γ---⎰⎰⎰作为练习,同学们试用Fourier 变换求解三维热传导方程Cauchy 问题. §4⋅2 波动方程Cauchy 问题4.2.1 一维波动方程Cauchy 问题考虑如下定解问题2(,), , 0 (2.1)(,0)(), (,0)(), . (2.2)tt xx t u a u f x t x t u x x u x x x ϕψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩20, , 0 (2.3)(,0)0, (,0)(), . (2.4)tt xx t u a u x t u x u x x x ψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 若记(2.3)—(2.4)的解为(,)(,)u x t M x t ψ=,则由叠加原理和齐次化原理可得(2.1)—(2.2)的解为0(,)(,)(,)(,)t f u x t M x t M x t M x t d tτϕψττ∂=++-∂⎰ (2.5)因此,只须求解定解问题(2.3)—(2.4).对(2.3)—(2.4)关于空间变量x 作Fourier 变换得2222ˆ(,)ˆ(,)0, 0ˆˆˆ(,0)0, (,0)().t d u t a u t t dt u uωωωωωψω⎧+=>⎪⎨⎪==⎩ 解之可得sin ˆˆ(,)() .a tut a ωωψωω= 记1, 2(;) 0, .x atax t x at ⎧<⎪Γ=⎨⎪≥⎩查Fourier 变换表或直接计算可得sin ˆ((;))()(,)a t F x t t a ωωωωΓ=Γ= 故有ˆˆˆ(,)()(,),ut t ωψωω=Γ 对上式取Fourier 逆变换并利用卷积公式得(,)()*Γ(,)u x t x x t ψ=()(,)x t d ψξξξ∞-∞=Γ-⎰1()2x atx at d aψξξ+-=⎰ . 利用(2.5)便得(2.1)—(2.2)的解为0(,)(,)(,)(,)t f u x t M x t M x t M x t d t τϕψττ∂=++-∂⎰11(())()22x at x at x at x at d d t a aϕξξψξξ++--∂=+∂⎰⎰ ()0()1(,)2tx a t x a t d f d a τττξτξ+---+⎰⎰[]11()()()22x at x atx at x at d a ϕϕψξξ+-=++-+⎰ ()0()1(,)2tx a t x a t d f d a τττξτξ+---+⎰⎰ (2.6)当0f ≡时,(2.6)称为一维波方程Cauchy 问题的达朗贝尔(D ’Alembert )公式.注1 在(2.4)中取()()x x ψδ=,则有(,)(;)u x t x t =Γ,即(;)x t Γ是如下定解问题20, , 0(,0)0, (,0)() .tt xx t u a u x t u x u x x x δ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 的解,称其为一维波动方程的基本解. 利用基本解(;)x t Γ,就可写出(2.1)—(2.2)的解(2.6)式. (;)x t Γ在(2.6)的表达式中也起到一个“基”的作用.4.2.2* 二维和三维波动方程Cauchy 问题下面,首先利用Fourier 变换求解三维波动方程Cauchy 问题,然后用降维法求出二维波动方程Cauchy 问题的解.考虑三维波动方程Cauchy 问题2333(,,,),(,,),0, (2.7)(,,,0)(,,),(,,), (2.8)(,,,0)(,,),(,,). (2.9)tt tu a u f x y z t x y z R t u x y z x y z x y z R u x y z x y z x y z R ϕψ⎧-∆=∈>⎪=∈⎨⎪=∈⎩为求解定解问题(2.7)—(2.9),先求出三维波动方程的基本解,即如下问题的解,23330, (,,), 0 (2.10)(,,,0)0, (,,) (2.11)(,,,0)()()(), (,,). tt t u a u x y z R t u x y z x y z R u x y z x y z x y z R δδδ-∆=∈>=∈=∈ (2.12)⎧⎪⎨⎪⎩记2222123123(,,) , ωωωωωωωω==++. 对定解问题(2.10)—(2.12)关于空间变量 作Fourier 变换得2222ˆ(,)ˆ(,)0, 0ˆˆ(,0)0, (,0) 1.t d u t a u t t dt uu ωωωωω⎧+=>⎪⎨⎪==⎩解之可得sin ||ˆ(,).||a tut a ωωω=故有123312331 ()1233()1233ˆ(,,,)()(,,,)1ˆ =(,)(2)sin 1 =(2)i x y z R i x y z R u x y z t F u x y z t u t e d d d a t e d d d a ωωωωωωωωωωπωωωωπω-++++=⎰⎰⎰⎰⎰⎰为计算上面积分,首先对上面积分作变量代换v A ωT T =,其中123(,,) , v v v v =A 为三阶正交矩阵. 选A 使得将(,,)x y z 变为(0,0,)r ,222 r x y z =++. 根据正交变换的保内积性可得,该变换将123 , x y z ωωωω++分别变为3 , v rv .故有33 1233sin 1(,,,)(2)i rv R a v t u x y z t e dv dv dv a v π=⎰⎰⎰,再利用球坐标变换123cos sin sin sin cos v v v ρθϕρθϕρϕ=⎧⎪=⎨⎪=⎩ 可得22 cos 3000cos 200 201sin (,,,)=sin (2) =sin (2) =sin() ()(2)i r i r i r i r a t u x y z t d d e d a i a t e d ar i a t e e d ar ππρϕπρϕρρρθρρϕϕπρρρπρρπ∞∞∞---⎰⎰⎰⎰⎰22sin() ()81( )()16i r i r i a t i a t i r ir i a t e e d ar e e e e d ar ρρρρρρρρπρπ∞--∞∞---∞=--=---⎰⎰. 注意到()(0)2()2()i i e d F e αραρωρπδωαπδα∞=-∞==-=⎰,221(,,,)( )()1612[()(())()()]161[()()]41().4 i a t i a t i r ir u x y z t e e e e d ar r at r at r at at r ar r at r at ar r at ar ρρρρρππδδδδπδδπδπ∞---∞=---=-⋅++-+----=--+=-⎰记1(,,,)()4x y z t r at arδπΓ=- (,,,)x y z t Γ即为三维波动方程的基本解.因此,当0==ϕf 时,(2.7)—(2.9)的解为33R R (,,,)(,,,) =(,,)(,,,)=(,,)(,,,)()=(,,). 4u x y z t M x y z t x y z x y z t x y z t d d d r at d d d arψψψξηζξηζξηζδψξηζξηζπ=*ΓΓ----⎰⎰⎰⎰⎰⎰其中222()()()r x y z ξηζ=-+-+-.对任一0t >, 记以点(,,)x y z 为心at 为半径的球面为(,,)r S x y z ,即3(,,){ (,,) }r S x y z R r at ξηζ=∈=. 将上面的积分化为累次积分并由δ函数的定义可得(,,)(,,)2(,,)2(,,,)(,,,)(,,)=()()4(,,)=4(,,) =41=(,,4rratS x y z r atS x y z S x y z u x y z t M x y z t r at ds dr ar ds ar ds a t a t ψψξηζδπψξηζπψξηζπψξηζπ∞==-⎰⎰⎰⎰⎰⎰⎰(,,)) . (2.13)at S x y z ds ⎰⎰最后,由叠加原理和齐次化原理便得(2.7)—(2.9)的解为22(,,)(,,)111(,,,)(,,)(,,)44at at S x y z S x y z u x y z t ds ds a t t a t ϕξηζψξηζππ⎛⎫∂=+ ⎪ ⎪∂⎝⎭⎰⎰⎰⎰211(,,,)4r f t d d d a r a ξηζξηζπΩ+-⎰⎰⎰ (2.14) 其中 (,,){ (,,) }a t B x y z r at ξηζΩ==<.(2.14)称为三维波动方程Cauchy 问题的克希霍夫(Kirchhoff )公式.利用Fourier 变换求二维波动方程的基本解比较难. 利用三维空间中已有的结果(2.13),下面用降维法求二维波动方程Cauchy 问题.考虑如下三维波动方程Cauchy 问题2330,(,,),0 (,,,0)0,(,,,0)(,),(,,) tt t u a u x y z R t u x y z u x y z x y x y z R ψ⎧-∆=∈>⎪⎨==∈⎪⎩(2.15)(2.16)对于定解问题(2.15)—(2.16),由于初始数据与z 无关,可推知解u 与z 也无关,故有zz u =0,即定解问题(2.15)-(2.16)其实是一个二维波动方程Cauchy 问题, 由(2.13)可得该问题的解为2(,,)2(,,)1(,,)(,)41=(,) (2.17)2at atS x y z S x y z u x y t ds a t ds a t ψξηπψξηπ+=⎰⎰⎰⎰其中22222(,,){(,,)|()()(),}atS x y z x y z a t z ξηζξηζξ+=-+-+-=≥. 对于上半球面(,,)atS x y z +直接计算得 2222221()() ()()ds d d atd d a t x y ζζξηξηξηξη∂∂=++∂∂=----将上式代入到(2.17)中便得222(,)(,,)(,,)1(,). (2.18)2at B x y u x y t x y t d d a a t rψξηξζπ=Γ=-⎰⎰其中22()()r x y ξη=---,(,){(,)|}at B x y r at ξη=<.和三维情形类似,由(2.18)可得二维波动方程Cauchy 问题2222(,,),(,),0 (2.19)(,,,0)(,),(,), (2.20)(,,,0)(,),(,). (2.21)tt tu a u f x y t x y R t u x y z x y x y R u x y z x y x y R ϕψ⎧-∆=∈>⎪=∈⎨⎪=∈⎩ 的解为222222(,)(,)1(,)1(,)(,,)22at at B x y B x y u x y t d d d d a t a a t r a t r ϕξηψξηξηξηππ∂=++∂--⎰⎰⎰⎰()222(,)1(,,)2()a t tB x y f d d d a a t r τξηττξηπτ---⎰⎰⎰(2.22) (2.22)称为二维波动方程Cauchy 问题的波以松(Poisson )公式.4.2.3 解的物理意义对一维波动方程Cauchy 问题,如果无外力作用,则解由D’Alembert 公式给出,即[]11(,)()()() .22x atx at u x t x at x at d aϕϕψξξ+-=++-+⎰ 将上式改写为(,)()() ,u x t f x at g x at =++-其中011()()() ,22x atf x at x at d a ϕψξξ++=++⎰11()()() .22x at g x at x at d aϕψξξ--=-+⎰ 记1(,)()u x t f x at =+,2(,)()u x t g x at =-,则12(,)(,)(,).u x t u x t u x t =+.首先考虑1(,)() .u x t f x at =+当0t =时1(,0)() .u x f x =在(,)x u 平面上画出函数()f x 的图形,则()f x at +的图形可通过()f x 的图形向左平移at 个单位长度而得. 随着t 的增加,()f x 的图形不断向左平移,移动速度为a ,故称1(,)u x t 为左传播波,a 为波速. 同样道理,2(,)()u x t g x at =-称为右传播波. D’Alembert 公式表明:弦线在t 时刻的振动是初始振动所产生的右传播波和左传播波的叠加.其次,从D’Alember t 公式还可看出:u 在(,)x t 的值(,)u x t 只与x 轴上区间[],x at x at -+上初始值有关,而与其它点的初始值无关. 这是由于波速为a ,在区间[],x at x at -+外的初始扰动在时刻t 还未传播到点x ,故称区间[],x at x at -+为点(,)x t 的依赖区间. 在(,)x t 平面上,过(,)x t 点分别作斜率为1a±的直线,两条直线在x 轴上所截得的区间便是[],x at x at -+(图2.1()a ).给定x 轴上的区间[]12,x x ,过点1(,0)x 作直线1x x at =+,过点2(,0)x 作直线2x x at =-,它们和x 轴构成了一个三角形区域(图2.1()b ).由于该区域内任一点的依赖区间都落在区间[]12,x x 内,因此,解在此三角形区域内的值完全由区间[]12,x x 上的初始值决定,而与此区间外的初始值无关,故称此三角形区域为区间[]12,x x 的决定区域. 同理,过点1(,0)x 作直线1x x at =-,过点2(,0)x 作直线2x x at =+,它们和x 轴构成一个梯形区域(图2.1()c ),该区域称为区间[]12,x x 的影响区域,它表示区间[]12,x x 上初始扰动对弦线振动的作用范围.t (x , t ) t t决定区域 影响区域x x x 0 x at - x a t + 0 1x 2x 0 1x 2x(a ) (b ) (c )图2.1由上面分析可得,波以常速a 沿两族直线x at c ±=向左﹑右两个方向传播,这是波动现象的一个基本特征. 直线x at c ±= 称为一维波动方程的特征线,它们在一维波动问题的研究中起着重要作用.当0f =时,对公式(2.14)和(2.22)进行分析,便可得到和上面类似的结论.对二维波动方程,一点(,,)x y t 的依赖区域是以(,)x y 为心,at 为半径的圆域;而对三维波动方程,一点(,,,)x y z t 的依赖区域是以(,,)x y z 为心,at 为半径的球面,而不是球形区域. 反映在波的传播过程中,平面波有前阵面而无后阵面,正像把一块石子扔在湖中,在湖面上激起层层浪花,这种现象称为波的弥漫现象;而空间波既有前阵面又有后阵面,正像人们听到声音,一会儿就消失了,这种现象称为空间波传播的无后效现象,此即Huygens 原理.§4⋅3 积分变换法举例在前二节中,利用Fourier 变换求出了热传导方程和波动方程Cauchy 问题的解. 下面再进一步举例,说明积分变换法在求解偏微分方程定解问题中的作用.例3.1 求解如下定解问题(,), , 0 (3.1)(,0)(), . (3.2)t x u au f x t x t u x x x ϕ+=-∞<<∞>⎧⎨=-∞<<∞⎩其中a 为实数.解 对(3.1)—(3.2)关于空间变量x 作Fourier 变换得ˆ(,)ˆˆ(,)(,), 0ˆˆ(,0)().du t ai u t f t t dtuωωωωωϕω⎧+=>⎪⎨⎪=⎩ 解之可得()0ˆˆˆ(,)()(,)tait ai t ut e f e d ωτωωϕωωττ---=+⎰ (3.3) 由于(())()ait F x at e ωδω--= ()((()))()ai t F x a t e τωδτω----=故(3.3)可表示为ˆˆˆˆˆ(,)()()()(,)(())()t u t x at f x a t d ωϕωδωωτδτωτ=-+--⎰对上式取Fourier 逆变换得0(,)()()(,)(()tu x t x x at f x x a t d ϕδτδττ=*-+*--⎰()()(,)(())t x at d d f x a t d ϕξδξξτξτδξτξ∞∞-∞-∞=--+---⎰⎰⎰()((),).t x at f x a t d ϕτττ=-+--⎰ 例3.2 求半平面上调和方程边值问题的有界解(,0)(), . (3.5)xx yy u x f x x ⎨=-∞<<∞⎩ 解 对(3.4)—(3.5)关于变量x 作Fourier 变换得222ˆ(,)ˆ()(,)0, 0ˆˆ(,0)().d u y i u y t dy u f ωωωωω⎧+=>⎪⎨⎪=⎩ 解之可得12ˆ(,)y y uy C e C e ωωω-=+ 由于u 有界,故20 .C =结合初始条件可得ˆˆ(,)() y u t f e ωωω-= (3.6) 直接求yeω-的Fourier 逆变换得11()()2y y ix F e x ee d ωωωωπ∞----∞=⎰1cos()y e x d ωωωπ∞-=⎰2201sin()cos()y x x y x e x y ωωωπ∞--=+ 221(,)yg x y x y π==+故(3.6)可表示为ˆˆˆ(,)() g(,)uy f y ωωω= 对上式取Fourier 逆变换得)))(,(*)((),(x y g f y x u ⋅⋅=() g(,)f x y d ξξξ∞-∞=-⎰221().()yf d x y ξξπξ∞-∞=-+⎰ 例3.3* 设有一单位长度均匀杆,侧面绝热,两端温度为零度.若初始温度为sin 2x π,求杆内的温度分布.解 设(,)u x t 为杆内温度分布,则u 满足如下定解问题(0,)(1,)0, 0 (3.8)(,0)sin 2, 0 1. (t xx u t u t t u x x x π==≥=≤≤ 3.9)⎪⎨⎪⎩对(3.7)—(3.9)关于时间变量t 作Laplace 变换,并记(,)u x t 的像函数为(,)u x s 可得222(,)(,)(,0)0(0,)(1,)0.d u x s su x s u x a dx u s u s ⎧--=⎪⎨⎪==⎩即2222(,)1(,)sin 2 (3.10) (0,)(1,)0 (3.11)d u x s s u x s x dx a a u s u s π⎧-=-⎪⎨⎪==⎩(3.10)是常系数二阶线性常微分方程,非齐次项为三角函数. 易得该方程 通解为1222sin 2(,)4s s x x aaxu x s C eC es a ππ-=+++利用边界条件(3.11)得10C =,20,C =故22sin 2(,)4xu x s s a ππ=+取Laplace 逆变换可得224(,)sin 2a tu x t e x ππ-=.例3.4* 求下面半无界弦振动问题有界的解2cos , 0, 0 (3.12) (,0)0, (,0)0, 0 (3.13)(0,)0, 0. tt xx t u a u t x t u x u x x u t t ρω-=>>==≥=≥ (3.14)⎧⎪⎨⎪⎩解 对(3.12)—(3.14)关于时间变量t 作Laplace 变换得222222(,)(,)(0,)0,d u x s s s u x s adx s u s u ρω⎧-=⎪+⎨⎪=⎩有界. 或者2222222(,)(1)(,)()(0,)0,d u x s s su x s dxa a s u s u ρω⎧--=⎪+⎨⎪=⎩有界. 解之可得1222(,)()s s x x aau x s C e C es s ρω-=+++由于u 有界,故20 .C =结合初始条件可得22(,)(1)()s x au x s es s ρω-=-+ (3.15)对(3.15)取Laplace 逆变换可得)()(),(221t s s t x u -⎪⎪⎭⎫ ⎝⎛+=ωρL )()(221t s s e x a s -⎪⎪⎪⎭⎫ ⎝⎛+--ωρL (3.16) 由于)()(221t s s -⎪⎪⎭⎫ ⎝⎛+ωρL =)()1(2221t s s s -⎪⎭⎫ ⎝⎛+-ωωρL =))(())(1(221212t s s t s --ωωρωρ+- L L 2222(1cos )sin 2tt ρρωωωω=-= (3.17) 利用Laplace 变换的延迟性质)()))(()((s f e s t u t f s τττ-=--L其中()u t 为阶跃函数. 取x aτ=得 )()(221t s s e x a s -⎪⎪⎪⎭⎫ ⎝⎛+-ωρL =)()()(221a x t u a x t s s ---⎪⎪⎭⎫ ⎝⎛+ωρL 22()2sin ()2xt x a u t aωρω-=- 22()2sin , 2 0, 0 . x t x a t a x t a ωρω⎧-⎪≥⎪=⎨⎪≤<⎪⎩(3.18)将(3.17)—(3.18)代入到(3.16)中便得222222sin sin () , 22(,)2sin , 0 . 2t x x t t a au x t t x t a ρωωωρωω⎧⎡⎤--≥⎪⎢⎥⎪⎣⎦=⎨⎪≤<⎪⎩ 注1 定解问题(3.7)—(3.9)也可用分离变量法求解. 一般而言,Laplace变换方法的求解过程比较繁琐,而分离变量法已成固定模式,求解过程相对简明.习 题 四1. 用Fourier 变换求解如下定解问题(1) 20, , 0, 2(,0)0, 2.t xx u a u x t A x u x x ⎧-=-∞<<∞>⎪>⎨⎧=⎨⎪<⎩⎩(2) 40, , 0, 1(,0) 0, 1t xx u u x t h x u x x -=-∞<<∞>⎧⎪⎧<⎨⎪=⎨⎪>⎪⎩⎩2*用Fourier 变换求解如下定解问题(1) 20, , 0 , 0(,0) 0, 0.t xx x u a u x t e x u x x -⎧-=-∞<<∞>⎪⎧>⎨=⎨⎪<⎩⎩(2) 2, , 0(,0)0, . t t xx u a u e x t u x x -⎧-=-∞<<∞>⎨=-∞<<∞⎩ 3. 用Fourier 变换求解如下定解问题(1) 2, , 0(,0)sin , .t t x u u xe x t u x x x -⎧+=-∞<<∞>⎨=-∞<<∞⎩(2) 23, , 0(,0)(), . t x u u u x t u x x x ϕ=+-∞<<∞>⎧⎨=-∞<<∞⎩4. 求解如下一维波动方程Cauchy 问题(1) sin , , 0(,0)0, (,0)0, . tt xx t u u t x x t u x u x x -=-∞<<∞>⎧⎨==-∞<<∞⎩(2) 22, , 0 1(,0)sin , (,0), . 1tt xx t u a u tx x t u x x u x x x ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪+⎩5*求解如下Cauchy 问题(1) 222()0, (,)R , 0(,,0), (,,0)=0, (,)R . tt xx yy t u a u u x y t u x y xy u x y x y ⎧-+=∈>⎪⎨=∈⎪⎩(2) 2222()0, (,)R , 0(,,0)0, (,,0)=, (,)R . tt xx yy t u a u u x y t u x y u x y x y x y ⎧-+=∈>⎪⎨=∈⎪⎩ (3) 2323()0, (,,)R , 0 (,,,0)0, (,,,0)=, (,,)R .tt xx yy zzt u a u u u x y z t u x y z u x y z x z x y z ⎧-++=∈>⎪⎨=∈⎪⎩6. 由三维波动方程Cauchy 问题解的公式,利用降维法求解如下问题20, , 0 (,0)0, (,0)(), .tt xx t u a u x t u x u x x x ψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩7. 考虑如下定解问题20, , 0(,0)(), (,0)(), .tt xx t u a u x t u x x u x x x ϕψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 设()x ϕ和()x ψ为直线R 上奇(偶,周期为T 的)函数,证明该问题的解(,)u x t 关于变量x 也是奇(偶,周期为T 的)函数. 对于一维热传导方程Cauchy 问题,类似结果是否成立?8*设()x ϕ和()x ψ在{0}x x ≥二阶连续可导,(0)(0)0ϕψ==,求解如下波动方程半无界问题20, 0, 0(0,)0, 0(,0)(), (,0)(), 0 . tt xx t u a u x t u t t u x x u x x x ϕψ⎧-=<<∞>⎪=≥⎨⎪==<<∞⎩如将该问题的边界条件换为 (0,)0, 0x u t t =≥,如何求解相应的定解问题?9.考虑如下定解问题000, , 0(), 0, .tt xx t t t u u x t u x u x ϕ==-=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩其中初始波形为如下锯齿波1, 12()3, 230, .x x x x x ϕ-<<⎧⎪=-<<⎨⎪⎩其它(1)分别画出1,2t =时刻的(,)u x t 的波形图.(2)如果将初始位移换为1()()()x x x ϕϕϕ=--,分别画出1,2t =时刻的(,)u x t 的波形图.10. 考虑如下定解问题030, , 00, (), . tt xx t t t u u x t u u x x ψ==-=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩其中2e , 13()0, .x x x ψ⎧<<⎪=⎨⎪⎩其它 试找出(,)u x t 恒为零的区域,又弦线上10x =-的点在那个时刻开始振动. 11. 考虑如下定解问题2200()0, (,), 00, (,), (,).tt xx yy t t t u u u x y R t u u x y x y R ψ==⎧-+=∈>⎪⎨==∈⎪⎩ 其中, (,)(,)0, .x y x y ψ∈Ω⎧=⎨⎩正值其它 若区域Ω为正方形{ (,) 1 1 , 1 1 }x y x y -<<-<<,试指出(,,10)u x y 恒为零的区域.12. 考虑如下定解问题3300()0, (,,), 00, (,,), (,,).tt xx yy zz t t t u u u u x y z R t u u x y z x y z R ψ==⎧-++=∈>⎪⎨==∈⎪⎩ 若(,,)x y z ψ除在球形域222{ (,,) (1) 1 }x y z x y z -++≤取正值外其它恒为零,试指出(,,,10)u x y z 恒为零的区域.13*求解下面定解问题21200, , 0(,0),0, .tt xx x t t u u u x t u x e u x -=-+=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩ 14*考虑下面定解问题20, , 0 (,0)cos , . t xx u a u x t u x x x ⎧-=-∞<<∞>⎨=-∞<<∞⎩求出该定解问题解的有限表达形式.[利用结果2240cos ,0]4b ax ae bxdx ea aπ-∞-=>⎰.15*考虑下面定解问题230, , 0 (,0), . t xx u a u x t u x x x ⎧-=-∞<<∞>⎪⎨=-∞<<∞⎪⎩求出该问题解的有限表达形式.16*利用误差函数求解下面定解问题20, , 0 (,0)(), . t xx u a u x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩其中, 0(), 0.A x xB x ϕ>⎧=⎨<⎩。
数学物理方程第四章 积分变换法(课堂PPT)
❖ 傅里叶变换建立R了信号时域与频域之间的关系,
频率是信号的物理本质之一。
6
❖ 设f(x)为[-π,π]上的有限信号,则f(x)的傅 里叶变换可简化为:
fˆ ( ) π f (x)eix dx π
❖ 对于只在有限区间,例如在上有定义的函数,可 采取延拓的方法,使其成为某种周期函数,而在 上,。然后再对作傅里叶级数展开,其级数和在 区间上代表f(x),由于f(x)在x=0和x=l无定义,因此 可以有无数种延拓方式,因而有无数种展开式, 它们在上均代表.有时,对函数在边界(区间的 端点)上的行为提出限制,即满足一定的边界条 件,这常常就决定了如何延拓。
第四章 积分变换法 傅立叶变换与拉普拉斯变换
数学物理方程
1
1777年以前,人们普遍采用多项式函数P(x)来对信 N 1
号f(x)进行表征:f (x) P(x) anxn。 n0 1777年,数学家Euler在研究天文学时发现某些函
数可以通过余弦函数之和来表达。1807年,法国科学
家傅里叶进一步提出周期为2π的函数f(x)可以表示为
( x ,t 0)
U ' (t; k) k 2a2U (t; k) F(t; k) U (t; k) |t0 0
其中 U (t; k) 为u(x,t)的傅里叶变换。为求解这个非齐次
e 常微分方程,用 k2a2t 遍乘方程各项 18
d [U (t; k)ek2a2t ] F (t; k)ek2a2t dt
19
❖ 交换积分次序
u(x,t) t
1
= 0
f ( , )[2
e e dk] k2a2 (t ) ik (x ) d d
引用积分公式
e2k2 ek dk =
数理方程参考答案4第四章 积分变换法
若 在 点连续,则
1
定义
设函数 f ( x) 在 (−∞, +∞) 上的任意有限区间上满足狄利克雷条件,在 (−∞, +∞) 上绝
对可积,则称广义积分
为
的傅里叶变换,或者称为 定义 称
的像函数。通常记为
,或
。
为
的傅里叶逆变换,或者称为 傅里叶变换及其逆变换的基本性质
的像原函数。记为
.
性质 1(线性性质) 傅里叶变换及其逆变换都是线性变换,即
其中 , 是任意常数。 性质 2(相似性质) 对于任意实常数 ,有 . 性质 3(位移性质)对于任意实常数 ,有 , 性质 4(微分性质)设 , 的傅里叶变换存在,则 . 一般地,若 , ,…, 的傅里叶变换存在,则 . 性质 5(乘多项式性质)设 的傅里叶变换存在,则
2
.
. 性质 6 (积分性质) . 性质 7 (对称性质) . 定义 于所有的 设函数 和 是 上定义的函数。 如果广义积分 对
2 ∂ 2u 2 ∂ u a − = 0 (−∞ < x < +∞, t > 0), ∂t 2 ∂x 2 ∂u u| ψ ( x). ( x), = = t =0 ϕ ∂t t =0
的解为
二维拉普拉斯方程的边值问题
∂ 2u ∂ 2u = 0 ( −∞ < x < +∞, y > 0), ∂x 2 + ∂ y2 u | = f ( x ), x =0 u = 0. |xlim |→+∞ 的解为
2
s2
例3 解
求函数 F ( p ) = 因为
p 的拉普拉斯逆变换 p − 2 p +5
数学积分变换法
1 a
F
p a
,
a 0.
6) 卷积性质 定义
f
g
x
x
0
f
x
t
g
t dt
则 L f g L f Lg
例 设 y yt 求解常微分方程的初值问题:
y''2 y'3y et y |t0 0, y'|t0 1 解 对 t 进行拉普拉斯变换, 设 yt Fp, 则
et 1 p 1
y' pFp y0 pF( p)
可解得
U ,t e2t
由于
F 1[e2t ]
1
x2
e 4t
2 t
即
F
1
x2
e 4t
e2t
2 t
则
U ,t F
1
x2
e 4t F[ ]F
1
x2 e 4t
2 t
2 t
u x,t F 1 U(,t)
F 1 F ( )F
1
x2 e 4t
)e1 4(t )
x2
de4( t
)
d
2 0 t 2 (t )
*
1
x2
e 4t
2 t
t
f ( x, )*
1
x2
e d 4(t )
0
2 (t )
傅立叶逆变换是一种把分析运算化为代数运算的 有效方法,但
1.傅立叶变换要求原象函数在R上绝对可积.,大 部分函数不能作傅立叶变换
.
而 u x,0 x 0
解: 则
作关于 x
F
ux,t
的U 2傅1,立tt叶e变x42t u换x。,et设e2t i
数学物理方程与特殊函数积分变换法
(3)求拉式逆变换
F (s)
2s 3 s2 9
练习
(1)求傅式变换
f
(t)
sin t
t
0
t
其他
解 F f (t) f (t)e jtdt F()
sin te jtdt
2i sin t sintdt
0
2sin 2 1 i
练习
(2)求傅式变换 f (t) e jw0t u(t)
例1 解定解问题
u a 2 2u ,
t
x 2
u(x,0) (x),
x ,t 0 x
解:利用傅立 叶变换的性质
dU (,t) a 2 2U (,t), t 0
dt
U (,0) 1,
U (,t) Cea22t
C 1
U (,t) ea22t
1
e e x2 2 2
0
L f1(t) f2(t) L f1(t) L f2(t)
L 1 F1(s) F2 (s) f1(t) f2 (t)
常见结果
L
ekt
s
1
k
L cos
kt
s2
s
k2
L
sin
kt
s2
k
k2
L 1 1
s
练习
(1)求傅式变换
f
(t)
sin t
t
0
t
其他
(2)求傅式变换 f (t) e jw0t u(t)
2 2
2
u
2
1
x2
e e 4 2t
a2 2t
2a t
u(x,t)
1
x2
e 4 2t
第四章.积分变换法---求解偏微分方程
记作:F [ f ( x)] = f (k ) ,即
F [ f ( x)] = f (k ) = ∫
f f(x): (k ) 的傅里叶逆变换
∞ −∞
f ( x) e −ikx dx
记作: f ( x) = F −1[ f (k )] ,即
1 F [ f (k )] = f ( x) = 2π
−1
∫
∞
9
可以证明: 如果定义在 (−∞, ∞) 的函数在任一有限区间上满足 狄利克莱条件,且绝对可积( ∫ | f ( x) |dx 有界),则在
−∞ ∞
f(x)的连续点处,傅里叶积分存在:
1 f ( x) = 2π ⎡∞ ⎤ ikx −ikξ ∫∞ ⎢−∫∞ f (ξ )e dξ ⎥ e dk − ⎣ ⎦
——对于发生了任意位移x 0 的函数,其傅里叶变换 − ikx 等于 f(x)的傅里叶变换乘以一相位因子 e 0 证明:由定义:
F [ f ( x − x0 )] = ∫ f ( x − x0 ) e −ikx dx
−∞ u = x − x0 ∞
=
∫
∞
−∞
f (u ) e −ik (u + x0 ) du
频率域 波矢域
e − ikx
↔Leabharlann 412.1 傅里叶变换 一、傅里叶级数和复数形式的傅里叶级数 一个以2l为周期的函数f(x),若在区间[-l, l]满足 狄利克莱条件:(1)连续或只有有限个第一类间断 点;(2)只有有限个极值点,则 f(x) 在[-l, l]上可展开 为傅里叶级数
a0 ∞ nπ x nπ x + bn sin ) f ( x) = + ∑ (an cos 2 n =1 l l
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
0
u
t
a
f
t
a
e
st
dt
0 f t aestdt 0
令 v t a ,t v a
dt dv
f
v esvadv
0
esa f t est dt 0
esa L f t
F s eas
结论
easF s Lut a f t a
表示Fs乘以eas后,相当于f t
在 t 軸向右平移了a距離。
f
(x)
a0 2
ak
k 1
cos kx
而当信号具有反对称性(奇)特征时,ak=0,
f (x)
a0 2
bk sin kx
k 1
❖ 在研究热传导方程的过程中,为了简化原问题, 傅里叶建议将热导方程从时间域变换到频率域,
为此他提出了著名的傅里叶变换的概念。信号
f(x)的傅里叶变换定义为:
❖
fˆ ( ) f (x)eix dx,i 1
若函数 f (x) 以 f (x 2l) f (x) 为周期,即
则可取三角函数族
1,cos x,cos 2 sin xl,sin 2 xl
x, ,
… …
cos sin
n
l
,x …
n x, …
l
l
l
作为基本函数族,将f (x) 展开为级数
f (x)
= a0 + (an
n1
cosn
l
x
+bncos
f
(t
)]
1
1 e sT
T f (t)estdt
0
卷积与卷积定理
t
卷积定义: f1(t) f2 (t) 0 f1( ) f2 (t )d
卷积定理: l[ f1(t) f2 (t)] F1(S )F2 (S )
注意到前面所给出的约定,即函数 f (t) 等价于函数
f (t)u(t),因此,这里所给出的卷积实际上与Fourier
2 j j
其中,F(s) 称为函数 f (t) 的像函数,f (t)称为 F(s) 的像原函数.
注1:函数f (t)的Laplace变换就是函数 f (t)u(t)et 的Fourier变换.
注2:由于Laplace变换只用到了函数f (t)在 t 0
的部分,为方便起见,在Laplace变换中所提
其中 F(ω,t)=
1
2
u(x,t)[eix ]*dx
❖ 1 用傅里叶级数法解决有界细杆的热传导问题
ut a2uxx 0
第二类齐次边界条件下的本征函数:cos n x
l (0,1,2,…),
u(x,t)=
n0
Tn
(t
)co
s
n
l
x
把这个级数代入泛定方程,
[Tn' (t)
n0
n2 2a2
其中 U (t; k) 为u(x,t)的傅里叶变换。为求解这个非齐次
e 常微分方程,用 k2a2t 遍乘方程各项
d [U (t; k)ek2a2t ] F (t; k)ek2a2t dt
❖ 对t积分一次,计及零初始值,
U (t; k)
e k 2a2t
=
t F ( ; k )ek2a2 d
0
= t f ( , )eik ek2a2 (t )d d 0
(a,b是常数)
线形性质: l[af (t) bg(t)] aF(s) bG(s)
相似性质: l[af (t)] 1 F ( s ) aa
延迟性质 : l[ f (t )u(t )] es F (s)
微分性质: l[ f (n) (t)] snF (s) sn1F (0) sn2F(0) L F (n1) (0)
2
e 4 2
u(x,t) t
=
f ( , )[
1
e ]d d
( x )2 4a2 (t
)
0
2a (t )
4.4 Laplace变换的定义和基本性质
❖ Laplace变换应用范围: Laplace变换方法广泛应用于求解非稳态 热传导问题,将对时间的偏导数消去。
❖ Laplace变换方法简单,但对变换后得到 的解进行反变换则相当复杂。
变换中的卷积是一致的.
❖阶梯函数的Laplace变换
u
t a
0,t 1,t
a a
a
0
Lut a
ut a est dt
0
0
a
1
e
st
dt
1 e st s
a
0
1
s e as
1 e as s
❖Laplace变换的移位特性
若a>0,L[f(t)]=F(s) 則L[u(t-a)f(t-a)]
n
l
x
)
an
1
nl
l f ( ) cos n d
l
l
bn
1 l
l l
f ( )sin n
l
d
其中
n
2 1
(n 0) (n 0)
周期函数f(x)可以理解为由正弦波(含余弦与正 弦函数)叠加而成,其中an,bn为叠加的权值,表 示信号在不同频率时刻的谱幅值大小。
显然,当信号具有对称性(偶)特征时,bk=0,
g(x)=
a0 +
(anco s
n1
n
l
x
bn
sin
n
l
x
)
❖ 在l→∞时的极限形式就是所要寻找的非 周期函数f(x)的傅里叶展开。
f(x)= 0 A() cos xd 0 B()sin xd
其中
1
A(ω)= f(ξ)sinωξdξ
B(ω)=
1
f(ξ)cosωξdξ
复数形式的傅里叶积分
f(x)= F(ω) ei x dω
零初始条件
T0 (0) 0
1 l
l
( )d
o
Tn (0)
n
2 l
l
( ) cos
n
d
o
l
❖Tn(t)的常微分方程在初始条件下的解:
Tn
(t)=
e [
n2 2 l2
a2
t
fn
n2 2a2
(t)e l2
t
dt
n
fn (t)dt]
u(x,t)=
{
e [
n2 2a2 l2
t
n0
fn
n2 2a2
l 1[F (n) (s)] (1)n t n f (t)
❖Laplace变换的性质
积分性质:
l
t 0
f
(t)dt
1 F(s) s
l
1
s
F
(s)ds
f (t) t
周期函数的像函数性质:设 f (t) 是 [0 , )内以T为
周期的函数,且 f (t) 在一个周期内逐段光滑,则
.
l[
(t)e l2
t
dt
n
fn
(t
)dt
]
}cos
n
l
x
4.3无界空间的有源导热问题
❖ 1.一维无源导热问题和基本解 ❖ 2. 一维热传导问题 ❖ 3.一维有源导热问题。
❖傅里叶变换法求解无界细杆的热传 导问题
ut a2uxx f (x,t) u |t0 0
( x ,t 0)
U ' (t; k) k 2a2U (t; k) F(t; k) U (t; k) |t0 0
到的函数一般均约定在 t 0的部分为零.
换句话说,函数 f (t) 等价于函数 f (t)u(t).
注3:像函数F (s)通常仅在复平面s上的某个区域内 存在,称此区域为存在域,它一般是一个右半平 面.当函数 f (t)只要不比某个指数函数增长得快时, 则它的Laplace变换一定存在,因此我们所接触 到的绝大多数函数的Laplace变换都是存在的.在 进行Laplace变换时,常常略去存在域.
❖ 如要求 f (0) f (l) 0
这时应延拓为奇的周期函数,因为
sin
n
l
x│x
0
=0,
sin n
l
x∣x
l
=0;
如要求 f ' (0) f ' (l) 0
这时应延拓为偶的周期函数,因为余弦级 数的和的导数在 x 0 和 xl 为零
❖ 对于函数u(x,t),-l<x<l,t≥0,展开为傅里叶级 数时,可将t视为参数,仅关于x展开为傅 里叶级数
u(x,t)=a0
(t)+
n1
(an
(t
)co
s
n
l
x
bn
(t
)
sin
n
l
x
)
其中展开系数不是常数,而是关于t的函数,
1
an (t) nl
l u( ,t) cos n d
l
l
1
bn (t) l
l u( ,t) sin n
l
l
d
4.2 傅里叶变换
❖ 一般说来,定义在区间(-∞<x<∞)上的函数 f(x)是非周期的,不能展开为傅里叶级数。为 了研究这样的函数的傅里叶展开问题,可试 将非周期函数f(x)看作是某个周期函数g(x)于 周期2l→∞时的极限情形。这样,g(x)的傅里 叶级数展开式
§1 Laplace变换的定义、性质
❖ Laplace变换所考虑的对象通常是定义在 [0 , )上的 实值函数 f (t)
Laplace(正)变换: