热传导方程

2018-11-14 22:23

2.1初边值问题的求解

?ut?a2uxx?0(t?0,0?x?l).........(2.1)?.........................(2.2)?t?0:u??(x).........初边值问题?

.........................(2.3)?x?0:u?0..........?.................(2.4)?x?l:ux?hu?0..........其中h为正常数。

用分离变量求解。令u(x,t)?X(x)T(t),这里X(x)和T(t)分别表示仅与x有关和仅与t有

T,即关的函数。把它代入方程,得到XT??aX??才成立。令此常数为-?,则有

2T?X???.这等式只在两边均等于常数时2aTXT???a2T?0................(2.5)

X????X?0.................(2.6)先考虑(2.6).根据边界条件(2.3)、(2.4),X(x)应当满足边界条件

X(0)?0,X?(l)?hX(l)?0...........(2.7)

对于边值问题(2.6)、(2、7),分类进行讨论。 ①当??0时,只有平凡解X?0; ②当??0时, X(x)?Acos?x?Bsin?x.........(2.8)

利用边界条件X?0??0,得A?0.于是(2.7)的第二个边界条件得到

?2.9? B(?cos?l?hsin?l)?0.................为使X(x)为非平凡解,?应满足

?cos?l?hsin?l?0..........(2.10)

即?应是下述超越方程的正解:

tan?l??令 v??h..........(2.11)

?l.................(2.12)

v.........(2.13) lh则(2.11)式变为 tanv??利用图解法或数值求解法可得出这个方程的根。方程(2.13)有可列无穷多个正根

1(k?)??vk?k?.因此,特征值问题(2.6)、(2.7)存在着无穷vk?0(k?1,2,....),满足

2多个固有值

?k?(vk2)(k?1,2......)...........(2.14) lvkx(k?1,2,.....)........(2.15) l及相应的固有函数

Xk?x??Bksin?kx?Bksin把???k代入方程(2.5),可得 Tk(t)?Cke?a2?kt(k?1,2,......).........(2.16)

于是得到一列可分离变量的特解 uk(x,t)?Ake?a2?ktsin?kx(k?1,2,...).........(2.17)

由于方程(2.1)及边界条件(2.3)、(2.4)都是齐次的,故可以利用叠加原理构造级数形式的解

u?x,t???Ake?a?ktsin?kx...............(2.18)

2?k?1以下来决定常数Ak,使(2.18)满足初始条件(2.2) 由(2.2),为使在t?0时u(x,t)取到初值??x?,应成立

??x???Aksin?kx...............(2.19)

k?1?为确定系数Ak,须先证明固有函数系?Xk??sin??kx在?0,l?上正交。

?设固有函数Xn和Xm分别对应不同的固有值?n和?m,即

????nXn?0,Xm??mXm?0. Xn以Xm和Xn分别乘上面第一式和第二式,得到

????nXmXn?0,XmXn

??XnXm??mXnXm?0.两式相减后在?0,l?上积分,有

?(X0lm??-Xm??Xn)dx?(?n??m)?XnXmdx?0. Xn0l再利用Xm和Xn都满足边界条件(2.3)(2.4),就得到

