微分方程数值解课程设计.doc

微分方程数值解课程设计 题一：重调和方程的二阶边值问题的求解 Possion方程第一边值问题是数学、物理和力学等领域中一类最常见的也是最重要的偏微分方程的求解问题，对于这些问题已有很多求解的方法，如五点差分格式。但在实际问题中，常会遇到如下重调和方程的二阶边值问题 （1） 其中是Laplace算子，是二维平面上的有限区域，是其光滑边界，a，b非负常数。它可以用来描述平板在不同荷载下的横向变形，以及stokes流体的流动等。问题（1）如果直接采用差商逼近导数离散可以得到经典的二阶精度的13点差分格式，但由于差分方程涉及的离散结点数较多，从而对邻近边界点处的计算带来困难。另一种方法是，采用适当的技巧将其分解成两个Possion方程求解。具体如下 假设若（若 分解方式不同，此处不讨论）， 令，那么边值问题（1）可转化为下列两个边值问题 （2） （3） Dirichlet边值问题（2）及（3）可以采用经典的五点差分格式逼近，从而求出问题（1）的数值解。 数值算列： 1. ，，使问题（1）存在精确解（边界条件由此解给出）。 2. ，，问题（1）存在精确解（边界条件由此解给出）。 计算不同步长选取下的数值解，且分别取不同的a，b值且使，如 ① ② ③ ④. 内容包括： (1) 格式推导、误差分析； 数值解与精确解的比较，不同步长选取下的误差比较(讨论是否收敛及收敛阶)。 数值计算收敛阶: 其中Error1,Error2分别是取步长h1和h2下计算所得误差值. 注: h1和h2同一个方向的两个步长, 另外一个方向的步长固定且尽可能小. 题二：热传导现象数值模拟 设一根长度为的均匀热导细杆，其两端固定温度源或绝热，并给其侧面外界热源，讨论细杆内温度变化规律。 温度函数满足如下热传导抛物型方程 （4） 不妨假设。 设初始时刻细杆上的温度为，并且细杆两端温度固定为0，即热传导方程（4）满足如下定解初边值条件： 此时方程存在精确解 设初始时刻细杆上的温度为，并且细杆左端温度固定为0，而右端绝热，即热传导方程（4）满足如下定解初边值条件： 设初始时刻细杆上的温度为，并且细杆左端绝热，而右端温度固定为0，即热传导方程（4）满足如下定解初边值条件： 用最简显格式求解上面的定解问题，并试讨论如下问题等： 各时刻的大概最高温度点； 各点温度随时间的变化规律； 各时刻温度随位置的变化规律； 2和3两种情况下有何关联等。 另外，内容还包括：格式推导（注意边界条件的处理！）、截断误差分析；并试讨论1的稳定性（选取不同的步长计算，观察计算结果是否稳定）； 如 ① ② ③ ④ ⑤ 课程设计须知： 严禁抄袭！若发现雷同，不区分抄与被抄，一并处理，成绩不及格或要求重做。 一、格式要求具体见附件1(课程设计报告模板) （1）题目 (恰当) （2）摘要：（150字左右） 摘要是对论文（设计报告）内容不加注释和评论的简短陈述，要求扼要说明研究工作的目的、主要方法、研究结果、结论、科学意义或应用价值等，是一篇具有独立性和完整性的短文。摘要中不宜使用公式、图表以及非公知公用的符号和术语， 不标注引用文献编号。 （3）关键词：四个（按重要性、频率排序，分号隔开） （4）引言：（背景、问题的提出等） （5）正文主体 正文基本内容包括（可进一步发挥）： 方法介绍、理论分析结论； 设计思路、算法步骤（或流程图）； 程序清单（函数文件、命令文件）（加上必要的注释）； 程序运行操作过程与输出结果（必须附上相应的截屏图，图中须有任务栏和命令历史窗口中的日期、时间， 为了行文的方便, 可以放在附件中）； 对计算过程与结果的分析（如误差分析，收敛性，稳定性，计算量，存储量，方法比较，等等）；各设计的优缺点（如特色、自己最满意之处、需改进的地方等）； 重点考虑方法的收敛精度，及稳定性 （6）心得体会、自我评价（优、良、中、及格、差之一）及其支持依据 （7）参考文献 (8）附录（可有可无） 其他要求 1. 课程设计期间，可以相互讨论，但是具体算法实现、数据分析和报告的撰写要求独立完成，严禁抄袭。 判断抄袭的参考标准：出现下列情形之一 无截屏图； 截屏图与他人相同； 需编写的M文件（含注释）与他人完全相同； 无“课程设计中的心得体会”或过于简短； “课程设计中的心得体会”不含“遇到的问题与困难及解决办法”； “课程设计中的心得体会” 与他人相同。 2. 课程设计提交时间： 按要求完成后，统一交给学习委员，于20周之前交至任课教师处，包括光盘和纸质的课程设计报告。过期不交作为缺考处理。 请将全班