傅里叶级数的应用:热流 [学习笔记]

本文将运用傅里叶变换来解决数学物理方法中的热流问题:
假设我们有一个圆环,初始温度分布为 f ( x ) ,求任意时刻圆环各处的温度情况


求解过程如下

U ( x , t ) 为圆环在 x 位置, t 时刻时的温度
显然 U ( x , 0 ) = f ( x )
为了便于计算,我们设绕圆环一周的总长度为 1 ,则我们有: f ( x + 1 ) = f ( x ) U ( x + 1 , t ) = U ( x , t ) ,即一圈后回到原来的位置
由于函数 U 是周期函数,因此我们可以写出它的傅里叶级数表达式:

U ( x ) = k = C k e 2 π i k x

由于函数 U x t 的函数,而 e 2 π i k x 部分展现了其周期性,因此,系数 C k 应当是 t 的函数,且与 x 无关(否则将失去周期性)。于是我们得到:

(1) U ( x , t ) = k = C k ( t ) e 2 π i k x

由数学物理方法知道,一维热方程(又称扩散方程):

U t = a U x x

其中, a 是由圆环的环境、材质等决定的常数,为了便于计算,我们将其设为 1 2 ,从而得到(注意此处的下标 t x x 表示求导)

(2) U t = 1 2 U x x

( 1 ) 式代入 ( 2 ) 式,我们有

U t = k = C k ( t ) e 2 π i k x

U x x = k = C k ( t ) ( 4 π 2 k 2 ) e 2 π i k x

k = C k ( t ) e 2 π i k x = k = C k ( t ) ( 2 π 2 k 2 ) e 2 π i k x

由于 { e 2 π i k x } 的正交性,等式左右两边的系数必然一一对应,于是我们有:

C k ( t ) = 2 π 2 k 2 C k ( t )

解该常微分方程,有:

d ( C k ) d t = 2 π 2 k 2 C k d ( C k ) C k = 2 π 2 k 2 d t 1 C k d ( C k ) = 2 π 2 k 2 d t ln C k = a 2 π 2 k 2 t a (3) C k ( t ) = C k ( 0 ) e 2 π 2 k 2 t

因而进一步地,我们需要计算 C k ( 0 )

t = 0 时,根据 ( 1 ) 式,有

U ( x , 0 ) = f ( x ) = k = C k ( 0 ) e 2 π i k x

考虑到 f ( x ) 也是周期函数,因此可以写出其傅里叶级数:

f ( x ) = k = f ^ ( k ) e 2 π i k x

对比两式可知, C k ( 0 ) 可以看作 f ( x ) 傅里叶展级数开式的第 k 项系数 f ^ ( k )
根据这个结论,以及 ( 3 ) 式,我们有:

U ( x , t ) = k = f ^ ( k ) e 2 π 2 k 2 t e 2 π i k x

从这个式中,我们可以知道:当 t U ( x , t ) 0
这是本题的解的第一种表达方式


我们继续将 f ^ ( k ) f ( x ) 表示,此处,为了区分于 U ( x , t ) 中的 ( x ) 变量,我们将 f ( x ) 改写为 f ( y ) ,从而有:

f ^ ( k ) = 0 1 e 2 π i k y f ( y ) d y

因此,我们可以最终写出 U ( x , t ) 的表达式:

U ( x , t ) = k = 0 1 e 2 π 2 k 2 t e 2 π i k x e 2 π i k y f ( y ) d y (4) = 0 1 k = ( e 2 π i k ( x y ) e 2 π 2 k 2 t ) f ( y ) d y

这是本题的解的第二种表达方式


此外,通过 ( 4 ) 式,我们可以引入卷积的概念
g ( x , t ) = k = e 2 π i k x e 2 π 2 k 2 t ,那么 U ( x , t ) 可以表示为:

U ( x , t ) = 0 1 g ( x y , t ) f ( y ) d y

这个式子的含义是: U ( x , t ) f ( x ) g ( x , t ) 的卷积
此处, g ( x , t ) 被称为热核函数(也称热方程格林函数或热方程基本解),在其他数学物理方法的题目中,还会看到许许多多不同的“核”函数,他们同样会以卷积的形式构成类似的解

猜你喜欢

转载自blog.csdn.net/qq_22943397/article/details/80303265