数学物理方程课件第二章1分离变量法
【精品课件】数学物理方程分离变量法
l
) 1,2,3,
)
特征值与 特征函数
2u t2
a2
2u x2
,
0xl,t 0
u(0,t) 0, u(l,t) 0,
t 0
u(x,0) (x),
u(x,0) (x),
0 x l
t
T''n(t)a2nl222Tn(t)0
X''(x)X(x)0
T''(t)a2T(t)0
X nn(xn)2 l2B 2nsi(n nn l1x ,2,3(,n )1 ,2,3 ,
驻 波 ： 两 列 反 向 行 进 的 同 频 率 的 波 形 形 成 驻 波 。 波 腹 ： 振 幅 最 大 的 点 ; 节 点 ： 振 幅 最 小 的 点
求方程的通解的步骤为：
(1)写出微分方程的特征方程 r2 0,
(2)求出特征根 r1 , r2，
(3)根据特征根的情况按下表写出所给微分方程 的通解。
特征根
通
解
yC1er1xC2er2x
y(C1C2x)
y ( C 1 c o sx C 2 s inx )
一 有界弦的自由振动
1 求两端固定的弦自由振动的规律
u(x,t)un(x,t) n1
n 1(CncosnlatD nsinnlat)sinnlx (n1,2,3, )
步骤3，其余的定解条件求出系数。
un 1(C nco n la st D nsin ln at)sin lnx
n
u(x,t)t 0u(x,0 )n 1C nsinl x(x)
X(x)AexBex
AB0
AB0 X0
Ae l Be l 0
《数学物理方程》第2章 分离变量法
1) 分离变量 偏→常 常 设定解问题有形如 u(x,t)=X(x)T(t) 的非零解
代入DE得 代入 得 X ( x)T ′′(t ) = a 2 X ′′( x)T (t )
u( 0, t ) = X ( 0 )T ( t ) = 0 代入BC得 代入 得 u( l , t ) = X ( l )T ( t ) = 0 X ′′( x) T ′′(t ) ∴ = 2 = λ(常数) X满足 X ′′( x ) + λ X ( x ) = 0(*) 满足 X ( x ) a T (t ) X ( 0) = X ( l ) = 0 T 满足 T ′′( t ) + a 2 λ T ( t ) = 0 (**)
双曲型方程－ 一 双曲型方程－有界弦的自由振动
utt = a 2uxx (t > 0, 0 < x < l ) 1 定解问题 u(0, t ) = 0, u(l , t ) = 0 (t ≥ 0) u( x,0) = ( x) ut ( x,0) = ψ( x) (0 ≤ x ≤ l )
2 求解步骤
第2章 分离变量法 章
椭型的混合问题和边值 混合问题和边值问题 适用于 双 抛 椭型的混合问题和边值问题 非齐DE非齐 非齐BC 非齐DE齐次 齐次BC 非齐 非齐 令u( x , t ) = v( x , t ) + w( x , t ) 非齐 齐次
utt = a 2 uxx + f ( x, t ) u(0, t ) = u1 (t ) u(l , t ) = u2 (t ) u( x,0) = ( x) ut ( x,0) = ψ( x)
2 n
注 由2.6节Sturm Liouville 理论 , 特征值问题的特征函
数理方程第二章分离变量法
分离变量法得到的解可能不唯一，有时需要额外的条件或参数才能 确定唯一解。
数值稳定性
分离变量法在数值实现时可能存在数值稳定性问题，如数值误差的 累积和扩散等，需要采取适当的措施进行控制和校正。
06
CATALOGUE
分离变量法的改进与拓展
改进方向一：提高求解精度
数值稳定性
通过改进数值算法，提高求解过程中数值的稳定性， 减少误差的传播和累积。
原理推导
01
首先，将偏微分方程中的多个变量分离出来，使方程变为一个 关于各个变量的常微分方程。
02
然后，对每个常微分方程分别求解，得到各个变量的解。
最后，将各个变量的解代回原偏微分方程，得到整个问题的解
03 。
原理应用
在物理学中，分离变量法广泛应用于求解具有多个独立变量的偏微分方程 ，如波动方程、热传导方程等。
高阶近似方法
研究高阶近似方法，以更精确地逼近真实解，提高求 解精度。
自适应步长控制
引入自适应步长控制策略，根据解的精度要求动态调 整步长，提高求解精度。
改进方向二：拓展应用范围
复杂边界条件
研究如何处理更复杂的边界条件，使得分离变 量法能够应用于更广泛的数理方程问题。
多维问题
将分离变量法拓展到多维问题，以解决更复杂 的数学模型。
04
CATALOGUE
分离变量法的实例
实例一：一维波动方程的分离变量法
总结词
通过将一维波动方程转化为常微 分方程，分离变量法能够简化求 解过程。
详细描述
一维波动方程是描述一维波动现 象的基本方程，通过分离变量法 ，我们可以将该方程转化为多个 常微分方程，从而逐个求解，得 到波动问题的解。
数学表达式
《分离变量法》课件
06
总结与展望
总结
内容回顾
详细梳理了分离变量法的基本概 念、应用场景、实施步骤和注意 事项，帮助学习者全面理解这一
方法。
案例分析
通过具体的案例分析，展示了分离 变量法在解决实际问题中的应用， 加深学习者对方法的理解和掌握。
互动问答
鼓励学习者在课程结束前提出疑问 ，并对常见问题进行了解答，有助 于巩固学习效果。
展望
新应用领域
实践应用建议
探讨分离变量法在未来可能的应用领 域，如人工智能、大数据分析等，为 学习者提供新的思路和方向。
为学习者提供将分离变量法应用于实 际问题的建议和指导，帮助他们更好 地实现学以致用。
方法改进
介绍分离变量法的最新研究进展和可 能的改进方向，激发学习者进一步探 索和研究。
谢谢您的聆听
02
分离变量法的原理
原理概述
分离变量法是一种求解偏微分方程的 方法，通过将多个变量分离，将复杂 的偏微分方程简化为一系列简单的常 微分方程，从而求解。
该方法适用于具有多个变量的偏微分 方程，特别是当各变量之间相互独立 时。
数学模型建立
首先，需要建立偏微分方程，并确定变量 的个数。
然后，通过适当的变换，将偏微分方程转 化为全微分方程。
求解过程
通过分离变量法，可以将 $u(x, t) = X(x) T(t)$，从而将波动方程 转化为 $X''(x) = -lambda X(x)$ 和 $T''(t) = -omega^2 T(t)$， 其中 $lambda$ 和 $omega$ 是常数。
应用实例二：化学反应动力学模型
总结词
描述化学反应速率
THANKS
《数学物理方程》分离变量法-精PPT共29页
文 家 。汉 族 ，东 晋 浔阳 柴桑 人 （今 江西 九江 ） 。曾 做过 几 年小 官， 后辞 官 回家 ，从 此 隐居 ，田 园生 活 是陶 渊明 诗 的主 要题 材， 相 关作 品有 《饮 酒 》 、 《 归 园 田 居 》 、 《 桃花 源 记 》 、 《 五 柳先 生 传 》 、 《 归 去来 兮 辞 》 等 。
1
0
、
倚
南
窗
以
寄
傲
，
审
ቤተ መጻሕፍቲ ባይዱ
容
膝
之
易
安
。
31、只有永远躺在泥坑里的人，才不会再掉进坑里。——黑格尔 32、希望的灯一旦熄灭，生活刹那间变成了一片黑暗。——普列姆昌德 33、希望是人生的乳母。——科策布 34、形成天才的决定因素应该是勤奋。——郭沫若 35、学到很多东西的诀窍，就是一下子不要学很多。——洛克
《数学物理方程》分离变量法-精
6
、
露
凝
无
游
氛
，
天
高
风
景
澈
。
7、翩翩新 来燕，双双入我庐 ，先巢故尚在，相 将还旧居。
8
、
吁
嗟
身
后
名
，
于
我
若
浮
烟
。
9、 陶渊 明（ 约 365年 —427年 ），字 元亮， （又 一说名 潜，字 渊明 ）号五 柳先生 ，私 谥“靖 节”， 东晋 末期南 朝宋初 期诗 人、文 学家、 辞赋 家、散
数学物理方程分离变量法研究生高校本科生PPT课件
l n , (n 1, 2,.....)
从而得到了固有值为
n
n2 2
l2
,
(n 1, 2,.....)
相应的固有函数为
(11) 固有值问题
Xn (x)
B sin
n
l
x
,
(n 1, 2,.....)
X ''(x) X (x) 0 (6)
(12)
X
(0)
0,
X
(l)
0
(10)
第12页/共97页
ut |t0
x)
Cn sin
n1
n
l
x
n a
(x) Dn
n1
l
sin
n
l
x
第14页/共97页
§2.2.1 齐次方程定解问题的解法
u
ut
|t0 ( |t0
x)
Cn
n1
sin
n
l
x
n a
(x) Dn
n1
l
sin
n
l
x
这两式正好是(x)和 (x)关于 sin
n
l
x的正弦展开。
根据Fourier级数展开法则(见下页附录)，便可得到
§2.2.1 齐次方程定解问题的解法
iii)求方程满足边界条件的特解。
设u(x,t) X (x)T (t) (4)
为了求出T (t)，把(11)式代入(7)式，得 T ''(t) a2T(t) 0 (7)
T
''(t)
n2 2a2
l2
T
(t)
0
其通解为一对共轭复根，即
n
数学物理方程第二章分离变量法
X ( x ) X ( x) 0
(2.6)
要 求非 零解 u( x , t ) 0, T ( t ) 0, X (0) X ( l ) 0
X ( x ) X ( x) 0 x (0) x (l ) 0
第二步 求固有（值）函数
通解形式
r1 r2
r1 r2
r α iβ
y( x) Ae
r1x
Be
r2 x
y( x) Axer1x Ber1x
y ( x) ex ( A cos βx B sin βx)
X ( x ) X ( x ) 0
(1) 0
X ( x) Ae
即 在 固 有 值 已 知 情 况求 下解 另 一 个 常 微 分 方的 程通 解
第 四 步 合 成 偏 微 分 方 程 的 通并 解求 满 足 初 始 条 件 的解 定
说 明 关 于 求 固 有 值 的题 问： 1)固 有 值 是 否 存 在 ？ 2)固 有 函 数 系 是 否 正 交 ？ 3)函 数 是 否 可 以 按 固 有 数 函系 展 开 ？
X ( x ) T ( t ) 2 X ( x ) a T (t ) (2.5)
X ( x )T (t ) a 2 X ( x )T (t )
u ( x , t ) X ( x )T ( t ) T (t ) a T (t ) 0 (2.7) x 由条件 u x (0, t ) X (0)T (t ) 0, u x ( l , t ) X (l )T (t ) 0,
nat nat nx un ( x, t ) X n ( x )Tn (t ) ( En cos Fn sin ) sin l l l
数学物理方程-第二章分离变量法
第二章 分离变量法分离变量法是求解偏微分方程定解问题最常用的方法之一，它和积分变换法一起统称为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）的解是唯一的.。
数学物理方法技巧分离变量法
目录
CONTENTS
• 引言 • 分离变量法的基本原理 • 分离变量法的具体应用 • 分离变量法的注意事项 • 分离变量法的优缺点 • 分离变量法的未来发展与展望
01
引言
分离变量法的定义
分离变量法是一种数学物理方法，用 于将多变量问题转化为多个单变量问 题，以便于求解。
它通过将偏微分方程转化为常微分方 程，或者将高阶微分方程转化为一系 列在解决具有多个相互独立变量的物理 问题时，如波动、热传导、流体动力 学等，分离变量法是非常有效的工具 。
它适用于具有周期性边界条件或对称 性边界条件的问题，如无限大区域、 周期性结构等。
THANKS
感谢您的观看
结合其他数值方法进行优化
1 2 3
混合方法
结合分离变量法和有限元法、有限差分法等其他 数值方法，形成混合方法，取长补短，提高求解 精度和效率。
自适应方法
结合分离变量法和自适应方法，根据问题特性和 求解需求，动态调整算法参数和求解精度，实现 高效求解。
无网格方法
结合分离变量法和无网格方法，克服传统数值方 法的网格依赖性，提高求解灵活性。
02
偏微分方程的解需要满足一定的边界条件和初始条件。
03
偏微分方程的解需要满足可解性条件，即解需要是有限个变 量的函数。
分离变量法的步骤
01
将偏微分方程转化为常微分方程。
02 对常微分方程进行求解，得到各个变量的解。
03 将各个变量的解组合起来，得到原偏微分方程的 解。
03
分离变量法的具体
应用
一维波动方程的分离变量法
为了确保数值解法的稳定性，可以采用多种方法，如增加计算步长、使用更精确的数值格式等。同时 ，也需要不断尝试和改进计算方法，以提高数值解法的稳定性和准确性。
数学物理方法分离变量法
数学物理方法分离变量法分离变量法是数学物理中常用的一种解微分方程的方法，它适用于一些特定形式的偏微分方程，能够将原方程分解成一系列简单的常微分方程，从而求得方程的解。
在物理学中，分离变量法常常用于描述热传导、波动、量子力学等问题的求解。
本文将介绍分离变量法的基本思想和应用，以及一些实际问题中的案例分析。
首先，我们来看一般形式的偏微分方程：\[F(x,y,u,\frac{\partial u}{\partial x},\frac{\partial u}{\partial y},\frac{\partial^2u}{\partial x^2},\frac{\partial^2 u}{\partial y^2},...) = 0\]其中，\(u = u(x,y)\) 是未知函数，\(F\) 是关于 \(x,y,u,\frac{\partial u}{\partial x},\frac{\partial u}{\partial y},\frac{\partial^2 u}{\partial x^2},\frac{\partial^2u}{\partial y^2},...\) 的已知函数。
我们的目标是求解这个偏微分方程，找到满足条件的 \(u\) 函数。
分离变量法的基本思想是假设未知函数 \(u(x,y)\) 可以表示为两个独立变量 \(x\) 和 \(y\) 的乘积形式，即 \(u(x,y) = X(x)Y(y)\)。
将这个形式代入原方程中，然后通过变量分离的方法，将方程化为两个关于 \(x\) 和 \(y\) 的常微分方程。
最后再对这两个方程分别进行积分，得到原偏微分方程的解。
下面我们通过一个具体的例子来说明分离变量法的应用。
考虑二维热传导方程：\[\frac{\partial u}{\partial t} = k\left(\frac{\partial^2 u}{\partial x^2} +\frac{\partial^2 u}{\partial y^2}\right)\]其中，\(u(x,y,t)\) 表示温度分布，\(k\) 是热传导系数。
数学物理方程- 2 分离变量法PPT课件
第 二
上式要恒成立，只有它们均等于同一个常数，记为 章
因此， 0 不是特征值。
数 学 物 理 方 程
N o 当 0 时，方程 X (x )X (x ) 0 的通解为 X ( x ) c 1 c o s x c 2 s inx
第 二 章 分 离
利用边界条件，确定常数
变
量
即有
Image c10,c2sinl0 sin l 0
法
所以
西安交通大学理学院
的各种各样的定解问题，其基本思想是把偏微分方程分解为几个常
No 数学物微分方程，其中有的常微分方程带有附加条件而构成本特征值问题第二
理 方
章
程
2.1 特征值问题
分 离 变
量
Image 2.1.1 矩阵的特征值问题
矩阵的特征值问题 设A为一n阶实矩阵，其特征值满足
法
一般来说，特征值和线性无关的特征向量不多于n个。任意n阶矩阵 都有n个线性无关的广义特征向量，以此n个线性无关的广义特征向
的一个基。将x，b按此基展开为
，则 变
量
法
等价于Image
或
西安交通大学理学院
由于 线性无关，比较系数有
数则
为原问题的解。
第
学 物 理 方 程
例2 设
No
二 章 分 离
求解非齐次常微分方程组
变 量
法
Image 其中
为已知向量函数，
解 和例1相似 ，将 开
西安交通大学理学院
按基
分别展
则（4）等价于 数
No 数 若原方程的系数均为常数，并分别用小写的
学
物
理
方 程
a,c,d,e, f
数理方程第讲
X(x)lX(x)0. (2.5) 6
再利用边界条件(2.2), 由于u(x,t)=X(x)T(t),
X(0)T(t)=0, X(l)T(t)=0. 但T(t)0, 如果T(t)=0, 这种解称为平凡解, 所 以
X(0)=X(l)=0
(2.6)
因此, 要求方程(2.1)满足条件(2.2)的变量分离
由于方程(2.13)与边界条件(2.14)都是齐次的,
所以 u (x,t)C n e- n 2 a 2 tsinnx (2 .2 2 )
n 1
仍满足方程与边界条件. 最后考虑u(x,t)能否
n
xd x
0
10
2
5n 3
3
(1 -
cos
n
)
0, 当 n为偶数 ,
4
5n 3
3
,
当
n 为奇数
.
23
因此, 所求的解为 u(x,t)
543n 0(2n1 1)3sin(2n1 01)xcos10(2n1)t
24
解题中常用到的积分表的内容:
xsin
axd
x
1 a2
sin
ax
-
1 a
x
cosax
x
(2.11)
16
u(x,t) un(x,t) n1
n1Cn
cosnat
l
Dn
sinnl at sinnl
x
(2.11)
将初始条件(2.3)代入上式得:
u(x,t)|t0u(x,0)n1Cnsinnlx(x)
u
t t0
Dn
n1
nasinn
ll
x(x)
17
复习高等数学中周期为2l的傅立叶级数: 如果周期为2l的周期函数f(x)为奇函数, 则有
数学物理方程第二章分离变量法
第五讲补充常微分方程求解相关知识。
第二章 分离变量法偏微分方程定解问题常用解法，分离变量法。
解常微分方程定解问题时，通常总是先求出微分方程的特解，由线性无关的特解叠加出通解，而后用定解条件定出叠加系数一阶线性偏微分方程的求解问题，基本方法也是转化为一阶线性常微分方程组的求解问题对于二阶以及更高阶的偏微分方程定解问题，情况有些不同：即使可以先求出通解，由于通解中含有待定函数，一般来说，很难直接根据定解条件定出，因此，通常的办法就是把它转化为常微分方程问题（第六讲）§2.1 有界弦的自由振动什么是分离变量法？使用分离变量法应具备那些条件？ 下面通过两端固定的弦的自由振动问题来说明。
定解问题：考虑长为l ，两端固定的弦的自由振动，其数理方程及定解条件为.0 ),(u ),(u 0,,0u ,0u 0, l,0 ,0t0022222l x x x t t x xu a t u t t l x x ≤≤==>==><<∂∂=∂∂====ψϕ分析：1. 方程和边界条件都是齐次的，求这样的问题可用叠加原理。
2. 我们知道，在解常微分方程定解问题时，通常总是先求出微分方程的特解，由线性无关的特解叠加出通解，而后用定解条件定出叠加系数。
启发：能否运用类似求常微分方程定解问题的方法求偏微分方程？也既是能否先找出满足齐次方程及齐次边界条件的足够多的特解，再用其作线性组合使其满足初始条件。
由分析，我们现在试求方程的变量分离形式：)()(),(t T x X t x u =的非零解。
将),(t x u 代入方程，可得)()()()()()()()(2''''''2''x T a x T x X x X t T x X a t T x X =⇒= 此式中，左端是关于x 的函数，右端是关于t 的函数。
因此，左端和右端相等，就必须等于一个与t x ,无关的常数。
数学物理第二章-分离变量法
例1 设 b Rn ，求解线性方程组 Ax b.
4
解 A的n个线性无关的特征向量{Ti}(1 i n) 可以作为 Rn
n
n
的一组基。将x，b按此基展开为 x xi Ti ,b bi Ti，则
Ax b 等价于 n
i1
i1
n
xi ATi bi Ti
i1
i1
或
n
n
xi iTi bi Ti
l n ,n 1
n
n
l
2
,n
1
所以，可得
11
Xn (x)
sin
n
l
x, n
1
因此,特征值问题（1）的解为
n
n
l
2
,n
1,
Xn (x)
sin
n
l
x, n
1.
注：
特征值问题是分离变量法的理论基础；
改变边界条件，相应的特征函数系也会改变；
Sturm-Liouville定理：特征函数系的正交性和完备性。
（3）导出 Tn (t)满足的方程，给出通解（傅里叶展开）；
（4）由初始条件确定通解系数.
注2 对齐次问题
u(x,t) 2 l(s)sin( n s)ds cos n a t sin n x
l0 n1
l
l
l
2
l
(s)sin( n
s)ds sin
0
xi0 i ,
n
f (t) fi T (t)6 i.
i 1
i 1
i 1
则原问题等价于 dx Ax f (t), x(0) x0
dt
T T n dxi
i1 dt
第二章-分离变量法-1
T = F (x )
0 ≤ x ≤L ，τ＝ 0 ＝
解：1．分离函数 ．
假定该问题的解可以分解成空间函数与 时间函数的乘积形式
T ( x,τ ) = X ( x )Γ(τ )
代入微分方程及定界条件，转化为 个常 代入微分方程及定界条件，转化为2个常 微分方程——分离方程 微分方程 分离方程
T ( x,τ ) = X ( x )Γ (τ )
上式所示的解既满足原导热问题的微分方程， 上式所示的解既满足原导热问题的微分方程，又满 足边界条件，但它不一定满足初始条件。因此， 足边界条件，但它不一定满足初始条件。因此，还 需将初始条件应用于上式。 需将初始条件应用于上式。
F ( x) = ∫
∞
β =0
C ( β )[β cos( β x) + H sin( β x)]dβ
数学描述： 数学描述：
h
1
初始时 T=F(x)
1 ∂T ( x,τ ) ∂ 2T ( x,τ ) = a ∂τ ∂x 2
x
0 ＜ x ＜ ∞，τ＞0 ， ＞
O
∂T λ − hT = 0 ∂x ∂x
x ＝0 ，τ＞0 ＞
半无限大物体的导热
T = F (x )
0 ≤x ≤L ，τ＝ 0 ＝
解：1．分离函数 ．
1 d 2 X ( x) 1 d Γ (τ ) = X ( x ) dx 2 a Γ (τ ) d τ
dΓ(τ ) + aβ 2 Γ(τ ) = 0 dτ
1 ∂T ( x,τ ) ∂ 2T ( x,τ ) = a ∂τ ∂x 2
= -β
2
d 2 X ( x) + β 2 X( x) = 0 dx 2
∂X + HX = 0 ∂x
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第2章分离变量法
第二章 分离变量法
一、有界弦的自由振动 二、有限长杆上的热传导 三、拉普拉斯方程的定解问题 四、非齐次方程的解法 五、非齐次边界条件的处理 六、关于二阶常微分方程特征值问题的一些结论
数学物理方程与特殊函数
第2章分离变量法
2u t2
a2
2u x2
,
0xl,t 0
其中： A nC n 2 D n 2
na nl
n a rc ta nC D n n
x=x0时： un(x0,t)A nsinn l x0cos(ntn)
t=t0时： un(x,t0)A ncos(nt0n)sinn l x (n1,2,3, )
sin n x
l
n
2 n
l
l
fn
n 2
na t
na t
T n ( t) C 'n c o sl D 'n s inl ( n 1 ,2 ,3 , ) n a n a n
u n ( x ,t ) ( C n c o slt D n s i n lt ) s i n lx( n 1 ,2 ,3 ,)
u(x,t)un(x,t)
n1
▪确定常数 Cn2 l 0l(x)sin nlxdx D nn 2a0 l(x)sin nlxdx
分离变量法可以求解具有齐次边界条件的齐次偏微分方程。
数学物理方程与特殊函数
第2章分离变量法
2 解的性质
un(x,t)(CncosnlatDnsinnlat)sinnlxAncos(ntn)sinnlx
数学物理方程与特殊函数
第2章分离变量法
X ''( x ) X ( x ) 0 X ( 0 ) 0 ,X ( l ) 0
特征（固有）值问题：含有待定常数常微分方程在一定条
件下的求解问题
特征（固有）值：使方程有非零解的常数值
特征（固有）函数：和特征值相对应的非零解
分情况讨论：
1） 0
X(x)AexBex
C
m
Cm2 l 0l(x)sin m lxdx
Cn2 l 0l(x)sin nlxdx D nn 2a0 l(x)sin nlxdx
数学物理方程与特殊函数
第2章分离变量法
u(2t0u2,t)a20,ux2u2(l,,t)0,
0xl,t 0 t 0
u(x,0)(x),u(xt,0) (x), 0xl
n 1(CncosnlatD nsinnlat)sinnlx (n1,2,3, )
数学物理方程与特殊函数
第2章分离变量法
un 1(C nco n la st D nsin ln at)sin lnx
u(x,t)t 0u(x,0 )n 1C nsinn l x(x) u(x t,t)t0n 1D nnlasinnlx(x)
•适用范围： 波动问题、热传导问题、稳定场问题等
数学物理方程与特殊函数
第2章分离变量法
一 有界弦的自由振动
1 求两端固定的弦自由振动的规律
2u a2 2u,
t2
x2
0xl,t 0
u(0,t) 0, u(l,t) 0,
t 0
u(x,0) (x),
u(x,0) (x),
0 x l
t
令 u(x,t)X(x)T(t)
带入方程： X(x)T''(t)a2X''(x)T(t)
令
X ''(x) T ''(t) X (x) a2T(t)
X ''( x ) X ( x ) 0T ''( t) a 2 T ( t) 0
带入边界条件 X (0 ) T ( t) 0 , X ( l) T ( t) 0
X(0)0, X(l)0
Байду номын сангаас
na 2l
v
fnn
na 2l
2l n
a
T
驻波法
数学物理方程与特殊函数
第2章分离变量法
例1：设有一根长为10个单位的弦，两端固定，初速为零，初
位移为(x)x(10x)1000 ，求弦作微小横向振动时的位移。
u(t20u2,t)10u4(1x2u02,t,)0,
0x10,t 0 t 0
u(x,0)x(11000x)0,u(xt,0)0, 0x10
AB0
AB0
Ae l Be l 0
X0
2） 0 X(x)AxB
AB0 X0
3） 0 令 2 ， 为非零实数 X (x ) A c o sx B s inx
A0 B sin l 0
n (n1,2,3, )
l
n2 2
l2
n
n22
l2
(n1,2,3, )
n
X n(x)B nsinl x (n1 ,2,3 , )
解： u(x,t)X(x)T(t)
u(0 ,t)X (0 )T (t)0 X(0)0
0 lsin 2n lx d x 0 l1 c o s2 2 n /ld x 2 l
0 ls i n n lx s i n m lx d x 1 2 0 l c o s n lm x c o s n lm x d x 0
l 0
(x )sim lnx d x0 ln 1 C nsin ln x sim lnx d x 2l
数学物理方程与特殊函数
第2章分离变量法
2u t2
a2
2u x2
,
0xl,t 0
u(0,t) 0, u(l,t) 0,
t 0
u(x,0) (x),
u(x,0) (x),
0 x l
t
T''n(t)a2nl222Tn(t)0
X''(x)X(x)0
T''(t)a2T(t)0
X nn(xn)2 l2B 2nsi(n nn l1x ,2,3(,n )1 ,2,3 ,
u(0,t) 0, u(l,t) 0,
t 0
u(x,0) (x),
u(x,0) (x),
0 x l
t
•基本思想：
首先求出具有变量分离形式且满足边界条件的特解，然后
由叠加原理作出这些解的线性组合，最后由其余的定解条件
确定叠加系数。
•特点： a.物理上由叠加原理作保证，数学上由解的唯一性作保证； b.把偏微分方程化为常微分方程来处理，使问题简单化。
▪分离变量 u(x,t)X(x)T(t) XX0 Ta2T0
▪▪▪求求求特 另 通征一解值个u和函特数n1征un函数n1T Xnn TnC n c nn 1o (C n n ln s a c /tl o 2n D la n sst Xin D nn (ln xa s)tin Bn lna st) is nnlin xlnx