当前位置：文档之家 › 微分方程课程设计

微分方程课程设计

  • 格式：doc
  • 大小：205.65 KB
  • 文档页数：15

湖南工程学院课程设计任务书课程名称微分方程数值解法课题常微分方程初值问题专业班级信息与计算科学班学生姓名学号指导老师审批任务书下达日期2011 年12 月20 日任务完成日期2011 年12 月30 日一、设计内容与设计要求1．设计内容：对课程《微分方程数值解法》中的常见算法进行综合设计或应用。

具体地，对型方程的差分解法进行研究，包括选取合理的差分网格、建立差分格式、求解代数方程组、考察差分格式的收敛性等问题。

2．设计要求：●课程设计报告正文内容a.问题的描述及算法设计；b.算法的流程图；c.算法的理论依据及其推导；d.相关的数值结果；e.数值计算结果的分析；f.附件。

●书写格式a．要求用A4纸打印成册b．正文格式：一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。

c．正文的内容:正文总字数要求在3000字左右（不含程序原代码）。

d．封面格式如下页。

●考核方式指导老师负责验收程序的运行结果，并结合学生的工作态度、实际动手能力、创新精神和设计报告等进行综合考评，并按优秀、良好、中等、及格和不及格五个等级给出每位同学的课程设计成绩。

具体考核标准包含以下几个部分：a．平时出勤（占10%）b．系统需求分析、功能设计、数据结构设计及程序总体结构合理与否（占10%）c．程序能否完整、准确地运行，个人能否独立、熟练地调试程序（占40%）d．设计报告（占30%）注意：不得抄袭他人的报告（或给他人抄袭），一旦发现，成绩为零分。

e．独立完成情况（占10%）。

课程验收要求a．判定算法设计的合理性，运行相关程序，获得正确的数值结果。

b．回答有关问题。

c．提交课程设计报告。

d．提交软盘（源程序、设计报告文档）。

e．依内容的创新程度，完善程序情况及对程序讲解情况打分。

二、进度安排1、班级：信息与计算科学专业2008级本科班2、主讲教师：杨继明3、辅导教师：杨继明4、时间安排：2011年12月20日—2011年12月30日每天8：30—11：30， 3：00—5：30目录一、问题的描述 (5)二、算法设计及流程图 (5)1.1算法设计 (5)1.2流程图 (6)三、算法的理论依据及其推导 (7)3.1 欧拉法理论推导 (4)3.2 龙格库塔法理论推导 (5)四、数值结果分析及分析 (9)五、总结 (12)六、附件 (13)常微分方程的初值问题一、问题的描述很多科学技术和工程问题常用微分方程的形式建立数学模型，因此微分方程的求解是很有意义的。

建立微分方程只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并加以检验。

如果能得到解析形式的解固然是便于分析和应用的，但是我们知道，虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些典型的方程，而对于绝大多数的微分方程问题，很难或者根本不可能得到它的解析解，实际问题终归结出来的微分方程主要靠数值解法。

因此，研究微分方程求解的数值方法是非常有意义的。

本文介绍了欧拉法、和四阶龙格_库塔方法二种单步法，通过matlab 的平台运行。

二、算法设计及流程图2.1算法设计常微分方程初值问题的数学模型是：求y ()y x =，使之满足0'()(,)()y x f x y x b y a y =≤≤⎧⎨=⎩，(a ),, (1)其中0y ，a ，b 是已知的常数，(,)f x y 是已知函数，且满足条件：(,)(,')'f x y f x y L y y -≤-，式中的L 是不依赖于y ，'y 的常数，称为利普希茨（Lipschitz ）常数。

由常微分方程理论知识可知，上述问题存在唯一解()y x 。

现在的目标就是计算区间[],a b 上等节点i x a i h =+处该未知函数的函数值()i y x ，其中()/h b a N =-，0,1,2,...,i N =。

为此，可以用求解常微分方程问题的单步法，即欧拉方法和龙格_库塔方法求解。

2.2 流程图2.2.1欧拉法流程图：程序代码见附录。

x0+h=>x1y0+h*f(x0,y0)=>y1n=1输出x1,y1n=1+nx1=> x0y1=> y0 结束n=N ?读入x0,y0,b,h开始计算N=fix((b-x0)/h)2.2.2 四阶龙格_库塔法流程图：程序代码见附录。