?-XmXn?)(?n??m)?XnXmdx?(XnXm?0 00ll由于?n??m,故得固有函数系的正交性:

?l0XnXmdx??sin?nxsin?mx?0,m?n........(2.20)

0lMk??sin0l2?kxdx??记?lsin2?kll??? 2222?k1?tan?kl4?klh?...........(2.21)222(h??k)1?cos2?kxdx02tan?kl1l?于是,在(2.19)两边乘以sin?kx,在进行积分,利用正交性(2.20)即可得

1Ak?Mk?????sin0L?n?d?.......(2.22)

将它代入(2.18)式,就得到初边值问题(2.1)-(2.4)的形式解为

u?x,t???1L?a2?kt????sin??d??esin?kx......(2.23) k?0k?1Mk?a2?kt?下面证明分离变量法得到的形式解是初值问题(2.1)-(2.4)的经典解。 表达式(2.23)中含有因子e?2k,因此对任意??0,当t??时,对任意p?0级数

p?a?t??ke均是一致收敛的。若?为有界函数及(2.21)式,可得 k?1?????sin?k?d??Ml及0l12?........(2.24) Mkl因而,由(2.23)表示的级数,当t?0时,关于x和t是无穷次可导的,并且求导与求和可以交换。由于级数的每一项都满足方程(2.1)及边界条件(2.3)、(2.4),从而(2.23)式表示的级数在t?0时确实满足方程及边界条件。为了保证当t?0时,对任意的x?[0,l],由(2.23)式给出的级数趋于初值??x?,还需要对??x?加上进一步的条件。例如在

??x??C2,??0??0,???l??h??l??0时,可以证明,由(2.23)式给出的级数确实是初边值

问题(2.1)-(2.4)的经典解。 2.1柯西问题 ⑴傅里叶变换及其性质

-?,?)上的函数,它在[?l,l]上有一阶连续导数,则在设f(x)是定义在( (?l,l)中f(x)可以展开为傅里叶级数a0?n?n?f(x)???(ancosx?bnsinx)..........(2.25)

2n?1ll并且

1ln?1ln?an??f(?)cos?d?,bn??f(?)sin?d?(n?0,1,2...).....(2.26)

l?lll?ll若

?f(x)在(??,?)???上绝对可积,则有

f(x)???10d??f???cos?(x??)d?........(2.27)

积分表达式(2.27)称为f(x)的傅里叶积分。

(2.27)式也可以写成复数形式。由于cos?(x??)是?的偶函数,sin?(x??)是?奇函数,可以将(2.27)式写成

?1?f(x)?d?f(?)[cos?(x??)?isin?(x??)]d?2???????

?1??d??f???ei?(x??)d?..................(2.28)?2?????若令g????就有f?x??称

??..(2.29) ?f???ed?..........i???12?????g???ei?xd?.........(2.30)

g(?)为f(x)的傅里叶变换,记为F[f];称f(x)为g(?)的傅里叶逆变换,记为F?1[g].当

f(x)在(??,?)上连续可导且绝对可积时,它的傅里叶变换存在,且其逆变换等于f(x).傅里叶变换的性质

① 傅里叶变

线

于 任

?,?以及函数f1,f2,成立F[?f1??f2]??F[f1]??F[f2]........(2.31)于f1(x)和f2(x)的傅里叶变换的乘积,②f1(x)和f2(x)的卷积的傅里叶变换等即 F[f1?f2]?F[f1]?F[f2].........(2.32)

③f1(x)和f2(x)乘积的傅里叶变换等于f1(x)和f2(x)的傅里叶变换的卷积乘以1,即 2?F[f1?f2]?1F[f1]?F[f2].........(2.33) 2?④如果f(x),f?(x)都是可以进行傅里叶变换的,而且当x??时,f(x)?0,则成立

F[f?(x)]?i?F[f(x)]................(2.34)

⑤如果f(x)及xf(x)都可以进行傅里叶变换,那么

F[?ixf(x)]?dF[f]..............(2.35) d?(2)齐次热传导方程柯西问题的求解

2??u2?u??a.................(2.36) ??t?x2?.........(2.37)?u(x,0)??(x).........关于x进行傅里叶变换,记

~(?,t),F[u(x,t)]?u~(?). F[?(x)]??在(2.36)两边关于x进行傅里叶变换,利用性质④,就得到

~du~..........??a2?2u....(2.38)

dt类似,由(2.37)式可得

~(?).........~(?,0)??u.....(2.39)

(2.28)、(2.29)是带参数?的常微分方程的柯西问题,它的解为

~(?)e?a2?2t..........~(?,t)??u...(2.40)

函数e?a2?2t的傅里叶变换为

F?1[e?a?t]??12?2212?????e?(a?t?i?x)d??x24a2t22????e?a2t(??ix2)2a2t

d??e利用复变函数的积分计算知

????e?a2t(??ix2)2a2td???e?at?d????22?1at????e?ydy?2?at

.?x24a2t所以F?1[e?a22?t]?12a?te,


热传导方程.doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:北师大版小学二年级数学上册全册单元测试卷

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: