数理方程-分离变量法(2).ppt

  1. 1、本文档共51页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
解: 其中 * * 例4 求下列热传导方程的定解问题 解法一:令 * * 解法二:令 特征值和特征函数分别为 代入边界条件 * * 则分离变量形式的特解 解之得: 其中: 定解问题的解: 【解】杆上温度满足下列泛定方程和定解条件 试探解 代入方程和边界条件得 固有值问题 【例题3】研究细杆导热问题,初始时刻杆的一端温度为零度, 另一端跟外界绝热,杆上初始温度为 ,试求无热源时细杆上 温度的变化。 和常微分方程 分析:方程与边界条件均为齐次,用分离变量法。 特征值为: 特征函数: 由初值条件知 下面求解 由叠加原理,得 得 从而下列问题 的解为 图形如下: (程序:my1) (a) 精确解图 (b) 瀑布图 图形如下: (程序:my1) 1 矩形区域上拉普拉斯方程 【例题3】散热片的横截面为矩形。它的一边 处于较高温度 , 边处于冷却介质中而保持较低的温度 , 其他两边 , 温度保持为零, 求解这横截面上的稳定温度分布  . 【解】先写出定解问题定解问题 方程齐次 这组边界条件齐次 用分离变量法 设形式解为: 代入上述泛定方程,得到 得到固有值问题 和常微分方程 得特征值: 特征函数: 而 于是有 叠加得 为确定叠加系数,将 代入非齐次边界条件 将等式右边展开为傅里叶正弦级数,并两边比较系数,得 联立求解得 故原问题的解为 小结:对矩形域上拉普拉斯方程,只要一组边界条件 是齐次的,则可使用分离变量法求解。 图形如下: (程序:my2) (a) 精确解图 (b) 瀑布图 小 结 分离变量法求解偏微分方程的基本步骤: ? 分离变量(齐次条件) ? 求解本征值 ? 求出所有特解,叠加出一般解 ? 利用本征函数正交性定出叠加系数 其中: 解为: 故所求的解为 练习:求定解问题的解 解为 第二章 分 离 变 量 法 齐次稳定场问题的求解 本章小结 作业: 备用知识回顾 齐次边界问题的求解 二阶常系数非齐次线性方程解法 解法 待定系数法。 方程 §2.1 分 离 变 量 研究两端固定弦的自由振动. 定解问题为: 特点:方程和边界条件都是线性齐次的. 引例: 有界弦的自由振动 设 且 不恒为零,代入方程得 ① 由 不恒为零,有: 这个式子的左端是x的函数, 右端是t的函数,何时恒等? ② …..…….. ③ 思考:先解哪一个方程? 则 ⑤ 特征值问题 故 而 正是Fourier正弦级数的基本函数族. 再求解T: 其解为 所以 两端固定弦的特征振动 特解,未必满足初始条件(2.3) 受叠加原理启发, 第三步:先求特解,再叠加求出通解 代入初始条件得: ………………...⑥ 将 展开为Fourier正弦级数,比较系数得 第四步: 利用特征函数的正交归一性确定待定系数 由叠加原理,如果上式右端的无穷级数是收敛的,并且关于 x,t 能逐项微分两次,则和式 u(x,t) 确实是问题(2.1)-(2.3)的解(古典解)。 其中 和 由上面给出。 如果(*)定义的函数 u(x,t)不具备古典解的要求,则 称为问题(2.1)-(2.3)的形式解。 则无穷级数解 上, ,且 条件:若在区间 为混合问题(2.1)-(2.3)的古典解, 其中 注1: 本书不讨论所求形式解是否满足古典解要求的条件,只要求得了形式解,就认为定解问题得到了解决。 注2: 用分离变量法求解定解问题的关键是确定特征函数和运用叠加原理,这些运算能够进行,是条件限制。 (1)第一个限制:变系数的二阶线性偏微分 方程并非总能实施变量分离 (2)第二个限制:二阶线性偏微分方程的解, 不一定是分离变量的乘积形式 ⑵弦上各点振幅 因点而异 在 处,振幅永远为0 二、解的物理意义 节点 腹点 角频率 初位相 在 处,振幅最大,为 n N ⑴弦上各点的角频率 和初位相 都相同,因而没 有波形的

文档评论(0)

dajuhyy + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档