三、算法的理论依据及其推导3.1欧拉算法的理论推导将微分方程离散化，用向前商1()()n n y x y x h+-代替微分()n y x ，代入（1）中的微分方程，可得1()()(,())1,2,3,...n n n n y x y x f x y x n h+-==，（）化简可得计算N=fix((b-x0)/h)n=1x0+h=>x1f(x0,y0)=>K1f(x0+h/2,y0+h/2*K1)=>K2 f(x0+h/2,y0+h/2*K2)=>K3 f(x0+h,y0+h*K3)=>K4y0+h/6*(K1+2*K2+2*K3+K4)=>y1输出x1,y1n=n+1x1=>x0 y1=>y0n=N ?结束开始读入x0,y0,b,h1()()(,())1,2,3,...n n n n y x y x f x y x h n +=+=，（）如果用n y 近似()n y x 代入上式便可得到1()n y x +的近似值1n y +，计算式为：1(,)1,2,3,...n n n n y y f x y h n +=+=，（）(2) 这样问题(1)的近似解课通过求解下面的差分初值问题：10(,)1,2,3,...()n n n n y y f x y h n y a y +=+=⎧⎨=⎩，（）(3) 得到，按(3)式由初值0y 可逐次求出12y ,,...y 。

欧拉方法就是用差分方程初值问题(3)的解来近似微分方程初值问题(1)的解。

即由公式(3)算出()n y x 的近似值。

这组公式求问题(1)的数值解就是著名的欧拉(Euler)公式。

欧拉法是数值求解微分方程最简单的方法，它特别适合快速编程，尽管精度不是很高。

Euler 方法的局部截断误差为11()()(,())n n n n n T y x y x hf x y x ++=--()()'()n n n y x h y x h yx =+-- 232''()/2()()nh y x h h οο=+≈即局部截断误差是2h 阶的。

而数值算法的精度定义为：若一种算法的局部截断误差为1()p h ο+，则称该算法具有p阶精度。

显然p越大，方法的精度越高。

上式说明，欧拉方法是一阶方法，因此它的精度不高。

Euler 方法的整体截断误差为()n n n e y x y =-((1)1)/((1)1)nn e hL hL T ≤+-+-()(1)/()()L b a ehL T h ο-≤-=即整体截断误差是h 阶的。

3.2 龙格_库塔方法理论推导由Lagrange 微分中值定理，11()()'()()()(,())n n n n n y x y x y x x y x hf y ξξξ++=+-=+记*(,())k hf y ξξ=，则得到*1()()n n y x y x k +=+这样，给出*k 的一种算法，就得到求解微分方程初值问题的一种计算公式。

四阶龙格_库塔法是用1k ，2k ，3k 和4k 的加权平均值来近似*k 。

最经典的四阶龙格_库塔公式为：121324311234(,)(,)22(,)22(,)y (22)2n n n n n nn n n n K f x y h hK f x y K h h K f x y K K f x h y hK hy K K K K +=⎧⎪⎪=++⎪⎪⎪=++⎨⎪=++⎪⎪⎪=++++⎪⎩四阶龙格_库塔法的误差估计局部截断误差为5()O h 。

值得指出的是，龙格_库塔方法的推导基于泰勒展开方法，因而它要求的解具有较好的光滑性质。

反之，如果解得光滑性差，那么，使用四阶龙格_库塔法方法求得的数值解，其精度可能反而不如梯形方法。

实际计算中，应当针对问题的具体特点选择合适的算法。

对欧拉法、梯形法和龙格_库塔法进行比较：在相同步长的情况下，欧拉法每步只计算一个函数值，四阶龙格_库塔法每步需计算四个函数值，就是说，四阶龙格_库塔法的计算量差不多是欧拉法的四倍，为了比较它们的计算精度，可以将欧拉法的步长取为h ,将四阶龙格_库塔法的步长取为4h 。

这样，如果用二种方法求解同一初值问题，则它们的计算量相当。

在计算量相当的条件下，比较它们的计算结果，就能够看出它们的精度差异。

