傅里叶分析应用于热传导问题

1、傅里叶分析应用于热传导问题(物理系郭素梅指导教师陆立柱)摘要傅里叶分析是一种重要的数学工具，本文综述了用傅里叶分析解决细杆的热传导问题，并进行了讨论。傅里叶分析包括傅里叶级数和傅里叶积分，用傅里叶级数法解决有界细杆的热传导问题，用含参数的傅里叶变换法解决无界细杆的热传导问题，比其它方法更系统，体现出一种数学与物理对应的美感。关键词傅里叶级数傅里叶积分傅里叶变换细杆的热传导问题引言1822年，傅里叶在研究热传导问题时，创造了傅里叶分析，随着时代的进步，这一数学工具被广泛地应用于信号分析、匹配滤波、图象处理等方面，掌握这种具有广泛用途和发展前景的工具是十分必要的.热传导是历来研究的热点，尤其是随着

2、计算机电子设备的高集成化发展，机器内发热部件和集成电路元件的发热量随之增加，传统的强制冷方式已不能达到理想效果，因此，热传导设计成了重要问题。万变不离其宗，为了更好地掌握傅里叶分析，为了更好地掌握热传导问题，本文就一维热传导问题对傅里叶分析作了全面详尽的论述。1 .傅里叶分析1.1 傅里叶级数傅里叶级数在应用上有以下优点1:能表示不连续的函数、周期函数，能对任意函数作调和分析。若函数f(x)以21为周期，即f(x21)=f(x)2(1.1.1)则可取三角函数族1,cos-,cos2x11n二xcos1sin,sinsin(1.1.2)作为基本函数族，将f(x)展开为级数3-'n二xn二

3、xf(x)=a0+%(ancos+bncos)n111(1.1.3)可以证明，函数族(1.1.2)是正交完备的。根据三角函数族的正交性，可求得(1.1.3)中的展开系数为其中an=2，f(")cos千d之bn=1/()sinnj-d(1.1.4)(n=0)(n=0)(1.1.3)称为周期函数f(x)的傅里叶级数展开式，其中的展开系数(1.1.4)称为傅里叶系数。关于傅里叶级数的U敛性问题，有Dirichlet定理4。若周期函数是奇函数，则由傅里叶系数计算公式(1.1.4)可见，及诸氏均等于零，展开式(1.1.3)为n二x，、f(x)=£bnsin,。.1巧)n：1l这叫做傅

4、里叶正弦级数。由于对称性，其展开系数为1ln二.bn=jJf(-)sin1d-(1.1.6)同理，若周期函数是偶函数，则二n二x_f(x)=a0+，ancos(1.1.7)nWl这叫做傅里叶余弦级数，其中，11n二,、an=Jf(-)cosd(1.1.8)二nll对于只在有限区间，例如在(0,l)上有定义的函数f(x)，可采取延拓的方法，使其成为某种周期函数g(x),而在(0,l)上，g(x)三f(x)。然后再对g(x)作傅里叶级数展开，其级数和在区间(0,l)上代表f(x),由于f(x)在x=0和x=l无定义，因此可以有无数种延拓方式，因而有无数种展开式，它们在(0,l)上均代表f(x).有

5、时，对函数f(x)在边界(区间的端点)上的行为提出限制，即满足一定的边界条件，这常常就决定了如何延拓。例如要求f(0)=f(l)=0这时应延拓为奇的周期函数，因为n二xn二x,sinxz0=0,sinIx=t=0;又如要求f(0)-f(l)=0这时应延拓为偶的周期函数，因为余弦级数的和的导数在x=0和x=l为零。对于函数u(x,t),-l<x<l,t>0,展开为傅里叶级数时，可将t视为参数，仅关于x展开为傅里叶级数u(x,t)=a0(t)+x(an(t)cosnxbn(t)sinnx)(1.1.9)n4ll其中的展开系数不是常数，而是关于t的函数，an(t)=-u(,t)co

6、s-n-d(1.1.10)、nll1 i.bn(t)=i.u(,t)sin1.2傅里叶积分一般说来，定义在区间(-8<x<oo)上的函数f(x)是非周期的，不能展开为傅里叶级数。为了研究这样的函数的傅里叶展开问题，我们采取如下办法：试将非周期函数f(x)看作是某个周期函数g(x)于周期218时的极限情形。这样，g(x)的傅里叶级数展开式n二x,.n二x、g(x)-a0+”(ancosbnsin)n=4ll在l-8时的极限形式就是所要寻找的非周期函数的傅里叶展开。仔细研究这一极限过程4,可以得到:f(x)=0A()cosxd：'，IB()sinxd(1.2.1)其中1A(co

7、)=jiB(co)=一冗QOf(8f(1.2.2)(1.2.1)右边的积分称为傅里叶积分,(1.2.1)称为非周期函数f(x)的傅里叶积分表达式。(1.2.2)称为f(x)的傅里叶变换式。对f(x)的条件，有傅里叶积分定理5。复数形式的傅里叶积分为f(x)=-QOF(co)eixdo(1.2.3)F(co)=2二二.ix*一f(x)eixdx(1.2.4)1.3含参数的傅里叶变换对于函数u(x,t),(-的傅里叶展开类似可得：oo<x<8,t>0),可将t视为参数,仅将x成为臼变量)则与一元函数f(x)其中u(x,t)=F(3,t)i-xedw(1.3.1)F(3,t)=6八

