第一章
§1 方程的导出。定解条件
1．细杆（或弹簧）受某种外界原因而产生纵向振动，以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明),(t x u 满足方程
()??
?
??????=??? ??????x u E x t u x t ρ 其中ρ为杆的密度，E 为杨氏模量。
证：在杆上任取一段，其中两端于静止时的坐标分别为 x 与
+x x ?。现在计算这段杆在时刻t 的相对伸长。在时刻t 这段杆两
端的坐标分别为：
),();,(t x x u x x t x u x ?++?++
其
相
对
伸
长
等
于
),()],([)],([t x x u x
x
t x u x t x x u x x x ?+=??-+-?++?+θ
令
0→?x ，取极限得在点x 的相对伸长为x u ),(t x 。由虎克
定律，张力),(t x T 等于
),()(),(t x u x E t x T x =
其中)(x E 是在点x 的杨氏模量。
设杆的横截面面积为),(x S 则作用在杆段),(x x x ?+两端的力分别为
x u x S x E )()(x u x x S x x E t x )()();,(?+?+).,(t x x ?+
于
是
得
运
动
方
程
tt
u x x s x ???)()(ρx
ESu t x =),(x x x x x ESu x x |)(|)(-?+?+
利用微分中值定理，消去x ?，再令0→?x 得
tt u x s x )()(ρx
??
=
x ESu () 若=)
(x s 常量，则得
22)(t u x ??ρ=))((x
u x E x ????
即得所证。
2．在杆纵向振动时，假设(1)端点固定，(2)端点自由，(3)
端点固定在弹性支承上，试分别导出这三种情况下所对应的边界条
件。
解：(1)杆的两端被固定在l x x ==,0两点则相应的边界条件
为
.0),(,0),0(==t l u t u
(2)若l
x =为自由端，则杆在
l
x =的张力
x
u x E t l T ??=)(),(|l
x =等于零，因此相应的边界条件为
x
u
??|l x ==0 同理，若
0=x 为自由端，则相应的边界条件为
x
u ??∣
00
==x
(3)若l x
=端固定在弹性支承上，而弹性支承固定于某
点，且该点离开原来位置的偏移由函数)(t v 给出，则在l x =端支
承的伸长为)(),(t v t l u -。由虎克定律有
x
u
E
??∣)](),([t v t l u k l x --== 其中k 为支承的刚度系数。由此得边界条件
)(
u x
u
σ+??∣
)
(t f l
x == 其中
E
k =
σ 特别地，若支承固定于一定点上，则,0)(=t v 得边界条件
)(
u x
u
σ+??∣0==l x 。 同理，若0=x
端固定在弹性支承上，则得边界条件
x u
E ??∣)](),0([0t v t u k x -==
即 )(u x
u
σ-??∣).(0t f x -= 3. 试证：圆锥形枢轴的纵振动方程为
2
222)1(])1[(t u
h x x u h x x E ??-=??-??ρ 其中h 为圆锥的高(如图1)
证：如图，不妨设枢轴底面的半径为1，则x 点处截面的半径l 为：
h
x l -
=1 所以截面积2)1()
(h
x
x s -=π。利用第1题，得
])1([)1()(222
2
x u
h x E x
t u h x x ??-??=??-ππρ
若E x E =)
(为常量，则得
22
22)1(])1[(t
u
h x x u h x x E
??-=??-??ρ 4. 绝对柔软逐条而均匀的弦线有一端固定，在它本身重力作用下，此线处于铅垂平衡位置，试导出此线的微小横振动方程。
解：如图2，设弦长为l ，弦的线密度为ρ，则x 点处的张力
)(x T 为
)()(x l g x T -=ρ
且)(x T 的方向总是沿着弦在x 点处的切线方向。仍以),(t x u 表示弦上各点在时刻
t
沿垂直于
x
轴方向的位移，取弦段
),,(x x x ?+则弦段两端张力在u 轴方向的投影分别为
)(sin ))(();(sin )(x x x x l g x x l g ?+?+--θρθρ
其中)(x θ表示)(x T 方向与x 轴的夹角
又 .
sin x u tg ??=
≈θθ 于是得运动方程 x u
x x l t u x
???+-=???)]
([2
2
ρ∣
x
u
x l g x x ??--?+]
[ρ∣g x ρ 利用微分中值定理，消去x ?，再令0→?x 得
])[(2
2x u
x l x g t u ??-??=??。
5. 验证
2
221),,(y x t t y x u --=
在锥
222y x t -->0中都满足波动方程
2
22222y u
x u t u ??+??=??证：函数
2
2
2
1),,(y
x t t y x u --=
在锥
2
22y x t -->0内对变量
t y x ,,有
二阶连续偏导数。且 t y x t t
u
?---=??-
23
222)(
2
52222
3
2222
2)
(3)
(t
y x t y x t t
u ?--+---=??-
-
)2()(2
2223
222y x t y x t ++?--=-
x y x t x
u
?--=??-
23
222)(
()()
222232222
23y
x t y x t x
u -
---+--=??(
)()
222
25
2222y x t y x t -+-
-=-
同理
()()222252222
22y x t y x t y
u
+---=??-
所
以
()()
.22
22
2225222222
2t
u y x t y x t y
u x
u ??=++--=??+
??- 即得所证。
6. 在单性杆纵振动时,若考虑摩阻的影响,并设摩阻力密度涵数(即单位质量所受的摩阻力)
与杆件在该点的速度大小成正比(比例系数设为b), 但方向相反,试导出这时位移函数所满足的微分方程.
解: 利用第1题的推导,由题意知此时尚须考虑杆段
()
x x x ?+,上所受的摩阻力.由题设,单位质量所受摩阻力为
t
u b ??-,故
()x x x ?+,上所受摩阻力为
()()t
u x
x s x p b ?????-
运动方程为:
()()()()t u x x s x b x x u ES t u ES t u
x x s x x x ????-??-??? ????=?????+ρρ22 利用微分中值定理，消去x ?,再令0→?x 得
()()()().22t u
x s x b x u ES x t u x s x ??-??? ??????=??ρ
ρ
若=)
(x s 常数，则得
()()t u
x b x u E x t
u x ??-??? ??????=??ρρ22
若
()()则得方程
令也是常量是常量,.,2ρ
ρρE
a E x E x =
== .2
2
222x
u a t u b t u ??=??+??
§3混合问题的分离变量法 1. 用分离变量法求下列问题的解：
(1)
???
?
?
?
???==<<-=??=??=??==0),(),0()0()1(,3sin 0
2
2
222t l u t u l x x x t u l x u x u a t u o
t t π
解：边界条件齐次的且是第一类的，令
)()(),(t T x X t x u =
得固有函数x l
n x X n
π
sin )(=，且 t l
an B t l an A t T n n n π
πsin cos )(+=，)2,1(Λ=n
于是
∑∞
=+=1
sin )sin cos (),(n n n x l n t l an B t l an A t x u πππ 今由始值确定常数
n A 及n B ，由始值得
∑∞==1sin 3sin n n x l n A l x π
π
∑
∞
==-1
sin )(n n x l n B l an x l x π
π 所以
,13=A ,0=n A 当3≠n
?-=l
n xdx l n x l x an B 0
sin )(2π
π
??? ?
?+???? ??+-=n x l n n l x l n x n l l an ππππ
π
sin cos 2
222)}
(4cos 2sin 2443
333222l
an l x
l n n l x l n n x l =--π
ππ
ππ因此所求解为
∑∞
=--+
=1
44
3
sin sin )1(143sin 3cos ),(n n l n t l an n a l x l t l a t x u π
ππππ
(2) ???
?
?
????=??==??==??-??0)0,(,)0,(0),(0),0(022222x t
u
x l h x u t l t u
t u x u a t u 解：边界条件齐次的，令
)()(),(t T x X t x u =
得：
??
?='==+''0)(,0)0(0l X X X X λ (1) 及
)2(0
2=+''X a T λ。
求问题(1)的非平凡解，分以下三种情形讨论。 ο1 0<λ时，方程的通解为
x
x
e C e
C x X λλ--
-+=21)(
由
0)0(=X 得021=+c c 由
0)(='l X 得021=-----
-l
l
e C e
C λλλλ
解以上方程组，得01
=C ，02=C ，故0<λ时得不到非零
解。
ο2 0=λ时，方程的通解为x c c x X 21)(+=
由边值
0)0(=X 得01=c ，再由0)(='l X 得02=c ，仍
得不到非零解。
ο30>λ时，方程的通解为
x c x c x X λλsin cos )(21+=
由
0)0(=X 得01=c ，再由0)(='l X 得
0cos 2=l c λλ
为了使02
≠c ，必须 0cos =l λ，于是
2
212?
?
?
??+==πλλl n n
)2,1,0(Λ=n
且
相
应
地
得
到
x l
n x X n π21
2sin
)(+=
)2,1,0(Λ=n
将λ代入方程(2)，解得
t a l
n B t a l n A t T n n n ππ21
2sin 212cos
)(+++=
)2,1,0(Λ=n
于
是
∑∞
=++++=0
21
2sin )212sin 212cos
(),(n n n x l
n t a l n B t a l n A t x u πππ再由始值得
???
????++=+=∑∑∞
=∞
=00
212sin 2120212sin n n n n x
l n B a l n x l n A x l h
πππ 容易验证?
??
???+x l n π212sin )2,1,0(Λ=n 构成区间],0[l 上的正交函数系：
?????=≠=++?n m l n
m xdx l n x l m l
当当2
0212sin 212sin 0ππ 利用?
??
???
+x l n π212sin
正交性，得 xdx l
n x l h l A l
n π21
2sin 20+=?
n
n l x l n x n l l h 2
222sin )12(2212cos )12(22??
?
?????
? ??++++-=πππ
n
n h )1()12(82
2-+=
π
0=n B
所
以
∑∞
=+++-=02
2212sin 212cos )12()1(8),(n n x l n t a l n n h
t x u πππ
2。设弹簧一端固定，一端在外力作用下作周期振动，此时定解问题归结为
???
????=??===??=??0)0,()0,(sin ),(,
0),0(2
2
222x t u x u t A t l u t u x u a
t
u ω
求解此问题。
解：边值条件是非齐次的，首先将边值条件齐次化，取
t x l
A
t x U ωsin ),(=，则),(t x U 满足
0),0(=t U ，t A t l U ωsin ),(=
令),(),()
,(t x v t x U t x u +=代入原定解问题，则),(t x v 满足 )1()0,(0)0,(0),(,0),0(sin 2
2
2
222
???
?
???-=??===+??=??x l A x t v x v t l v t v t x l A x v a t v ωωω
)
,(t x v 满足第一类齐次边界条件，其相应固有函数为
x l
n x X n π
sin
)(=，)2,1,0(Λ=n 故设
)2(sin )(),(1
∑∞
==n n x
l n t T t x v π
将方程中非齐次项
t x l
A ωω
sin 2
及初始条件中
x l
A ω
-
按?
??
?
??x l n π
sin
展成级数，得 ∑∞
==1
2sin )(sin n n x l n t f t x l A π
ωω 其中
?=l
n xdx l
n t x l A l t f 02sin sin 2)(π
ωω
l
x l n n l x l n x n l t l A 022222sin cos sin 2???
???+-=ππππωω x l
A t n A n ω
ωπω--=+sin )1(212
x l
n n n π
ψsin
1
∑∞
== 其
中
n
l
n n A xdx l n x l A l )1(2sin 202-=-=?π
ωπωψ 将
(2)
代
入
问题(1)
，得
)
(t T n 满足
???
???
?-='=-=
??
?
??+''+n
n n n n n
n A T T t n A t T l an t T )1(2)0(,0)0(sin )1(2)()(12
2
π
ω
ωπωπ
解
方
程
，
得
通
解
212)(sin )1(2sin cos )(π?
π?ππ-
?-++=+l
an n A t l an B t l an A t T n n n n 由始值，得0=n A
2
2222231)(2)1(}))((2)1(2)1{(1an a A l an n l A n A an B n n n n ?π??ππ?π?π--=----=+所以 ∑∞
=--=12
2sin )
()(2)1({),(n n t l an l an al A t x v π
?π?
x l n t n l an l A n π?π
?π?sin }sin 1)()(2)1(22221?--++
x l n t n l t l an a l an l A n π
?π?π?π?sin }sin sin {)
()()1(212
22∑∞
=---= 因此所求解为
∑∞
=--+=12
22
)()()1(2sin ),(n l an l A t x l A t x u ?π?? x l
n t nt l t l an a π??πsin }sin sin {-? 3．用分离变量法求下面问题的解
???
?
?
?
???===??=+??=??====0||0||00022
222l x x t t u u t u u bshx x u a t u 解：边界条件是齐次的，相应的固有函数为
),2,1(sin
)(Λ==n x l
n x X n π
设
∑∞
==1
sin )(),(n n x l n t T t x u π
将非次项bshx 按}{sin
x l
n π
展开级数，得 ∑∞
==1
sin
)(n n x l
n t f bshx π
其
中
shl bn l
n xdx l n shx l b t f n l
n πππ2)1(sin 2)(2221
0+-==+? 将 ∑∞
==1
sin
)(),(n n x l
n t T t x u π
代入原定解问题，得)(t T n 满足
????
?='=+-=+''+0
)0(,0)0(2)1()()()(2
221
2n n n n n T T shl l n bn t T l an t T πππ 方程的通解为
shl l n bn an l t l an B t l an A t T n n n n 12222)1(2)(sin cos )(+-+?++=πππππ 由0)0(=n T ，得：shl l
n bn an l A n n 12222)1(2)(+-+-=ππ
π 由0)0(='n T ，得0=n B
所
以
)cos 1()1(2)1(
)(12222t l
an shl l n bn an t T n n ππππ--+=+ 所求解为
∑∞
=+-+-=1
2
22122sin )cos 1()()1(2),(n n x l n t l an l n n shl a bl t x u π
πππ 4．用分离变量法求下面问题的解：
????
?
????=??===>??=??+??====0|,|0
||)0(20002
2
222t t l x x t u x l h u u u b x u a t u b t u 解：方程和边界条件都是齐次的。令
)()(),(t T x X t x u =
代入方程及边界条件，得
λ-==+X X T
a bT T "
2
'"2
0)()0(==l X X
由此得边值问题
?
?
?===+0)()0(0
"l X X X X λ 因此得固有值2
?
?
?
??==l n n πλλ，相应的固有函数为
Λ,2,1,sin
)(==n x l
n x X n π
又)(t T 满足方程
022'"=++T a bT T λ
将n λλ
=代入，相应的)(t T 记作)(t T n ，得)(t T n 满足
022
'
"=??
? ??++T l an bT T n n π
一般言之，b 很小，即阻尼很小，故通常有
Λ
,2,1,2
2
=??
? ?? 故得通解 )sin cos ()(t B t A e t T n n n n bt n ωω+=- 其中 2 2 b l an n -?? ? ??=πω 所以 x l n t B t A e t x u n n n n n bt πωωsin )sin cos (),(1 +=∑∞ =- 再由始值，得 ??? ????+-==∑∑∞ =∞=x l n B bA x l n A x l h n n n n n n πωπsin )(0sin 11 所以 102)1(2sin 2+-==?n l n n h xdx l n x l h A ππ 1)1(2+-= = n n n n n n bh A b B πωω 所求解为 .sin )sin (cos )1(2),(1 1x l n t b t n e h t x u n n n n n bt π ωωωπ +-= ∑∞ =+- 数学物理方程小结 第七章 数学物理定解问题 数学物理定解问题包含两个部分：数学物理方程（即泛定方程）和定解条件。 §7.1数学物理方程的导出 一般方法： 第一确定所要研究的物理量u ,第二 分析体系中的任意一个小的部分与邻近部分的相互作用,根据物理规律, 抓住主要矛盾, 忽略次要矛盾。（在数学上为忽略高级小量.）第三 然后再把物理量u 随时间，空间的变为通过数学算式表示出来, 此表示式即为数学物理方程。 （一） 三类典型的数学物理方程 （1）波动方程： 0 :),(:),(:22222222==??-??=?-??→f 当无外力时t x f x u a t u 一维t r f u a t u 三维 此方程 适用于各类波动问题。（特别是微小振动情况.） （2）输运方程： 0 :).(:),(:2222==??-??=?-??→f 无外源时t x f x u a t u 一维t r f u a t u 三维 此方程 适用于热传导问题、扩散问题。 （3）Laplace 方程： . 0(:0 :).程时泊松方程退化拉氏方f f u 泊松方程u 拉氏方程t r ==?=?→ 稳定的温度和浓度分布适用的数学物理方程为Laplace 方程, 静电势u 在电荷密度为零处也满足Laplace 方程 。 §7.2定解条件 定解条件包含初始条件与边界条件。 （1） 初始条件的个数等于方程中对时间最高次导数 的次数。例如波动方程应有二个初始条件, 一般 选初始位移u （x,o ）和初始速度u t （x,0）。而输 运方程只有一个初始条件选为初始分布u （x,o ）, 而Laplace 方程没有初始条件。 （2） 三类边界条件 第一类边界条件： u( r ,t)|Σ = f （1） 第二类边界条件： u n |Σ = f （2） 第三类边界条件： ( u+Hu n )|Σ= f （3） 其中H 为常数. 7.3 二阶线性偏微分方程分类 判别式 , ,0,,0,,0221121222112122211212抛物型a a a 椭圆型a a a 双曲型a a a =-=?<-=?>-=? 波动方程是双曲型的,输运方程为抛物型的,而拉普拉斯方程为椭圆型的. 数学物理方程第三版答案 第一章． 波动方程 §1 方程的导出。定解条件 1．细杆（或弹簧）受某种外界原因而产生纵向振动，以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明),(t x u 满足方程 ()?? ? ??????=??? ??????x u E x t u x t ρ 其中ρ为杆的密度，E 为杨氏模量。 证：在杆上任取一段，其中两端于静止时的坐标分别为 x 与+x x ?。现在计算这段杆在时刻t 的相对伸长。在时刻t 这段杆两端的坐标分别为： ),();,(t x x u x x t x u x ?++?++ 其相对伸长等于 ),()],([)],([t x x u x x t x u x t x x u x x x ?+=??-+-?++?+θ 令 0→?x ，取极限得在点x 的相对伸长为x u ),(t x 。由虎克定律，张力),(t x T 等于 ),()(),(t x u x E t x T x = 其中)(x E 是在点x 的杨氏模量。 设杆的横截面面积为),(x S 则作用在杆段),(x x x ?+两端的力分别为 x u x S x E )()(x u x x S x x E t x )()();,(?+?+).,(t x x ?+ 于是得运动方程 tt u x x s x ???)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-?+?+ 利用微分中值定理，消去x ?，再令0→?x 得 tt u x s x )()(ρx ?? = x ESu () 若=)(x s 常量，则得 22)(t u x ??ρ=))((x u x E x ???? 竭诚为您提供优质文档/双击可除数学物理方法学习心得 篇一：数学物理方程的感想 数学物理方程的感想 通过对数学物理方程一学期的学习，我深深的感受到数学的伟大与博大精深。 当应用数学发展到一定高度时，就会变得越来越难懂，越来越抽象，没有多少实际的例子来说明；物理正好也要利用数学来进行解释和公式推导，所以就出现了数学物理方法。刚开始到结束这门课程都成了我的一大问题。很难理解它的真正意义（含义），做题不致从何入手，学起来越来越费劲。让我很是绞尽脑汁。 后来由于老师耐心的指导与帮助下我开始有了点理解。用数学物理方法来解释一些物理现象，列出微分方程，当然这些微分方程是以物理的理论列出来的，如果不借助于物理方法，数学也没有什么好办法来用于教学和实践，而物理的理论也借助于数学方法来列出方程，解出未知的参数。这就是数学物理方法的根本实质所在。真正要学好数学物理方程 不仅要数学好物理也不能够太差。 接下来我想先对数学物理方程做一个简单的介绍与解 释说明。数学物理方程——描述许多自然现象的数学形式都可以是偏微分方程式 特别是很多重要的物理力学及工程过程的基本规律的 数学描述都是偏微分方程，例如流体力学、电磁学的基本定律都是如此。这些反映物理及工程过程的规律的偏微分方程人们对偏微分方程的研究，从微分学产生后不久就开始了。例如，18世纪初期及对弦线的横向振动研究，其后，对热传导理论的研究，以及和对流体力学、对位函数的研究，都获得相应的数学物理方程信其有效的解法。到19世纪中叶，进一步从个别方程的深入研究逐渐形成了偏微分的一般理论，如方程的分类、特征理论等，这便是经典的偏微分方程理论的范畴。 然而到了20世纪随着科学技术的不断发展，在科学实践中提出了数学物理方程的新问题，电子计算机的出现为数学物理方程的研究成果提供了强有力的实现手段。又因为数学的其他分支（如泛函分析、拓扑学、群论、微分几何等等）也有了迅速发 展，为深入研究偏微分方程提供了有力的工具。因而，20世纪关于数学物理方程的研究有了前所未有的发展，这些发展呈如下特点和趋势： 书本个人总结： 由于物理学，力学和工程技术等方面的许多问题都可以归结为偏微分方程的定解问题，而在数学物理方程这门课上，我们的主要任务便是求解这些定解问题，也就是说在已经列出的方程与定解条件之后，怎样去求既满足方程又满足定解条件的解。 而我们的常用的解决偏微分方程的方法的统一思路是将一个偏微分方程的求解设法转化成一个常微分方程问题的求解。 而我们在学习过程中接触到的常用方法有：分离变量法，行波法，积分变换法和拉普拉斯方程的格林函数法 第二章： 本章主要介绍了分离变量法，介绍了有界弦的自由振动，有限长杆上的热传导，圆域内的二维拉普拉斯方程的定解问题等泛定方程和边界条件都是齐次的偏微分方程的求解，还介绍了非齐次方程的解法，非齐次边界条件的处理等等。 A ． 其中泛定方程和边界条件都是齐次的偏微分方程的求解步骤，取有界弦的自由振动的方程求解作为例子，定解问题为： 第一步：分离变量 目标：分离变量形式的非零解)()(),(t T x X t x u = 结果：函数)(x X 满足的常微分方程和边界条件以及)(t T 满足的常微分方程 条件：偏微分方程和边界条件都是齐次的 第二步：求解本征值问题 利用0)()(''=+x X x X λ和边界条件0)0(=X 和0)(=l X 求出本征值和本函数： 本征值： 本征函数： 第三步：求特解，并叠加出一般解 ? ??????====<<>??=??) ()0,(),()0,(,0),(),0(0 ,0 ,22222x x u x x u t L u t u L x t x u a t u t ψ?0 )(2 )(''=+t T a t T λ ,3,2,1 2)(==n l n n πλx l n πsin (x)X n =x l n at l n D at l n C t x u n n n πππsin )cos sin (),(1∑∞ =+= 第 二 章 热 传 导 方 程 §1 热传导方程及其定解问题的提 1. 一均匀细杆直径为l ，假设它在同一截面上的温度是相同的，杆的表面和周围介质发生热交换，服从于规律 dsdt u u k dQ )(11-= 又假设杆的密度为ρ，比热为c ，热传导系数为k ，试导出此时温度u 满足的方程。 解：引坐标系：以杆的对称轴为x 轴，此时杆为温度),(t x u u =。记杆的截面面积4 2 l π为S 。 由假设，在任意时刻t 到t t ?+内流入截面坐标为x 到x x ?+一小段细杆的热量为 t x s x u k t s x u k t s x u k dQ x x x x ????=???-???=?+221 杆表面和周围介质发生热交换，可看作一个“被动”的热源。由假设，在时刻t 到t t ?+在截面为 x 到x x ?+一小段中产生的热量为 ()()t x s u u l k t x l u u k dQ ??-- =??--=11 1124π 又在时刻t 到t t ?+在截面为x 到x x ?+这一小段内由于温度变化所需的热量为 ()()[]t x s t u c x s t x u t t x u c dQ t ????=?-?+=ρρ,,3 由热量守恒原理得： ()t x s u u l k t x s x u k t x s t u c x t ??-- ????=????11 2 24ρ 消去t x s ??，再令0→?x ，0→?t 得精确的关系： ()11 224u u l k x u k t u c -- ??=??ρ 或 ()()11 22 2112244u u l c k x u a u u l c k x u c k t u --??=--??=??ρρρ 其中 ρ c k a =2 2. 试直接推导扩散过程所满足的微分方程。 解：在扩散介质中任取一闭曲面s ，其包围的区域 为Ω，则从时刻1t 到2t 流入此闭曲面的溶质，由dsdt n u D dM ??-=，其中D 为扩散系数，得 ?????= 2 1 t t s dsdt n u D M 浓度由u 变到2u 所需之溶质为 ()()[]???????????ΩΩΩ ??=??=-=2 12 1121,,,,,,t t t t dvdt t u C dtdv t u C dxdydz t z y x u t z y x u C M 两者应该相等，由奥、高公式得： ????????Ω Ω??==????????? ??????+???? ??????+??? ??????=2 12 11t t t t dvdt t u C M dvdt z u D z y u D y x u D x M 其中C 叫做孔积系数=孔隙体积。一般情形1=C 。由于21,,t t Ω的任意性即得方程： ?? ? ??????+???? ??????+??? ??????=??z u D z y u D y x u D x t u C 3. 砼(混凝土)内部储藏着热量，称为水化热，在它浇筑后逐渐放出，放热速度和它所储藏的 水化热成正比。以()t Q 表示它在单位体积中所储的热量，0Q 为初始时刻所储的热量，则 Q dt dQ β-=，其中β为常数。又假设砼的比热为c ，密度为ρ，热传导系数为k ，求它在浇后温度u 满足的方程。 解： 可将水化热视为一热源。由Q dt dQ β-=及00Q Q t ==得()t e Q t Q β-=0。由假设，放 热速度为 t e Q ββ-0 它就是单位时间所产生的热量，因此，由原书71页，(1.7)式得 ??? ? ??-=+??? ? ????+??+??=??-ρρββc k a e c Q z u y u x u a t u t 20222222 2 4. 设一均匀的导线处在周围为常数温度0u 的介质中，试证:在常电流作用下导线的温度满足微分方程 ()2201224.0ρω ρωρc r i u u c P k x u c k t u +--??=?? 其中i 及r 分别表示导体的电流强度及电阻系数，表示横截面的周长，ω表示横截面面积，而k 表示导线对于介质的热交换系数。 解：问题可视为有热源的杆的热传导问题。因此由原71页(1.7)及(1.8)式知方程取形式为 数学物理方程第二版答案 第一章． 波动方程 §1 方程的导出。定解条件 4. 绝对柔软逐条而均匀的弦线有一端固定，在它本身重力作用下，此线处于铅垂平衡位置，试导出此线的微小横振动方程。 解：如图2，设弦长为l ，弦的线密度为ρ，则x 点处的张力)(x T 为 )()(x l g x T -=ρ 且)(x T 的方向总是沿着弦在x 点处的切线方向。仍以),(t x u 表示弦上各点在时刻t 沿垂直于x 轴方向的位移，取弦段),,(x x x ?+则弦段两端张力在u 轴方向的投影分别为 )(sin ))(();(sin )(x x x x l g x x l g ?+?+--θρθρ 其中)(x θ表示)(x T 方向与x 轴的夹角 又 . sin x u tg ??=≈θθ 于是得运动方程 x u x x l t u x ???+-=???)]([22ρ∣x u x l g x x ??--?+][ρ∣g x ρ 利用微分中值定理，消去x ?，再令0→?x 得 ])[(2 2x u x l x g t u ??-??=??。 5. 验证 2 221),,(y x t t y x u --= 在锥2 22y x t -->0中都满足波动方程 222222y u x u t u ??+??=??证：函数2221),,(y x t t y x u --=在锥2 22y x t -->0内对变量t y x ,,有 二阶连续偏导数。且 t y x t t u ?---=??- 23 222)( 22 52222 3 2222 2 ) (3) (t y x t y x t t u ?--+---=??- - 第一章． 波动方程 §1 方程的导出。定解条件 1．细杆（或弹簧）受某种外界原因而产生纵向振动，以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明),(t x u 满足方程 其中ρ为杆的密度，E 为杨氏模量。 证：在杆上任取一段，其中两端于静止时的坐标分别为 x 与+x x ?。现在计算这段杆在时刻t 的相对伸长。在时刻t 这段杆两端的坐标分别为： 其相对伸长等于 ),()],([)],([t x x u x x t x u x t x x u x x x ?+=??-+-?++?+θ 令 0→?x ，取极限得在点x 的相对伸长为x u ),(t x 。由虎克定律，张力),(t x T 等于 其中)(x E 是在点x 的杨氏模量。 设杆的横截面面积为),(x S 则作用在杆段),(x x x ?+两端的力分别为 于是得运动方程 tt u x x s x ???)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-?+?+ 利用微分中值定理，消去x ?，再令0→?x 得 若=)(x s 常量，则得 22)(t u x ??ρ=))((x u x E x ???? 即得所证。 2．在杆纵向振动时，假设(1)端点固定，(2)端点自由，(3)端点固定在弹性支承上，试分别导出这三种情况下所对应的边界条件。 解：(1)杆的两端被固定在l x x ==,0两点则相应的边界条件为 (2)若l x =为自由端，则杆在l x =的张力x u x E t l T ??=) (),(|l x =等于零，因此相应的边界条件为 x u ??|l x ==0 同理，若0=x 为自由端，则相应的边界条件为 x u ??∣00==x (3)若l x =端固定在弹性支承上，而弹性支承固定于某点，且该点离开原来位置的 偏移由函数)(t v 给出，则在l x =端支承的伸长为)(),(t v t l u -。由虎克定律有 x u E ??∣)](),([t v t l u k l x --== 数学物理方程小结 第七章 数学物理定解问题 数学物理定解问题包含两个部分：数学物理方程（即泛定方程）和定解条件。 §7.1数学物理方程的导出 一般方法： 第一确定所要研究的物理量u ,第二 分析体系中的任意一个小的部分与邻近部分的相互作用,根据物理规律, 抓住主要矛盾, 忽略次要矛盾。（在数学上为忽略高级小量.）第三 然后再把物理量u 随时间，空间的变为通过数学算式表示出来, 此表示式即为数学物理方程。 （一） 三类典型的数学物理方程 （1）波动方程： 0 :) ,(:) ,(:22 2222 22==??-??=?-??→f 当无外力时t x f x u a t u 一维t r f u a t u 三维 此方程 适用于各类波动问题。（特别是微小振动情况.） （2）输运方程： 0 :).(:) ,(:2 2 2 2 ==??-??=?-??→f 无外源时t x f x u a t u 一维t r f u a t u 三维 此方程 适用于热传导问题、扩散问题。 （3）Laplace 方程： . 0(:0 :) .程时泊松方程退化拉氏方f f u 泊松方程u 拉氏方程t r ==?=?→ 稳定的温度和浓度分布适用的数学物理方程为Laplace 方程, 静电势u 在电荷密度为零处也满足Laplace 方程 。 §7.2定解条件 定解条件包含初始条件与边界条件。 （1） 初始条件的个数等于方程中对时间最高次导数的次数。 例如波动方程应有二个初始条件, 一般选初始位移u （x,o ）和初始速度u t （x,0）。而输运方程只有一个初始条件选为初始分布u （x,o ）,而Laplace 方程没有初始条件。 （2） 三类边界条件 第一类边界条件： u( r ,t)|Σ = f （1） 第二类边界条件： u n |Σ = f （2） 第三类边界条件： ( u+Hu n )|Σ= f （3） 其中H 为常数. 7.3 二阶线性偏微分方程分类 判别式 , ,0,,0, ,022112 1222112 12 22112 12抛物型a a a 椭圆型a a a 双曲型a a a =-=?<-=?>-=? 波动方程是双曲型的,输运方程为抛物型的,而拉普拉斯方程为椭圆型的. 数学物理方程有感绝对牛 人写的 The document was prepared on January 2, 2021 书本个人总结： 由于物理学，力学和工程技术等方面的许多问题都可以归结为偏微分方程的定解问题，而在数学物理方程这门课上，我们的主要任务便是求解这些定解问题，也就是说在已经列出的方程与定解条件之后，怎样去求既满足方程又满足定解条件的解。 而我们的常用的解决偏微分方程的方法的统一思路是将一个偏微分方程的求解设法转化成一个常微分方程问题的求解。 而我们在学习过程中接触到的常用方法有：分离变量法，行波法，积分变换法和拉普拉斯方程的格林函数法 第二章： 本章主要介绍了分离变量法，介绍了有界弦的自由振动，有限长杆上的热传导，圆域内的二维拉普拉斯方程的定解问题等泛定方程和边界条件都是齐次的偏微分方程的求解，还介绍了非齐次方程的解法，非齐次边界条件的处理等等。 A ． 其中泛定方程和边界条件都是齐次的偏微分方程的求解步骤，取有界弦的自由振动的方程求解作为例子，定解问题为： 第一步：分离变量 目标：分离变量形式的非零解)()(),(t T x X t x u = 结果：函数)(x X 满足的常微分方程和边界条件以及)(t T 满足的常微分方 程 条件：偏微分方程和边界条件都是齐次的 第二步：求解本征值问题 利用0)()(''=+x X x X λ和边界条件0)0(=X 和0)(=l X 求出本征值和本函数： 本征值： 本征函数： 第三步：求特解，并叠加出一般解 这样的特解都满足齐次偏微分方程和齐次边界条件 第四步：利用本征函数正交性定叠加系数 总结：通过以上例子我们可以得出分离变量的一般方法，总的来说可以分成四步： 一． 首先将偏微分方程的定解问题通过分离变量转化为常微分方程的定解问题。 二． 确定特征值和特征函数。由于特征值是要经过叠加的，所以用来确定特征函数的方程与条件，当函数经过叠加之后仍旧要满足。当边界条件是齐次时，求特征函数就是求一个常分方程满足零边界条件的非零解。 三． 定出特征值和特征函数后，再解其他的常微分方程，把得到的解与特征函数乘起来成为Un(x,t). 四． 最后为了使解满足其余的定解条件，需要把U 叠加起来成为级数形式，叠加出一般解，再利用本征函数的正交性定叠加系数。 0)(2)(''=+t T a t T λ ,3,2,1 2)(==n l n n πλ 成都理工大学 《数学物理方程》模拟试题 一、填空题（3分?10=30分） 1.说明物理现象初始状态的条件叫（ ），说明边界上的约束情况的条件叫（ ），二者统称为 （ ）. 2.三维热传导齐次方程的一般形式是：（ ） . 3 .在平面极坐标系下，拉普拉斯方程算符为 （ ） . 4.边界条件 f u n u S =+??)(σ是第（ ）类边界条件，其中S 为边 界. 5.设函数),(t x u 的傅立叶变换式为),(t U ω，则方程22 222x u a t u ??=??的傅立叶变换为 （ ） . 6.由贝塞尔函数的递推公式有 =)(0x J dx d （ ） . 7.根据勒让德多项式的表达式有)(3 1)(3 202x P x P += （ ）. 8.计算积分 =? -dx x P 2 1 1 2)]([（ ） . 9.勒让德多项式)(1x P 的微分表达式为（ ） . 10.二维拉普拉斯方程的基本解是（ ） . 二、试用分离变量法求以下定解问题（30分）： 1.??? ? ? ????<<=??===>< 2.???? ? ?? ??===><<<+??=??====20,0,8,00,20,16200202 2 2 22x t u t x x u t u t t x x u u u 三、用达朗贝尔公式求解下列一维波动方程的初值问题（10分） ?? ???=??=>+∞<<-∞+??=??==0 ,2sin 0,,cos 0022 2 22t t t u x u t x x x u a t u 四、用积分变换法求解下列定解问题（10分）： ??? ? ???=+=>>=???==, 1, 10,0,1002y x u y u y x y x u 五、利用贝赛尔函数的递推公式证明下式（10分）： )(1)()(' 0' '02x J x x J x J -= 六、在半径为1的球内求调和函数u ，使它在球面上满足 θ21cos ==r u ，即所提问题归结为以下定解问题（10分）: 理工大学数学系 第一章：偏微分方程的基本概念 偏微分方程的一般形式：2211 (,,,,,,)0n u u u F x u x x x ???=???L L 其中12(,,...,)n x x x x =是自变量，12()(,,...,)n u x u x x x =是未知函数 偏微分方程的分类：线性PDE 和非线性PDE ，其中非线性PDE 又分为半线性PDE ，拟线 性PDE 和完全非线性PDE 。 二阶线性PDE 的分类（两个自变量情形）： 2221112222220u u u u u a a a a b cu x x y y x y ?????+++++=?????? （一般形式 记为 PDE （1）） 目的：可以通过自变量的非奇异变换来化简方程的主部，从而据此分类 (,) (,)x y x y ξξηη=?? =? 非奇异 0x y x y ξξηη≠ 根据复合求导公式最终可得到： 22211122222 20u u u u u A A A A B Cu ξξηηξη ?????+++++=??????其中： 22111112221211 122222221112 22()2()()()2()A a a a x x y y A a a a x x x y x y y y A a a a x x y y ξξξξξηξηηξξη ηηηη ?????=++???????????????=+++?????????? ?????=++?????? 考虑22111222( )2()0z z z z a a a x x y y ????++=????如果能找到两个相互独立的解 (,)z x y φ= (,)z x y ψ= 那么就做变换(,) (,)x y x y ξφηψ=?? =? 从而有11220A A == 在这里要用到下面两个引理： 引理1：假设(,)z x y φ=是方程22111222( )2()0z z z z a a a x x y y ????++=???? (1)的特解，则关主部 数学物理方程习题解 习题一 1，验证下面两个函数： (,)(,)sin x u x y u x y e y == 都是方程 0xx yy u u += 的解。 证明：（1 ）(,)u x y = 因为322 2 22 2222 2222 22 322 222 2222 2222 222222 222222 1 1()22 () 2()()11()22()2()()0()() x xx y yy xx yy x u x x y x y x y x x x y u x y x y y u y x y x y x y y y y x u x y x y x y y x u u x y x y =-? ?=- +++-?-=-=++=-??=-+++-?-=-=++--+=+=++ 所以(,)u x y =是方程0xx yy u u +=的解。 （2）(,)sin x u x y e y = 因为 sin ,sin cos ,sin x x x xx x x y yy u y e u y e u e y u e y =?=?=?=-? 所以 sin sin 0x x xx yy u u e y e y +=-= (,)sin x u x y e y =是方程0xx yy u u +=的解。 2，证明：()()u f x g y =满足方程 0xy x y uu u u -= 其中f 和g 都是任意的二次可微函数。 证明：因为 ()()u f x g y = 所以 ()(),()()()() ()()()()()()()()0 x y xy xy x y u g y f x u f x g y u f x g y uu u u f x g y f x g y g y f x f x g y ''=?=?''=?''''-=?-??= 得证。 3， 已知解的形式为(,)()u x y f x y λ=+，其中λ是一个待定的常数，求方程 430xx xy yy u u u -+= 的通解。 解：令x y ξλ=+则(,)()u x y f ξ= 所以2 (),()x xx u f u f ξλξλ'''=?=? (),(),()xy y yy u f u f u f λξξξ'''''=?== 将上式带入原方程得2 (43)()0f λλξ''-+= 因为f 是一个具有二阶连续可导的任意函数，所以2 -430 λλ+=从而12 =3,1λλ=， 故1122(,)(3),(,)()u x y f x y u x y f x y =+=+都是原方程的解，12,f f 为任意的二阶可微函数，根据迭加原理有 12(,)(3)()u x y f x y f x y =+++为通解。 4，试导出均匀等截面的弹性杆作微小纵振动的运动方程（略去空气的阻力和杆的重量）。 解：弹性杆的假设，垂直于杆的每一个截面上的每一点受力与位移的情形都是相 同的，取杆的左端截面的形心为原点，杆轴为x 轴。在杆上任意截取位于 [,]x x x +?的一段微元，杆的截面积为s ，由材料力学可知，微元两端处的相对伸长（应 变）分别是 (,)u x t x ??与(,)u x x t x ?+??，又由胡克定律，微元两端面受杆的截去部分的拉力分别为()(,)u SE x x t x ??与()(,)u SE x x x x t x ?+?+??，因此微元受杆的截去部分的作用力的合力为：()(,)()(,)u u SE x x x x t SE x x t x x ??+?+?-?? 数学物理方程学习总结 四年前匡老师作为我的高数老师走进我的大学生活，如今作为一名研究生，很荣幸又能跟着匡老师学习数学。我本科主修土木工程专业，现在学的是岩石力学专业，主要是跟着导师从事一些关于应力波的研究，所以数学物理方程这门课成了我的必修课。 数学物理方程研究的主要对象是从物理学中提出来的一些偏微分方程。这些方程中的自变量和函数有着鲜明的物理意义，有些问题的解可以通过实验给出，这给偏微分方程的研究指明了方向，同时由于物理学上的需求，就诞生了专门研究有物理意义的偏微分方程的解法。 本学期数学物理方程起初学习了拉普拉斯和傅立叶变换概念、性质以及卷积定理，了解其在微分方程求解中的应用，并着重介绍了Γ函数和β函数的性质以及其两者的关系。然后介绍了三大经典方程的建立和定解条件（泊松方程与拉普拉斯方程都是描述恒稳场状态，与初始状态无关，所以不提初始条件）的提出和表示。第四章和第五章分别详细的讲了分离变量法、行波法和积分变换法在求解经典方程中的应用，主要针对求解热传导方程和波动方程。三种方法有时候可以通用但有时候还是有区别，分离变量法主要用来求解有限区域内定解问题；行波法是一种针对无界域的一维波动方程的求解方法；积分变换法主要是求解一个无界域上不受方程类型限制的方法。第六章主要讲述用格林函数法求解拉普拉斯方程，伊始提出两种拉普拉斯方程的边值问题（狄氏内问题、狄氏外问题、牛曼内问题、牛曼外问题），然后介绍几种格林函数的取得，最后简介求解狄氏问题。最后三章分别介绍几个特殊类型的常微分方程（贝塞尔方程和勒让德方程）的引入和他们性质和求解。数学物理方程概括起来就是使用四种方法求解三种经典方程，介绍求解过程中产生的两种特殊函数的一门学科。 作为数理方程的学习者，本人觉得它确实是一门比较难的课程，真正的难点却并不是只有数理方程课程本身，而是对以前高等数学学过的知识的理解与记忆的加深。所以，我觉得想学好这门课程，不仅要把时间放在对相关内容的巩固、复习上，还得多做课本上的例题、习题。 无限长弦的一般强迫振动定解问题 200(,)(,0)() () tt xx t t t u a u f x t x R t u x u x ?ψ==?=+∈>? =?? =? 解()()().() .0()1 11(,)(,)222x at t x a t x at x a t u x t x at x at d f d d a a ττ??ψξξατατ++----??=++-+ +??????? ???? 三维空间的自由振动的波动方程定解问题 ()22 22222220001,,,,0(,,) (,,)t t u u u a x y z t t x y z u x y z u x y z t ??==???????=++-∞<<+∞>? ????????? =????=??? 在球坐标变换 sin cos sin sin (0,02,0)cos x r y r r z r θ?θ??πθπθ=?? =≤<+∞≤≤≤≤??=? L 21()1 () (,)44M M at r S S M M u M t dS dS a t r a r ?ψππ??''?=+??????????? 乙 (r=at) 221()1() (,)44M M at at S S M M u M t dS dS a t t a t ?ψππ??''?=+??????? ???? 乙无界三维空间自由振动的泊松公式 ()sin cos ()sin sin (02,0)()cos x x at y y at z z at θ?θ??πθπθ'=+?? '=+≤≤≤≤??'=+? L 2()sin dS at d d θθ?= 二维空间的自由振动的波动方程定解问题 ()22 2222200,,,0(,)(,)t t u u u a x y t t x y u u x y x y t ?ψ==??????=+-∞<<+∞>? ???????? ?? ==??? 22000011(,,)22at at u x y t a t a ππθθππ?????= +????????? ???? 傅立叶变换 数学物理方程期末考试试题及答案 一、求解方程（15分） ?????===-=+=-. )()(0002x u x u u a u at x at x xx tt ψ? 其中)0()0(ψ?=。 解：设? ??+=-at x at x ηξ＝则方程变为： 0=ξηu ，)()(at x G at x F u ++-=(8’)由边值条件可得： )()0()2(),()2()0(x G x F x x G F ψ?=+=+ 由)0()0(ψ?=即得： )0()2 ()2( ),(?ψ?--++=at x at x t x u 。 二、利用变量分离法求解方程。（15分） ?????==≥==∈=-====)(,)(, 0,0,),(,00002x u x u t u u Q t x u a u t t t l x x xx tt ψ? 其中l x ≤≤0。0>a 为常数 解：设)()(t T x X u =代于方程得： 0''=+X X λ，0''2=+T a T λ(8’) x C x C X λλsin cos 21+=，at C at C T λλsin cos 21+= 由边值条件得： 21)( ,0l n C πλ== l x n at A at B u n n n πλλsin )sin cos (1 +=∑∞= ?=l n dx l x n x l B 0sin )(2π?，?=l n dx l x n x an A 0sin )(2πψπ 三．证明方程02=--cu u a u xx t )0(≥c 具有狄利克雷边界条件的初边值问题解的唯一性与 稳定性. (15分) 证明：设u e v ct -=代入方程： ?? ???====-=).(),(),(),0()(02102t g t l v t g t v x v v a v t xx t ? 设21,v v 都是方程的解设21v v v -=代入方程得： ?? ???====-=0),(,),0(0002t l v t v v v a v t xx t 由极值原理得0=v 唯一性得证。(8’)由 ≤-21v v ετ≤-2 1v v ，稳定性得证由u e v ct -=知u 的唯一性稳定性 得证。 四．求解二维调和方程在半平面上的狄利克雷问题(15分). ,0,0>=++=?z u u u u zz yy xx ).(0x f u z == 解：设),,(ζηξp 是上半平面内一点，在该点放置单位点电荷，其对称点 ),,(?ηξ-p 格林函数： 222)()()(141 ),,,(?ηξπ ηξ-+-+--=z y x y x G 222)()()(141 ?ηξπ++-+-+z y x 数学物理方程试题（一） 一、填空题（每小题5分，共20分） 1.长为π的两端固定的弦的自由振动，如果初始位移为x x 2sin ，初始速度为 x 2cos 。则其定解条件是 2. 方程 03=??-??x u t u 的通解为 3.已知边值问题???===+0 )()0(0 )()('"πλX X x X x X ，则其固有函数)(x X n = 4.方程0)(222'"2=-++y n x xy y x α的通解为 二.单项选择题（每小题5分，共15分） 1. 拉普拉斯方程02222=??+??y u x u 的一个解是（ ） （A ）xy e y x u x sin ),(= （B ）22),(y x y x u += （C ）2 21),(y x y x u += （D ）22ln ),(y x y x u += 2. 一细杆中每点都在发散热量，其热流密度为),(t x F ,热传导系数为k ，侧面绝热,体密度为ρ,比热为c ，则热传导方程是 ( ) （A ）ρc t x F x u a t u ),(222 22+??=?? （B ）ρc t x F x u a t u ),(222+??=?? (C ) ρc t x u x F a t F ),(22222+??=?? (D) ρc t x u x F a t F ),(22 2+??=?? (其中ρc k a =2) 3. 理想传输线上电压问题??? ??? ?=??=??=??=x aA t u x A x u x u a t u t ωωωsin ,cos )0,(0 2 2 222 ( 其中C L a 1 2 = )的解为（ ） （A ）)(cos ),(at x A t x u +=ω （B ）t a x A t x u ωωcos cos ),(= （C ）t a x A t x u ωωsin cos ),(= （D ）)(cos ),(t a x A t x u -=ω 平流方程.这个方程是： (,)(,)(,)u u x t x t x t t t γκ??=-+?? (1) 其中γ是一个正常数，κ是(,)x t 的函数.这个方程可以用来模拟艾滋病的传播、流体动力学，以及其他有关空气或水中的物质流动问题.作为具体的应用，考察风的运动，假设风沿着某个方向运动，比如x 轴的正方向，其速度为每秒γ米.设在原点处有一工厂，风从工厂处携带一些污染物.令(,)u x t 表示t 时刻x 处的污染物的（线）密度（每米的颗粒个数）.再设污染颗粒以正比于(,)u x t 的速度掉下来，比例常数是0r >，则u 满足方程(1)，(,)(,)x t ru x t γ=-. (a)试证方程 (,)(,)(,)u u x t x t ru x t t t γ??=--?? 的所有解都具有形式(,)()rt u x t e f x t γ-=-. (b)令M 是初始时刻空气中的总颗粒数.试证：在时刻0t >空气中的总颗粒数是rt e M -. -------------------------------------------------------------------------- 对于[,]x x x + 段t 到t t + 时间段内的颗粒物数有以下关系 (,)(,)(,)(,)(,)x x x t t x x t t t x t u x t t u x t dx ku x t ku x x t dt ru x t dtdx +++++-?=-+-??? (,)(,)()(,)x x t t t t x x x x t t x t t x x t u x t u x t dtdx k dx dt ru x t dtdx t x ++++++??∴=--???????? 由于,,x t u 任意性知 (,)(,)[(,)]0x x t t x t u x t u x t k ru x t dtdx t x ++??++=???? (,)(,)(,)u x t u x t k ru x t t x ??∴=--?? ---------------------------------------------------------------------------- 推导过程 k 为风速/m s ，(,)u x t 为t 时刻x 处的污染的（线）物密度（每米的颗粒个数），则(,)ku x t 为单位时间内x 处的颗粒数 ---------------------------------------------------------------------------------------------- 研究t 到t t + 段通过[,]t t t + 短的颗粒物。(,)u x t 为t 时刻x 处的颗粒物（线）密度 (,)x x x u x t dx +? 为t 时刻x 到x x + 段内的颗粒物 ---------------------------------------------------------------------------------------- +(,)x x x u x t t dx ?+??为t t +?的间段x 到x x +?段内的颗粒物 从而[(,)(,)]x x x u x t t u x t dx +?+?-?表示(),t t t +?时间段内 (),x x x +?段颗粒物的变化量 则 (,)t t t ku x t dt +??为(),t t t +?时间段内吹过x 处的颗粒物 则(,)t t t ku x x t dt +?+??为(),t t t +?时间段内吹过x x +?处的颗粒物 污染颗粒以正比于(,)u x t 的速度下落比例常数为0r > 则下落速度为(,)(,)r x t ru x t =- 从而在x x +?段内t 到t t +?的间段下落的颗粒物总数 为x x t t x t rudxdt +?+??? 由于总颗粒物的变化量=吹入(),x x x +?的颗粒物-吹出 的颗粒物-下落的颗粒物 第一章 §1 方程的导出。定解条件 1．细杆（或弹簧）受某种外界原因而产生纵向振动，以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明),(t x u 满足方程 ()?? ? ??????=??? ??????x u E x t u x t ρ 其中ρ为杆的密度，E 为杨氏模量。 证：在杆上任取一段，其中两端于静止时的坐标分别为 x 与 +x x ?。现在计算这段杆在时刻t 的相对伸长。在时刻t 这段杆两 端的坐标分别为： ),();,(t x x u x x t x u x ?++?++ 其 相 对 伸 长 等 于 ),()],([)],([t x x u x x t x u x t x x u x x x ?+=??-+-?++?+θ 令 0→?x ，取极限得在点x 的相对伸长为x u ),(t x 。由虎克 定律，张力),(t x T 等于 ),()(),(t x u x E t x T x = 其中)(x E 是在点x 的杨氏模量。 设杆的横截面面积为),(x S 则作用在杆段),(x x x ?+两端的力分别为 x u x S x E )()(x u x x S x x E t x )()();,(?+?+).,(t x x ?+ 于 是 得 运 动 方 程 tt u x x s x ???)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-?+?+ 利用微分中值定理，消去x ?，再令0→?x 得 tt u x s x )()(ρx ?? = x ESu () 若=) (x s 常量，则得 22)(t u x ??ρ=))((x u x E x ???? 即得所证。 2．在杆纵向振动时，假设(1)端点固定，(2)端点自由，(3) 端点固定在弹性支承上，试分别导出这三种情况下所对应的边界条 件。 解：(1)杆的两端被固定在l x x ==,0两点则相应的边界条件 为 .0),(,0),0(==t l u t u (2)若l x =为自由端，则杆在 l x =的张力 x u x E t l T ??=)(),(|l x =等于零，因此相应的边界条件为 x u ??|l x ==0 同理，若 0=x 为自由端，则相应的边界条件为 x u ??∣ 00 ==x (3)若l x =端固定在弹性支承上，而弹性支承固定于某 点，且该点离开原来位置的偏移由函数)(t v 给出，则在l x =端支 承的伸长为)(),(t v t l u -。由虎克定律有 x u E ??∣)](),([t v t l u k l x --== 其中k 为支承的刚度系数。由此得边界条件 )( u x u σ+??∣ ) (t f l x == 其中 E k = σ 特别地，若支承固定于一定点上，则,0)(=t v 得边界条件 )( u x u σ+??∣0==l x 。 同理，若0=x 端固定在弹性支承上，则得边界条件 x u E ??∣)](),0([0t v t u k x -== 即 )(u x u σ-??∣).(0t f x -= 3. 试证：圆锥形枢轴的纵振动方程为 2 222)1(])1[(t u h x x u h x x E ??-=??-??ρ 其中h 为圆锥的高(如图1) 证：如图，不妨设枢轴底面的半径为1，则x 点处截面的半径l 为：数学物理方程小结85856
数学物理方程第三版第一章答案(全)
数学物理方法学习心得
数学物理方程有感
数学物理方程谷超豪版第二章课后答案
数学物理方程第二版答案
数学物理方程 答案 谷超豪
数学物理方程小结
数学物理方程有感绝对牛人写的
数学物理方程与特殊函数-模拟试题及参考答案
数学物理方程总结材料
数学物理方程习题解答案
数学物理方程学习总结
数学物理方程公式总结-14页文档资料
数学物理方程期末考试试题及答案
数学物理方程考试试题及解答
数学物理方程报告
数学物理方程第一章答案
相关主题
文本预览