四、数值结果及分析考虑下列常微分方程的初值问题：4.1欧拉方法.1)0(,5'=-=y y y采用欧拉法进行求解 分别取步长为h=0.1 数值格式为：采用Matlab 程序设计语言编程实现该问题的数值求解。

在Matlab 运行得出的结果是： R =0 1.0000 0.1000 0.8979 0.2000 0.8273 0.3000 0.7726 0.4000 0.7406 0.5000 0.7129 0.6000 0.6956 0.7000 0.6932 0.8000 0.6896 0.9000 0.7038 1.0000 0.73424.2四阶龙格库塔方法⎩⎨⎧=≤≤-=1)0(10'y x y x y 采用四阶龙格库塔方法进行求解。

数值格式为：10=y),(1i i y t f k =)2,2(12k h y h t f k i i ++= )2,2(23k h y h t f k i i ++=),(34hk y h t f k i i ++=)22(643211k k k k h y y i i ++++=+其中h 为区间步长。

采用Matlab 程序设计语言编程实现该问题的数值求解。

在Matlab 命令窗口中输入命令： R= rk4(@f,0,1,1,10) 运行程序得到实验结果为： R =0 1.0000 0.1000 0.9097 0.2000 0.8375 0.3000 0.7816 0.4000 0.7406 0.5000 0.7131 0.6000 0.6976 0.7000 0.6932 0.8000 0.6987 0.9000 0.7131 1.0000 0.7358五、总结通过两周的课程设计，进一步了解了常微分方程的初值问题的解法，通过用欧拉法、四阶龙格_库塔法解决，只要输入相关给出的初始值和步长就可以相应的近似解的方法。

完成这次课程设计，我选用的运行软件是matlab，这次课程设计是对matlab功能的进一步熟悉和了解，尤其在函数的运用方面提高了很多。

合集下载

常微分方程课程设计论文

常微分方程课程设计论文

常微分方程课程设计论文一、教学目标本课程的教学目标是使学生掌握常微分方程的基本概念、方法和应用。

通过本课程的学习，学生应能理解并熟练运用常微分方程解决实际问题，具备一定的数学建模能力。

具体来说，知识目标包括：1.掌握常微分方程的定义、解的概念和性质；2.熟悉一阶、二阶线性微分方程的求解方法；3.了解常微分方程在自然科学和工程技术中的应用。

技能目标包括：1.能够熟练地求解一阶、二阶线性微分方程；2.能够运用常微分方程进行简单的数学建模；3.能够运用计算机软件辅助求解常微分方程。

情感态度价值观目标包括：1.培养学生的逻辑思维能力和科学精神；2.增强学生对数学应用价值的认识，提高学习兴趣；3.培养学生团队协作和自主学习能力。

二、教学内容根据教学目标，本课程的教学内容主要包括：1.常微分方程的基本概念，如解、通解、特解等；2.一阶微分方程的求解方法，如可分离变量法、齐次方程法、伯努利方程法等；3.二阶线性微分方程的求解方法，如常系数方程、变系数方程、线性非齐次方程等；4.常微分方程的应用，如物理、生物学、经济学等领域的问题。

三、教学方法为了达到教学目标，本课程将采用以下教学方法：1.讲授法：系统地传授常微分方程的基本概念、方法和应用；2.讨论法：学生分组讨论，培养学生的思考能力和团队协作精神；3.案例分析法：通过分析实际问题，引导学生运用常微分方程进行数学建模；4.实验法：利用计算机软件，让学生亲自动手求解实际问题，提高实际操作能力。

四、教学资源为了支持教学内容和教学方法的实施，本课程将准备以下教学资源：1.教材：《常微分方程》；2.参考书：相关领域的学术论文、专著等；3.多媒体资料：教学PPT、视频讲座等；4.实验设备：计算机、数学软件等。

五、教学评估本课程的评估方式包括平时表现、作业、考试等，以全面客观地评价学生的学习成果。

平时表现主要考察学生的课堂参与、提问、讨论等，占总评的20%；作业包括练习题和数学建模项目，占总评的30%；考试包括期中考试和期末考试，占总评的50%。

微分方程基础及应用课程设计

微分方程基础及应用课程设计

微分方程基础及应用课程设计1. 课程概述微分方程是数学中的一门重要的学科，它在自然科学和工程技术等领域中都有着广泛的应用。

本课程介绍微分方程的基本概念、基本解法和常见应用。

主要内容包括：一阶微分方程、二阶及高阶微分方程、常系数线性微分方程、解析解和数值解等内容。

2. 教学目标2.1 知识目标1.掌握一阶微分方程的解法和应用；2.掌握二阶及高阶微分方程的解法和应用；3.掌握常系数线性微分方程的解法和应用；4.了解解析解和数值解的基本研究方法。

2.2 能力目标1.能够运用微分方程解决实际问题；2.能够对数学模型进行合理的建立和求解；3.能够开展科研和工程设计中的微分方程相关的基础性研究。

3. 教学大纲3.1 课程内容1.一阶微分方程1.可分离变量的微分方程2.齐次线性微分方程3.Bernoulli方程4.恰当微分方程5.变量可分离微分方程6.常系数线性微分方程7.应用2.二阶及高阶微分方程1.齐次线性微分方程2.非齐次线性微分方程3.应用3.常系数线性微分方程1.一阶常系数微分方程2.二阶及高阶常系数微分方程3.应用4.解析解和数值解1.Taylor公式2.Euler法3.Runge-Kutta法4.应用3.2 教学进度和安排章节内容学时1.1 可分离变量的微分方程 21.2 齐次线性微分方程 21.3 Bernoulli方程 21.4 恰当微分方程 21.5 变量可分离微分方程 21.6 常系数线性微分方程 2章节内容学时1.7 应用 22.1 齐次线性微分方程 22.2 非齐次线性微分方程 22.3 应用 23.1 一阶常系数微分方程 23.2 二阶及高阶常系数微分方程 23.3 应用 24.1 Taylor公式 24.2 Euler法 24.3 Runge-Kutta法 24.4 应用 24. 教学方法本课程采用理论讲授、课堂讨论、案例分析和课程实践相结合的教学方法，注重培养学生的综合运用能力和分析问题的能力。

高中数学教案微分方程

高中数学教案微分方程

高中数学教案微分方程微分方程教案高中数学教案摘要：本教案主要介绍微分方程的基本概念、求解方法和应用，并设计了相关的教学活动和练习。

教学目标：1.了解微分方程的概念与分类，并理解微分方程的意义。

2.能够运用常微分方程的解法，解决简单的微分方程问题。

3.了解微分方程在实际问题中的应用，并能够将数学知识与实际问题相结合，解决实际问题。

教学重点：1.微分方程的概念与分类。

2.常微分方程的解法。

3.微分方程在实际问题中的应用。

教学难点：1.应用题中的问题分析和建立微分方程的能力。

2.求解复杂微分方程的能力。

教学准备：1.教师：PPT课件、教案、多媒体设备。

2.学生：教材、笔记本、计算器。

教学过程：一、导入（约5分钟）教师通过给学生出示一些实际问题，引发学生对微分方程的思考，激发学生的学习兴趣。

例如：一辆汽车在某段路程上的速度是多少？一杯冷水从什么温度下降到什么温度需要多长时间？二、知识讲解（约25分钟）1.微分方程的概念与分类（10分钟）教师结合多媒体展示，详细介绍微分方程的定义和分类，包括常微分方程和偏微分方程的区别，以及一阶、二阶微分方程等。

2.常微分方程的解法（15分钟）教师重点讲解常微分方程的解法，包括分离变量法、齐次方程法、一阶线性方程法等，并通过示例演示每种方法的具体步骤和应用场景。

三、教学活动（约15分钟）1.小组讨论（10分钟）将学生分成小组，让他们根据所学的知识，自行应用解题方法解决教师提供的实际问题。

鼓励学生自主思考、合作探讨，培养学生的问题解决能力和团队精神。

2.展示与总结（5分钟）请每个小组派代表展示解题过程和结果，并让其他小组评价和提问。

教师及时纠正错误，总结解题思路和方法。

四、知识拓展（约20分钟）教师通过讲解微分方程在实际问题中的应用，如放射性衰变问题、人口增长问题等，引导学生将数学知识与实际问题相结合，培养学生应用数学解决实际问题的思维能力。

五、教学总结（约5分钟）教师对本节课的重点内容进行总结，并展示一些典型的习题，让学生巩固所学知识。

常微分方程课程设计

常微分方程课程设计

常微分方程课程设计1. 背景介绍常微分方程是数学中的一门重要的分支，其应用广泛，尤其在工程、物理、经济等学科中有着重要的地位。

常微分方程的研究对象是描述性能或物理现象的数学关系，其解法方法多种多样，如常系数线性微分方程、变系数线性微分方程、常系数非线性微分方程等。

针对这些不同类型的微分方程，研究者们提出了丰富的解法，如变量分离法、常数变易法、欧拉公式等。

2. 目的和意义常微分方程的解法非常实用，是解决多种实际问题的有效工具。

因此，本次课程设计旨在探究不同类型的常微分方程，并结合实例进行分析和解决问题，加深学生对常微分方程的理解和应用。

通过本次课程设计，学生可以很好地掌握常微分方程的基本概念及其解法方法，增强其在实际问题中应用数学知识解决实际问题的能力，提高其数学思维能力，帮助学生打牢数学基础。

3. 设计内容3.1 基本概念常微分方程的基本概念是本课程设计的基础。

在教学过程中，需要向学生简要介绍基本概念，包括常微分方程的定义、解、常微分方程的分类等。

3.2 常系数非齐次线性微分方程常系数非齐次线性微分方程是常微分方程的重要类型之一，解决实际问题时常常需要采用非齐次线性微分方程的解法。

因此，在本次课程设计中，需要详细介绍常系数非齐次线性微分方程的定义、特点、解法等，并在此基础上，结合实例进行分析和应用。

3.3 高阶常微分方程高阶常微分方程是常微分方程的重要类型之一，也是常微分方程的拓展，其解决了微分方程中仅含一自变量的情况，具有广泛的应用场景。

因此，在本次课程设计中，需要向学生详细介绍高阶常微分方程的定义、特点、解法等，并结合实例进行分析和应用。

3.4 常系数线性微分方程组常系数线性微分方程组是常微分方程的一种特殊情况，也是一种常见的微分方程类型，其解法多样，具有一定的难度。

因此，在本次课程设计中，需要向学生详细介绍常系数线性微分方程组的定义、特点、解法等，并结合实例进行分析和应用。

4. 总结本次课程设计中，我们主要介绍了常微分方程的各种类型和解法方法，可以帮助学生更好地理解和掌握常微分方程的基本概念和解法方法。