8、*uu(x,t)ec<xdx-CO(1.3.2)(1.3.1)是u(x,t)傅里叶积分表达式，(1.3.2)是u(x,t)的傅里叶变换式。2.细杆的热传导问题由于温度不均匀，热量从温度高的地方向温度低的地方转移，这种现象叫做热传导。在细杆的热传导问题中研究的是温度在一维空间中的分布和在时间中的变化u(x,t)。应用热传导定理和能量守恒定律，可导出6可导出热传导方程：2utauxx=0(无热源、汇)2ut-auxx=f(x,t)(有热源、汇)还需初始条件u(x,t)|t40=：(x)和三类边界条件口：第一类u(x,t)|x%=W(t)第二类ux(x,t)|xx="(t)第三类u(

9、x,t)|x=xo+Hux(x,t)|x=xo=w(t)这样构成完整的一维热传导问题8。根据空间变量的范围可分为以下两种细杆的热传导问题。2.1有界细杆的热传导问题这里仅选第二类边界条件作讨论，构成220,、ut-auxx=f(x,t)ux|x=0=0ux|xd：0口匕=*(x)2.2无界细杆的热传导问题2ut-auxx=f(x,t)u|t王：=0(0：二x；l,t0)(2.1.1)(-二：x:二二,t0)',7(2.2.1)对半无界细杆的热传导问题，根据边界条件延拓到无界，转化为无界细杆的定解问题。对第一类齐次边界条件的定解问题(x>0,t>0)2ut-auxx=f(x,

10、t)Ulx-0u|t'二：(x)作奇延拓2Ut-aUxx-f(x,t)(x)U|t'-U(-x)对第二类边界条件2Ut-aUxx=f(x,t)(x0)(x:二0)(x>0,t>0)作偶延拓Ux|xz0-0U|t/=：(x)2Ut-aUxx=f(x,t).=(x)U|t'-V(-x)(x0)(x:二0)3.傅里叶分析应用于细杆的热传导问题3.1用傅里叶级数法解决有界细杆的热传导问题傅里叶级数法是直接求解非齐次方程的定解问题。对问题(2.1.1),把所求解U(x,t)本身展开为傅里叶级数，基本函数族应是相应齐次方程q-a2Uxx=0在第二类齐次边界条件下的本征函

11、数cosn(0,1,2,)，这样试把所求解展开为傅里叶余弦级数lU(x,t)=COzn=0nxTn(t)cosl(3.1.1)把这个级数代入泛定方程,8'Tn(t)n=0n2n2a2nHx(3.1.2)+ATn(t)cosT=f(x,t)方程左边是傅里叶余弦级数，提示我们把方程右边也展开为傅里叶余弦级数，得到:oO”Tnn旦222-n二an二xn二x2-Tncos-='、fn(t)cos-l1nz0l(3.1.3)其中fn(t)为f(x,t)的傅里叶余弦级数的第微分方程n个傅里叶系数。比较两边的系数，分离出Tn(t)的常222(3.1.4)(3.1.5)一'n刀a-Tn

