数学物理方法
解题策略:一定要先看清方程属于哪一类(时间二阶导波动方程、时间一阶导热传导方程、跟时间项无关的泊松方程),给的哪一类边界条件。
行波法(直接求通解,无界问题)
注意:方程中都是二阶偏导才能用达氏公式。
用行波法求解无界弦的强迫振动
无界弦的自由振动问题:即范定方程为齐次方程,采用行波法(又叫达朗贝尔公式、特征线法)。
对特征线法的初步理解:先对范定方程因式分解,然后引入两个新变量作为原范定方程中变量的变量,利用链式法则求解出利用两新变量表示的原范定方程的变量,然后再用原范定方程的两变量表示引入的新变量,实现回代。
对于纯强迫振动问题:(1)先将有源问题按冲量原理化为无源问题;(2)再利用达氏公式求解。
无界:是指空间变量的无界。
行波法的核心思想:先求通解再求特解。所以并不一定得无界才能利用这种思想求解,而是利用特征线法的思想进行变换求解。
一维一般强迫振动的解(叠加原理):
考点:
- 叠加原理
- 达朗贝尔公式(方程齐次)
- 冲量原理(边界齐次)
特别注意:用达式公式时一定是在两类初始条件下才能用,一定要先看清给的何种形式的条件,才能判断是否能用达式公式。
一维纯强迫振动的解(冲量原理)
用行波法求解半无界弦的自由振动
自我提示:自由运动即方程是齐次的。
公式依旧是利用达朗贝尔公式,关键是需要对
x
−
a
t
x-at
x−at的大小进行判断。当
x
−
a
t
x-at
x−at大于等于0时,解与无界弦的自由运动相同;当
x
−
a
t
x-at
x−at小于0时,解要进行相应的变号。
注意:这是由于
u
(
0
,
t
)
=
0
u(0,t)=0
u(0,t)=0。当
u
(
0
,
t
)
=
f
(
t
)
u(0,t)=f(t)
u(0,t)=f(t)时,当
x
−
a
t
x-at
x−at大于等于0时(1)公式不变;当
x
−
a
t
x-at
x−at小于0时,(2)式加上
f
(
t
−
x
/
a
)
f(t-x/a)
f(t−x/a)。根据P175习题9.2的解来猜测的。
分离变量法(直接求特解,有界问题)
思路:两端固定的弦形成驻波,故可用驻波法(即分离变量法)求解。
设原方程的特解为
u
(
x
,
t
)
=
X
(
x
)
T
(
t
)
u(x,t)=X(x)T(t)
u(x,t)=X(x)T(t),代入原方程将范定方程转换为两个分别关于
X
X
X和
T
T
T的方程。同时,将边界条件也转换为关于
X
X
X或
T
T
T的边界条件。
第一步:解本征值问题(先解关于
X
X
X的方程)
注意:分离变量法中求和时
n
n
n的取值从
0
0
0还是
1
1
1开始取决于求出的本征函数,本征函数为正弦函数则从
1
1
1开始,本征函数为余弦函数则从
0
0
0开始。
本征值问题:求齐次方程带有齐次边界条件的本征值和本征函数的问题
二阶齐次常微分方程采用特征线法求解。
二维圆盘的狄氏问题求解
狄氏问题为第一类边界条件
考点:
- 坐标变换
- 分离变量法
- 本征函数法
- 傅里叶级数
- 欧拉型公式的解
- 引入周期性边界条件和有界性边界条件
傅里叶级数(分离变量法求系数时基本都会用到)
设以
2
l
2l
2l为周期的函数
f
(
x
)
f(x)
f(x)在
[
−
l
,
l
]
[-l,l]
[−l,l]上可以展开为傅里叶级数。(二维圆盘狄氏问题采用分离变量法求解本征值问题时引入了周期性边界条件,所以最后得到本征函数形式与傅氏级数形式是相同的)。
分离变量法求解圆盘狄氏问题
注意:欧拉型公式的每一导数都是连续降阶的。
电像法求解圆盘狄氏问题
在物理上,格林函数
G
(
M
,
M
0
)
G(M,M_0)
G(M,M0)表示位于
M
0
M_0
M0点的点源,在一定边界条件下在
M
M
M点产生的场。
注意:
r
r
r与
ρ
0
\rho_0
ρ0和
a
a
a都有关
格林函数
求四分之一平面的狄氏问题
傅氏变换与拉氏变换
性质对比
常用变换式
一阶与二阶常微分方程的解(变换时常涉及)
一阶常微分方程的解
一阶齐次与非齐次常微分方程的解
二阶常微分方程的解
范定方程的解
二阶齐次常微分方程用特征线法,非齐次考虑为 P m ( x ) e λ P_m(x)e^\lambda Pm(x)eλ形式时,当非齐次项中 λ \lambda λ与特征线法解出的特征值不相等时,设一个特解为 c x + b cx+b cx+b代入方程求解出系数 c c c和 b b b。
初值问题的解(拉氏变换)
解析延拓
对于半无界问题,如 u ( 0 , t ) = f ( t ) , t > 0 u(0,t)=f(t),t>0 u(0,t)=f(t),t>0:先将 f ( t ) f(t) f(t)延拓为 f ( t ) H ( t ) f(t)H(t) f(t)H(t),并将 x x x看成参数来求解。 H ( t ) H(t) H(t)为阶跃函数 ,是连接半无界问题的桥梁。
每天增加一点常识
调和级数是发散的。