十章行波法和分离变量法本征值问题课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《十章行波法和分离变量法本征值问题课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 行波 分离 变量 问题 课件
- 资源描述:
-
1、驻波驻波:u=u1+u2=A cos(t-kx)+A cos(t+kx)=2A 波腹波腹波节波节2.557.51012.515-1-0.50.51u1u2 驻波的数学表示驻波的数学表示:u(x,t)=X(x)T(t)。这是一个变量可分离函数这是一个变量可分离函数。驻波的数学描述驻波的数学描述:为简便起见研究为简便起见研究两列反向两列反向行进的同频率的波的叠加行进的同频率的波的叠加 u1=A cos(t-kx),u2=A cos(t+kx)弦在弦在平衡位置的振动可用函数平衡位置的振动可用函数 T(t)描述;描述;振幅随位置的变化可用函数振幅随位置的变化可用函数 X(x)表示表示弦横振动方程:弦横
2、振动方程:20,ttxxux ta ux t边界条件:边界条件:0()tux 0()ttux 00(,)xu x t0(,)x lu x t初始条件:初始条件:由于两端固定,解应是驻波解,由于两端固定,解应是驻波解,求解求解两端固定两端固定的弦长为的弦长为l的的均匀细弦均匀细弦的横振动。的横振动。写出定解问题:写出定解问题:2.557.51012.515-1-0.50.51(,)()(u x tXT tx(0,0)xl t2TXXa T2()()()()TtXxX xa T t 代入波动方程代入波动方程移项整理得变量分离等式:移项整理得变量分离等式:20ttxxua u(,)()()u x t
3、X x T t把把20XTa XT得得:得出得出:200Ta TXX 0|0 xX|0 x lX()()0X l T t 5(,)()()u x tX x T t把把代入边界条件代入边界条件00(,)xu x t0(,)x lu x t本征问题本征问题:XXxX xX l 00,0,0(0)()0XT t 0 12()xxX xC eC e 00()X120CC0()X l 120llC eC e 120CC0 12()X xC xC20C 120C lC120CC000|0|0 xx lxx lXXXXXXXX 0 00()X0()X l 边界条件:边界条件:边界条件:边界条件:0 12()
4、cossinX xxCCx10C 22200000sin,sinsin,lCCllC C2非零解非零解222nnl 1 2 3,n 2()sinnnxXxCl Xn(x)称称00()X0()X l ii12()xxX xC eC e边界条件:边界条件:200,sinCl 12()xxX xC eC e n只能取分列特定值只能取分列特定值(正整数正整数)-)-称称 222nnnlT tTt n=1,2,31,2,3 解方程解方程 02 TaT nnn aTtTl 2()()0()cossinnn an aT tAtBtll 与与n对应的对应的本征解本征解un(x,t):令令C2A=An C2B=
5、BnA、B 是积分常数是积分常数本征振动的线性叠加本征振动的线性叠加.上式正好是随位置变化的傅里叶正弦级数上式正好是随位置变化的傅里叶正弦级数.1(,)(cossin)sinnnnn atn atn xu x tABlll1 2 3,n(,)()()(cossin)sinnnnnnn atn atn xux tXx T tABlll1 2 3,n 解解一般解一般解 对应于本征解对应于本征解Xn(x)的的含时函数解:含时函数解:1(,)(cossin)sinnnnn atn atn xu x tABlll0txu()初始条件:初始条件:1nnxn xAl sin()02()sinlnnAdll
6、0ttxu ()1nnn an xBxll s n(i)02()sinlnnBdn al (,)nux t是驻波,是两端固定弦的是驻波,是两端固定弦的相邻节点之间距离等于半波长相邻节点之间距离等于半波长 2ln波长波长=节点数节点数 n+1 n+1,位置位置 即即:x=kl/n(k=0,1,2,n)n xkl (4 4)、用初始条件确定一般解的系数(傅立叶展开)、用初始条件确定一般解的系数(傅立叶展开 )(1 1)、将齐次偏微分方程求解问题转化为若干常微分方程)、将齐次偏微分方程求解问题转化为若干常微分方程定解问题。定解问题。(2 2)、常微分方程与齐次边界)、常微分方程与齐次边界(或周期性或
7、周期性)条件构成本征值条件构成本征值问题得出本征值与本征函数问题得出本征值与本征函数(3 3)、将本征函数)、将本征函数(满足边界条件满足边界条件)叠加成无穷级数,给叠加成无穷级数,给出一般解出一般解10本征频率本征频率lnavlannn22,n=1=1时时,1la 基频基频基波基波(决定了音调)(决定了音调)n 1 1 时时lann 谐谐频频谐波(决定了音色)谐波(决定了音色)例例10.1:长为长为l的均匀细的均匀细杆,侧面绝热。杆的杆,侧面绝热。杆的x=0端温度保持为端温度保持为零度,另一端(零度,另一端(x=l)按牛顿冷却定律与外界进行热交换,设)按牛顿冷却定律与外界进行热交换,设外界温
8、度恒为零度。已知外界温度恒为零度。已知杆的初始温度分布为杆的初始温度分布为f(x)。试求杆上。试求杆上温度的变化(设热交换系数为温度的变化(设热交换系数为h=b/k,其中其中b为传热系数和为传热系数和k为导为导热系数热系数)。)。22,(/)txxux ta ux tak c 00(,)xu x t;第一类边界条件第一类边界条件第三类边界条件第三类边界条件0(,)(,)/xx lux tbu x tk泛定方程和泛定方程和边界条件皆边界条件皆是齐次的是齐次的,可以应用分离变量法可以应用分离变量法 00()tufxxl解:解:记杆上温度为记杆上温度为u(x,t),写出热传导方程写出热传导方程(如图
9、如图)初始条件:初始条件:边界条件:边界条件:0lx000u(,)u x t00/u x l1 1 分离变量分离变量:分离变量分离变量:(,)()()u x tX x T t20XTa XT20Ta T 0XX 12()cossinX xCxCx20000 0,txxxxxltua uuuhuufx12定解问题转化为定解问题转化为:即关于时间与位置的二个方程即关于时间与位置的二个方程:0()00()()XT t 0()()()()Xl T thX l T t2 2 X(x)方程和边界条件构成本征问题方程和边界条件构成本征问题,得解得解:00()X0()由边界条件:由边界条件:10C lh ta
10、n00()X0()()XlhX l得出:得出:220()sinX xCxC 上式应用到第二个边界条件:上式应用到第二个边界条件:引入量纲为引入量纲为1 1的量的量:1,lhltan 2221 2 3,nnnl 13前式改写为前式改写为:此方程是一个超越方程,只能用数值图解法求解(如图)此方程是一个超越方程,只能用数值图解法求解(如图)得本征值:得本征值:相应本征函数:相应本征函数:2nnXxCx()sintanlh 21 2()sin(,)nnXxCxn 1 2 3,n 220nnnTa T;22()exp;nnnT tBa t 221nnnnu x tBxa t(,)sinexp1420Ta
11、 T :改写为改写为此方程的解此方程的解本征解本征解222(,)()()sinexpnnBC Bnnnnnnux tXx T tBxa t000()()sin();nntnu xf xBxxl 2002002022dddd2dsin/sinsin/sincossin/llnnnllnnlnnBfllflflh 由初始条件由初始条件()()t ou xf x利用三角函数的正交性,等式两边同乘正弦函数,即有利用三角函数的正交性,等式两边同乘正弦函数,即有 22021d2lnnnnnxa tfu x tllh sinexpsin(,)cos220(,)sinexpnnnnu x tBxa t例例10
12、.2*:边长为边长为l1,l2的矩形薄板,两的矩形薄板,两板板面不透热,它的一边面不透热,它的一边y=l2 2为绝热,其为绝热,其余三边保持温度为零(见图)余三边保持温度为零(见图).设板的设板的的初始温度分布为的初始温度分布为f(x,y)。试求板内的。试求板内的温度变化。温度变化。21200txxyyux y taux y tux y txlyl,(,)10(,)x lu x t;第一类边界条件第一类边界条件第二类边界条件第二类边界条件泛定方程和泛定方程和边界条件皆边界条件皆是齐次的是齐次的,可以应用分离变量法可以应用分离变量法0(,),tu x y tfx y解:解:记板内任一点记板内任一
13、点(x,y)的温度为的温度为u(x,y,t),满足定解问题满足定解问题(如图如图)初始条件:初始条件:边界条件:边界条件:0l1l2绝热绝热00(,)yu x t20yy lu00(,)xu x t1 1 分离变量和本征问题分离变量和本征问题:分离变量分离变量:(,)()()()u x y tX x Y y T t 20XYTtaXYXYT t12200000 00 0,txxyyxxlyyyltuauuuuuuufx y17定解问题转化为定解问题转化为:等式两边除等式两边除:且且则则 相应边界条件:相应边界条件:2()()()a X x Y y T t得出:得出:2()()TtYY yXXa
14、 T tx 0()XXXX x 令令0()YYY yY 0()TtTTtT 12000000(),(),(),()XX lYY l 初始条件:初始条件:0()()(),X x Y y Tfx y0XX 180()2 2 X(x)和和Y(y)的常微分的常微分方程和边界条件构成本征问题方程和边界条件构成本征问题:对于对于X(x)由前面讨论:由前面讨论:00()X10()X l0YY 00()Y20()Y l3 3 本征值和本征函数本征值和本征函数:2221nnl 1 2 3,n 21()sinnnnxXxCl 0 12()cossinmmY yyCyC120()sinmmY yCyC 222222
15、00000cos,coscos,mmmlCCllC Cm2非零解非零解2222210 1 2 32(,)mmml 22212()()sinmmmyYyCl 00()Y20()Y lii12()xxmmY yC eCe边界条件:边界条件:2200,cosmCl 对于对于Y(y)类似前面讨论:类似前面讨论:1 2 30 1 2 3,nm2220;nmnmnmTaT 222222212212()()exp;mnmnnmTtBatll2222222011212212122,()()(,)sinsinexpnmmnn xmynmu x tCatllll2020TaT :改写为:改写为此方程的解此方程的解
16、本征解本征解2222222221212212122(,)()()()()()sinsinexpnmnmnmnmnmnmCCCBnmux tXx Yy Ttn xmynmCatllll0112212,()(,)sinsinnmmnn xmyf x yCll12001 212421d d2(),sinsinllnmn xmyCfx yx yl lll 由初始条件由初始条件t ou x y tf x y(,)(,)利用三角函数的正交性,等式两边同乘正弦函数,即有利用三角函数的正交性,等式两边同乘正弦函数,即有 10.620000ttxxxx ltttua uutuv tuxux,(),(),(),(
17、),000 xx lwu x tx tw x tw(,)(,)(,)目目标标x tA t xB t(,)()()0000()()xxxwutB tB tt设:设:1 1、一般处理方法、一般处理方法齐次方程齐次方程非零初值非零初值A(t)与与 B(t)是是待定参数,由边界条件确定:待定参数,由边界条件确定:解:解:设计一个设计一个新函数新函数,构成齐次边界条件的定解问题。构成齐次边界条件的定解问题。0()()()()()()x lx lx lwuv tA t lB tv ttA tl tx,令函数令函数 使使u(x,t)构成下列关系构成下列关系 (,)()v ttx txtl 即有:即有:即:即
18、:2000,0,0(),(),0(),(),0ttxxxx ltttua uxl tut uv ttux uxxl现在定解问题转换为:现在定解问题转换为:220000000ttxxxxttxx lttttttwa wawwwxwx,(),(),200000ttxxxx ltttxtvtwa wtlwwwxv ttx ltwxv ttx lt ,()/(),()/(),000 xxx lwu x tx tw x tw(,)(,)(,)目目标标(,)()()x tA t xB t 000000 xxxwuuB tB tu 仍设:仍设:则:则:解:定解问题解:定解问题000/()()/xx lxx
19、lx lwuqkA tA tqk 2000000txxxxx ltua uuu uqkuu,/,.例例10.310.3 长为长为l、侧面绝热的均匀细杆的导热问题。它的侧面绝热的均匀细杆的导热问题。它的x=0端端保持恒温保持恒温u0,另一端(另一端(x=l)有热流量为)有热流量为q0的常热流进入。设杆的常热流进入。设杆的初始温度分布也是的初始温度分布也是u0,求杆上的温度分布。求杆上的温度分布。所以所以00(,)qx txuk 将定解问题分解为将定解问题分解为(x,t)和和w(x,t)两个定解问题之和:两个定解问题之和:现在定解问题转换为:现在定解问题转换为:020000000000txxxxl
20、ttxtwa wwwuuq xwq xkuk,(/),/1 1 分离变量分离变量:分离变量分离变量:w x tX x T t(,)()()20XTa XT20Ta T 0XX 定解问题转化为定解问题转化为:即关于时间与位置的二个方程即关于时间与位置的二个方程:0()00()()XT t 0Xl T t()()00 xX0 x lX0XX 12()cossinX xCxCx270()2 2 X(x)方程和边界条件构成本征问题方程和边界条件构成本征问题,得解得解:00 xX0()由边界条件:由边界条件:10C 00 xX0 x lX220()sinX xCxC 上式应用到第二个边界条件:上式应用到
21、第二个边界条件:解简化为:解简化为:0 x lX22200000lCCllCcos,coscos,C2非零解非零解22222 310 12nnnl ,()2212nnxXxCl()()sin 200Cl,cos 0 1 2 3n ,220nnnTaT;2222212nnnT tBatl()()exp;22220212122nnnxnw x tCatll()()(,)sinexp2820TaT :改写为:改写为此方程的解此方程的解本征解本征解22222212122nnnnnCC Bnwx tXx T tnxnCatll(,)()()()()sinexp2121202201218221,esinn
22、natlnnxq lw x tkln 得:得:212120002201218221,esinnnatlnnxqq lu x txukkln 本定解问题的解本定解问题的解00212nnq x knxCl()sin/0022008221121d2221llnq lnnCllklnqk ()()sinsin 由初始条件由初始条件0t oqtkxxw(,/)利用三角函数的正交性,等式两边同乘正弦函数,即有利用三角函数的正交性,等式两边同乘正弦函数,即有0221002222022810 2 42118812128121211 3 521nq lnknq lq lnkq lknnnkn,()sin,222
23、20212122nnnxnw x tCatll()()(,)sinexpu x tx tw x t(,)(,)(,),0 xxyyux yux y0,0yu x y0(,)0,xu x y0(,)x au x yu例例10.410.430 求半带形区域(求半带形区域(00 xa,y0)内的静电势内的静电势,已知边界,已知边界x=和和y=0=0上的电势都是零上的电势都是零,而且边界而且边界x=a上的电势保势为上的电势保势为u0 0(常数)。(常数)。解:由于本题是解:由于本题是无源无源静电场问题静电场问题,所以定解问题可写为:,所以定解问题可写为:为了齐次化为了齐次化x=0=0和和x=a处的边界
24、条件,设处的边界条件,设000 xx aw x yu x yx yw x yw x y,(,)(,)(,),目目标标x yA y xB y(,)()()00000 xxxwuB yB yo 仍设:仍设:则:则:000 x ax ax awuuA y aA yua()()/所以所以0ux yxa(,)现在定解问题转换为:现在定解问题转换为:000000 xxyyxx aywwu x awww,/,利用分离变量法和上式第一行与第二行方程得出利用分离变量法和上式第一行与第二行方程得出:1n yn yaannnn xw x yA eB ea(,)sin yw是有限值这里这里An n和和Bn n是待定常
25、数。但是待定常数。但不能确定不能确定An和和Bn。注意到注意到u|y y时应有限,所以相应地有时应有限,所以相应地有:0201 2 311 2 3nnnuAnBnn,0001221n ynanuuun xu x yxeana,sin 结合:结合:u|y y有限有限和和00ywu x a/得出得出定解问题解为:定解问题解为:u x yx yw x y(,)(,)(,)20ttxxua u00tu00ttu0(,)0 xu x t(,)sinx lu x tAt(,)()sinx tA xt 由于弦在由于弦在x=l端端受迫作受迫作Asint的的振动振动,它一定存在它一定存在(x,t),满足齐次方程
展开阅读全文