12、-2Tn=fn(t)又把(3.1.1)代入初始条件，得：Tn(0)COS=：(X)=；COS-n卫ln=0l由于基Tn(t)的非其中Q为5(x)的傅里叶余弦级数的第n个傅里叶系数。(3.1.5)式两边都是傅里叶余弦级数,本函数族cos-的正交性，等式两边对应同一基本函数的傅里叶系数必然相等，于是得l零初始条件1,.To(O)=°<,J()dTn(0)=：n=j.O()COS十d(3.1.7)Tn(t)的常微分方程(求解9)在初始条件(3.1.7)下的解是n2壬2tn2二2a2t2t22tTn(t)=e,fn(t)edtn-.fn(t)dt(3.1.8)这样所求解是2_222_2

13、2二2atn2at"Xu(x,t)=ve'fn(t)e'dtn一.fn(t)dtcos丁n=0l(3.1.9)可以证明(3.1.9)是存在且唯一的10.3.2用傅里叶变换法求解无界细杆的热传导问题对问题(2.2.1)应用含参数的傅里叶变换，即用不着遍乘方程及定解条件各项，并对空间变数积分(时间变数视作参数)，原来的定解问题变成(3.2.1)U'(t;k)k2a2U(t;k)=F(t;k)U(t;k)|«=02.2x其中U(t;k)为u(x,t)的傅里叶变换。为求解这个非齐次常微方程，用e遍乘方程各项,dk2a2tk2a2t-U(t;k)ekat=F(

14、t;k)edt对t积分一次，计及零初始值，,22t22U(t;k)=e*LF&k)edzt.22.=J；fffieee(tT)d：di进行傅里叶逆变换，t2-2,u(x,t)=f(,)ee*a(t)dd.?eikxdk2-：0-二交换积分次序u(x,t)=j丛匕e*2a2(Jeik(x/dkdM引用积分公式-.2二e=2k2e%k=e«-0°a可得结果(xFt二一.12(t).u(x,t)=!.f(,)e-)dd.(3.2.2)12；2aj：(t-)可以验证11(3.2.2)确实符合(2.2.1).有热源或热汇的热传导问题，即泛定方程是齐次的，求解更容易。4.讨论4

15、.1 一维热传导问题方法和结论的推广用傅里叶分析法解决细杆的热传导问题，以及得到的结论均可推广到二维、三维空间，用到的理论基础是二、三重傅里叶级数13和二14、三15重傅里叶变换，求解过程与一维类似。4.2 傅里叶分析应用于其它定解问题用傅里叶分析法求解热传导问题时，只是对所求解进行了傅里叶展开或变换，并未对方程限制，常见的其它定解问题13:振动问题，扩散问题等均可用傅里叶分析法。参考文献1近藤次郎等.微分方程付里叶分析.于溶渤译者沈阳：辽宁人民出版社，19812董延.级数.上海：上海科学技术出版社，19823周肇锡.积分变换.国防工业出版社，19824梁昆淼.数学物理方法（第三版）.北京：高

16、等教育出版社，19985管平等.数学物理方法.北京：高等教育出版社，20016郭敦仁.数学物理方法.北京：人民教育出版社，19657陆全康等.数学物理方法自学辅导.上海：上海科学技术出版社，19898杨应辰徐明聪.数学物理方法与特殊函数.国防工业出版社，19809四川大学数学系高等数学教研室.高等数学（第三版）.北京：高等教育出版社，199510TynMyint-U.数学物理中的偏微分方程.徐元钟译.上海：上海科学技术出版社，198311陈庆益.数学物理方程.人民教育出版社，197912陆立柱.数学物理方法.太原：山西高校联合出版社，199313周祥龙.数学物理方程.浙江大学出版社，19911

17、4孙仲康.快速傅里叶变换及其应用.北京：人民邮电出版社，198215C.哈普尔.数学物理引论.肖布森译.北京：科学出版社，1989FourieranalysisapplicationtoHeat-conductionquestion(DepartmentofPhysicsGuoSumeiDirectorLuLizhu)AbstractFourieranalysisisanimportantMathematicstool.ThethesisappliesFourieranalysistoone-dimensionalspaceHeat-conductionquestion,andhaveadiscussion.FourieranalysisconsistsofFourierserieswhichcansolvelimitedone-dimensionalspaceHeat-conductionquestionandFourierintegralwhichcansolveinfiniteone-dimensionalspaceHeat-conductionquestio