数学物理方法 课件 09傅里叶级数展开法

目录第9章

傅里叶级数展开法9.1强迫振动的定解问题9.2有源热传导的定解问题9.3泊松方程的定解问题9.4非齐次边界的处理第2篇

数学物理方程9.1强迫振动的定解问题9.1强迫振动的定解问题我们以一个长为l的细杆的振动方程为例，不过这时细杆受到外界强迫力

f(x，t)的作用。此外，我们还假设细杆的两端是自由振动的，即对应于第二类齐次边界条件。这样，细杆振动的泛定方程和定解条件为由上一章的讨论可知，一维齐次方程在齐次边界条件下的本征解为三角函数族，它具有完备正交性。因此，本章求解非齐次泛定方程的基本思想是:将非齐次泛定方程的解u(x，t)按照对应的齐次方程在齐次边界条件下的本征解

Xn(x)进行傅里叶级数展开，即其中展开系数Tn(t)是随时间变化的，由初始条件(9.1-3)来确定。可以看到，非齐次振动方程(9.1-1)所对应的齐次方程的本征解问题为9.1强迫振动的定解问题不难验证，它的本征解为这样可以将非齐次方程(9.1-1)的通解表示为为了确定出展开系数Tn(t)，我们还需要把方程(9.1-1)右边的非齐次项及初始条件(9.1-3)按这个本征函数来展开，即9.1强迫振动的定解问题其中展开系数为9.1强迫振动的定解问题分别把式(9.1-5)-(9.1-8)代入泛定方程(9.1-1)及初始条件(9.1-3)，则可以得到其中kn=nπ/l。这样通过上述级数展开，我们就把原来的一个非齐次二阶偏微分方程的求解问题转化成一个非齐次二阶常微分方程的求解问题。接下来，我们采用第六章介绍的拉普拉斯变换来求解方程(9.1-11)。设与

Tn(t)对应的像函数为

Tn(p)，则借助于拉普拉斯变换，可以把方程(9.1-11)转化为如下代数方程其中

fn(p)是函数

fn(t)对应的像函数。将初始条件(9.1-12)代入方程(9.1-13)，则可以得到然后再进行拉普拉斯反演，可以得到9.1强迫振动的定解问题最后，我们得到非齐次振动方程(9.1-1)的一般解为可以看到，一旦给定非齐次项f(x，t)及初始位移ϕ(x)和初始速度ψ(x)的函数形式，就可以确定出展开系数fn(t)，ϕn

及ψn

的形式，进而可以确定出非齐次振动方程的解u(x，t)。上面介绍的求解非齐次方程的解法是一种较为普遍的方法，可以用来求解任意形式的线性非齐次方程，只是对于不同形式的齐次边界条件，本征函数和本征值有所不同。9.1强迫振动的定解问题求解如下细杆强迫振动的定解问题其中A，ω

及a

均为常数。

9.1强迫振动的定解问题将以上结果代入式(9.1-17)，则可以得到

9.2有源热传导的定解问题9.2有源热传导的定解问题对于有源热传导问题，也可以采用上节介绍的傅里叶级数展开法进行处理。考虑如下细杆的热传导问题

并代入方程(9.2-1)，则可以得到其中9.2有源热传导的定解问题对于初始温度ϕ(x)也做类似地展开，可以得到其中采用拉普拉斯变换方法来求解非齐次常微分方程(9.2-5)，可以得到将初始条件(9.2-7)代入上式，并进行反演，则有最后，可以得到非齐次热传导方程(9.2-1)的一般解为这样，一旦给定了非齐次项

f(x，t)和初始温度ϕ(x)的形式，非齐次热传导方程的解就完全被确定。9.2有源热传导的定解问题在上述非齐次热传导问题中，如果初始温度为零，即ϕ(x)=0，以及非齐次项为

f(x，t)=Asin(ωt)，求细杆上的温度分布。解：由于ϕ(x)=0及

f(x，t)=Asin(ωt)，则有这样有可见在Tn(t)的表示中，第一部分为弛豫项，随着t→∞，很快衰减为零;第二项为外界热源维持的强迫项，随外界热源的变化而变化。将式(9.2-11)代入式(9.2-10)，即可以得到细杆上任意时刻的温度场分布u(x，t)。9.2有源热传导的定解问题下面讨论利用傅里叶级数展开法求解二维有源热传导问题。考虑一个长和宽分别为a

和b的矩形薄板，在外界热源作用下，对应的定解问题为

将上式代入方程(9.2-12)，并利用三角函数的正交性，则可以得到9.2有源热传导的定解问题其中对初始条件(9.2-14)，也可以做类似地展开，有其中利用拉普拉斯变换，很容易得到常微分方程(9.2-16)的解为将式(9.2-19)代入式(9.2-15)，即可以得到二维非齐次热传导方程(9.2-12)的一般解。9.3泊松方程的定解问题9.3泊松方程的定解问题我们再讨论一下二维泊松方程的定解问题。考虑一个长和宽分别为a和b的矩形区域，静电势u(x，y)满足如下泊松方程及边界条件其中

