数学物理方程-第3章-
数学物理方程(谷超豪)-第三、四章 课后习题答案
第三章调和方程§1建立方程定解条件1．设)(),,,(21r f x x x u n = )(221n x x r ++=是n 维调和函数（即满足方程022212=∂∂++∂∂nx ux u），试证明221)(-+=n rc c r f )2(≠n rInc c r f 1)(21+=)2(=n 其中21,c c 为常数。
证：)(r f u =,rx r f x rr f x u i i i ⋅=∂∂⋅=∂∂)()(''32''22"22)(1)()(r x r f r r f rx r f x ui i i ⋅-⋅+⋅=∂∂312''212"122)()()(rx r f r nr f rx r f x uni i ni i ni i∑∑∑===⋅-⋅+⋅=∂∂)(1)('"r f rn r f -+=即方程0=∆u 化为0)(1)('"=-+r f rn r f rn r f r f 1)()('"--=所以)1(1')(--=n r A r f 若2≠n ，积分得1212)(c r n A r f n ++-=+-即2≠n ，则221)(-+=n r c c r f 若2=n ，则rA r f 1')(=故Inr A c r f 11)(+=即2=n ，则rInc c r f 1)(21+=2．证明拉普拉斯算子在球面坐标),,(ϕθr 下，可以写成sin 1)(sin sin 1(12222222=∂∂⋅+∂∂∂∂⋅+∂∂∂∂⋅=∆ϕθθθθθur u r r u r r r u 证：球坐标),,(ϕθr 与直角坐标),,(z y x 的关系：ϕθcos sin r x =，ϕθsin sin r y =，θcos r z =（1）222222z u yu xu u ∂∂+∂∂+∂∂=∆为作变量的置换，首先令θρsin r =，则变换（1）可分作两步进行ϕρcos =x ，ϕρsin =y （2）θρsin r =，θcos r z =（3）由（2）⎪⎪⎭⎪⎪⎬⎫∂∂+-∂∂=∂∂∂∂+∂∂=∂∂)cos ()sin (sin cos ϕρϕρϕϕϕρy ux u u y u x u u 由此解出⎪⎭⎪⎪⎬⎫⋅∂∂+∂∂=∂∂⋅∂∂-∂∂=∂∂ρϕϕϕρρϕϕϕρcos sin sin cos u u y u u u x u （4）再微分一次，并利用以上关系，得)sin cos (22ρϕϕϕρ⋅∂∂-∂∂∂∂=∂∂u u x xu)sin cos (sin )sin cos (cos ρϕϕϕρϕρϕρϕϕϕρρϕ⋅∂∂-∂∂∂∂⋅-⋅∂∂-∂∂∂∂=u u u u +∂∂⋅+∂∂∂⋅-∂∂=22222222sin cos sin 2cos ϕρϕϕρρϕϕρϕuu u ρρϕϕρϕϕ∂∂⋅+∂∂⋅+u u 22sin cos sin 2cos sin (22ρϕϕϕρ⋅∂∂+∂∂∂∂=∂∂u u y yu)cos sin (cos )cos sin (sin ρϕϕϕρϕρϕρϕϕϕρρϕ⋅∂∂+∂∂∂∂++⋅∂∂+∂∂∂∂=u u u u ρρϕϕρϕϕϕρϕϕρρϕϕρ∂∂⋅+∂∂⋅--∂∂⋅+∂∂∂+∂∂=u u uu u2222222222cos cos sin 2cos cos sin 2sin 所以ρρϕρρ∂∂⋅+∂∂⋅+∂∂=∂∂+∂∂uu u yu xu 11222222222（5）ρρϕρρ∂∂⋅+∂∂⋅+∂∂+∂∂=∂∂+∂∂+∂∂uuz uu z u y u x u112222222222222再用（3）式，变换2222zu u ∂∂+∂∂ρ。
数学物理方程第三章_行波法和积分变换法
[x − at , x + at ] 上的值,而与其他点上的初始条件无关,这个区间称为点 (x, t ) 的依赖区间,
它是过 ( x, t ) 点分别作斜率为 ±
1 的直线与 x 轴相交所截得的区间,如图 3-2 所示. a
(x,t0)
y
x O x-at0 x+at0
图 3-1
初 始 时 刻 t = 0 时 , 取 x 轴 上 的 一 个 区 间 [x1 , x 2 ] , 过 点 x1 作 斜 率 为
同理可得
2 ∂ 2u ∂ 2u ∂ 2u ⎤ 2⎡∂ u = + a + 2 ⎢ 2 ∂ξ∂η ∂η 2 ⎥ ∂t 2 ⎣ ∂ξ ⎦
将其代入式(3.1.1),得
∂ 2u =0 ∂ξ∂η
对 ξ 积分,得
∂u = f (η ) ∂η
对此式再关于η 积分,得
u = ∫ f (η )dη + f1 (ξ ) = f1 (ξ ) + f 2 (η )
第三章 行波法与积分变换法 本章我们介绍两个常用的解题方法：行波法和积分变换法。行波法只用于求解无界区 域上的波动方程定解问题， 积分变换法不受方程类型的限制， 一般应用于无界区域的定界问 题，有时也应用于有界域的定解问题.
3.1 达朗贝尔公式及波的传播 在求解常微分方程的特解时,一般先求出方程的通解,然后利用所给的定解条件去解出 通解中含有的任意常数,最后得到了满足所给条件的特解.这个想法能否推广到求解偏微分方 程的过程中呢?一般情况下,随着自变量个数的增加,偏微分方程的通解非常难求,并且偏微分 方程的通解一般都含有任意函数,这种任意函数很难由定解条件确定为具体的函数.所以在求 解数学物理方程时,主要采用通过分析各类具体的定解问题,直接求出符合定解条件的特解的 方法.但事情没有绝对的,在有些情况下,我们可以先求出含任意函数的通解,然后根据定解条 件确定出符合要求的特解.本节我们研究一维波动方程的求解,就采用这种方式. 3.1.1 达朗贝尔公式 如果我们所考察的弦无限长,或者我们只研究弦振动刚一开始的阶段,且距弦的边界较远 的一段,此时可以认为弦的边界,对此端振动的弦不产生影响.这样,定解问题就归结为如下形 式
数学物理方程-3
其中ϕ(x, y, z) 和 ψ (x, y, z) 均为已知函数。
u
3-3 高维波动方程的初值问题
平均值法：不考虑函数 平均值法：不考虑函数 u(x, y, z, t) 本 身，而是研究u(x, y, z, t)在以点 M(x, y, z) 为球心，以r 为球心，以r为半径的球面上的平 均值 u ，当暂时选定 M(x, y, z) 后， u 就是关于r 就是关于r，t的函数。当我们很方 便地求出 u (r, t) 后，令 r →0 则 u(r, t) →u(x, y, z, t) ，问题就得到了 解决。
第3章 行波法与积分变换法
原柯西问题的通解为 u = f1 (x + at) + f2 (x − at) 初始条件代入其中，有 ϕ(x) = f1 (x) + f 2 (x) ′ ψ (x) = af1′(x) − af 2 (x) 无界弦振动的柯西问题的解（达朗贝尔解 无界弦振动的柯西问题的解（达朗贝尔解 ） 1 1 x+at 为： u(x, t) = [ϕ(x + at) +ϕ(x − at)] + ∫ ψ (ξ )dξ
3-2 延拓法求解半无限长振动问题
延拓后的定解问题：
2 ∂2v 2 ∂ v + F(x, t) (−∞ < x < +∞, t > 0) 2 =a 2 ∂x ∂t ∂v v(x,0) = Φ(x), |t=0 = Ψ(x) ∂t v(0, t) = 0
x >0 ϕ(x), Φ(x) = −ϕ(−x), x < 0
x >0 ψ (x), Ψ(x) = −ψ (−x), x < 0
x >0 f (x, t), F(x, t) = − f (−x, t), x < 0
数学物理方程课后参考答案第三章
解：令
又 故取 则 满足调和方程
即
代入原定解问题，得 满足
用分离变量法零解 ，得
.
所以
再由另一对边值得
所以 .
得
最后得
8.举例与说明在二维调和方程的狄利克莱外问题，如对解 不加在无穷远处为有界的限制，那末定解问题的解以不是唯一的。
是区域 中的调和函数（无穷远点除外）.
如果区域 为球面K以外的无界区域，则函数u 在 中除去原点O外是调和的，函数 称为函数 的凯尔文（Kelvin）变换。
证明：只需证明 满足 。
=
=
代入 的表达式，有
=
=
若u在包含原点O的有界区域内处处式调和的即 ,则除无穷远点（O的反演点）外， 即除 点外v是调和的。若u在无界域 上是调和的，则除去O点外，v也是调和的。证毕。
且矩阵( )是正定的，即
由于矩阵( )是非正定的，故 可以写成 的线性齐次式的平方和，即
=
所以
于是
因此在 点
与 在 点满足方程是矛盾的，故 不能在 内部达到正的最大值。
7．证明第6题中讨论的椭圆形方程第一边值问题的唯一性与稳定性。
证：唯一性。只须证明方程在齐次边值条件只的零解。
设 在 内满足方程，在 边界 上 。因 在 上连续，故 是有界的，
第三章调和方程
§1建立方程定解条件
1．设 是n维调和函数（即满足方程
），试证明
其中 为常数。
证： ,
即方程 化为
所以
若 ，积分得
即 ，则
若 ，则 故
即 ，则
2．证明拉普拉斯算子在球面坐标 下，可以写成
数学物理方程(谷超豪)第三章调和方程习题解答
∆u
=
1 r2
⋅
∂ ∂r
(r 2
∂u ) ∂r
+
r2
1 sin θ
⋅
∂ ∂θ
(sin θ
∂u ∂θ
)
+
r2
1 sin
2
θ
⋅
∂2u ∂ϕ 2
=0
证：球坐标 (r,θ ,ϕ) 与直角坐标 (x, y, z) 的关系：
x = r sinθ cosϕ ， y = r sin θ sin ϕ ， z = r cosθ
f
(r)
=
−
A1 n+
2
r −n+2
+
c1
即 n ≠ 2 ，则
f
(r)
=
c1
+
c2 r n−2
若 n = 2 ，则 即 n = 2 ，则
f ' (r) = A1 故 f (r) = c1 + A1Inr r
f (r) = c1 + c2 In 1 r
2． 证明拉普拉斯算子在球面坐标 (r,θ ,ϕ) 下，可以写成
⋅
∂u ∂ρ
（5）
∂ 2u ∂x 2
+
∂2u ∂y 2
+
∂2u ∂z 2
=
∂2u ∂ρ 2
+
∂2u ∂z 2
+
1 ρ2
⋅
∂2u ∂ϕ 2
+
1 ρ
⋅
∂u ∂ρ
∂2u 再用（3）式，变换 ∂ρ 2
+
∂ 2u ∂z 2
。这又可以直接利用（5）式，得
∂2u ∂ρ 2
数学物理方程课件第三章行波法与积分变换法
U (,0)
a 2 2U (, t), (), dU (,0)
dt
(),
t0
U (,t) Acosat Bsin at
U (,0) A ()
B () a
U (,t) () cos at () sin at
a
f(x ) F ()e j
x
f()d
F ()
0
j
数学物理方程与特殊函数
u(x,t) 1 (x at) (x at) 1
xat
( )d
t2
2a xat
t
P( x, t )
依赖区间
x
x at x at
x x1 at
x x2 at
决定区域
x1
x2
x
t
x x1 at
影响区域
x1
x2
x x2 at
x at C 特征线 x at x at 特征变换
第3章行波法与积分变换法
补充作业： 解定解问题
4
2u t 2
25
2u x2
,
u(
x,
0)
sin
x,
u ( x, t
0)
3x,
y 0, x x
数学物理方程与特殊函数
第3章行波法与积分变换法
二 积分变换法
1 傅立叶变换法
傅立叶变换的定义
U (, t) u(x, t)e jxdx
数学物理方程与特殊函数
第3章行波法与积分变换法
u(x,t) 1 (x at) (x at) 1
xat
( )d
2
2a xat
5 达朗贝尔公式的应用
utt
a
u |t0
数学物理方程----3-Bessel-函数PPT优秀课件
贝 塞
erx (r2 + pr + q) = 0 .
尔
由于erx 0，因此，只要 r 满足方程
函 数
r2 + pr + q = 0，
⑤
即 r 是上述一元二次方程的根时，y = erx 就是
④式的解. 方程⑤称为方程④的特征方程. 特征方
程根称为特征根.
西安交通大学 数学与统计学院
例 1 求方程 y - 2y - 3y = 0 的通解.
西安交通大学 数学与统计学院
定理 4 设二阶线性非齐次方程为
数
y + p(x)y + q(x)y = f1 (x) + f2 (x)， ① 第
学
三
物
理 且y1*与y2*分 别 是
方
程
y + p(x)y + q(x)y = f1 (x)，
章 贝 ②塞
尔
和
函
y + p(x)y + q(x)y = f2 (x)
数 学
y = C1 y1 + C2 y2
第 三
物理是该方程的通解，其中
方
C1,
C2为任意常数.
程
章 贝
塞
尔
函
数
西安交通大学 数学与统计学院
定理 3 如果函数 y* 是线性非齐次方程的一个
特解，Y 是该方程所对应的线性齐次方程的通解，则
数
第
学 物
y = Y + y*，
理 方
是线性非齐次方程的通解.
程
数
解 该方程的特征方程为 r2 - 4r + 4 = 0，它 有 第
学 物
重根
数学物理方程第三章练习题
2012-10-3 3 / 69
建立方程、定解条件
∂2u ∂x2i
=
x2i r2
f
′′(r)
+
( 1 r
−
x2i r3
)
f
′(r),
(i = 1, 2, . . . , n)
将上式代入调和方程得
f
′′(r)
+
n
−
1 f
′(r)
=
0,
r
即
f ′′(r) f ′(r)
=
−n
− r
1.
对上式两边积分即得结论.
πx a
,
u(x, b)
=
0.
齐海涛 (SDU)
数学物理方程
2012-10-3 17 / 69
建立方程、定解条件
.E.xample 1.6
用分离变量法求解由下述调和方程的第一边值问题所描述的矩形平板 (0 ≤ x ≤ a, 0 ≤ y ≤ b) 上的稳定温度分布:
.
uxx + uyy = 0,
u(0, y) = u(a, y) = 0,
,
∂r ∂R
=
sin θ,
∂θ ∂R
=
cos θ . r
由 (1.2) 及 (1.3) 知
(1.3)
∂2u ∂z2
=
cos2
θ
∂2u ∂r2
+
sin2 r2
θ
∂2u ∂θ2
+
sin2 r
θ
∂u ∂r
+
sin 2θ r2
∂u ∂θ
−
sin 2θ r
∂2u ∂r∂θ
,
数学物理方程第三章 行波法
(1.7)
f 1 ( ), f 2 ( )是待定的任意二次连续 可微函数 .
u( x,0) f 1 ( x) f 2 ( x) ( x)
u ( x) - x (1.2) t t 0
x
1 f ( x ) f ( x ) ( )d 1 2 u( x ,0) a0 af 1( x ) af 2( x ) ( x ) t x 1 1 x f 1 ( x ) ( x ) ( )d 1 2 2a 0 f 1 ( x ) f 2 ( x ) ( )d
第三章 行波法
• • • • • 主要内容 掌握一维弦振动的解 掌握类比方法求三维、二维问题的解 了解偏微分方程的分类 会求偏微分方程的特征线
§1 弦振动的初值问题
无限长均匀细杆的振动问题，就可以表达成如下形式
2 2u 2 u x , t 0 (1.1) 2 a 2 x t u ( x ) - x (1.2) u(0, x ) ( x ), t t 0
(1.8)
（1.8）称作达朗贝尔公式。这种求解方法也称达朗贝尔解法或行波法（特征线法）。 这种方法对一般的偏微分方程来说是十分困难的。因此只适合波动方程定解问题的求解。
1.2 达朗贝尔公式的物理意义
达朗贝尔公式的物理意义
u( x, t ) f 1 ( x at) f 2 ( x at) (1.7) 先考察 u2 ( x, t ) f 2 ( x at) 的意义
2 2u u 2 a 波 动 方 程 ， 双 曲 型 方 .程 2 2 t x 2 u u 2 a 热 传 导 方 程 ， 抛 物 型程 方. 2 t x
数学物理方程第三章练习题
∂u ∂r
−
sin θ r
∂u ∂θ
,
∂u ∂R
=
sin
θ
∂u ∂r
+
cos θ r
∂u ∂θ
.
R2 + z2 = r2,
tan θ
=
R z
,
(1.1) (1.2)
齐海涛 (SDU)
数学物理方程
2012-10-3 6 / 69
建立方程、定解条件
故有
∂r ∂z
=
cos θ,
∂θ ∂z
=
−
sin r
θ
H1
=
√( ∂x )2 ∂q1
( ∂y )2 + ∂q1
+
(
∂z ∂q1
)2 ,
H2
=
√( ∂x )2 ∂q2
( ∂y )2 + ∂q2
+
(
∂z ∂q2
)2 ,
H3
=
√( ∂x )2 ∂q3
( ∂y )2 + ∂q3
+
(
∂z ∂q3
)2 ,
齐海涛 (SDU)
数学物理方程
2012-10-3
8 / 69
数学物理方程
2012-10-3 2 / 69
1. 建立方程、定解条件 2. 格林公式及其应用 3. 格林函数 4. 强极值原理、第二边值问题解的唯一性
齐海涛 (SDU)
数学物理方程
2012-10-3 3 / 69
建立方程、定解条件
.E.xample 1.1
√
设 u(x1, . . . , xn) = f(r) (其中 r = x21 + · · · + x2n ) 是 n 维调和函数, 试证明
行波法与积分变换法——数学物理方程
1 3 u f1 3 x f1 f2 x f 2 1 3 f1 0 f2 0 C
其解中得f1 , ff21是3两x个二94x次2连34续C可微函数.
于是原方程 f的1 通x 解 为14 x 2
4
4
3.1 一维波动方程的达朗贝尔公式
例 求方程 u x x 2 s in xu x y c o s2x u y y c o sx u y 0
的一般解. 解 特征方程为
d y 2 2 s in x d x d y c o s 2x d x 2 0
dy sinx1 dx
rat 1( )d , r at 0
at)
u
(r,t)
(r
at
)
0
(r
at
)
(at
r
)
0
(at
r
)
2r
1 2ar
atr
atr 1( )d , r at 0
3.2 三维波动方程的泊松公式
二. 一般情况
令
u(r,
t)
f1 x f2 x x ……………①
u t| t 0 a f ' 1 ( x a 0 ) a f 2 ' ( x a 0 )
a '1 x f a '2 x fx ……………②
由第二式得
f1xf2xa10xdC.............③
进一步有
2tu 2 a22ru rr2u 20 2(tr2u)a22(rr2u)0
数学物理方程学习指导书第3章经典方程建立及定解条件
第3章经典方程的成立和定解条件在议论数学物理方程的解法从前，我们第一要弄清楚数学物理方程所研究的问题应当如何提，为此，我们从双方面来议论，一方面要将一个详细的物理、力学等自然科学识题化为数学识题，即成立描绘某种物理过程的微分方程——数学物理方程，称此方程为泛定方程；另一方面要把一个特定的物理现象自己所拥有的详细条件用数学形式表达出来，即列出相应的初始条件和界限条件，二者合称为定解条件.定解条件提出详细的物理问题，泛定方程提供解决问题的依照，作为一个整体称之为定解问题.3．1 经典方程的成立在本节，我们将经过几个不一样的物理模型推导出数学物理方程中三种典型的方程，这些方程组成我们的主要研究对象.经典方程的导出步骤：（1）确立出所要研究的是哪一个物理量u；（2）用数学的“微元法”从所研究的系统中切割出一小部分，再依据相应的物理（力学）规律剖析周边部分和这个小部分间的作用（抓住主要作用，略去次要要素，即高等数学中的抓主部，略去高阶无量小），这类互相作用在一个短的时间间隔是如何影响物理量u3）把这类关系用数学算式（方程）表达出来，经化简整理就是所需求的数学物理方程.例1弦的振动弦的振动问题，固然是一个古典问题，但对于初学者仍旧拥有必定的启迪性.设有一根平均柔嫩的细弦，均衡时沿直线拉紧，并且除受不随时间而变的张力作用及弦自己的重力外，不受外力影响，下边研究弦的细小横向振动，即假定所有运动出此刻一个平面上，并且弦上的点沿垂直于x轴的方向运动（图3-1）.图3-1设弦上拥有横坐标为x的点，在时辰t时的地点为M，位移NM记作u.明显，在振动过程中位移u是变量x与t的函数u(x,t).此刻来成立位移u知足的方程.我们把弦上点的运动先看作小弧段的运动，而后再考虑小弧段趋于零的极限状况.在弦上任取一弧段MM，其长为ds，设是弦的线密度，弧段MM两头所受的张力记作T，T，此刻考虑孤段MM在t 时辰的受力状况，用牛顿运动定律，作用于弧段上任一方向上的力的总和等于这段孤的质量乘以该方向上的加快度.在x 轴方向弧段受力的总和为TcosTcos ，因为弦只作横向振动，所以TcosTcos0.（）假如弦的振动很小，并且在振动过程中弦上的切线倾角也很小，即0,0，则由24cos14!2!可知，当为无量小量时，cos 与1的差量是的高阶无量小量，能够略去不计，所以当0, 0时cos1,cos 1代入(3.1)式，即可近似获得TT .在u 方向弧段受力的总和为Tsin Tsingds ，此中是单位弧段的质量，gds 是弧段MM 的重力.又因当0， 0时 sin1 tg tgu(x,t)，tg 2xsin ' tg 'u(xdx,t)，x2ds 1u(x,t)dxdx,x2且小弧段在时辰 t 沿u 方向运动的加快度为u(x,t)，小弧段的质量为t 2TsinTsin2u(x,t)gdsds2t或gds ，所以（）Tu(xdx,t)u(x,t)gds2u(x,t)dx, x xt 2上式左侧方括号内的部分是因为x 产生dx 的变化而惹起的u(x,t)的改变量，可用微x分取代，即u(xdx,t)u(x,t)xu(x,t)dx2u(x,t)dx, xxxx 2于是T2u(x,t) gdx2u(x,t)dxx 2x 2或T2u(x,t) 2u(x,t) g.x2t22一般说来，张力较大时弧振动速度变化很快，即u要比g 大得多，所以又能够把g 略去.t 2经过这样逐渐略去一些次要的量，抓住主要的量，最后得出 u(x,t)应近似地知足方程2ua 22u（）t 2x 2这里的a2T.式(3.3)称为一维颠簸方程.假如在振动过程中，弦上此外还遇到一个与弦的振动方向平行的外力，且假定单位长 度所受外力的 F(x,t)，明显，在这里(3.1)及(3.2)分别为TcosTcos 0,FdsTsinTsin2u gdsds2.t利用上边的推导方法并略去弦自己的重量，可得弦的逼迫振动方程为2u22uf(x,t),（）’t2x2此中f(x,t)1F(x,t).方程(3.3)与(3.3)’的差异在于(3.3)’的右端多了一个与未知函数u 没关的项f(x,t)，这个项称为自由项，包括有非零自由项的方程称为非齐次方程，自由项恒等于零的方程称为齐次方程.(3.3)为齐次一维颠簸方程， (3.3)’为非齐次一维颠簸方程 .例2 传输线方程对于直流电或低频的沟通电，电路的基尔霍夫定律指出同一支路中电流相等 .但对于较高频次的电流（指频次还没有高到能明显地幅射电磁波的状况） ，电路中导线的自感和电容的效应不行忽视，因此同一支路中电流未必相等.现考虑一来一往的高频传输线，它被看作拥有散布参数的导体（图3-2）.在拥有散布参数的导体中，电流经过的状况，能够用电流强度与电压v来描绘，此处i与v都是x,t的函数，记作i(x,t)与v(x,t)，以R，L，C，G分别表示以下参数：R——每一回路单位的趾串连电阻，L——每一回路单位的串连电感，C——每单位长度的分路电容，G——每单位长度的分路电导.依据基尔霍夫第二定律，在长度为x的传输线中，电压降应等于电动势之和，即v(v v)R xii Lx.t而vx, vx故上式可写成v Ri L i.x T此外，由基尔霍夫第必定律，流入节点x的电流应等于流出该节点的电流，即i(i i)C x iGxv, t或i C vGv.x t将方程(3.4)与(3.5)归并，即得i与v应近似地知足以下方程组i C v Gv0,x tv iRi0.Ltx （）（（（（（（（（）为了确立函数i与v，将方程(3.5)对x微分，同时在方程(3.4)两头乘以C后再对t微分，并把两个结果相减，即得2i G v LC2i RC i x2x t20,t将(3.4)中的v代入上式，得x2i2i(RC i（）x2LG2GL)GRi,t t这就是电流i近似知足的微分方程，采纳近似的方法从（）与（）中消去i可得电压v近似知足的方程2v LG2v(RC GL)v GRv,（）x2t2t方程(3.6)或(3.7)称为传输线方程.依据不一样的详细状况，对参数R,L,C,G作不一样的假定，就能够获得传输线方程的各种特别形式.比如，在高频传输的状况下，电导与电阻所产生的效应能够忽视不计，也就是说可令G R0，此时方程(3.6)与(3.7)可简化为2i12ix2LC t2,2v12vt2LC t2.这两个方程称为高频传输线方程.若令a21这两个方程与(3.3)完整同样.因而可知，同一个方程能够用来描绘不一样的LC物理现象，一维颠簸方程不过颠簸方程中最简单的状况，在流体力学、声学及电磁场理论中，还要研究高维的颠簸方程.例3电磁场方程从物理学我们知道，电磁场的特征能够用电场强度E与磁场强度H以及电感觉强度D 与磁感觉强度B来描绘，联系这些量的麦克斯韦(Maxwell)方程组为rotH J D（）,trotE B,(3.9) tdivB 0,（）divD.(3.11)此中J为传导电流的体密度，为电荷的体密度.这组方程还一定与下述场的物质方程D eE,（）B H,（）J E,（）相联立，此中是介质的介电常数，是导磁率，为导电率，我们假定介质是平均并且是各向同性的，此时，，均为常数.方程(3.8)与(3.9)都同时包括有E与H，从中消去一个变量，就能够获得对于另一个变量的微分方程，比如先消去H，在(3,8)式两头求旋度并利用(3.12)与（）得rotrotH rotE rotE,t将(3.9)与(3.13)代入得rot rotH2H Ht2,t而rotrotH grad div2H,且divH1divB0,所以最后获得H所知足的方程为2H2HH t2t;同理，若消去H即得E所知足的方程2E2E E.t2t假如介质不导电(0)，则上边两个方程简化为2H12H,（）t22E12E,（）2t(3.15)与(3.16)称为三维颠簸方程.若将三维颠簸方程以标量函数的形式表示出来，则可写成2ua 2ua22u2u2u（t2x2y2z2,）此中a21,u是E或H的随意一个重量.从方程(3.11)与(3.12)还能够推导出静电场的电位所知足的微分方程.事实上，以(3.12)代入(3.11)得divD div E divE,而电场强度E与电位u之间存在关系E gradu,所以可得div(gradu)或2u,（）这个非齐次方程称为泊松（Poisson）方程.假如静电场是无源的，即0,则(3.18)变为2u 0，（）这个方程称为拉普拉斯（Laplace）方程.例4热传导方程一块热的物体，假如体内每一点的温度不全同样，则在温度较高的点处的热量就要向温度较低的点处流动，这类现象就是热传导.在工程技术上有很多传热问题都要归纳为求物体内温度的散布，此刻我们来推导传热过程中温度所知足的微分方程，与上例近似，我们不是先议论一点处的温度，而应当先考虑一个地区的温度.为此，在物体中任取一闭曲面S，它所包围的地区记作V（图3-3）.假定在时辰t，地区V内点M(x,y,z)处的温度为u(x,y,z,t)，n为曲面元素S的外法向（从V内指向V外）.图3-3由传热学可知，在t,t t时间内，从S流入地区V的热量与时间t，面积S，以及沿曲面的法线方向的温度变化率三者的乘积成正比，即Q kuk(gradu)n St Stnk(gradu)S t.此中k称为物体的热传导系数，当物体为平均导热体时，k为常数.于是，从时辰t1到时辰t2，经过曲面S流入地区V的所有热量为Q1t2dSdt.kgradut1S流入的热量使V内温度发生了变化，在△t时间内地区V内各点温度从u(x,y,z,t)变化到u(x,y,z,t+△t)，则在△t内V内温度高升所需要的热量为c[u(x,y,z,t t)u(x,y,z,t)]dVVu(x,y,z,t)ct tdV.V进而从时辰t1到时辰t2，因为温度高升所汲取的热量为t2uQ2cdV dt,t1tV此中c为物体的比热，为物体的密度，对平均物体来说，它们都是常数.因为热量守恒，流入的热量应等于物体温度高升所需汲取的热量，即t2t2ukgradudSdt c dVdt.t1t1tS V此式左端的由面积分中S是关闭曲面，能够利用奥-高公式将它化为三重积分，即kgradudS kdiv(gradu)dV k2udV,S V V 所以有t22udVdt t2udVdt.（）k ct1t1tV V因为时间间隔t,tt及地区V都是随意取的，并且被积函数是连续的，所以(3.20)式左右恒等的条件是它们的被积函数恒等，即u a22ua22u2u2u ,（）tx 2y 2z 2此中a 2k .方程(3.21)称为三维热传导方程.c若物体内有热源，其强度为F(x,y,z)，则相应的热传导方程为ua 22u2u2uf(x,y,z,t),tx 2y 2z 2此中fF .c作为特例，假如所考虑的物体是一根细杆 （或一块薄板），或许即便不是细杆（或薄板）而此中的温度u 只与x,t （或x,y,t ）相关，则方程(3.21)就变为一维热传导方程2u a 2u2;tx或二维热传导方程u a 22u2utx 22.y假如我们考虑稳恒温度场，即在热传导方程中物体的温度趋于某种均衡状态，这时温度u 已与时间t 没关，所以u 0,此时方程(3.21)就变为拉普拉斯方程(3.19).因而可知稳恒t温度场内的温度 u 也知足拉普拉斯方程 .在研究气体或液体的扩散过程时，若扩散系数是常数，则所得的扩散方程与热传导方程完整同样.3．2 初始条件与界限条件上边所议论的是如何将过程的物理规律用数学式子表达出来.除此之外，我们还需要把 详细条件也用数学形式表达出来， 这是因为任何一个详细的物理现象都是处在特定条件之下 的.比如弦振动问题，上节所推导出来的方程是全部柔嫩平均的弦作细小横向振动的共同规 律，在推导这个方程时没有考虑到弦在初始时辰物状态以及弦所受的拘束状况.假如我们不 是平常地研究弦的振动，必然就要考虑到弦所拥有的特定条件.因为任何一个详细振动现象 老是在某时辰的振动状态和此时辰从前的状态相关，进而就与初始时辰的状态相关.此外， 弦的两头所受的拘束也会影响弦的振动，端点所处的物理条件不一样会产生不一样的影响，因此弦的振动也不一样 .所以对弦振动问题来说，除了成立振动方程之外，还需列出它的详细条件对热传导方程，拉普拉斯方程也是这样.提出的条件应当恰好能够说明某一详细物理现象的初始状态以及界限上的拘束状况，.用以说明系统的初始状态的条件称为初始条件.用以说明界限上的拘束状况的条件称为界限条件.下边详细说明初始条件和界限条件的表达形式，先谈初始条件，对于弦振动问题来说，初始条件就是弦在开始时辰的位移及速度，若以(x),(x)分别表示初位移和初速度，则初始条件能够表达为u t0(x)u（）t(x) t0而对热传导方程来说，初始条件是指在开始时辰物体温度的散布状况，若以(M)表示t 0时物体内任一点M处的温度，则热传导方程的初始条件就是u(M,t)t0(M).(3.23)泊松方程与拉普拉斯方程都是描绘稳恒状态的，与初始状态无头，所以不提初始条件.再谈界限条件.假如界限条件直接给出了未知函数u(M,t)在界限S上的值，以s表示界限S上的动点，则这样的界限条件可表为u(M,t)MS(s,t),或简写成u S.（）这类界限条件称为第一类界限条件，此中(s,t)表示在界限S上给定的已知函数.比如，在杆的导热问题中，若在端点x a处温度保持为常数u0，这时在端点x a的界限条件为u xa u0.若在端点x a处温度随时间的变化规律f(t)为已知，在这点的界限条件为uxaf(t).又如在弦振动问题中，若弦的某端点x a是固定的，则在该点的位移为零，即uxa0.以上都是第一类界限条件的例子.总之，第一类界限条件直接给出了未知函数u(M,t)在边界S上的值但在很多状况下，界限上的物理条件其实不可以用第一类界限条件来描绘.比如，在杆的导热问题中，若杆的一端xa绝热，那么绝热这个条件就不可以直接给出杆的端点处的温度变化.因为从杆外经过杆端流入杆内的热量为kuSt(此中t为时间间隔，S为杆nxa的截面积，n为杆在端点x a处的外法向，若x a是杆的左端点，n的正向与x轴正向相反，则u u,若x a是杆的右端点，则n的正向与x轴正向同样，则uu), n x n x所以绝热这个条件能够表达为k uSt0, nxa即u0.nxa若在单位时间内经过x a端单位面积流入杆内的热量是t的已知函数f(t)，则这个条件可表示为k u f(t).nxa弦在对于弦振动问题来说，假如弦在x a处沿位移方向的张力（参照x a处是自由的，即沿着位移方向不受外力，中例1的推导）为则此时Tu0,nx a即u0.xx a总之，有时界限条件一定表达为u（）(s,t).n S的形式，此中u.表示函数沿界限外法向的变化率，这类界限条件称为第二类界限条件n除了上述两类界限条件外，有时还会碰到其余形式的界限条件.比如在杆的导热热问题中，若杆在某个端点x a自由冷却，那么自由冷却这个条件就是K uH(u1u xa), nxa(此中u1为四周介质的温度)即uu1h kuh.n xa H这是因为在单位时间内从四周介质传到杆的x a端单位面积上的热量与介质和杆端的温度差成正比，而在单位时间内经过u xa端单位面积传向杆内的热量与n考取例4）.成正比（参xa对于有界杆(0 x l)，若两头都是自由冷却，则在x l处，上述条件可表为uu h u1;nx i在x0处，这个条件可表为uu h u1.nx 0一般地，这类界限条件的形式为u(s,t).（）uhn s这样的界限条件称为第三类界限条件.无论哪一种界限条件，假如它的数学表达式中的右端自由项恒为零，则这类界限条件称为齐次的.定解问题的提法前方两节我们推导了三种不一样种类的偏微分方程并议论了与它们相应的初始条件与边界条件的表达方式.因为这些方程中出现的未知函数的偏导数的最高阶都是二阶，并且它们对于未知函数及其各阶偏导数来说都是线性的，所以这类方程称为二阶线性偏微分方程*）1 .在工程技术上二介线性偏微分方程碰到最多.假如一个函数拥有所需要的各阶连续编导数，并且代入某偏微方程中能使该方程变为恒等式，则此函数称为该方程的解.因为每一个物理过程都处在特定的条件之下，所以我们的任务是要求出合适初始条件和界限条件的解.初始条件和界限条件都称为定解条件.求一个偏微方程知足定解条件的解的问题称为定解问题.只有初始条件，没有界限条件的定解问题称为始值问题（或柯西问题）；而没有初始条件，只有界限条件的定解问题称为边值问题；既有初始条件也有界限条件的定解问题称为混合问题.一个定解问题提得能否切合实质状况，自然一定靠实质来证明，但是从数学角度来看，能够从三方面加以查验. 1）解的存在性，即看所结出来的定解问题能否有解；2）解的独一性，即看能否只有一个解；3）解的稳固性，即看当定解条件有细小改动时，解能否相应地只有细小的改动，假如*）二阶线性编微分方程能够按它们的二阶导数的系数的代数性质进行分类，在§中所推导的颠簸方程属于双曲型,拉普拉斯（或泊松）方程属于椭圆型，热传导方程属于抛物型，对于二阶线性偏微分方程的分类方法,读者可参阅复旦大学数学系编《数学物理方程》（第二版，上海科学技术第一版社第一版）第一章§5.确立这样，此解便称定的，否所得的解就无用价 .因定解条件往常是利用方法得的，因此所获得的果， 有必定的差，假如所以而解的很大， 那么种解然不可以切合客的要求 .假如一个定解存在独一且定的解，此称适定的，在此后中我把着眼点放在定解的解法上， 而极少它的适定性， 是因定解的适定 性常常十分困，而本所的定解都是古典的，它的适定性都是了然的 .习题一1. l 的平均杆，面，一端温度零，另一端有恒定流 q 入（即位内通位截面流入的量q ），杆的初始温度散布是x(ix)，写出相的定解.22. l 的弦两头固定，开始在xc 遇到冲量的作用，写出相的定解.有一平均杆，只需杆中任一小段有向位移或速度，必致段的或伸，种仲开去，就有波沿着杆播，推杆的振方程 .4.一平均杆原l ，一端固定，另一端沿杆的方向被拉e 而静止，忽然松手任其振，成立振方程与定解条件.若F(z),G(z)是随意二可微函数，uF(xat)G(xat)足方程2ua2 2ut 2x 26.若函数u 1(x,t),u 2(x,t), ,u n (x,t),⋯均性次方程2up2uq u r ux 2t 2xt的解，此中p,q,r 不过x,t 的函数，并且数uu k (x,t)收，并x,t 能够行两次k 1逐微分，求数uuk(x,t)足原方程（个叫做性次方程的叠加原理）.k1。
《数学物理方程》第3章 行波法与积分变换法
2 2 dS 1 dd
M Cat : ( x ) 2 ( y) 2 ( at ) 2 at dd
(at ) 2 ( x ) 2 ( y ) 2
sin( x at ) sin( x at ) 1 x at cos d sin( x at ) 解 u( x , t ) 2 2 x at
u( x ,0) f ( x ) g ( x ) 3 x 2 u ( x ,0) 1 f ( x ) g ( x ) 0 1 f ( x ) g ( x ) C y 3 3 9 2 3 2 解 出f ( x ) x C , g( x ) x C 4 4
2 2
2a t
wtt a 2 w xx ( t , x ) 设w( x , t , )是 的解 w( x , , ) 0, wt ( x , , ) f ( x , ) t x a( t ) 1 t f (,)d 则 u( x, t ) w( x, t , )d d 0 x a ( t ) 2a 0 2
x at uxx u 2u u 变换 u 0 2 x at utt a (u 2u u ) u h( )d g() f ( ) g() 方程通解 u( x , t ) f ( x at ) g( x at )
M Sat
0
0
积分中x y z是常数
2 ( x at sin cos ) ( y at sin sin ) ( z at cos ) 2 d ( at ) sin d 0 t 0 4a at
数学物理方程-第3章-2012
对于以密度ρ(x,y,z)分布在区域Ω上的质量而言，根据叠加原理，它所产生的 总引力位势为
( x, y, z) G ( , , )
ddd ( x ) ( y ) ( z )
2 2 2
(1.3)
通过直接计算可以验证，φ(x,y,z)在Ω外满足调和方程
§3-2 格林公式及其应用
1.格林（Green）公式 高等数学中的高斯公式如下
(
在上式中，令 P u
v v v ，于是有 ,Q u ,R u x y z 2 v 2 v 2 v u v u v u v [u( x 2 y 2 z 2 ) x x y y z z ]d
n
内流问题的求解在边界Г 内部进行，例如气罐或管道内流动。外流问题 需要在边界Г 外部的无限大区域内求解，例如翼型或飞机的绕流问题。 从直观认识来看，对于外问题，前述的定解条件下外问题的解并不唯一。 以二维翼型流动为例，仅有法向无穿透条件是不够的，还需要在无穷远 处施加来流条件(速度大小和迎角等)。 （这个问题课本上也有举例）
*
1 u(M 0 ) 4
1 1 u(M ) (u(M ) n ( rM M ) rM M n )dSM 0 0
(2.6)
0 1 1 u ( M ) (u ( M ) ( ) )dSM 2u ( M 0 ) n rM 0 M rM 0 M n 4u ( M ) 0
对于泊松方程Δu=F ，也有类似公式
M0在Ω外 M0在Г上 M0在Ω内
1 1 1 u(M ) F (M ) u(M 0 ) (u(M ) ( ) )dSM d M 4 n rM 0 M rM 0 M n rM 0 M
数学物理方法_第3章 二阶线性常微分方程的幂级数解法本征值问题
y ( x) 2 1a2 3 2a3 x (k 2)( k 1)ak 2 x k
把以上结果代入方程，比较系 数得 2 2
2 1a2 a0 0, 3 2a3 a1 0, 4 3a4 2 a2 0, 5 4a5 2 a3 0,
(2k 1)!
a1.
于是方程的级数解为
1 1 1 y( x) a0 1 ( x)2 ( x) 4 (1) k ( x) 2 k 4! (2k )! 2! 2 k 1 a1 1 1 3 5 k ( x) x ( x) ( x) (1) 3! 5! (2k 1)! a a0 cos x 1 sin x.
n 1
n 1
cn1 (n 1)( x x0 )n ,
n 0
可将式（3.1.4）写成
c
n 0 n n
n ( n 2)( n 1)( x x ) [ ( k 1) a c ]( x x ) n2 0 n k k 1 0 n n 0 k 0
y( x) an x n
n 0
(3.3.2)
于是
y( x) nan x
n 1 n 1
(k 1)ak 1 x k ,
k 0
y( x) n(n 1)an x
n2
2 k 0
n2
(k 2)(k 1)ak 2 x k ,
(1 l )(l 2) 3 (3 l )(1 l )(l 2)(l 4) 5 y1 ( x) x x x 3! 5! (2k 1 l )(2k 3 l ) (1 l )(l 2)(l 4) (l 2k ) 2k 1 x (2k 1)!
数理方程第二版(谷超豪)答案第一章-第三章
3.利用传播波法，求解波动方程的特征问题（又称古尔沙问题）
(ϕ (0) = ψ (0) )
数学物理方程答案
解：u(x,t)=F(x-at)+G(x+at) 令 x-at=0 得 令 x+at=0 得 所以
T ( x) = ρg (l − x)
且 T ( x) 的方向总是沿着弦在 x 点处的切线方向。仍以 u ( x, t ) 表示弦上各点在时刻 t 沿垂直 于 x 轴方向的位移，取弦段 ( x, x + ∆x), 则弦段两端张力在 u 轴方向的投影分别为
ρg (l − x) sin θ ( x); ρg (l − ( x + ∆x)) sin θ ( x + ∆x)
t有
G(x+at) ≡ 常数.
即对任何 x, G(x) ≡ C 0 又 G（x）=
1 1 x C ϕ ( x) + ∫ ψ (α )dα − 2 2 a x0 2a
所以 ϕ ( x),ψ ( x) 应满足
ϕ ( x) +
或
1 x ψ (α )dα = C1 （常数） a ∫x0 1 ϕ ' (x)+ ψ ( x) =0 a
+
x + at 1 (h − α )ψ ( α )dα . 2a (h − x) ∫x − at
即为初值问题的解散。 齐次波动方程初值问题的解仅由右传 ２． 问初始条件 ϕ ( x) 与ψ ( x) 满足怎样的条件时， 播波组成？ 解：波动方程的通解为 u=F(x-at)+G(x+at) 其中 F，G 由初始条件 ϕ ( x) 与ψ ( x) 决定。初值问题的解仅由右传播组成，必须且只须对 于任何 x,
数学物理方程与特征函数-03
0 x l,t 0 t0 0 xl
t
T T 0
a2
0 X 0 (x) B0 n n2 2 / l 2 , n 1,2,3,
n
X n (x) Bn cos l x
0
T 0 T0 (t) D0t C0 C0 u0 B0C0 C0
0 xl X (l) 0
X (0) A B 0, X (l) Ael Bel 0
AB0
X (x) 0
0
X 0 X (x) Ax B A 0 X 0 (x) B0
2 0 X 2 X 0 X (x) Asin x B cosx
2 l
l
n
x cos
0
l
xdx 2l
n
(1)2
l
C0l
C0
dx
0
l 0
C0
n1
Cn
cos
n
l
xdx
l
xdx
0
C0
1 l
l xdx l
0
2
u(x,0) n n
t
Dn
n1
la
cos
l
x0
Dn 0
C0
n1
Cn
cos n
la
t
Dn
sin
n
la
t
cos
n
l
x
u(ux2xu2(,0x0,)ta)2x,t2uu2u(,(lx,xt,)0)0,0,
t
0 x
t0
数学物理方程学习指导书第3章经典方程的建立和定解条件
第3章经典方程的建立和定解条件在讨论数学物理方程的解法以前，我们首先要弄清楚数学物理方程所研究的问题应该怎样提，为此，我们从两方面来讨论，一方面要将一个具体的物理、力学等自然科学问题化为数学问题，即建立描述某种物理过程的微分方程一一数学物理方程，称此方程为泛定方程；另一方面要把一个特定的物理现象本身所具有的具体条件用数学形式表达出来，即列出相应的初始条件和边界条件，两者合称为定解条件•定解条件提出具体的物理问题，泛定方程提供解决问题的依据，作为一个整体称之为定解问题3. 1经典方程的建立在本节，我们将通过几个不同的物理模型推导出数学物理方程中三种典型的方程，这些方程构成我们的主要研究对象•经典方程的导出步骤：(1)确定出所要研究的是哪一个物理量u ；(2)用数学的“微元法”从所研究的系统中分割出一小部分，再根据相应的物理(力学)规律分析邻近部分和这个小部分间的作用(抓住主要作用，略去次要因素，即高等数学中的抓主部，略去高阶无穷小)，这种相互作用在一个短的时间间隔是如何影响物理量u(3)把这种关系用数学算式(方程)表达出来，经化简整理就是所需求的数学物理方程•例1弦的振动弦的振动问题，虽然是一个古典问题，但对于初学者仍然具有一定的启发性设有一根均匀柔软的细弦，平衡时沿直线拉紧，而且除受不随时间而变的张力作用及弦本身的重力外，不受外力影响，下面研究弦的微小横向振动，即假定全部运动出现在一个平面上，而且弦上的点沿垂直于x轴的方向运动(图3-1).图3-1设弦上具有横坐标为X的点，在时刻t时的位置为M,位移NM记作u .显然，在振动过程中位移u 是变量X与t的函数u(x,t) •现在来建立位移u满足的方程我们把弦上点的运动先看作小弧段的运动，然后再考虑小弧段趋于零的极限情况•在弦上任取一弧段MM •，其长为ds,设是弦的线密度，弧段MM •两端所受的张力记作T , T ,现在考虑孤段MM在t时刻的受力情况，用牛顿运动定律，作用于弧段上任一方向上的力的总和等于这段孤的质量乘以该方向上的加速度•在x轴方向弧段受力的总和为_T cos T cos，由于弦只作横向振动，所以T cos - T cos = 0 .( 3.1)如果弦的振动很小，并且在振动过程中弦上的切线倾角也很小，即:-:-0^ 0，则由2 4“ a a _ _cos -1 -2! 4!可知，当:为无穷小量时，cos与1的差量是:的高阶无穷小量，可以略去不计，因此当在u方向弧段受力的总和为-T sin驀"T Sin:「：gds，其中「是单位弧段的质量,-gds是弧段MM •的重力又因当：一0，厂0时上式左边方括号内的部分是由于x产生dx的变化而引起的u(x,t)的改变量，可用微代入(3.1)式，便可近似得到COS：1,cos： 1T T .tg：sin：J i +tg2«tg，皿)，ex=tg：：：u(x dx,t)：u(x，t)lx dx.且小弧段在时刻t沿u方向运动的加速度为2:u(x,t)-t2小弧段的质量为'gds，所以-T sin 二"T sin : 一：gds」ds：：2u(x, t):t2(3.2)T/^-晋」gds -：2u(x,t)吐2dx,:x分代替，即:x£u(x+dx,t) £u(x,t) ◎「初匕,灯［故 $u(x,t )d x—'ex 1 —x,2 2T j u(x,t) j u(x,t)2 ------- --- --------- + 、x 2:xT”g dx 「'5.x 2t 2 g. 般说来，张力较大时弧振动速度变化很快，即u 要比g 大得多，所以又可以把 g 略去•:t经过这样逐步略去一些次要的量，抓住主要的量, 最后得出u( x, t)应近似地满足方程-2:U 2 a 2 一 t x■2u (3.3)T这里的a.式(3.3)称为一维波动方程.如果在振动过程中，弦上另外还受到一个与弦的振动方向平行的外力，且假定单位长 度所受外力的F(x,t)，显然，在这里(3.1)及(3.2)分别为T cos ； " -T cos ：二 0, o2 u Fds-Tsin 二■ T sin ： -，gds …ds —— ct利用上面的推导方法并略去弦本身的重量，可得弦的强迫振动方程为-2 -.2:U 2： U 亨=口 —-1: Xf(x,t),(3.3)'1其中 f(x,t) F(x,t).方程(3.3)与(3.3)'的差别在于(3.3)'的右端多了一个与未知函数 U 无关的项f(x,t),这个项称为自由项，包含有非零自由项的方程称为非齐次方程，自由项恒等于零的方程称为齐次方程.(3.3)为齐次一维波动方程，(3.3)'为非齐次一维波动方程. 例2传输线方程对于直流电或低频的交流电，电路的基尔霍夫定律指出同一支路中电流相等 .但对于较高频率的电流(指频率还没有高到能显著地幅射电磁波的情况)，电路中导线的自感和电容的效应不可忽略，因而同一支路中电流未必相等现考虑一来一往的高频传输线，它被当作具有分布参数的导体(图3-2).在具有分布参数的导体中，电流通过的情况，可以用与电压v来描述，此处i与v都是x,t的函数，记作i(x,t)与v(x,t)，以R，L，C，G分别表示下列参数：R――每一回路单位的趾串联电阻，L——每一回路单位的串联电感，C――每单位长度的分路电容，G――每单位长度的分路电导•根据基尔霍夫第二定律，在长度为x的传输线中，电压降应等于电动势之和，即v -(v ：-v) = R x i L x而v x,ex故上式可写成Ri - L -(3.4)x 汀另外，由基尔霍夫第一定律，流入节点x的电流应等于流出该节点的电流，即=(i ：=i) C x G x v,(3.5) 将方程(3.4)与(3.5)合并，即得i与v 应近似地满足如下方程组这两个方程称为高频传输线方程.X 2=2:V;t 2LC ft 2 , 1；：2V1若令a 2这两个方程与(3.3)完全相同.由此可见，同一个方程可以用来描述不同的LC物理现象，一维波动方程只是波动方程中最简单的情况， 在流体力学、声学及电磁场理论中，还要研究高维的波动方程. 例3电磁场方程从物理学我们知道，电磁场的特性可以用电场强度 E 与磁场强度H 以及电感应强度 D与磁感应强度B 来描述，联系这些量的麦克斯韦(Maxwell)方程组为rotH = J, ct(3.8)cB耐百,(3.9)divB 二 0,(3.10)为了确定函数i 与v ，将方程(3.5)对x 微分，同时在方程(3.4)两端乘以C 后再对t 微分, 并把两个结果相减，即得将(3.4)中的兰代入上式，得CX这就是电流i 近似满足的微分方程，采用类似的方法从( 3.4)与(3.5)中消去i 可得电压V近似满足的方程方程(3.6)或(3.7)称为传输线方程.R, L,C,G 作不同的假定，就可以得到传输线方程的各种特殊形式•例如，在高频传输的情况下，电导与电阻所产生的效应可以忽略不计，也就是 说可令G = R = 0，此时方程(3.6)与(3.7)可简化为；：2i 1；：2i2・X 2iG 兰 一 LC ；x【RC 5, ct :t 2：2i2・i —二 LG X ;t 2 (RC GL)厂 GRi, ct(3.6).2轨LGxV(RC GL)二 GRv, .t :t(3.7)根据不同的具体情况，对参数divD 「.(3.11)将(3.9)与(3.13)代入得而 rot rotH 程为rot rotHrot rotHrotE 「rotE, :t21二 grad div - ' H ,且 divH divB = 0,所以最后得到H 所满足的方同理，若消去 H 即得E 所满足的方程2E如果介质不导电(卞-0)，则上面两个方程简化为2H 1 k 2Hr H , .t 二(3.15)(3.16)(3.15)与(3.佝称为三维波动方程.若将三维波动方程以标量函数的形式表示出来，则可写成其中J 为传导电流的体密度，「为电荷的体密度.这组方程还必须与下述场的物质方程D =eE,(3.13) (1.14)相联立，其中；是介质的介电常数， 」是导磁率，匚为导电率，我们假定介质是均匀而且是各向同性的，此时；，」，；「均为常数.方程(3.8)与(3.9)都同时包含有 E 与H ，从中消去一个变量，就可以得到关于另一个变 H ，在(3,8)式两端求旋度并利用(3.12)与(3.14)得(3.12) 量的微分方程，例如先消去(3.17)2 2 22- 2『召u d u & u= aJ = a — +一 +一2 , ^c x cy 氐丿 其中a_W ，u 是E 或H 的任意一个分量从方程(3.11)与(3.12)还可以推导出静电场的电位所满足的微分方程 •事实上，以(3.12)代入(3.11)得divD = div ；E = ； div E = ■,而电场强度E 与电位u 之间存在关系E = -gradu,所以可得Pdiv(gradu)= -一z或'、u = - — ,( 3.18)z这个非齐次方程称为泊松(Poisson )方程.如果静电场是无源的，即?=0,则(3.18)变成(3.19)这个方程称为拉普拉斯(Lap lace )方程. 例4热传导方程一块热的物体，如果体内每一点的温度不全一样，则在温度较高的点处的热量就要向 温度较低的点处流动，这种现象就是热传导.在工程技术上有许多传热问题都要归结为求物体内温度的分布，现在我们来推导传热过程中温度所满足的微分方程， 与上例类似，我们不是先讨论一点处的温度，而应该先考虑一个区域的温度.为此，在物体中任取一闭曲面S ，它所包围的区域记作 V (图3-3).假设在时刻t ，区域V 内点M(x,y,z)处的温度为u(x,y,z,t)，n 为曲面元素 S 的外法向(从V 内指向V 夕卜).图3-3:t2由传热学可知，在lt,^ t ］时间内，从 S 流入区域V 的热量与时间At ，面积 S ,以及沿曲面的法线方向的温度变化率三者的乘积成正比，即-k(grad u ) S t.其中k 称为物体的热传导系数，当物体为均匀导热体时，k 为常数.于是，从时刻t i 到时刻t 2，通过曲面S 流入区域V 的全部热量为t 2_Q i = tkgrad u1一 S流入的热量使 V 内温度发生了变化，在厶t 时间内区域V 内各点温度从u(x,y,z,t)变化到 u(x,y,z,t+ △ t)，则在△ t 内V 内温度升高所需要的热量为I I I c [u(x, y,z,t t)-u(x,y,z,t)]dVV从而从时刻t i 到时刻t 2，由于温度升高所吸收的热量为tJ ■ u二 c dV dt,飞 _ v ：t其中C 为物体的比热， '为物体的密度，对均匀物体来说，它们都是常数•由于热量守恒，流入的热量应等于物体温度升高所需吸收的热量，即上2 I.t 2. .：uJ t IJJkgradudS dt= J t | JJJ c°〒dV dt.ti.S _t^V醴—此式左端的由面积分中 S 是封闭曲面，可以利用奥 -高公式将它化为三重积分，即11 kgradudS 二 kdiv(gradu)dVSV因此有由于时间间隔't,r 11及区域V 都是任意取的，并且被积函数是连续的，所以 (3.20)式左 右恒等的条件是它们的被积函数恒等，即= Q= k S = t =cnk(grad udS dt.fu(x, y,z,t)dttdV.Q 2二 k 2udV,V：k 2udVti_ Vt2u IF t 」川 cP^dV dt. V ：t(3.20)3(3.21 ).u a :t或二维热传导方程如果我们考虑稳恒温度场，即在热传导方程中物体的温度趋于某种平衡状态，这时温 度u 已与时间t 无关，所以 —=0,此时方程(3.21)就变成拉普拉斯方程(3.19).由此可见稳恒Ct温度场内的温度 U 也满足拉普拉斯方程•在研究气体或液体的扩散过程时，若扩散系数是常数，则所得的扩散方程与热传导方 程完全相同•3. 2初始条件与边界条件上面所讨论的是如何将过程的物理规律用数学式子表达出来•除此以外，我们还需要把具体条件也用数学形式表达出来，这是因为任何一个具体的物理现象都是处在特定条件之下 的.例如弦振动问题，上节所推导出来的方程是一切柔软均匀的弦作微小横向振动的共同规 律，在推导这个方程时没有考虑到弦在初始时刻物状态以及弦所受的约束情况 •如果我们不是泛泛地研究弦的振动，势必就要考虑到弦所具有的特定条件 .因为任何一个具体振动现象总是在某时刻的振动状态和此时刻以前的状态有关，从而就与初始时刻的状态有关•另外，弦的两端所受的约束也会影响弦的振动， 端点所处的物理条件不同会产生不同的影响， 因而弦的振动也不同.所以对弦振动问题来说，除了建立振动方程以外，还需列出它的具体条件 对热传导方程，拉普拉斯方程也是如此.提出的条件应该恰恰能够说明某一具体物理现象的初始状态以及边界上的约束情况，其中a 2 = 2u=a 2；t2 2 2\ d u d u d u r + r + r cy cz厂■方程(3.21)称为三维热传导方程.若物体内有热源，其强度为F (x, y,z)，则相应的热传导方程为:u 2 a :t2 2 2u u u +r 2" 2x : yf(x, y,z,t),如果所考虑的物体是一根细杆 作为特例， 而其中的温度u 只与x,t (或x, y,t )有关,(或一块薄板)，或者即使不是细杆 (或薄板) 则方程(3.21)就变成一维热传导方程:U 2a：t-2；u-2:y 2用以说明系统的初始状态的条件称为初始条件 •用以说明边界上的约束情况的条件称为边界条件•下面具体说明初始条件和边界条件的表达形式，先谈初始条件，对于弦振动问题来说， 初始条件就是弦在开始时刻的位移及速度，若以 「(x),t (X )分别表示初位移和初速度，则 初始条件可以表达为t」= (x)-(X )t =0而对热传导方程来说，初始条件是指在开始时刻物体温度的分布情况，若以 「(M)表 示t=0时物体内任一点 M 处的温度，则热传导方程的初始条件就是u(M,t )tn 「(M).泊松方程与拉普拉斯方程都是描述稳恒状态的，与初始状态无头，所以不提初始条件 再谈边界条件•如果边界条件直接给出了未知函数 u(M ,t)在边界S 上的值，以s 表示边界S 上的动点，则这样的边界条件可表为u(M,t )M ：s 「(St),或简写成(3.24)这种边界条件称为第一类边界条件，其中(s,t)表示在边界S 上给定的已知函数•例如，在杆的导热问题中，若在端点 x = a 处温度保持为常数 u 0，这时在端点x = a 的边界条件为uxn =u 0・若在端点X =a 处温度随时间的变化规律 f (t)为已知，在这点的边界条件为ux 「f (t)・又如在弦振动问题中，若弦的某端点X = a 是固定的，则在该点的位移为零，即Uxn=0・以上都是第一类边界条件的例子 •总之，第一类边界条件直接给出了未知函数u(M ,t)在边界S 上的值但在许多情况下，边界上的物理条件并不能用第一类边界条件来描述 .例如，在杆的导热问题中，若杆的一端x 二a 绝热，那末绝热这个条件就不能直接给出杆的端点处的温度变cu 丄丄化.由于从杆外通过杆端流入杆内的热量为k —— 心S^t (其中A t 为时间间隔，AS 为杆cn(3.22)(3.23)的截面积，n 为杆在端点x=a 处的外法向，若 x=a 是杆的左端点，n 的正向与x 轴正向PuQuQu Qu相反，贝厂■,若x=a 是杆的右端点，则n 的正向与x 轴正向相同，贝厂 '), cnex cn ex所以绝热这个条件可以表达为0( Ik —也S 也 t = 0,cn xy即二 0.x =a若在单位时间内通过 x=a 端单位面积流入杆内的热量是 t 的已知函数f(t)，则这个条件可表示为对于弦振动问题来说， 如果弦在x=a 处是自由的，即沿着位移方向不受外力， 则此时弦在x=a 处沿位移方向的张力(参照 3.1中例1的推导) 为的形式，其中—表示函数沿边界外法向的变化率，这种边界条件称为第二类边界条件•cn除了上述两类边界条件外，有时还会遇到其他形式的边界条件 •例如在杆的导热热问题中，若杆在某个端点 x = a 自由冷却，那末自由冷却这个条件就是cuK 石=H (5 -ux =ax=a ),(其中U i 为周围介质的温度)即(加(k )u + h ———:u 1 h -1州」 X HI H 丿=f(t).x^a总之，有时边界条件必须表达为=0,(3.25)对于有界杆（0乞x 乞I ），若两端都是自由冷却，则在X = I 处，上述条件可表为-U 在x = 0处，这个条件可表为这是由于在单位时间内从周围介质传到杆的 X = a 端单位面积上的热量与介质和杆端的温度差成正比，而在单位时间内通过 考3.1中例4）.x = a 端单位面积传向杆内的热量与x=acn 丿X -i 成正比般地，这种边界条件的形式为(3.26)这样的边界条件称为第三类边界条件•不论哪一种边界条件，如果它的数学表达式中的右端自由项恒为零，则这种边界条件称为齐次的.3.3定解问题的提法前面两节我们推导了三种不同类型的偏微分方程并讨论了与它们相应的初始条件与边界条件的表达方式•由于这些方程中出现的未知函数的偏导数的最高阶都是二阶，而且它们对于未知函数及其各阶偏导数来说都是线性的，所以这种方程称为二阶线性偏微分方程*）1在工程技术上二介线性偏微分方程遇到最多•如果一个函数具有所需要的各阶连续编导数，并且代入某偏微方程中能使该方程变成恒等式，则此函数称为该方程的解•由于每一个物理过程都处在特定的条件之下，所以我们的任务是要求出适合初始条件和边界条件的解•初始条件和边界条件都称为定解条件•求一个偏微方程满足定解条件的解的问题称为定解问题只有初始条件，没有边界条件的定解问题称为始值问题（或柯西问题）；而没有初始条件，只有边界条件的定解问题称为边值问题；既有初始条件也有边界条件的定解问题称为混合问题•一个定解问题提得是否符合实际情况，当然必须靠实际来证实，然而从数学角度来看，可以从三方面加以检验•1）解的存在性，即看所结出来的定解问题是否有解；2）解的唯一性，即看是否只有一个解；3）解的稳定性，即看当定解条件有微小变动时，解是否相应地只有微小的变动，如果*）二阶线性编微分方程可以按它们的二阶导数的系数的代数性质进行分类，在§ 1・1中所推导的波动方程属于双曲型，拉普拉斯（或泊松）方程属于椭圆型，热传导方程属于抛物型，关于二阶线性偏微分方程的分类方法，读者可参阅复旦大学数学系编《数学物理方程》（第二版，上海科学技术岀版社岀版）第一章§5.确定如此，此解便称为稳定的，否则所得的解就无实用价值•因为定解条件通常总是利用实验方法获得的，因而所得到的结果，总有一定的误差，如果因此而解的变动很大，那末这种解显然不能符合客观实际的要求•如果一个定解问题存在唯一且稳定的解，则此问题称为适定的，在以后讨论中我们把着眼点放在讨论定解问题的解法上，而很少讨论它的适定性，这是因为讨论定解问题的适定性往往十分困难，而本书所讨论的定解问题都是古典的，它们的适定性都是经过证明了的.习题一1.长为I的均匀杆，侧面绝缘，一端温度为零，另一端有恒定热流q进入(即单位时间内通过单位截面积流入的热量为q ),杆的初始温度分布是x— x)，试写出相应的定解问题22.长为I的弦两端固定，开始时在x = C受到冲量的作用，试写出相应的定解问题3.有一均匀杆，只要杆中任一小段有纵向位移或速度，必导致邻段的压缩或伸长，这种仲缩传开去，就有纵波沿着杆传播，试推导杆的纵振动方程4.一均匀杆原长为I，一端固定，另一端沿杆的轴线方向被拉长e而静止，突然放手任其振动，试建立振动方程与定解条件.5.若F(z),G(z)是任意二阶可微函数，验证u = F (x at) G(x -at)满足方程-2 -2:-U 2： U2—a 2.t x6.若函数U i(x,t), U2(x, t)JM, U n(x, t)，…均为线性齐次方程-2 -2:u : u-~2 P 2~x :t的解，其中p,q, r只是x, t的函数，而且级数u k(x,t)收敛，并对x,t可以进行两次kz!逐项微分，求证级数u=：£u k(x,t)满足原方程(这个结论叫做线性齐次方程的叠加原理)k m。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
• 物理背景：用于描述稳定或平衡的物理现象
§3-1 方程的建立及其定解条件
调和方程，又称拉普拉斯(Laplace)方程，其三维形式为
2 u 2 u 2u u 2 2 2 0 x y z
这个方程相应的非齐次方程，称为泊松(Poisson)方程，即
(1.1)
1.方程的导出 数学史上导致调和方程的一个著名实例来自牛顿万有引力。根据万有 引力定律，位于(x0，y0，z0)处质量为M的质点对位于(x，y，z)处具有单位 质量的质点的引力,其大小等于GM/r2，而作用方向沿着这两点的连线，指 向(x0，y0，z0)点，其中r为两点之间的距离。写为向量形式，即为
因此，对于狄利克雷或诺依曼外问题而言，还需要在无穷远处对解添加一 定的限制条件。在三维情况下，一般要求解在无穷远处的极限为零（或者 说极限为某个特定的值），即
lim u( x, y, z ) 0
r
r x2 y 2 z 2
(1.1) (1.2)
要在空间的某个区域中确定方程(1.1)和(1.2)的解，还必须附加一些定解条 件。现在这两个方程中并未出现时间变量，因此它们的解与时间无关，所 以在定解条件中只有边界条件，其定解问题是一种边值问题。 与前面的波动方程和热传导方程类似，对方程(1.1)和(1.2)也可以提出三种 类型的边界条件。本次课程只研究第一及第二边值问题。 1)第一边值问题(狄利克雷条件)： 2)第二边值问题(诺依曼条件)：
r ( x x0 ) 2 ( y y0 ) 2 ( z z0 ) 2
显然引力场函数是位势函数φ的梯度
GM ( x, y , z ) r
F grad i j k x y z
除了允许相差一个任意常数外，位势函数是任意确定的。
并注意到G的任意性，可得
divE = 4πρ
(1.4)
又由库仑定律可知，静电场是有势的，即存在静电位势u=u(x,y,z) ，使
u u u E gradu ( i j k) x y z
Δu＝－4πρ
于是得到静电位势u满足以下的泊松方程
特别地，当某区域内没有电荷存在时，此区域内的静电位势满足调和方程。
0
还可以进一步验证，)在Ω内满足泊 松方程
4
另一个例子是静电场的电位势。设空间有一电荷密度为ρ(x,y,z)的静电场， 在此电场内任取一个封闭曲面Σ包围的区域G ，由静电学知，通过Σ向外的 电通量等于G中总电量的4π倍，即成立
E nds 4 dxdydz
G
其中E为电场强度矢量，而n为Σ上的单位外法线向量。
利用格林公式
E ndS [ E
x
cos(n, x) E y cos(n, y ) E z cos(n, z )]dS
E x E y E z ( )d divEd Ed x y z G G G
对于以密度ρ(x,y,z)分布在区域Ω上的质量而言，根据叠加原理，它所产生的 总引力位势为
( x, y, z) G ( , , )
ddd ( x ) ( y ) ( z )
2 2 2
(1.3)
通过直接计算可以验证，φ(x,y,z)在Ω 外满足调和方程
n
0
法向无穿透条件
内流问题的求解在边界Г 内部进行，例如气罐或管道内流动。外流问题需 要在边界Г 外部的无限大区域内求解，例如翼型或飞机的绕流问题。 从直 观认识来看，对于外问题，前述的定解条件下外问题的解并不唯一。以二 维翼型流动为例，仅有法向无穿透条件是不够的，还需要在无穷远处施加 来流条件(速度大小和迎角等)。 （这个问题课本上也有举例）
F ( x, y , z )
r r
GM x x0 y y0 z z0 ( , , ) F(x,y,z)称为引力场函数 2 r r r r r r r x x0 y y0 z z0 i j k i j k x y z r r r
2u 2u 2u u 2 2 2 f ( x, y, z ) x y z
(1.2)
这类方程在力学、物理学问题中经常遇到。前面两章推导的波动方程和热 传导方程如果去掉了时间导数项，那么方程就可以转化为泊松方程或调和 方程。流体力学中的速度势和流函数都满足调和方程；静电场中的电位势 满足泊松方程。
u ( x, y, z ) ( x , y , z ) g ( x, y, z ) (1.5)
u n
( x , y , z )
g ( x, y , z )
(1.6)
习惯思维中，上述定解问题都认为是在有界区域考虑的。也就是说在某光 滑的闭曲面Г的内部寻找满足边界条件的调和函数。 但在实际运用中，常 常会遇到一些无界区域的问题。例如：要确定一个热源物体外部的稳定温 度场。这种情况下，需要在闭曲面Г的外部寻找满足边界条件的调和函数。 为了显示区别，我们把前一种定解问题称为狄利克雷内问题和诺依曼内问 题，把后一类定解问题称为狄利克雷外问题和诺依曼外问题。 流体力学的内流问题和外流问题就是上述问题的典型代表。考虑不可压无 粘势流，其速度势在流动区域内满足拉普拉斯方程，且在物面边界Г上有
调和函数定义：我们把具有关于空间变量的二阶连续偏导数，且满足调和 方程的函数称为调和函数。复变函数中涉及的只是二元函数。
2.定解条件和定解问题
2 u 2 u 2u u 2 2 2 0 x y z 2u 2u 2u u 2 2 2 f ( x, y, z ) x y z