第二章 数学物理定解问题
什么是定解问题
§1.2 什么是定解问题1. 定解问题定解问题是根据已知物理规律求解特定物理过程的数学条件,它由泛定方程和定解条件两个部分组成,泛定方程也称为数学物理方程。
2. 泛定方程泛定方程是待解物理过程所遵循的物理规律的数学表达式,具体表现为某物理量关于时间和空间变量的偏微分方程,同一类物理过程遵循相同的物理规律,因此泛定方程反映一类物理过程的共性。
方程中物理量对时间变量的偏微分项反映物理过程的因果关联。
方程中物理量对空间变量的偏微分项反映物理过程的内部作用,或内在关联。
例1. 质点运动状态的演化问题在质点动力学问题中常求质点的运动轨迹,一旦求出运动轨迹,则一切与质点运动有关的物理量(如动能、动量、角动量等)都可求出。
质点的运动状态是由质点的位矢和动量完全确定,求质点运动轨迹的方法就是求解质点的运动状态随时间演变的过程,即由前一时刻的位矢和动量推算出下一时刻位矢和动量,从物理上看前后二时刻质点的运动状态的联系为dt t p m t r dt t r t r dt t r )(1)()()()(K K K K K +=+=+, dt t F t p dt t p t p dt t p )()()()()(K K K K K +=+=+ 因此,只要知道质点的受力情况就能由前一时刻的运动状态求出下一时刻的运动状态,这样的推演过程就是求解常微分方程F t r m K K =)(满足初始条件“0000)(,)(v t r r t r K K K K ==”的解。
§1.3 定解条件。
一、初始条件初始条件描述特定物理过程的起因,就t 这个自变数而言,如果泛定方程中物理量u 对t 最高阶偏导数是n 阶偏导数n n tu ∂∂,则要确定具体的定解问题,需要n 个初始条件。
例1:均匀细杆的导热问题满足的泛定方程为02=−xx t u a u ,则要确定具体的导热问题的解只需一个初始条件:)(0x u t ϕ==,即要已知初始温度分布。
第2章定解问题
第2章定解问题第2章定解问题1、何谓数理⽅程?按其描绘的物理过程,它可分为哪⼏类?2、何谓定解问题?它分为哪⼏类?试写出⼀维波动⽅程的Cauchy问题的数学表⽰。
3、何谓定解条件?它包括哪些内容?4、何谓边界条件?它分为哪⼏类?⼀个边界需⽤⼏个边界条件来描述?5、⽤数理⽅程来研究物理问题需要经历哪⼏个步骤?6、在静电场问题中,由介电常数分别为和的两种介质组成的系统的交界⾯S 处的衔接条件有⼏个?应如何表⽰?7、如何导出物理模型的数理⽅程?在推导弦的横振动⽅程时采⽤了哪些近似?由⼩⾓度近似我们得到什么结论?8、热传导⽅程的扩散⽅程有何共同和不同之处?9、在杆的纵振动问题中,若端⾃由,这个边界条件如何写?你能从Hooke定律出发证明吗?10、在杆的导热问题中,若端绝热,这个边界条件该如何写?你能从⼀物理定律出发证明吗?11、在热传导问题中,若热源密度不随时间⽽变化,则热传导⽅程会发⽣怎样的变化?12、在弦的横振动问题中,若弦受到了⼀与速度成正⽐的阻⼒,该阻⼒对于弦的振动问题是否起到了源的作⽤?若受到了⼀与位移成正⽐的回复⼒呢?第3章⾏波法1、⾏波法的解题要领是什么?它适合⽤来求解哪⼀类定解问题?为什么?2、⼀维波动⽅程的通解为什么含有两个任意函数?他们各个有怎样的形式和怎样的物理意义?靠什么确定他们的具体函数形式?3、公式是⽤⾏波法求解弦的横振动问题时推得的,能否⽤公式求解如下定解问题?请说明原因?4、能否⽤公式求解如下定解问题?5、能否⽤⾏波法求解如下定解问题?6、你能否根据直⾓坐标系中的导出球坐标中球对称情况下的的表达式请记住这个结论:7、何谓平均值法?你能通过引⼊球⾯的平均值,将三维的波动⽅程化为关于平均值的⼀维⽅程吗?8、在Poisson 公式中,?若已知9、对于定解问题除了可⽤Poisson 公式求解外?你能否有其他的求解法?10、在弦的横振动⽅程单位质量的弦所受的外⼒,若将则怎样的物理含意?它的量纲是什么?11、冲量原理的精神是什么?12、你能否⽤纯强迫振动的解来求解定解问题13、试述推迟势的物理意义,在推迟势中,若,且局限于⼀单位球内,则其中的体积分该如何计算?14、对于定解问题按下述⽅法进⾏求解是否正确?为什么?令使由公式可求得⽽显然,所满⾜的定解问题的解为所以,原定理问题的解为第4章分离变量法1、分离变量法的物理背景是什么?为什么能将未知函数表⽰为单元函数的乘积?2、分离变量法适于求解哪些定解问题?能⽤分离变量法求解⽆界问题吗?4、分离变量法有哪⼏个求解步骤?其中最关键的是哪⼀步?5、何谓本征值问题?以下两个定解问题是否构成本征值问题?(1)(2)6、仿照上章⽤冲量原理求解⽆界弦的纯迫振动的思想和⽅法,你能否写出⽤冲量原理求有界弦的纯强迫振动的公式?7、在将边界条件齐次化时,为什么通常可选辅助函数为X的⼀次式,⽽当问题的两个端点均有第⼆类边界条件时,必须选辅助数为X的⼆次式?8、在⽤分离变量法求解圆的Dirichlet问题时,需要将边界条件齐次化吗?为什么?9、在⽤分离变量法求解下述问题时,是否需将边界条件齐次化?如何齐次化?10、在柱坐标和极坐标中对分离变量,所得到的的⽅程为…其后为什么要注明…?它是怎样得来的?11、在扇形区域中,⽤分离变量法求Dirichlet问题应选择什么坐标系?所得到的的⽅程仍是…吗?为什么?12、在⽤分离变量法求解定解问题时,应如何选择坐标系?能在直⾓坐标系中求解吗?5章特殊函数>> 1)勒让德多项式1、⽅程是什么⽅程?你能写出它在中的⼀有限解吗?2、试述Legendre⽅程本征值问题的提法,其本征值、本征函数是什么?3、你能证明吗?你能由和之值算出吗?4、Legendre多项式的母函数是什么?何谓母函数法?它有哪些⽤途?5、Legendre多项式的归⼀化因⼦是什么?模是什么?你能得到⼀正交归⼀的Legendre多项式吗?6、积分和之值分别是多少?和7、你能将⽤Legendre多项式表⽰吗?8、你能否⽤关系式导出递推公式9、在球坐标系中,在轴对称的情况下,△u=0的变量分离形式的解是什么?在球内的解是什么?在球外的解呢?10、什么是缔合Legendre函数?它是否⼀定是多项式?为什么?11、试述缔合Legendre⽅程本征值问题的提法,其本征值和本征函数是什么?12、缔合Legendre函数的模和归⼀化因⼦是什么?13、是否等同于?与有何关系?你能否由的正交归⼀性导出的正交归⼀性?15、何谓球函数⽅程?它满⾜下列条件的特解是什么?16、独⽴的l阶球函数共有多少个?17、你能⽤两种不同的形式,写出在球坐标系中,在⾮轴对称的情况下△u=0的解吗?它们对于球内和球外的具体情况,⼜分别是怎样的呢?2)贝塞⽿函数1、⽅程叫什么⽅程?你能写出它的⼀有限解吗?2、何谓Bessel函数的零点?它与Bessel⽅程的何种本征值问题有关?有什么样的关系?3、Bessel函数的母函数是什么?当v不为整数时有⽆母函数?为什么?4、你能利⽤Bessel函数的母函数关系式推导出Bessel函数的递推公式吗?5、Bessel函数有⽆微分表达式?若有,试写出;若⽆,说明为什么?6、什么是三类柱函数?它们是否均满⾜Bessel⽅程?它们互相的关系是怎样的?7、第⼆、三类柱函数是否也满⾜Bessel函数递推公式?为什么?8、9、10、Bessel⽅程的通解是什么?其有限解是什么?11、什么是虚宗量的Bessel⽅程?它经过什么样的代换可变成Bessel⽅程?由此你能推得虚宗量的Bessel ⽅程的⼀个特解吗?12、什么是虚宗量的Bessel函数和虚宗量的Neumann函数?虚宗量Bessel⽅程的通解是什么?13、你能完整地写出在柱坐标中对分离变量后所得到的在柱体内的分离变量形式的解吗?14、⽅程在柱坐标系下分离变量,在什么样的边界条件下会出现虚宗量Bessel⽅程?虚宗量的Bessel⽅程是否会构成本征值问题?15、球Bessel⽅程是什么样的情况下出现的?它与半整数的Bessel⽅程有什么关系?你能理解式给出的⼏个函数是球Bessel ⽅程的特解吗?16、试述球Bessel⽅程本征值问题的提法,其本征值和本征函数是什么?17、你能写出在球坐标系中对所得到的分离变量形式的解吗?第6章积分变换法1、何谓积分变换法?他的解题步骤是怎样的?2、Fourier变换的定义是什么?它的存在条件是什么?你能由周期函数的Fourier级数⽽导出⾮周期函数的Fourier积分从⽽引⼊Fourier变换吗?3、试求函数的Fourier变换(a>0),你能利⽤Fourier变换的某些性质求出和吗?其中,a为常数,t为参变量。
数学物理方程_定解问题
根据边界条件确定任意函数 f:
令 故
规定,当
时
4、定解问题是一个整体
达朗贝尔公式的求解过程,与大家熟知的常 微分方程的求结果成完全类似。
但遗憾的是,绝大多数偏微分方程很难求出 通解;即是求出通解,用定解条件确定其中待 定函数往往更为困难。这说明,达朗贝尔公式 不适用于普遍的数学物理定解问题的求解?
(7.1.8)
称式(7.1.8)为弦的自由振动方程。
(2) 如果在弦的单位长度上还有横向外力 作用,则式(7.1.8)应该改写为
(7.1.9)
式中
称为力密度,为 时刻作用于
处单位质量上的横向外力
式(7.1.9)称为弦的受迫振动方程.
2、 均匀杆的纵振动
一根杆,只要其中任一小段做纵向移动,必然使 它的邻段压缩或伸长,这邻段的压缩或伸长又使 它自己的邻段压缩或伸长。这样,任一小段的纵 振动必然传播到整个杆,这种振动的传播是纵波.
泊松方程和拉普拉斯方程的定解条件不包含初始条件, 而只有边界条件. 边界条件分为三类:
1、在边界上直接给定未知函数 , 即为第一类边界条件.
2、在边界上给定未知函数导数的值,即为第二类边界条件.
3、在边界上给定未知函数和它的导数的某种线性组合, 即第三类边界条件.
第一、二、三类边界条件可以统一地写成
第二类边界条件 规定了所研究的物理量在边界外法线方向上方向导数 的数值
u n
x0 , y0 ,z0
f (x0, y0, z0,t)
(7.2.3)
第三类边界条件 规定了所研究的物理量及其外法向导数的线性组合在 边界上的数值
(7.2.4)
其中 是时间 的已知函数, 为常系数.
7.2.2 泊松方程和拉普拉斯方程的定解条件
数学物理方程第一章、第二章习题全解
18
数学物理方程与特殊函数导教·导学·导考
2δρ ut ( x , 0 ) = k ( c - δ≤ x ≤ c + δ) 在这个小段外,初速度仍为零, 我们想得到的是 x = c 处受到冲 击的初速度 , 所 以 最后 还 要 令 δ→ 0。此 外 , 弦是 没 有 初 位 移的 , 即 u( x, 0) = 0 , 于是初始条件为
3. 有一均匀杆 , 只要杆中任一小段有纵向位移或速度 , 必导致 邻段的压缩或伸长, 这种伸缩传开去, 就有纵波沿着杆传播, 试推导 杆的纵振动方程。
解 如图 1 9 所示, 取杆
长方向为 x 轴正向, 垂直于杆长
方向的 各截 面 均 用 它 的 平 衡 位 置 x 标记 , 在时刻 t, 此截面相对
u( x, 0) = 0 0,
ut ( x , 0 ) = δkρ,
| x - c| >δ | x - c | ≤ δ (δ→ 0)
所以定解问题为
utt - a2 uxx = 0
u(0 , t) = u( l, t) = 0 u( x, 0) = 0 , ut ( x , 0 ) =
0, | x - c| > δ δkρ, | x - c | ≤ δ (δ→ 0 )
16
数学物理方程与特殊函数导教·导学·导考
第一章 课后习题全解
1 .4 习题全解
1. 长为 l 的均匀杆 , 侧面绝缘 , 一端温度为零 , 另一端有恒定热
流 q进入 ( 即单位时间内通过单位截面积流入的热量为 q) , 杆的初始
温度分布是 x( l 2
x) ,试写出相应的定解问题。
解 见图 1 8, 该问题是一维热传导方程, 初始条件题中已给
u x
定解问题复习
选择合 适的坐 标系
非齐次方程, 齐次边界条件
特解法
非齐次方程, 齐次定解条件 特征函数法
齐次方程,齐 次边界条件 分离变量法
齐次边界非齐次初始条件下非齐 次方程的解法:
齐次定解条件非齐次方程的解:
方程类型 波动方程 Tn 通解
输运方程
齐次边界非齐次初始条件非齐次方 程的解:
方程类型 Tn
波动方程
解:
定解问题
令 U (x, t) = v (x , t) + w (x, t) ,代入定解问题
即 定解问题
定解问题
定解问题的特解为:
例题
有一长为 l ,侧面绝热而初始温度为零度的均匀细杆, 它的一端保持温度始终为零度,而另一端温度随时间直 线上升,求杆的温度分布。
设杆长方向为 x 轴,x = l 端保持温度始终为零度, x = 0 端 温度随时间直线上升,比例系数为常数 c ,则定解问题为:
2
3
研究细杆导热问题,初始时刻杆的一端温度为零度, 另一端 跟外界绝热,杆上初始温度为 ,试求无热源时细杆上 温度的变化。
长为 l ,两端固定的弦,在单位长度上受横向力 g(x) sinwx 的 作用下做小振动,已知弦的初始位移 和 速度分别为j (x) 和 f (x) ,求其横振动的规律。 有一长为 l ,侧面绝热而初始温度为零度的均匀细杆,它的 一端保持温度始终为零度,而另一端温度随时间直线上升, 求杆的温度分布。
解
令 U (x, t) = v (x , t) + w (x, t) ,代入定解问题
视 v(x, t) 为原方程的特解,考虑到非齐次边界条件,取
将 v(x, t) 代入原定解问题的边界条件,得
v x 0 B ct v x l Al B 0 A ct l
数学物理方程-福州大学-江飞-2.1热传导方程及其定解问题的导出
n
k
u n
k1 u
u1 uΒιβλιοθήκη k k1u nu1
一般形式:u
u n (x,y,z)
g(x, y, z,t)
或
u
u n
g
泛定方程:u t
a2
2u x2
2u y2
2u z2
f
柯 西
初始条件 u
a2u f
问 题
t0
初 边
u g
值边
热管道 1D : ut a2uxx f
t1
则有热源的热传导方程为 ut a2u f a2u F / c .
2. 扩散方程的导出
扩散物从浓度高流向浓度低
* Nerst扩散定律
在该点的扩散系数
扩散物在无穷小时段dt内沿法线方向流过一个无穷
小面积dS的质量dm与扩散物浓度沿曲面dS法线
方向的方向导数N 成正比,即
n
t1,t2
由能量守恒:Q流入 Q吸收
t2 k(x, y, z) udSdt
t1
n
N-L公式及交换下积分次
c(x, y, z)(x, y, z)[u(x, y, z,t2) u(x, y序, z,t1)]dxdydz
t2 t1
ctudxdydzdt
利用高维N-L积分公式,
左端 t2 k(x, y, z) udSdt
dm D(x, y, z) N dSdt
n
因此类似热方程推导:
t2 D(x, y, z) NdSdt
t1
n
(N(x, y, z,t2) N(x, y, z,t1))dxdydz
tN(x, y, z,t) x DxN x DyN x DzN
第二章定解问题
T (x,t) T (x) (t)
对x受力分析,由牛顿第二定律 得
T2 cos 2 T1 cos
T2 sin 2 T1 sin 1 F(x 1x,t)x (x)utt (x 2x, t) ( 2 1)
注意到
sin
tan tan2
ux 1 ux2
ux
sin 1 ux (x, t)
§2.1 引言
一、数学物理方程简介:
数学物理方程是指从物理问题中导 出的反映客观物理量在各个空间、时刻 之间相互制约关系的一些偏微分方程。 方程可以分为线性和非线性方程。
偏微分方程的基本概念:
u u u
mu
F(x1, x2,
, xn,u, x1 , x2 ,
,, xn
, x1m1x2m2
xnmn ) 0
(2)物理问题的数学抽象:
1)由于弦是“细长”的,所以 (x,t) t
忽略重力
2)由于弦“绷紧”于AB两点,这说明弦中各相邻部分之间有 拉力即“张力”作用;由于弦是“柔软”的,所以相邻小段张 力总是弦线的切线方向;
3)由于弦作“微小”的横向振动,故相邻点沿振动方向位移的 差别很小,即
u | ux || x | 1 无穷小量
若 f 0
称为弦的自由振动,振动过程中不受外力。
utt a2uxx
齐次波动方程
事实上,除了以上一维波动方程,像薄膜振动(二维),电 磁场方程(三维)等,均属于波动方程:
utt a2u f (x, y, t)
uxx
uyy
2u x2
2u y 2
utt a22u a2 (uxx uyy uzz )
§2.2 三类数理方程的导出
一、弦的横振动方程(波动方程的建立)
数学物理方程 2-3章课后部分习题答案 李明奇主编 电子科技大学出版社
数学物理方程 李明奇主编 电子科技大学出版社2-3章部分习题答案习题2.14.一根长为L 、截面面积为1的均匀细杆,其x=0端固定,以槌水平击其x=L 端,使之获得冲量I 。
试写出定解问题。
解:由Newton 定律: tt x x Sdxu t x YSu t dx x SYu ρ=-+),(),(,其中,Y 为杨氏模量,S 为均匀细杆的横截面积,x u 为相对伸长率。
化简之后,可以得到定解问题为:⎪⎪⎩⎪⎪⎨⎧-==========)(|,0|0|,0|)/(0002L x Iu u u u u a u Y u t t t L x x x xx xx tt δρρ。
习题2.23.设物体表面的绝对温度为u ,它向外辐射出去的热量,按斯特凡-波尔兹曼定律正比于4u ,即dSdt ku dQ 4=,设物体与周围介质之间,只有热辐射而无热传导,周围介质的绝对温度为已知函数),,,(t z y x ϕ,。
试写出边界条件。
解:由Fourier 热传导实验定律dSdt nuk dQ ∂∂-=1,其中1k 称为热传导系数。
可得dSdt u k dSdt nuk )(441ϕ-=∂∂-,即可得边界条件:)(441ϕ--=∂∂u k k nus。
习题2.34.由静电场Gauss 定理⎰⎰⎰⎰⎰⋅=⋅VsdV dS E ρε01,求证:0ερ=⋅∇E ,并由此导出静电势u 所满足的Poisson 方程。
证明:⎰⎰⎰⎰⎰⎰⎰⎰⋅=⋅=⋅VVsdV dV divE dS E ρε01,所以可以得到:0ερ=divE 。
由E divE ⋅∇=与u E -∇=,可得静电势u 所满足的Poisson 方程:2ερ-=∇u 。
习题2.42.求下列方程的通解:(2):;032=-+yy xy xx u u u (5):;031616=++yy xy xx u u u解:(2):特征方程:03)(2)(2=--dx dy dx dy解得:1-=dx dy 和3=dxdy。
数学物理方程课后参考答案第二章
第 二 章 热 传 导 方 程§1 热传导方程及其定解问题的提1. 一均匀细杆直径为l ,假设它在同一截面上的温度是相同的,杆的表面和周围介质发生热交换,服从于规律dsdt u u k dQ )(11-= 又假设杆的密度为ρ,比热为c ,热传导系数为k ,试导出此时温度u 满足的方程。
解:引坐标系:以杆的对称轴为x 轴,此时杆为温度),(t x u u =。
记杆的截面面积42l π为S 。
由假设,在任意时刻t 到t t ∆+内流入截面坐标为x 到x x ∆+一小段细杆的热量为t x s xuk 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 lkt x l u u k dQ ∆∆--=∆∆--=111124π又在时刻t 到t t ∆+在截面为x 到x x ∆+这一小段内由于温度变化所需的热量为()()[]t x s tuc x s t x u t t x u c dQ t ∆∆∂∂=∆-∆+=ρρ,,3由热量守恒原理得:()t x s u u lk t x s x uk t x s t u c x t ∆∆--∆∆∂∂=∆∆∂∂11224ρ消去t x s ∆∆,再令0→∆x ,0→∆t 得精确的关系:()11224u u l kxu k t u c --∂∂=∂∂ρ或 ()()11222112244u u l c k xu a u u l c k x u c k t u --∂∂=--∂∂=∂∂ρρρ 其中 ρc k a =22. 试直接推导扩散过程所满足的微分方程。
解:在扩散介质中任取一闭曲面s ,其包围的区域 为Ω,则从时刻1t 到2t 流入此闭曲面的溶质,由dsdt nuDdM ∂∂-=,其中D 为扩散系数,得 ⎰⎰⎰∂∂=21t t sdsdt nuDM 浓度由u 变到2u 所需之溶质为()()[]⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰ΩΩΩ∂∂=∂∂=-=2121121,,,,,,t t tt dvdt t uC dtdv t u C dxdydz t z y x u t z y x u C M两者应该相等,由奥、高公式得:⎰⎰⎰⎰⎰⎰⎰⎰ΩΩ∂∂==⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂=21211t t t t dvdt t uC M dvdt z uD z y u D y x u D x M 其中C 叫做孔积系数=孔隙体积。
数学物理方法第2章
§3.1 (1+1)维齐次方程的分离变量法
一、有界弦的自由振动 由第2章的讨论可知,讨论两端固定弦的自由振动规律问题可以归结 为求解下列定解问题:
2u 2u a 2 2 (0 x l , t 0), 2 x t u x 0 0, u x l 0 t 0 , u ( x), u x 0 x l . t 0 t t 0
§3.4
非齐次边界条件的处理
前面所讨论的定解问题,无论方程是齐次的还是 非齐次的,边界条件都是齐次的。如果遇到非齐次边界条件 的情况,应该如何处理?总的原则是设法将边界条件化成齐 次的。具体地说,就是取一个适当的未知函数之间的代换, 使对新的未知函数,边界条件是齐次的。现在仍以一维波动 方程的定解问题为例,说明选取代换的方法。
物理学、力学、工程科学甚至经济和社会科学中等许 多问题都可以归结为偏微分方程的定解问题。第二章中我 们讨论了怎样将一个物理问题表达为定解问题,这一章以 及以下几章的任务是怎样去求解这些定解问题,也就是说 在已经列出方程和定解条件之后,怎样去求既满足方程又 满足定解条件的解. 从微积分学得知,在计算诸如多元函数的微分和积分 (重积分等)时总是把它们转化为单元函数的相应问题来 解决,与此类似,求解偏微分方程的定解问题也可以设法 把它们转化为常微分方程的定解问题来求解。分离变量法 就是这样一种常用的转化方法。在这一章中,我们将通过 一些实例,讨论分离变量法及其应用。
2.1.1 2.1.2
(2.1.3)
这个问题的特点是,偏微分方程是线性齐次的,边界条件也是齐次的。 求解这样的问题可以运用叠加原理。我们知道,在求解常系数齐次常微分方 程的初值问题时,是在先不考虑初始条件的情况下,求出满足方程的足够多 的特解,再利用叠加原理做出这些特解的线性组合,构成方程的通解,然后 利用初始条件来确定通解中的任意常数,得到初值问题的特解。这就启发我 们要求解定解问题(2.1.1)——(2.1.3),须首先寻求齐次方程(2.1.1) 满足边界条件(2.1.2)的足够多的具有简单形式(变量被分离的形式)的 特解,再利用它们做线性组合,得到方程满足边界条件的一般解,再使这个 一般解满足初始条件(2.1.3)。
数学物理方程第一章、第二章习题全解
u x
d
x,
因此小段( x, x + d x) 的伸长( 压缩 ) 为 ud x, 其相对 伸长 (压 缩) 为 x
u x
,
即
x 点处的应变为
u x
(
x,
t)
。若 略
去垂
直杆 长方
向
的形
变
,
根
据
Hooke 定律 , 应力与应变 u 成正比 , 即 x
P=
E
u x
比例系数 E 称为杆的杨氏模量,故所求的纵振动方程为
18
数学物理方程与特殊函数导教·导学·导考
2δρ ut ( x , 0 ) = k ( c - δ≤ x ≤ c + δ) 在这个小段外,初速度仍为零, 我们想得到的是 x = c 处受到冲 击的初速度 , 所 以 最后 还 要 令 δ→ 0。此 外 , 弦是 没 有 初 位 移的 , 即 u( x, 0) = 0 , 于是初始条件为
h c
x
l
h -
c(
l
-
x)
(0 ≤ x ≤ c) ( c < x ≤ l)
ut ( x, 0) = ψ( x ) = 0
则 u( x, t) 是下列定解问题的解 :
utt - a2 uxx = 0
( 0 < x < l, t > 0)
u( x, 0) = φ( x ) , ut ( x, 0 ) = ψ( x )
nπl x
= φ( x )
得
∫ An =
2
l
φ(
x) sin
nπx d x
=
l0
l
∫ ∫ 2 l h xsin nπxd x + 2 l h ( l - x ) si n nπxd x =
数学物理方程-第二章分离变量法
第二章 分离变量法分离变量法是求解偏微分方程定解问题最常用的方法之一,它和积分变换法一起统称为Fourier 方法. 分离变量法的本质是把偏微分方程定解问题通过变量分离,转化为一个所谓的特征值问题和一个常微分方程的定解问题,并把原定解问题的解表示成按特征函数展开的级数形式. 本章介绍两个自变量的分离变量法,更多变量的情形放在其他章节中专门讨论.§21 特征值问题⋅2.1.1 矩阵特征值问题在线性代数中,我们已学过线性变换的特征值问题. 设为一阶实矩阵,A n 可视为到自身的线性变换。
该变换的特征值问题(eigenvalue problem )A n R 即是求方程:,,n Ax x x R λ=∈(1.1)的非零解,其中为待定常数. 如果对某个,问题(1.1)有非零解C λ∈λ,则就称为矩阵的特征值(eigenvalue),相应的称为矩阵n x R λ∈λA n x R λ∈的特征向量(eigenvector). 一般来讲,特征值问题(1.1)有不多于个相A n 异的特征值和线性无关的特征向量. 但可证明: 任一阶矩阵都有个线性无n n 关的广义特征向量,以此个线性无关的广义特征向量作为的一组新基,矩n n R 阵就能够化为标准型.Jordan 若为一阶实对称矩阵,在线性代数中有一个重要结果,即存在一个正A n 交矩阵使得T , 1T AT D -=(1.2)其中diag 为实对角阵. 设,为矩阵的第列D =12(,,...,)n λλλ12[ ... ]n T T T T =i T T i 向量,则式(1.2)可写为如下形式(1)i n ≤≤ ,1212 [ ... ][ ... ]n n A T T T T T T D =或, 1.i i i A T T i n λ=≤≤(1.3)上式说明,正交矩阵的每一列都是实对称矩阵的特征向量,并且这T A 个特征向量是相互正交的. 由于此结论在一定意义下具有普遍性,我们以定n 理的形式给出.定理1.1 设为一阶实对称矩阵,考虑以下特征值问题A n ,,n Ax x x R λ=∈则的所有特征值为实数,且存在个特征向量,它们是相互正交的A n ,1i T i n ≤≤(正交性orthogonality ),可做为的一组基(完备性completeness ).n R 特征值问题在线性问题求解中具有重要的意义,下面举例说明之.为简单起见,在下面两个例子中取为阶非奇异实矩阵,故的所有特A n A 征值非零,并且假设有个线性无关的特征向量 相应的特征值为A n ,i T ., 1i i n λ≤≤例1.1 设,求解线性方程组 .n b R ∈Ax b =解 由于向量组线性无关,故可做为的一组基. 将按此{1}i T i n ≤≤n R ,x b 组基分别展开为,则等价于11 ,nni i i i i i x x T b bT ====∑∑Ax b =,11nni ii ii i x AT bT ===∑∑或,11nni i ii ii i x T bT λ===∑∑比较上式两边的系数可得i T ,1, 1i i i x b i n λ-=≤≤便是原问题的解.12( ... )n x x x x T =例1.2 设,. 求解非齐次常微0n x R ∈12()((),(),...,()), 0n n f t x t x t x t R t T =∈>分方程组, 0(), (0)dxAx f t x x dt=+=(1.4)其中 . '''12((),(),...,()),0n dx x t x t x t t dtT =>解 类似于上例,将按基分别展开为0,,()x x f t {1}i T i n ≤≤ .0111, , ()()nn n i i i ii i i i i x x T x x T f t f t T ======∑∑∑则(1.4)等价于,0111()() +(), (0), 1n n ni i i i i i i i i i i dx t T x t AT f t T x x i n dt =====≤≤∑∑∑或,011()(()()), (0),1nni i i i i i i i i i dx t T x t f t T x x i n dt λ===+=≤≤∑∑比较上式两边的系数可得i T . 0()()(), (0), 1i i i i i i dx t x t f t x x i n dtλ=+=≤≤(1.5)(1.5)是个一阶线性方程的初始值问题,很容易求出其解.请同学们给出解n 的具体表达式.(),1i x t i n ≤≤2.1.2 一个二阶线性微分算子的特征值问题在这一小节,我们讨论在本章常用的一些特征值问题. 代替上节的有限维线性空间和阶实对称矩阵,在这儿要用到线性空间的某个子空间n R n A [0,]C l 和该子空间上的二阶线性微分算子. 一般地取H A在满足齐次边界条件.2{()[0,]()H X x C l X x =∈0,x l =}(1.6)下面我们讨论二阶线性微分算子的特征值问题. 先取边界条件为22d A dx=-,设是的特征函数,即且满足(0)0,()0X X l ==()X x H ∈A ()0X x ≠.()()AX x X x λ=此问题等价于是下面问题的非零解()X x "()()0, 0(0)()0 .X x X x x l X X l λ⎧+=<<⎨==⎩(1.7)(1.7)便是二阶线性微分算子的特征值问题,即要找出所有使22d A dx=-得该问题有非零解的. 下面求解特征值问题(1.7).λ首先证明要使(1.7)具有非零解,必须非负.λ设是相应于的一个非零解,用乘(1.7)中的方程,并在)(x X λ)(x X 上积分得[]l ,0,0)()()()("=+x X x X x X x X λ,0)()()( 0 2 0 "=+⎰⎰dx x X dx x X x X llλ.0)())(()()( 0 2 0 2'0'=+-⎰⎰dx x X dx x X x X x X lll λ由于,故有0)()0(==l X X ,2'2 0()(())llX x dx X x dx λ=⎰⎰.'22 0(())()0llX x dxX x dx λ=≥⎰⎰(1.8)当时,方程的通解为. 利用边界条件0λ=0)()("=+x X x X λ12()X x c c x =+可得,即. 因此,不是特征值.0)()0(==l X X 120c c ==()0X x =0λ=当时,方程的通解为0λ>0)()("=+x X x X λ. (1.9x C x C x X λλsin cos )(21+=)利用边界条件确定常数如下0)()0(==l X X 21,C C , ,10C =l C l C λλsin cos 021+=或.0sin 2=l C λ由于要求(1.7)中齐次微分方程的非零解,故不能为零. 故有2C .0sin =l λ,从而有0> , ,πλn l =1n ≥, .2)(ln n πλ=1n ≥将代入到(1.8)中,并略去任意非零常数得n C C λ,,212C , .x ln x X n πsin)(=1n ≥故特征值问题(1.7)的解为, , 2(l n n πλ=x ln x X n πsin )(=1n ≥(1.10)注1 特征值问题是分离变量法的理论基础. 上面已求出特征值问题(1.7)的解为. 在高等数学中知道,在一定条件下区间{ sin 1 }n x n lπ≥的任一函数可按特征函数系展开为Fourier 级数. 换言[0 , ]l { sin 1 }n x n lπ≥之,特征函数系是区间上满足一定条件的函数所成无穷维空间的一组基,{ sin 1 }n x n lπ≥[0 , ]l 而且还是该空间上的一组正交基,即有. 特征函0sinsin 0 , ln m x n m l lππ=≠⎰数系的这两个根本性质:正交性和完备性(基),和定理1.1{ sin1 }n x n lπ≥有限维空间中相应结论很相似,只是现在的特征值和特征函数是无穷个. 另n R 外,若改变(1.7)中的边界条件,其相应的特征值和特征函数也会有所变化.如将边界条件变为,则特征值和特征函数分别为(0)0,'()0X X l ==. 2(21)(21)(),()sin ,022n n n n X x x n l lππλ++==≥该特征函数系也具有和特征函数系类似(21){ sin1 }2n x n l π+≥{ sin 1 }n x n lπ≥的性质,既正交性和完备性.此类问题的一般结果便是著名的Sturm—Liouville定理,有兴趣的同学可参阅参考文献.[1][4]-将以上的结果以定理的形式给出.定理1.2 考虑二阶线性微分算子的特征值问题[1],[4]22d A dx=- "()()()()0 , 0 ,(0)0,()0 .k m X x X x x l X X l λ⎧+=<<⎪⎨==⎪⎩(1.11)其中. 则该问题的特征值非负,且满足0,1k m ≤≤.120......n λλλ≤<<<<→∞相应的特征函数系在上是相互正交的. 且对于任一在区间上1{()}n n X x ≥[0,]l [0,]l 分段光滑的函数,可按特征函数系展开为如下的级数()f x 1{()}n n X x ≥Fourier ,1()()n n n f x f X x ∞==∑其中系数为Fourier .20()(), 1()l nn lnf x Xx dxf n Xx dx =≥⎰⎰为后面需要,下面再求解二阶线性微分算子带有周期边界条件的22d A dx=-特征值问题. 在偏微分方程教材中,习惯上用表示周期函数,即考虑下面()θΦ二阶线性微分算子的周期边值问题22d A dx=- "()()0, () (2), .θλθθθπθθ⎧Φ+Φ=-∞<<+∞⎨Φ=Φ+-∞<<+∞⎩(1.12)可证(1.12)和以下问题等价"''()()0, 02(0) (2), (0) (2).θλθθπππ⎧Φ+Φ=≤≤⎪⎨Φ=ΦΦ=Φ⎪⎩(1.13)和(1.8)的证明相似易得(1.13)中的特征值.当时,0≥λ0λ=, 由周期边界条件可得. 所以为特征函数.12()c c θθΦ=+20c =0()1θΦ=当时,方程通解为0λ>,θλθλθsin cos )(21c c +=Φ求导得.'()c c θΦ=-+由周期边界条件可得112cos(2sin(2c c c c c c ππ⎧=+⎪⎨=-+⎪⎩或1212[1cos(2sin(20sin(2[1cos(20.c c c c ππ⎧--=⎪⎨+-=⎪⎩(1.14)由于要求非零解,故不能同时为零. 因此,齐次方程组(1.14)的系数矩12,c c 阵行列式必为零,即 .解之可得1cos(20-=,2n n =λ()cos sin .n n n c n d n θθθΦ=+此时对每个正特征值,特征函数有二个,既,. 总结所得2n n =λθn cos θn sin 结果为如下定理.定理1.3 考虑二阶线性微分算子带有周期边界条件的特征值问22d A d θ=-题"''()()0, 02(0) (2), (0) (2).θλθθπππ⎧Φ+Φ=≤≤⎪⎨Φ=ΦΦ=Φ⎪⎩则该问题的特征值和特征函数分别为,.00,λ=0()1;θΦ=2n n =λ(){cos ,sin }, 1n n n n θθθΦ=≥§22 分离变量法⋅本节结合具体定解问题的求解来介绍分离变量法(method of separation of variables ). 所举例子仅限于一维弦振动方程,一维热传导方程混合问题以及平面上一些特殊区域上的位势方程边值问题. 对高维问题的处理放在其它章节中介绍.以下多数例子均假定定解问题带有齐次边界条件. 否则,可利用边界条件齐次化方法转化之. 我们以弦振动方程的一个定解问题为例介绍分离变量法.2.2.1 弦振动方程定解问题例2.1求解两端固定弦振动方程的混合问题2(,), 0, 0 (2.1)(0,)0, (,)0, 0 (2.2)(,0)(), (,0)(),0. tt xx t u a u f x t x l t u t u l t t u x x u x x x l ϕψ-=<<>==≥==≤≤ (2.3)⎧⎪⎨⎪⎩解 分四步求解.第一步 导出并求解特征值问题. 即由齐次方程和齐次边界条件,利用变量分离法导出该定解问题的特征值问题并求解.令,并代入到齐次方程中得)()(),(t T x X t x u =,0)()()()(''2''=-t T x X a x X t T 或.''''2()()()()X x T t X x a T t =上式左端是的函数而右端是的函数,要二者相等,只能等于同一常数.x t 令此常数为-,则有λ , ,λ-=)()("x X x X "2()()T t a T t λ=-上面的第一个方程为.0)()("=+x X x X λ利用齐次边界条件(2.2),并结合得0)(≠t T .0)()0(==l X X 由此便得该定解问题的特征值问题为"()()0, 0(0)()0.X x X x x l X X l λ⎧+=<<⎨==⎩其解为特征值:特征函数: 2() , 1 ;n n n lπλ=≥()sin, 1 .n n X x x n lπ=≥第二步 正交分解过程. 即将初值和自由项按特征函数系展成{}1()n n X x ≥Fourier 级数,并将也用特征函数表出.),(t x u {}1()n n X x ≥ ,11()()sinn n n n n n x X x x lπϕϕϕ∞∞====∑∑(2.4), 11()()sinn n n n n n x X x x lπψψψ∞∞====∑∑(2.5), 11(,)()()()sinn n n n n n f x t f t X x f t x lπ∞∞====∑∑(2.6)(2.711(,)()()()sinn n n n n n u x t T t X x T t x lπ∞∞====∑∑)这里,和分别为,和的Fourier 系数,具体表示如n ϕn ψ)(t f n )(x ϕ)(x ψ),(t x f 下,02()sin l n n d l l πϕϕααα=⎰,02()sin l n n d l l πψψααα=⎰,02()(,)sin l n n f t f t d l lπααα=⎰而为待定函数.)(t T n 第三步 待定系数法. 即先将和的Fourier 级数代入到(2.1)),(t x f ),(t x u 中,导出关于满足的常微分方程. 再利用初值条件(2.3)得出满足)(t T n )(t T n 的初始条件.假设(2.7)中的级数可逐项求导,并将(2.6)和(2.7)代入到(2.1)中得,"2"111()()()()()()nnnnn n n n n T t Xx aT t Xx f t X x ∞∞∞===-=∑∑∑,"2111()()()(())()()nnn nnn n n n n T t Xx aT t Xx f t X x λ∞∞∞===--=∑∑∑ . (2.8"211(()())()()()nn n n n n n n T t a T t X x f t X x λ∞∞==+=∑∑)由于Fourier 展式是唯一的,比较(2.8)两端系数得)(x X n(2.9"2()()(), 1.n n n n T t a T t f t n λ+=≥)在(2.7)中令并结合(2.4)得0=t (2.10()(0)()()n n n n n n x T X x X x ϕϕ∞∞====∑∑)比较(2.10)两端系数得)(x X n(0), 1.n n T n ϕ=≥(2.11)类似地可得'(0), 1.n n T n ψ=≥(2.12)结合(2.9),(2.11)和(2.12)便得出关于满足的二阶常系数非齐)(t T n (1)n ≥次方程初始值问题"2'()()(), 0(0), (0).n n n n n n n n T t a T t f t t T T λϕψ⎧+=>⎪⎨==⎪⎩(2.13)第四步 求解关于的定解问题(2.13),并将其结果代入到(2.7)中)(t T n 即可.为简单起见,我们设. 将代入到(2.13)中可得方程的通()0,1n f t n =≥n λ解为, t lan d t l a n c t T n n n ππsin cos)(+=利用初始条件确定常数如下,n n c d.'(0), (0)n n n n nn aT c T d lπϕψ====故有. ()cossin n n n l n a n a T t t t l n a lψππϕπ=+最后将上式代入到(2.7)中便得定解问题(2.1)—(2.3)的解为12(,)()sin cos sin l n n n a n u x t d t xlll lπππϕααα∞==∑⎰ (2.14)012()sin sin sin l n n n a n d t x n a l l l πππψαααπ∞=+∑⎰注1 利用分离变量法求解(2.1)—(2.3),需要假设在(2.7)中可通过无穷求和号逐项求导. 而通过号求导要对无穷级数加某些条件,在这里就∑∑不做专门讨论了. 今后遇到此类问题,我们均假设一切运算是可行的,即对求解过程只作形式上的推导而不考虑对问题应加什么条件. 通常称这样得出的解为形式解. 验证形式解是否为真解的问题,属于偏微分方程正则性理论的范围. 一般地讲,偏微分方程定解问题的解大多数是以无穷级数或含参变量积分形式给出的. 对这两类函数可微性的研究需要较深的数学知识,也有一定的难度,有兴趣的同学可查阅参考文献和. 我们约定:本书只求定解问题的形式解.[1][2]注2 当时,由(2.14)可以看出:两端固定弦振动的解是许多(,)0f x t =简单振动的叠加,当时,对任意的(,)()sinn n n u x t T t x l π=(11)k klx x k n n==≤≤-时刻,,即在振动的过程中有个点永远保持不动,所t (,)0n k u x t =(,)n u x t (1)n +以称这样的振动为驻波,而称为该驻波的节点.显然当k x 时,在这些点上振幅最大,称这些点为驻波的21(11)2k x l k n n+=≤≤-sin 1x =腹点. 因此,求特征函数实际上就是求由偏微分方程及边界条件所构定的系统所固有的一切驻波. 利用由系统本身所确定的简单振动来表示一些复杂的振动,便是分类变量法求解波动问题的物理解释.注3 例2.1的求解方法也叫特征函数法(eigenfunction method ),现已成为固定模式,也具有普适性. 初学者似乎会感到有些繁琐,但随着进一步的学习,同学们就会熟练掌握这一方法. 特征函数法的关键之处是求解偏微分方程定解问题相应的特征值问题,而基本思想就是笛卡尔(Descartes )坐标系的思想.如在三维空间中,每个向量可由基的线性组合表出,两个向量3R {,,}i j k 111222 , a i b j c k a i b j c kαβ=++=++相等当且仅当在基下两个向量的坐标相等. 既.{,,}i j k121212 , , a a b b c c ===与此相类似,在例2.1求解中也是比较方程或初始条件两边的系数而得()n X x 到(2.13). 与三维空间相比较,例2.1中特征函数系相当3R { sin1 }n x n lπ≥于3R 中的基,而也就相当于上面的,即定解问题的解{,,}i j k{ T () 1 }n t n ≥111{,,}a b c 关于基函数的坐标. 因此,在具有可数基的无穷维空间中,特{ sin1 }n x n lπ≥征函数法也称为待定系数法.例2.2 设有一均匀细弦,其线密度为. 若端为自由端,端固ρ0x =x l =定.初始速度和初始位移分别为零,并受到垂直于弦线的外力作用,其单位长度所受外力为. 求此弦的振动. sin t ω 解 所求定解问题为(2.1521 sin , 0, 0(0,)0, (,)0, 0(,0)0, (,0)0, 0.tt xx x t u a u t x l t u t u l t t u x u x x l ρω-⎧-=<<>⎪==≥⎨⎪==≤≤⎩)利用特征函数法求解该问题.情形1 非共振问题,即.22, 0n a n ωλ≠≥ 该定解问题的特征值问题为(2.16)"'()()0, 0(0)0, ()0.X x X x x l X X l λ⎧+=<<⎪⎨==⎪⎩其解为, , 2(21)()2n n l πλ+=(21)()cos 2n n X x x lπ+=0n ≥将按特征函数展开成Fourier 级数得1sin t ρω-{}0)(≥n n x X , (2.17)11sin ()()n n n t f t X x ωρ∞==∑.021214()sin sin sin sin 2(21)l n n n f t t d t f t l l n ωπααωωρπρ+===+⎰令(,)()()n n n u x t T t X x ∞==∑(2.18)完全类似例2.1的求解过程可得,对于任意满足下面问题0, ()n n T t ≥(2.19"2'()()sin , 0(0)0, (0)0.n n n n n n T t a T t f t t T T λω⎧+=>⎪⎨==⎪⎩)初值问题(2.19)中齐次方程的通解为,12()cos sin n T t c c =+而非齐次方程的一个特解为.22()sin nn n f T t t a ωλω=-因此,(2.19)的通解为. 1222()cos sin sin nn n f T t c c t a ωλω=++-(2.20)由初始条件可确定出120, c c ==最后将所得到的代入到(2.18)中便得(2.15)的解.()n T t 情形2 共振问题,即存在某个 使得.0,n ≥22n a ωλ=不妨假设.此时,在情形1中求解所得到的不变.220a ωλ={ T () 1 }n t n ≥当时,要求解以下问题0n = "2000'00()()sin , 0(0)0, (0)0.T t T t f t t T T ωω⎧+=>⎪⎨==⎪⎩(2.21)(2.21)中齐次方程通解为.012()cos sin T t c t c t ωω=+为求得非齐次方程的一个特解,要将(2.21)中方程的自由項换为,而求0i t f e ω以下问题的一个特解"2000()().i t T t T t f e ωω+=令并代入到上面非齐次方程中可得 ,故有()i t T t Ate ω=02f iA ω=-,00()sin cos 22f t f tT t t i t ωωωω=-取其虚部便得(2.21)中方程的一个特解为. 00()Im(())cos 2f tT t T t t ωω==-结合以上所得结果便可得到(2.21)中方程的通解为,0012()cos sin cos 2f tT t c t c t t ωωωω=+-由初始条件确定出 ,由此可得01220, 2fc c ω==.0002()sin cos 22f f tT t t t ωωωω=-将代入到(2.18)中便得在共振条件下(2.15)的解为()n T t 000102112(,)()()()()()()(sin cos )cos ()()222 (,)(,) .n n n n n n n n n u x t T t X x T t X x T t X x f f t t t x T t X x l u x t u x t πωωωω∞=∞=∞===+=-+=+∑∑∑可以证明: 是有界的. 而在的表达式中取 ,则2(,)u x t 1(,)u x t 2k k t πω=中的基本波函数的振幅当逐渐变大时将趋于无穷大,最1(,)u x t cos2x lπ0()k T t k 终要导致弦线在某一时刻断裂,这种现象在物理上称为共振. 注意到在上面求解过程中我们取周期外力的频率等于系统的第一固有频率ω波函数分量上发生共振. 一般地讲,当周期外力的频率很接近或等于系统的ω某个固有频率时,系统都会有共振现象发生,即弦线上一些点的振幅将随着时间的增大而不断变大,导致弦线在某一时刻断裂.2.2.2 热传导方程定解问题例2.3 求解下面热方程定解问题(2.2220, 0, 0 (0,), (,)sin , 0(,0)0, 0.t xx x u a u x l t u t u u l t t t u x x l ω⎧=<<>⎪==≥⎨⎪=≤≤⎩)解 利用特征函数法求解(2.22).首先将边界条件齐次化,取,并令,则0(,)sin w x t u x t ω=+w u v -=(2.22)转化为(2.2320cos , 0, 0 (0,)0, (,)0, 0(,0), 0.t xx x v a v x t x l t v t v l t t v x u x l ωω⎧-=-<<>⎪==≥⎨⎪=-≤≤⎩)利用分离变量法可得(2.23)的特征值问题为"()()0, 0(0)0, '()0.X x X x x l X X l λ⎧+=<<⎨==⎩特征值和特征函数分别为,2(21)()2n n lπλ+=0≥n .(21)()sin 2n n X x x lπ+=0≥n 将,按特征函数展成Fourier 级数(,)cos f x t x t ωω=-0)(u x -=ϕ{}0)(≥n n x X 得, (2.24)cos ()()n n n x t f t X x ωω∞=-=∑,02(21)()(1)cos sin cos 2l n n n f t t d f t l lπωαωααω+=-=⎰其中. 1228(1)(12)n n l f n ωπ+-=+ , (2.25)00n n n u X ϕ∞=-=∑其中.00042(21)()sin 2(12)l n u n u d l l n πϕααπ-+=-=+⎰令(2.26)(,)()(), n n n v x t T x X x ∞==∑并将(2.26)代入到(2.23)中的方程得,'2"()()()()cos ()nnnnn n n n n T t Xx aT t Xx f tX x ω∞∞∞===-=∑∑∑.'2(()())()cos ()nn nnn n n n T t a T t Xx f tX x λω∞∞==+=∑∑在(2.26)中令并结合(2.25)得0=t .()(0)()()n n n n n n x T X x X x ϕϕ∞∞====∑∑比较上面两式中特征函数的系数便得()n X x(2.27'2()()cos , 0(0).n n n n n n T t a T t f t t T λωϕ⎧+=>⎪⎨=⎪⎩)(2.27)是一阶常系数常微分方程初值问题.齐次方程通解为.t a n n Ce t T λ2)(-=令,并利用待定系数法求特解可得()cos sin n T t A t B t ωω=+ ,2242242()cos sin n n nn n na f f T t t t a a λωωωωλωλ=+++故有(2.2822242242()cos sin n a tn n nn n na f f T t Cet t a a λλωωωωλωλ-=++++)在上式中代得0t =,2242n nn na f C a λϕωλ=++ . 2242n nn na f C a λϕωλ=-+最后将(2.28)代入到(2.26)中便得(2.23)的解为.0(21)(,)()sin2n n n v x t T t x lπ∞=+=∑故(2.21)的解为),(),(),(t x w t x v t x u +=0 (,)sin v x t u x t ω=++其中由(2.28)给出. )(t T n2.2.3 平面上位势方程边值问题考虑矩形域上Poisson 方程边值问题1212(,), , (,)(), (,)(), (,)(), (,)(), .xx yy u u f x y a x b c y d u a y g y u b y g y c y d u x c f x u x d f x a x b +=<<<<⎧⎪==≤≤⎨⎪==≤≤⎩(2.29)我们假设或. 否则,利用边界条件齐次化方法0)()(21==x f x f 0)()(21==y g y g 化非齐次边界条件为齐次边界条件. 当然,也可以利用叠加原理将(2.29)分解为二个问题,其中一个关于具有齐次边界条件,而另一个关于具有齐次边x y 界条件.例2.4 求解Dirichlet 问题(2.300, 02, 0 1 (0,)0, (2,)0, 01(,0)1, (,1)(1), 0 2.xx yy u u x y u y u y y u x u x x x x +=<<<<⎧⎪==≤≤⎨⎪==-≤≤⎩)解 令并将其代入到(2.29)中齐次方程得)()(),(y Y x X y x u =,0)()()()(""=+y Y x X y Y x X ,λ-=-=)()()()(""y Y y Y x X x X (2.31"()()0, 0 2(0)0, (2)0.X x X x x X X λ⎧+=<<⎨==⎩)0)()("=-y Y y Y λ(2.32)(2.31)便是(2.30)的特征值问题,其解为, , .2)2(πλn n =x n x X n 2sin)(π=1≥n 将代入到(2.32)中得n λ ,0)()("=-y Y y Y n λ(2.33)该方程有两个线性无关解,. 由于,也是(2.33)的y n e2πy n e2π-2n shy π2n ch y π解且线性无关,故(2.33)通解为.y n ch d y n shc y Y n n n 22)(ππ+=令(2.34)11(,)()()()sin 222n n n n n n n n n u x y X x Y y c shy d ch y x πππ∞∞====+∑∑则满足(2.30)中方程和关于的齐次边界条件. 利用关于的边界条),(y x u x y 件可如下确定,,n c n d ,∑∞==12sin1n n x n d π . (2.35))1(1(22sin12220n n n d n d --=⨯=⎰πααπ),x n n ch d n shc x x n n n ∑∞=+=-12sin )22()1(πππ . 22))1(1(22)1(416)1(163322ππππππn sh n chn n sh n n c n nnn -------=(2.36)故(2.30)解为(2.371(,)()sin ,222n n n n n n u x y c shy d ch y x πππ∞==+∑)其中,由(2.36)和(2.35)确定.n c n d 对于圆域,扇形域和圆环域上的Poisson 方程边值问题,求解方法和矩形域上的定解问题无本质区别,只是在此时要利用极坐标.同学们自己可验证:令,作自变量变换,则有θρcos =x θρsin =y .θθρρρρρu u u u u yy xx 211++=+令,将其代入到极坐标下的Laplace 方程中得)()(),(θρθρΦ=R u 222330216(1)164(1)(1)sin ,2222n nn n n n n n c sh d ch d n ππππααααπ----+=-=⎰,"'"211()()()()()()0R R R ρθρθρθρρΦ+Φ+Φ=,"'"211(()())()()()0R R R ρρθρθρρ+Φ+Φ=,"'"21()()()1()()R R R ρρθρλθρρ+Φ=-=-Φ故有, (2.380)()("=Φ+Φθλθ). (2.390)()()('"2=-+ρλρρρρR R R )方程(2.38)结合一定的边界条件便得相应定解问题的特征值问题,而(2.39)是欧拉(Euler )方程. 对(2.39)作自变量变换可得s e =ρ , ,s e =ρρln =s ,'1s dR dR ds R d ds d ρρρ==.2222'''2222211()ss s d R d R ds dR d s R R d ds d ds d ρρρρρ=+=-将以上各式代入到(2.39)得. (2.40''0ss R R λ-=)例2.5 求下面扇形域上Dirichlet 问题(2.4122220, 0, 0, 4(,0)0, 0 2(0,)0, 0 2 (,), 4. xx yy u u x y x y u x x u y y u x y xy x y ⎧+=>>+<⎪=≤≤⎪⎨=≤≤⎪⎪=+=⎩)的有界解.解 令,作自变量变换,(2.41)转化为θρcos =x θρsin =y(2.42)2110, 0, 0 2 2(,0)0, (,0, 022(2,)2sin 2, 0.2u u u u u u ρρρθθπθρρρπρρρπθθθ⎧++=<<<<⎪⎪⎪==≤≤⎨⎪⎪=≤≤⎪⎩令代入到(2.42)中的方程,并结合边界条件可得)()(),(θρθρΦ=R u"()()0, 0<</2(0)0, (/2)0.θλθθππ⎧Φ+Φ=⎨Φ=Φ=⎩(2.43). (2.440)()()('"2=-+ρλρρρρR R R )(2.43)便是(2.42)的特征值问题.求解特征值问题(2.43)可得, , .224)2/(n n n ==ππλθθn n 2sin )(=Φ1≥n 将代入到(2.44)中,并令作自变量变换可得n λs e =ρ,"240ss R n R -=.2222()ns ns n n n n n n n R c e d e c d ρρρ--=+=+由于是求(2.42)的有界解,故有,即. 从而有∞<)0(R 0=n d .n n n c R 2)(ρρ= 上面求出的对每个都满足(2.42)中的方程和齐(,)()()n n n u R ρθρθ=Φ1n ≥次边界条件,由叠加原理得, (2.45∑∑∞=∞==Φ=1212sin )()(),(n n n n n n n c R u θρθρθρ)也满足(2.42)中的方程和齐次边界条件.为使(2.42)中的非齐次边界条件得以满足,在(2.45)中令得(2,)2sin u θθ=2ρ= ,212sin 22sin 2n n n c n θθ∞==∑(2.46)比较上式两边特征函数的系数得θθn n 2sin )(=Φ , .112c =1)( 0≠=n c n 将,代入到(2.45)中便得(2.42)的解为1c 1)(≠n c n . θρθρ2sin 21),(2=u 例2.6 求解圆域上Dirichlet 问题2110, 0, 02(,)(), 02.u u u a u a ρρρθθρθπρρθϕθθπ⎧++=<<≤<⎪⎨⎪=≤≤⎩(2.47)解 圆域上的函数相当于关于变量具有周期. 令(,)u ρθθ2π并代入到(2.46)中的方程可得)()(),(θρθρΦ=R u(2.48"()()0() (2).θλθθπθ⎧Φ+Φ=⎨Φ=Φ+⎩). (2.490)()()('"2=-+ρλρρρρR R R )(2.48)是定解问题(2.47)的特征值问题. 由定理1.3知(2.48)的解为.2, ()cos sin , 0n n n n n c n d n n λθθθ=Φ=+≥将代入到(2.49)中可得(要利用自然边界条件)n λ(0,)u θ<∞,,00)(c R =ρn n n c R ρρ=)(1≥n 利用叠加原理可得(2.47)的如下形式解.∑∞=++=10)sin cos (),(n n n n n d n c c u θθρθρ(2.50)根据边界条件得)(),(θϕθ=a u ,01()(cos sin )n n n n c a c n d n ϕθθθ∞==++∑其中,2001()2c d πϕττπ=⎰,⎰=πτττϕπ20cos )(1d n a c n n .⎰=πτττϕπ20sin )(1d n a d n n 将以上各式代入到(2.50)中便得(2.47)的解为2 2 0 0111(,)()()(()cos cos 2n n u d n d n a ππρρθϕττϕτττθππ∞==+∑⎰⎰ .)sin sin )(12 0 ⎰+πθτττϕπn d n (2.51)注4 利用等式可将(2.51)化为如下形)Re()(cos 1)(1∑∑∞=-∞==-n in n n n e c n c τθτθ式(2.522222201()()(,),22cos()a u d a a πρϕτρθτπρρθτ-=+--⎰)式(2.52)称为圆域上调和函数的Poisson 公式. 在后面学习中还将用其它方法导出它. 注5 在例2.5和例2.6中,如果方程中自由项不为零,若),(θρf 特殊,可用函数代换将自由项化为零而转化齐次方程. 对于一般的),(θρf ,要利用特征函数方法求解.),(θρf 注6 上面例2.3—例2.6几个定解问题的求解思想和主要过程,是伟大的数学家和物理学家Fourier 给出的,详细内容见参考文献. 在这部著名论著[5]中,Fourier 首次利用偏微分方程来研究热问题,并系统地介绍了分离变量法的基本思想和主要步骤. 结合本节所举例子,请同学们小结一下在本章所学过的特征值问题,二阶常系数非齐次常微分方程和欧拉方程的求解方法. 习 题 二1. 设有如下定解问题2(,), 0, 0 (0,)0, (,)0, 0(,0)(), (,0)(), 0.tt xx x t u a u f x t x l t u t u l t t u x x u x x x l ϕψ⎧-=<<>⎪==≥⎨⎪==≤≤⎩利用分离变量法导出该定解问题的特征值问题并求解.2.求解下列特征值问题 (1) "''()()0, 0 (0)()0.X x X x x l X X l λ⎧+=<<⎪⎨==⎪⎩ (2) "()()0, 1 1 (1)0,(1)0X x X x x X X λ⎧+=-<<⎨-==⎩ (3) "()()0, 0 '(0)0, ()0.X x X x x l X X l λ⎧+=<<⎨==⎩ (4) "()()0, 02 (0)(2), '(0)'(2).X x X x x l X X l X X l λ⎧+=<<⎨==⎩3 考虑下面特征值问题*"()()0, 0 (0)0, '()()0.X x X x x l X X l X l λ⎧+=<<⎨=+=⎩(1)证明一切特征值0.λ>(2)证明不同的特征值对应的特征函数是正交的.(3)求出所有的特征值和相应的特征函数.4. 设在区间一阶连续可导且 考虑如下特(),()p x q x [0,]l ()0,()0.p x q x >≥征值问题[()()]()()(), 0 (0)0, ()0.d d p x X x q x X x X x x l dx dx X X l λ⎧-+=<<⎪⎨⎪==⎩(1)证明一切特征值0.λ≥(2)证明不同的特征值对应的特征函数是正交的.5.求解下列弦振动方程的定解问题(1)20, 0<, 0(0,)0, (,)0, 0(,0), (,0)0, 0.tt xx x x t u a u x l t u t u l t t u x x u x x l ⎧-=<>⎪==≥⎨⎪==≤≤⎩ (2) 20, 0<, 0(0,)0, (,)0, 035(,0)sin , (,0)sin , 0.22tt xx x t u a u x l t u t u l t t u x x u x x x l l l ππ⎧⎪-=<>⎪==≥⎨⎪⎪==≤≤⎩(3) 240, 0<1, 0(0,)0, (1,)0, 0(,0), (,0)0, 0 1.tt xx t u u u x t u t u t t u x x x u x x ⎧-+=<>⎪==≥⎨⎪=-=≤≤⎩(4) 242sin , 0<, 0(0,)0, (,)0, 0(,0)0, (,0)0, 0.tt xx x x t u u u x x t u t u t t u x u x x πππ⎧--=<>⎪==≥⎨⎪==≤≤⎩(5) 22, 0, 0 (0,) (,)0, 0(,0)0, (,0), 0.tt xx x t u a u x l t u t u l t t u x u x A x l ⎧-=<<>⎪==≥⎨⎪==≤≤⎩6.求解下列热传导方程的定解问题(1) 2cos , 0<, 02(0,)1, (,), 0(,0)0, 0<.t xx x x u a u x t u t u t t u x x ππππ⎧-=<>⎪⎪==≥⎨⎪=<⎪⎩(2) 22, 0<1, 0(0,)0, (1,)0, 0(,0)sin , 0< 1.t xx x u a u u x t u t u t t u x x x π⎧-=<>⎪==≥⎨⎪=<⎩(3) 220, 0<, 0(0,)0, (,)0, 0(,0)(), 0.t xx u a u b u x l t u t u l t t u x x x l ϕ⎧-+=<>⎪==≥⎨⎪=≤≤⎩(4) 2, 0, 0 (0,)0, (,)0, 0(,0)1, 0.t xx x x u a u xt x l t u t u l t t u x x l ⎧-=<<>⎪==≥⎨⎪=≤≤⎩7. 求解下面位势方程定解问题(1) , 0, 0 (,0)0, (,)0, 0(0,)0, (,), 0.xx yy y y u u x x a y b u x u x b x a u y u a y Ay y b +=<<<<⎧⎪==≤≤⎨⎪==≤≤⎩(2)22220, 0, , 4 (,0)0, 02, (,)0, 0(,), 4.xx yy u u y x y x y u x x u x x x u x y x y x y ⎧+=>>+<⎪⎪=≤≤=≤≤⎨⎪=++=⎪⎩(3) 22220, 4 (,)1, 4.xx yy u u x y u x y x x y ⎧+=+<⎪⎨=++=⎪⎩(4) 222222, 1< 4 (,)0, 1 (,), 4.xx yy u u xy x y u x y x y u x y x y x y ⎧+=+<⎪⎪=+=⎨⎪=++=⎪⎩8 设在区间的Fourier 展开式为 *()x ϕ[0,]l 1()sin ,k k k x x c l πϕ∞==∑(6.1)其部分和为 求解或证明以下结果.1()sin ,n n k k k x S x c l π==∑(1)设,求.()[0,]x C l ϕ∈20[()()]l n x S x dx ϕ-⎰(2)证明下面贝塞尔(Bessel )不等式 22012().l k k c x dx l ϕ∞=≤∑⎰(6.2)(3)设,的二阶导数的Fourier 展开式为2()[0,]x C l ϕ∈()x ϕ1''()sin ,n n n x x d l πϕ∞==∑如果 ,利用分部积分法证明(0)()0l ϕϕ==2, 1,n n d An c n =≥(6.3)其中为正常数.A (4)利用(6.2)和(6.3)证明(6.1)中的三角级数在区间上一致[0,]l 收敛,并且可以逐項求导.9 考虑如下定解问题* 2, 0, 0 (0,)0, (,)0, 0(,0)(), 0.t xx x x u a u x l t u t u l t t u x x x l ϕ⎧=<<>⎪==≥⎨⎪=≤≤⎩(1)给出该定解问题的物理解释.(2)当经过充分长的时间后,导热杆上的温度分布如何?(,)u x t (3)求极限.lim (,)t u x t →+∞10 考虑如下定解问题*2, 0, 0 (0,), (,), 0(,0)(), 0.t xx x u a u x l t u t A u l t B t u x x x l ϕ⎧=<<>⎪==≥⎨⎪=≤≤⎩(1)给出该定解问题的物理解释.(2)求极限.lim (,)t u x t →+∞11 考虑下面定解问题 *20, 0<, 0(0,)(,)0, 0(,0), (,0)0, 0.tt xx t t u u u u x t u t u t t u x x u x x πππ-++=<>⎧⎪==≥⎨⎪==≤≤⎩(1)解释该定解问题方程中各项的物理意义.(2)推导出问题的特征值问题并求解.(3)写出该问题解的待定表示式并求出表达式中第一特征函数的系数.12 考虑下面定解问题 * (,), 0<, 0(0,)(,)0, 0(,0)(), (,0)(), 0.tt xx x x t u u f x t x t u t u t t u x x u x x x ππϕψπ-=<>⎧⎪==≥⎨⎪==≤≤⎩(12.1)(1)写出该定解问题的特征值和特征函数 ,(),0.n n X x n λ≥(2)如果,而,求解该定解问题.()0,()0x x ϕψ==(,)f x t t =(3)如果,证明 ,下面等式(,)0f x t =0τ∀>,222200[(,)(,)][()()]l l t x x u x u x dx x x dx ττψϕ+=+⎰⎰(12.2)成立,解释该等式的物理意义.(4)证明(12.1)的解是唯一的.。
分离变量法——数学物理定解问题
分离变量法是求解偏微分方程最基本和 常用的方法。
理论依据:线性方程的叠加原理和 Sturm-Liouville 理论。 基本思想:将偏微分方程的求解化为对 常微分方程的求解
2.1 有界弦的自由振动
研究两端固定均匀的自由振动. 定解问题为:
2 2u u 2 0, 0 x l 2 a 2 t x t0 u x 0 0, u x l 0, u u t 0 ( x ), ( x ), 0 x l t t 0
④ 成立 X (0) 0, X ( l ) 0
2.1 有界弦的自由振动
则
X '' X 0 ⑤ X ( 0 ) 0, X ( l ) 0
参数
特征值问题 称为特征值.
函数X(x)称为特征函数 分三种情形讨论特征值问题的求解
2. 1 有界弦的自由振动
则无穷级数解 n at n at n x u( x , t ) ( An cos l Bn sin l ) sin l 为如下混合问题的解
n1
utt a 2 uxx 0 0 xl u xl 0 u x 0 0 0 xl u t 0 ( x ) u 0 xl t t 0 ( x )
特征方程 r 2 pr q 0
p 4q r1, 2 , 2 (1) 有两个不相等的实根 ( p 2 4q 0) r 1 , r2
特征根
2
p
两个线性无关的特解
得齐次方程的通解为
y1 e ,
r1 x
y2 e ,
r2 x
y C1e
r1 x
数学物理方程谷超豪版第二章课后答案.doc
第二章热传导方程§ 1热传导方程及其定解问题的提1. 一均匀细杆直径为 l ,假设它在同一截面上的温度是相同的,杆的表面和周围介质发生热交换,服从于规律dQ k 1(u u 1 )dsdt又假设杆的密度为,比热为 c ,热传导系数为 k ,试导出此时温度 u 满足的方程。
解:引坐标系:以杆的对称轴为x 轴,此时杆为温度u u( x,t) 。
记杆的截面面积 l 2为 S 。
t 到 tt 内流入截面坐标为 x 到 xx 一小段细杆的热量为 4由假设,在任意时刻dQu s t k u2u s x tkxs t k1x x x xx 2 xt 到 tt 在截面为杆表面和周围介质发生热交换,可看作一个“被动”的热源。
由假设,在时刻x 到 xx 一小段中产生的热量为4k 1dQ2k 1 u u l x tu u s x t1l1又在时刻 t 到 tt 在截面为 x 到 xx 这一小段内由于温度变化所需的热量为dQc u x,tt u x,t s x c u s x t由热量守恒原理得:3t tcu s x t k2us x t4k 1u u s x tt tx2 xl1消去 sx t ,再令x 0 , t 2 u 0 得精确的关系:cuk 4k 1 u ut x 2 l1u k 2u 4ka 22 u4k或t cx2c 1u u 1x2c 1u u 1ll其中a2kc2. 试直接推导扩散过程所满足的微分方程。
解:在扩散介质中任取一闭曲面s ,其包围的区域 为 ,则从时刻 t 1 到 t 2 流入此闭曲面的溶 质,由 dMDudsdt ,其中 D 为扩散系数,得nt 2D udsdtMt 1 snt 2t 2C udvdtM 1C u x, y, z, t 2 u x, y, z, t 1 dxdydzCudtdvt 1tt 1t两者应该相等,由奥、高公式得:t 2uuut 2C udvdtMD D D dvdt M 1t 1xx y y z zt 1t其中 C 叫做孔积系数 =孔隙体积。
定解问题讲解
Mathematical Methods for Physics第二篇数学物理方程Mathematical Equations for Physics要想探索自然界的奥秘就得解微分方程。
-牛顿中心:将物理问题翻译成数学语言 目的:1、如何用数理方程研究物理问题2、如何导出方程3、能正确写出定解问题§ 6.1 引言Introduction第六章 定解问题Mathematical Problem1、数学物理方程概念:数学物理方程是指从物理、工程问题中,导出的反映客观物理量在各个地点、时刻之间相互制约关系的一些偏微分方程。
数学物理方程 ♣ 线性方程♦♥ 非线性方程一、数理方程简介:§ 6.1 引言一、数理方程简介§ 6.1 引言ttu =a2⊗u +fut=D⊗u +f2、数理方程的产生和发展:(1)十八世纪初期(2)十九世纪中期三类数学物理方程:波动方程u -波动,a-波速,f-与源有关的函数输运方程u -浓度,D-系数,f -与源有关的已知量泊松方程h-与源有关的已知量,u-表示稳定物理量+fxx2Taylor :utt=a u⊗u =-h一、数理方程简介:§ 6.1 引言a u2、数理方程的产生和发展:(3)十九世纪末到二十世纪初高阶方程(梁的横振动):utt= 2xxxxf ( x, t )非线性方程KdV:ut+σuux+uxxx= 0∂ψh2schro&-dinger:i h∂t=-Δψ2μ+U(r)ψ+1、写出定解问题♣ 泛定方程:数理方程(一般规律)♦♥ 定解条件:初始、边界、衔接条件(个性)如:y '(t) - 4 y = 0♣y ' -4y = 0 -泛定方程♠y(0) = 0 ↔ y = C e 2t+ C e -2t♦ ← -定解条件 12-通解♠♥y '( 0) = 4↑♦1、写出定解问题2、求解:求解方法: 行波法、分离变量法、积分变换法、格林函数法、保角变换法、复变函数法、变分法 ♣ 物理意义3、分析解答:♠♠ ♣存在 ♠♥ 适定性 ♦唯一♠♥稳定数学物理方法物理(内容)桥梁数学(成果)、数理方法的特点三 § 6.1 引言。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图2-2(a)
图2-2(b)
(1)由x=c点的初始位移,求出其它点的初始位移。设。由于[0,c]
段是直线,设,由x=0时=0和x=c时=h定出A=h/c和B=0。同样,
[c,l]段也是直线,设,由x=c时,=h和x=l时,=0,定出和,综上得
另外,由题意,显然。 (2)先求初始位移,设拉力与张力T平衡时点x=c得位移为h,由 ,所以
, 而x点的应变则是
, 图2-1 弹性杆的微小纵振动 由于振动是微小的(不超过杆的弹性限度),由胡克定律有
。 设杆的横截面为S(设为常数),则由牛顿第二定律,段的运动方程是
其中常数满足,并利用了胡克定律式,而且将函数在处展开为泰勒级数 并取了前两项。以除上式的两端后, 令取极限,得到
记
则方程最后变为
2、长为的均匀杆,侧面绝缘一端温度为零,另一端有恒定热流进入 (即单位时间内通过单位截面积流入的热量为),杆的初始温度分布 是,试写出相应的定解问题。
3、一均匀杆的原长为,一端固定,另一端沿杆的轴线方向拉长而静 止,突然放手任其振动,试建立振动方程与定解条件。
4、把下列方程化为标准型。 5、 若是两个任意二次连续可微函数,验证 满足方程
的热量守恒问题(设杆的侧面是绝热的),有
其中S为杆的截面积,常数满足,将(1)代入,并令取极限,得到 (3)
类似地,对于小段,有 , 满足。将(1)代入,并令取极限,得到 (4)
(3) 所谓牛顿冷却定律是指,在介质边界面S上,热流强度地外法向分量是 同介质温
与周围介质温度之差成正比,即 (5)
其中的比例常数b>0称为热交换系数。在x=l处,的方向就是x的正方向,所 以(5)成
计算得原方程为
即
(3) 由 ,所以特征方程为
(Ⅰ)若,则,该方程为双曲型的。
特征线为和,或写为及。
令
(1)
(2)
原方程化为
(3)
将减式得, 所以 ,代入得标准型为
(Ⅱ)若,则,这时方程为椭圆型的。 其特征方程为 特征线为 及。 于是令 计算可得原方程化为
或
即
。
(4) 因为,所以该方程是椭圆型的,其特征方程为
x=l端
为绝热(即热流为零);(2)两端均有热流流入;(3)两端以牛顿冷却定
律与周围介质(其
温度恒为)进行热交换。
解:(1) 显然。设热流强度为,由热学中的Fourier定律,热流强
与介质的温度的关系是
,
(1)
是x方向的单位矢量),将(1)代入,记得
ห้องสมุดไป่ตู้
(2)
(2)设x=0处的热流强度为,x=i处的热流强度为。考虑杆的小段
性的,故可设任何方向上单位长度上所受的张力T也是常数。在非边界的膜 上任取一小块
到,现在来分析此小块横向上的受力和运动情况: 沿平行于r的方向,设张力与平衡位置的夹角为,有 在r边上受力 ,
边上受力 , 沿平行于的方向,设张力与平衡位置的夹角为,有
边上受力 , 边上受力 , 在小振动近似下, 由此并根据牛顿第二定律可
,
(4)
,
(5)
。
(6)
由方程(5)、(6)消去p后所得得方程再与方程(4)消去,即得
,
(7)
,由方程(3)、(7)和,得到
,
(8)
由于流体是无旋的,即 , 必存在速度势,使得
(9)
这样,求三个未知函数归结为求一个未知函数u。方程(5)和(6)消去p
后,将
(9)代入,有
(10)
方程(4)与(10)消去s后,将(9)代入,得到
(3)
其中k和是常数,它们由流体的性质决定。
原则上五个方程(1)-(3)((1)是矢量方程,它相当于三个标
量方程)可以确定五个未知函数和p,然而它们是非线性的,其求解已
超出本课程的范围。这里我们设法将它们线性化。
设平衡时流体的密度为,引入稠密度以代替。因为声波在传播过程
中,都是小量,略去高阶小量后,方程(1)-(3)成为
故定解问题为
2 解:杆的初始温度分布是,即有初始条件:
由杆的一端温度为零,有边界条件:
杆的另一端有恒定热流,即:
故定解问题为:
3 解:以表示小段的质心位移,设为杆的横截面积,为杆的密 度,是在点处所受的力。
由牛顿第二运动定律,有 当时, ,,又因为 ,故有
令 ,可得振动方程为 由题意放手时即是振动的初始时刻,此时杆振动的速度为零,即 而端拉离平衡位置使整个杆伸长了,故整个系统的初始位移为 再看边界条件,一端固定即该端没有位移,有 另一端由于放手任其振动时未受外力,有 故振动方程与定解条件为:
4 解 (1) 因为,所以该方程是椭圆型的,其特征方程为
特征线为 和。故可令,
为计算方便,又令
原方程化为 改变自变量的记号得:
(2) 由于 (Ⅰ)若,则,该方程为双曲型。 其特征方程为 ,和 特征线为 和。故可令,, 所以原方程化为
(Ⅱ) 若,则,该方程为椭圆型。 其特征方程为 ,和 特征线为 和。故可令,,为计算方便,又令
。在x=0段,的方向与x的正方向相反,所以有 。
例6 试将方程
化成标准型。 解:因为,所以该方程是双曲型的,其特征方程为
特征线为 和。令, 。 例7 讨论方程 的类型,并化为标准型. 解:所以
故所给方程为双曲型的,其特征方程为 即 解之有 特征曲线为
令 则原方程化为
1、长为的弦两端固定,开始时在受冲量的作用,试写出相应的定解问 题。
, 由此得,利用(1)小题的结果,得到
同样有。 (3)由题意有。在点x=c,由冲量定理,有,即是
弦的质量密度);在其它点处。因此。 说明: 也可以用极限的形式写出。将冲量I看成均匀分布在小段
为小量)内,则在此小段上单位长度受到的冲量是,于是,由冲 量定理,得
在计算结果最后取。
导出长为l的杆的一维热传导问题下的边界条件:(1)x=0端保持恒温,
第二章 数学物理定解问题
2.2.1 例题分析
2.2 基础训练
例1 导出均匀弹性杆的微小纵振动方程,设杆的弹性模量(杆伸长 单位长度所需要的力)为E,质量密度为,作用于杆的外力密度为 F(x,t)
解:取x轴沿杆的轴线方向,以u(x,t)表示x点,t时刻的纵向位移。 使用微元法,考虑杆上的一小段的运动情况。以记杆上x点,t时刻的应 力(杆在伸缩过程中各点相互之间单位截面上的作用力),其方向沿x 轴,现在求杆上x点,t时刻的应变(相对伸长)。如图(2。1)所示, A’B’表AB段(平衡位置)在t时刻所处的位置,则AB段的相对伸长是
特征线为 和。故令, 原方程变为
, 或
即: 5 验证:由题意有 及 及
可得 6证明:由题意只有验证 是方程的解即可。将代入方程(1)左边
式中,并利用条件是方程(1)的解,可得 证毕。
6、验证线性奇次方程的叠加原理。即若 均是线性二阶齐次方程 (1)
的解,其中都只是的函数,而且级数收敛,其中为任意常数,并对可以 逐次微分两次,求证仍是原方程的解。
2.2.3 解答与提示 1 解:设弦的两端为:,由题意有 弦的振动方程为
定解条件为:
在点,由冲量定理有 ,即 ;在其它处,或写成
注意:也可用极限的形式写出:将冲量看成均匀分布在小段 (是无穷 小量),则在此小段上单位长度受到的冲量是,于是,由冲量定理得
这就是杆的纵振动方程,也是一维波动方程。
例2 导出无旋流体中声波的传播方程。
解:运用规律法,直接从现有的理想流体动力学方程和热力学的物
态方程出发导出声波方程。
设流体密度为,空间中一点(x,y,z)的流体元的速度为, 压强为,而连续
性方程是
;
(1)
理想流体的动力学方程是
;
(2)
此外,声波传播过程是绝热过程,它的物态方程是
(11)
方程(7)、(8)和(11)都是声波的传播方程。如果都不
随t变化,则方程(7)、(8)和(11)都变成Laplace方程:
,
(12)
,
(13)
,
(14)
这些方程给出的是稳定场方程。
例3 均匀、各向同性的弹性圆膜,沿圆周固定,试列出膜的横振动方程和
边界条件。
解:设圆膜的质量面密度为为常数(因为圆模是均匀的),又因为圆膜是 弹性各向同
除方程两边,并且令取极限,得 , ,
。 边界条件为: 。 对弦的横振动问题导出下列情况下的初始条件:
(1)弦的两端点x=0和x=l固定, 用手将弦上的点x=c (0<c<l)拉 开使之与平衡位置的偏离为h(图2。2(a),并设h<<l),然后放手;
(2)弦的两端点x=0和x=l固定,用横向力拉弦上的点x=c (0<c<l), 达到平衡后放手(图2。2(b)); (3)弦的两端点x=0和x=l固定,以槌击弦上的点x=c (0<c<l) 使之获得冲量