f(x，y)为电荷源的空间分布。考虑到齐次边界条件(9.3-2)，首先可以将函数u(x，y)展开成如下形式其中系数Yn(y)待定。将式(9.3-4)代入方程(9.3-1)及边界条件(9.3-3)，则有9.3泊松方程的定解问题其中再考虑齐次边界条件(9.3-6)，可以进一步地把Yn(y)展开成如下级数形式其中Cnm

为待定系数。把式(9.3-8)代入方程(9.3-5)，则可以确定出展开系数Cnm

为其中9.3泊松方程的定解问题最后可以得到泊松方程的一般解为这样，一旦知道了非齐次项

f(x，y)的具体形式，就可以利用上述方法确定出泊松方程的解。在上面的讨论中，我们假定了所有的边界条件都是齐次的。实际上，对于直角坐标系中的二维泊松方程，只要其中的一对边界条件是齐次的，而另外一对边界条件是非齐次的，就可以利用傅里叶级数展开法进行求解。例如，在x=0及x=a处的边界条件是齐次的，如式(9.3-2);而在y=0及y=b处的边界条件是非齐次的，为这时仍然可以将泊松方程(9.3-1)的解写成式(9.3-4)的形式，其中系数Yn(y)也仍然满足方程(9.3-5)，不过方程(9.3-5)对应的边界条件不再是式(9.3-6)，而是原则上讲，利用常微分方程的求解方法，可以得到方程(9.3-5)在非齐次边界条件下的一般解，进而可以确定出泊松方程的一般解。下面再来讨论圆形区域内泊松方程的定解问题。在平面极坐标系中，齐次边界条件下泊松方程的定解问题为9.3泊松方程的定解问题其中a

为圆的半径。考虑到周期性边界条件u(r，φ)=u(r，φ+2π)，可以将方程(9.3-14)的解按照如下复数形式的傅里叶级数展开将式(9.3-15)代入方程(9.3-14)，则可以得到展开系数An(r)满足的方程为其中对应的齐次边界条件则变为另外，泊松方程在圆心处(r=0)的解应有限，因此有9.3泊松方程的定解问题方程(9.3-16)是一个二阶非齐次常微分方程，原则上讲可以求出它在边界条件(9.3-18)和(9.3-19)下的解。

而当n≠±1时，方程(9.3-16)变为欧拉方程，其解为An(r)=dnrn+enr-n，考虑边界条件(9.3-19)和(9.3-18)，则有dn=en=0，即

An(r)=0(n≠±1)。

最后，可以得到方程(9.3-14)的解为可以看到在平面极坐标系中，当函数

f(r，φ)形式比较复杂时，很难得到方程(9.3-16)的解析解。这说明，当所考虑的坐标系不是直角坐标系时，采用这种傅里叶级数展开法求解泛定方程的解要受到一定的限制。9.4非齐次边界的处理9.4非齐次边界的处理在前面的讨论中，无论方程是齐次的还是非齐次的，我们都假定边界条件是齐次的。那么在非齐次边界条件下，如何确定泛定方程的定解问题呢?由于所讨论的泛定方程的定解问题都是线性的，因此可以采用叠加原理把边界齐次化，其基本思路是:选择一个合适的辅助函数v(x，t)，且令使得关于函数w(x，t)的定解问题具有齐次边界条件。首先我们以细杆的自由振动为例来进行讨论，其定解问题如下其中u1(t)和u2(t)是时间变量的任意函数，所对应的边界条件为第一类非齐次边界条件。辅助函数v(x，t)的选取所要遵循的基本原则是:在保证能够使得边界条件齐次化的前提下，使得v(x，t)的形式最为简单。

对于上述第一类非齐次边界条件，可以选取v(x，t)是空间变量x

的线性函数，即9.4非齐次边界的处理将u(x，t)=v(x，t)+w(x，t)代入定解问题式(9.4-2)-式(9.4-4)，则可以得到关于w(x，t)的定解问题其中可以看到，经过上述处理后，原来的非齐次边界条件变成了齐次边界条件，原来的齐次方程变成了非齐次方程。这样，我们就可以利用

§9.1节介绍的方法来求解该非齐次方程的定解问题。9.4非齐次边界的处理还需要说明两点:(1)尽管我们是以齐次泛定方程(9.4-2)为例来讨论的，但如果在定解问题中，不仅边界是非齐次的，泛定方程也是非齐次的，我们仍然可以按照上面的方法来把边界齐次化。(2)辅助函数的形式依赖于边界条件的类型，如对于第二类非齐次边界条件可以选取辅助函数v(x，t)的形式为而对于“混合”边界条件则对应的辅助函数分别为9.4非齐次边界的处理求解如下定解问题其中A和ω

为常数。解：在这种情况下，可以选取辅助函数为并令u(x，t)=v(x，t)+w(x，t)，则可以把原来的定解问题转化为9.4非齐次边界的处理

则可以得到其中ωn=nπa/l及9.4非齐次边界的处理再利用拉普拉斯变换法来求解常微分方程(9.4-26)，可以得到最后就得到