微分方程建模(同名437)课件

微分方程建模1微分方程建模1

当我们描述实际对象的某些特性随时间（空间）而演变的过程、分析它的变化规律、预测它的未来形态、研究它的控制手段时。通常要建立对象的动态模型。

在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题。2感谢你的观看2019年8月21当我们描述实际对象的某些特性随时间（空间）而一、数学建模的基本思维过程

转化实际问题

1、对要讨论的问题所涉及的重要特征进行合理的数学刻画（转化），即用数学语言对问题涉及到的重要特征进行表述.2、寻求的实际问题的文字叙述，利用一些原则或定律，将其转化为数学描述。解数学问题用数学工具求解得到的数学问题。3感谢你的观看2019年8月21一、数学建模的基本思维过程转化实际问题2、寻求的二、微分方程模型

涉及“改变”、“变化”、“增加”、“减少”、“衰变”、“边际”、“速度”、“运动”、“追赶”、“逃跑”、、、等等词语的确定性连续问题。b、微分方程建模的基本手段微元法等a、微分方程建模的对象4感谢你的观看2019年8月21二、微分方程模型涉及“改变”、“变化

1、寻找改变量一般说来微分方程问题都遵循这样的文字等式变化率（微商）=单位增加量--单位减少量等式通常是利用已有的原则或定律。c、微分方程建模的基本规则2、对问题中的特征进行数学刻画3、用微元法建立微分方程；4、确定微分方程的定解条件（初边值条件）；5、求解或讨论方程（数值解或定性理论）6、模型和结果的讨论与分析。5感谢你的观看2019年8月211、寻找改变量一般说来微分方程问题都遵循这样的文字对论文的评价主要以“假设的合理性、建模的创造性、结果的正确性和文字表述的清晰性”为标准。所以，在论文中应努力反映出这些特点。三、

数学建模论文的撰写方法6感谢你的观看2019年8月21对论文的评价主要以“假设的合理性、建模的创造性、结果的正确性论文格式及组成题目摘要，关键词问题重述模型假设分析与建立数学模型数学模型的求解模型检验(总结与评价)模型推广参考文献附录（若有）正文7感谢你的观看2019年8月21论文格式及组成题目模型假设正文7感谢你的观看2019年8月21、题目论文题目是一篇论文给出的涉及论文范围及水平的第一个重要信息。要求简短精练、高度概括、准确得体、恰如其分。既要准确表达论文内容，恰当反映所研究的范围和深度，又要尽可能概括、精练。8感谢你的观看2019年8月211、题目论文题目是一篇论文给出的涉及论文范围及水平的第一个重2、摘要、关键词摘要是论文内容不加注释和评论的简短陈述，其作用是使读者不阅读论文全文既能获得必要的信息。在数学建模论文中，摘要是非常重要的一部分。数学建模论文的摘要应包含以下内容：所研究的实际问题、建立的模型、求解模型的方法、获得的基本结果以及对模型的检验或推广。论文摘要需要概括、简练的语言反映这些内容，尤其要突出论文的优点，如巧妙的建模方法、快速有效的算法、合理的推广等。一般科技论文的摘要要求不列举例证，不出现图、表和数学公式，不自我评价，且字数200以内。前几年，全国大学生数学建模竞赛要求摘要字数应在300字以内。但从2001年开始，为了提高论文评选效率，要求将论文第一页全用作摘要，对字数已无明确限制。在摘要中也可适当出现反映结果数学公式。9感谢你的观看2019年8月212、摘要、关键词摘要是论文内容不加注释和评论的简短陈述，其作3、问题重述数学建模比赛要求解决给定的问题，所以论文中应叙述给定问题。撰写这部分内容时，不要照抄原题，应把握问题的实质，再用较精练的语言叙述问题。10感谢你的观看2019年8月213、问题重述10感谢你的观看2019年8月214、模型假设建模时，要根据问题的特征和建模目的，抓住问题的本质，忽略次要因素，对问题进行必要的简化，做出一些合理的假设。模型假设部分要求用精练、准确的语言列出问题中所给出的假设，以及为了解决问题所做的必要、合理的假设。假设作得不合理或太简单，会导致错误的或无用的模型；假设作得过分详尽，试图把复杂对象的众多因素都考虑进去，会使工作很难或无法继续下去，因此常常需要在合理与简化之间作出恰当的折中。11感谢你的观看2019年8月214、模型假设建模时，要根据问题的特征和建模目的，5、分析与建立模型根据假设，用数学的语言、符号描述对象的内在规律，得到一个数学结构。建模时应尽量采用简单的数学工具，使建立的模型易于被人理解。在撰写这一部分时，对所用的变量、符号、计量单位应作解释，特定的变量和参数应在整篇文章保持一致。为使模型易懂，可借助于适当的图形、表格来描述问题或数据。12感谢你的观看2019年8月215、分析与建立模型12感谢你的观看2019年8月216、模型求解使用各种数学方法或软件包求解数学模型。此部分应包括求解过程的公式推导、算法步骤及计算结果。为求解而编写的计算机程序应放在附录部分。有时需要对求解结果进行数学上的分析，如结果的误差分析、模型对数据的稳定性或灵敏度分析等。13感谢你的观看2019年8月216、模型求解13感谢你的观看2019年8月217、模型检验把求解和分析结果翻译回到实际问题，与实际的现象、数据比较，检验模型的合理性和适用性。如果结果与实际不符，问题常出在模型假设上，应该修改、补充假设，重新建模。这一步对于模型是否真的有用十分关键。14感谢你的观看2019年8月217、模型检验14感谢你的观看2019年8月218、模型推广模型的评价不要流于空泛，需对模型的意义、可信程度、精度等可能问题，需要认真地思考和讨论。将该问题的模型推广到解决更多的类似问题，或讨论给出该模型的更一般情况下的解法，或指出可能的深化、推广及进一步研究的建议。15感谢你的观看2019年8月218、模型推广15感谢你的观看2019年8月219、参考文献在正文中提及或直接引用的材料或原始数据，应注明出处，并将相应的出版物列举在考文献中。需标明出版物名称、页码、著者姓名、出版日期、出版单位等。参考文献按正文中的引用次序列出，其中书籍的表述方式为：[编号]作者，书名[M]，出版地：出版社，出版年。参考文献中期刊杂志论文的表述方式为：[编号]作者，论文名[J]，杂志名，出版年，卷期号：起止页码。参考文献中网上资源的表述方式为：[编号]作者，资源标题，网址，访问时间（年月日）。16感谢你的观看2019年8月219、参考文献在正文中提及或直接引用的材料或原始数据，应注明出10、附录附录是正文的补充，与正文有关而又不便于编入正文的内容都收集在这里。包括：计算机程序、比较重要但数据量较大的中间结果等。为便于阅读，应在源程序中加入足够的注释和说明语句。17感谢你的观看2019年8月2110、附录附录是正文的补充，与正文有关而又不便于编入正文的内18最简单的数学模型之一——“航行问题”用x表示船速，y表示水速，列出方程：答：船速每小时20千米/小时.甲乙两地相距750千米，船从甲到乙顺水航行需30小时，从乙到甲逆水航行需50小时，问船的速度是多少?x=20y=5求解感谢你的观看2019年8月2118最简单的数学模型之一——“航行问题”用x表示船速，y19航行问题建立数学模型的基本步骤

作出简化假设（船速、水速为常数）；

用符号表示有关量（x,y表示船速和水速）；

用物理定律（匀速运动的距离等于速度乘以时间）列出数学式子（二元一次方程）；

求解得到数学解答（x=20,y=5）；

回答原问题（船速每小时20千米/小时）。感谢你的观看2019年8月2119航行问题建立数学模型的基本步骤作出简化假设（船速、水速四、导弹跟踪问题1、实验目的本实验主要涉及常微分方程。通过实验复习微分方程的建模和求解；介绍两种微分方程的数值方法：Euler法和改进的Euler法；还介绍了仿真方法。203、数学模型微分方程建模的方法主要是依据守恒律来建立等量关系式。对于这个问题，寻求等量关系是比较简单的。设坐标系如图3.1所示，取导弹基地为原点O(0,0)，x轴指向正东方,y轴指向正北方。2、实际问题某军的一导弹基地发现正北方向120km处海面上有敌舰一艘以90km/h的速度向正东方向行驶。该基地立即发射导弹跟踪追击敌艇，导弹速度为450km/h，自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌舰？感谢你的观看2019年8月21四、导弹跟踪问题1、实验目的203、数学模型2、实际问题感

当t=0时，导弹位于点O，敌艇位于点A(0,H)，其中H=120（km）。21设导弹在t时刻的位置为P(x(t),y(t))，由题意,感谢你的观看2019年8月21当t=0时，导弹位于点O，敌艇位于点A(0,H方程(3.1),(3.3)连同初值条件构成了一个关于时间变量t的一阶微分方程组的初值问题。其中另外在t时刻,敌艇位置为,其中。由于导弹轨迹的切线方向必须指向敌艇，即直线PM的方向就是导弹轨迹上点P的切线方向,故有或写为22两边对t求导为了寻求x与y的关系,要设法消去变量t,由式(3.2)得感谢你的观看2019年8月21方程(3.1),(3.3)连同初值条件其中22两边对t求导即有把式(3.1)写为代入上式,就得到轨迹方程.这是一个二阶非线性微分方程,加上初值条件,则初值问题23感谢你的观看2019年8月21即有23感谢你的观看2019年8月21就是导弹轨迹的数学模型。值得注意的是，前面导出的一阶微分方程组(3.1)，(3.3)和(3.4)实际上已经是一个数学模型了，不过多一个变量（或说参数）而已。24则（3.5）化为一阶可分离变量方程即易得由初值条件(3.7)即4、解析方法方程（3.5）可以降阶.令感谢你的观看2019年8月21就是导弹轨迹的数学模型。值得注意的是，前面导出的一阶微分方程注意到上式可改写为于是有这样我们又得到一个可分离变量方程积分得利用于是导弹轨迹方程为25感谢你的观看2019年8月21注意到上式可改写为25感谢你的观看2019年8月21设导弹击中敌艇于B(L,H),以Y=H代入(3.9)式,得而导弹击中敌艇的时刻将数据代入(3.10),(3.11)式,得265、数值方法

将初值问题(3.5)~(3.7)化为一阶微分方程组感谢你的观看2019年8月21设导弹击中敌艇于B(L,H),以Y=H代入(3.9)式取自变量y的步长为,于是得节点相应点上的x的值和p的值记为显然,有初值条件我们将介绍两种近似算法来进行数值处理.27设导弹到达处的时刻为那么得到计算的迭代格式ⅠEuler方法

Euler方法十分简单,就是用差商代替微商,即以代之以而代之以这样有感谢你的观看2019年8月21取自变量y的步长为,于于是表3.1是取n=4时的计算结果,读者可以用来检验程序或应用软件的正确性.

表3.1

kykxkpk000013000.052601.50.123905.00.22412011.50.42此时28感谢你的观看2019年8月21于是28感谢你的观看2019年8月21表3.2是对于不同的n值所对应的计算结果.显然,n越大(即h越小),结果就越精确.

表3.2n4812244896120240L11.5215.9617.9720.5522.2523.3323.5824.15T0.1280.1770.2000.2280.2470.2590.2620.268

29注意，由问题（3.1）,(3.3),(3.4)消去t推导出问题（3.5）∼(3.7)是较为巧妙和偶然的.一般而言，一个微分方程组未必能消去一些变元而减少方程的个数。那么，我们能否直接对初值问题(3.1),(3.3),(3.4)进行数值处理呢？答案是肯定的。实际上，只要由方程(3.1),(3.3)解出和的表达式，这样问题变为感谢你的观看2019年8月2129注意，由问题（3.1）,(3.3),(3.4)消去取时间步长对应时导弹轨迹上点的坐标为则Euler格式为当计算到即停止,于是,30感谢你的观看2019年8月21取时间步长对应表3.3和表3.4分别列出了取步长为0.1和0.05时的计算结果:

表3.3ktkxkyk

00.00.000000.0000010.10.0000045.0000020.25.3615489.6794630.322.67495131.21553此时取31感谢你的观看2019年8月21表3.3和表3.4分别列出了取步长为0.

表3.4ktkxkyk

10.050.0000022.5000020.101.0373642.9760730.153.4120567.3504140.207.6461589.4484350.2514.86790110.7579660.3029.19480128.10702此时取表3.5是对应不同的,用Euler法所得相应的步长推进次数n和计算结果.

表3.50.10.050.0050.001n3656278L22.6749529.1948025.6673125.04935T0.251940.324390.285190.2783332感谢你的观看2019年8月21

Euler方法较为简单,但也较为粗糙,对形式较复杂的微分方程更易有较大的误差.人们设计了不少更精确的近似算法,这里我们介绍其中的一种,进一步研究可看参考书.33其中其中而改进的Euler迭代格式则是Ⅱ.改进的Euler方法(预报---校正法)

以一维情况为例,对问题Euler迭代格式是感谢你的观看2019年8月21Euler方法较为简单,但也较为粗糙,对形式由积分表达式的几何意义看,右边为下方的曲边梯形,从图3.2我们可以看出Euler法是用矩形来代替曲边梯形,而改进的Euler法则是用梯形来代替曲边梯形.

对问题(3.18)~(3.20),我们写出相应的改进Euler迭代格式34感谢你的观看2019年8月21由积分表达式34感谢你的观看2019年8月2135感谢你的观看2019年8月2135感谢你的观看2019年8月21表3.6和表3.7分别列出了取步长为0.1和0.05时的计算结果:

表3.6ktkxkyk

00.00.000000.0000010.12.6807744.8397320.212.5752488.2867930.322.07242130.25569此时取表3.7ktkxkyk

10.050.5186822.4880420.102.1059644.9219530.155.0982167.2021340.2010.1609689.0690650.2519.65646108.9789860.3024.24089130.9903036感谢你的观看2019年8月21表3.6和表3.7分别列出了取步长为0.1此时取表3.8是对应不同的,用改进的Euler法所得相应的步长推进次数n和计算结果.

表3.80.10.050.0050.001n3656278L27.0724224.2408925.1355224.98112T0.300800.269340.279280.27757

图3.3画出了导弹轨迹由解析式所给出的精确曲线以及由Euler法和改进的Euler法进行数值计算所给出的近似曲线.37感谢你的观看2019年8月21此时取图3.3画出了导弹轨迹由解析式所给出的精确曲线以6、仿真方法

如果建立微分方程很困难，或者微分方程很复杂而较难作出数值处理，常常可以用仿真方法。所谓仿真方法，顾名思义，指的是模仿真实行为和过程的方法。在这个具体问题中，就是一步步地模拟导弹追踪敌艇的实际过程。而计算机仿真，则是在计算机上通过相应的程序和软件来实现对事件运行的实际过程的模拟。设导弹和敌艇在初始时刻（即t=0时）分别位于P0(0,0)和M0(0,H)。此时，导弹指向M0。而在t=时，导弹的位置P1(x1,y1)，其中，敌艇的位置则为这时导弹沿P1M1

方向飞行，P1M1的倾角为在t=2时，导弹的位置为P(x,y)，其中38感谢你的观看2019年8月216、仿真方法38感谢你的观看2019年8月21此时敌艇位置为,导弹沿P2M2方向飞行(见图3.4).以此方式,一般地,设时,导弹位置为敌艇的位置则为导弹将沿PkMk方向飞行,那么,PkMk的倾角为39感谢你的观看2019年8月21此时敌艇位置为,导弹沿从而时,导弹位置为,其中而敌艇位置为仍然可以如前那样,当时,仿真停止;或者事先给定误差界,当时,仿真停止,这时对于我们用仿真迭代格式(3.33)~(3.36)进行计算,结果与Euler迭代格式的结果完全一致(见表3.3~3.5).这两种迭代格式实际上确实是相同的,建议读者自己验证一下.

值得注意的是,在仿真方法中,我们根本没有用到微分方程组(3.18)~(3.20),却得到了它的一种离散形式,这是十分有意思的.40感谢你的观看2019年8月21从而时,导弹位置为五、实验任务

1.应用数学软件或编制计算程序对问题(3.12)~(3.14)进行数值计算，先运用Euler法，与表3.2以及表3.3的数据比较，并以更小的步长计算结果；再用改进的Euler法计算（步长与Euler法相同）。

2．在本实验介绍的计算过程中，我们是计算到即停止，然后取

,这样做法可能会有不小的误差。有时甚至会出现整体步长改小而结果却未必能改进的情况。由于Euler法或改进的Euler法的计算格式中每一步值的取得仅仅依赖上一步的值，因此在计算过程中改变步长是可行的，即当计算到而y远大于H时，可缩小步长（例如为原来的十分之一）以xy作为新起点继续进行迭代。试用这种变步长方法来改进在任务１中得到的结果。3．如果当基地发射导弹的同时，敌艇立即由仪器发觉。假定敌艇为一高速快艇，它即刻一135km/h的速度与导弹方向垂直的方向逃逸，问导弹何时何地击中快艇？试建立数学模型并求解。4、如果敌艇以135km/h的速度与导弹方向成固定夹角的方向逃逸，问导弹何时何地击中敌艇？试建立数学模型。并选择若干特殊角度进行计算。5、对问题5的结果，你发现敌艇与导弹方向成何夹角逃逸才好？从结论中你又能得到些什么看法？

41感谢你的观看2019年8月21五、实验任务41感谢你的观看2019年8月21微分方程建模42微分方程建模1

当我们描述实际对象的某些特性随时间（空间）而演变的过程、分析它的变化规律、预测它的未来形态、研究它的控制手段时。通常要建立对象的动态模型。

在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题。43感谢你的观看2019年8月21当我们描述实际对象的某些特性随时间（空间）而一、数学建模的基本思维过程

转化实际问题

1、对要讨论的问题所涉及的重要特征进行合理的数学刻画（转化），即用数学语言对问题涉及到的重要特征进行表述.2、寻求的实际问题的文字叙述，利用一些原则或定律，将其转化为数学描述。解数学问题用数学工具求解得到的数学问题。44感谢你的观看2019年8月21一、数学建模的基本思维过程转化实际问题2、寻求的二、微分方程模型

涉及“改变”、“变化”、“增加”、“减少”、“衰变”、“边际”、“速度”、“运动”、“追赶”、“逃跑”、、、等等词语的确定性连续问题。b、微分方程建模的基本手段微元法等a、微分方程建模的对象45感谢你的观看2019年8月21二、微分方程模型涉及“改变”、“变化

1、寻找改变量一般说来微分方程问题都遵循这样的文字等式变化率（微商）=单位增加量--单位减少量等式通常是利用已有的原则或定律。c、微分方程建模的基本规则2、对问题中的特征进行数学刻画3、用微元法建立微分方程；4、确定微分方程的定解条件（初边值条件）；5、求解或讨论方程（数值解或定性理论）6、模型和结果的讨论与分析。46感谢你的观看2019年8月211、寻找改变量一般说来微分方程问题都遵循这样的文字对论文的评价主要以“假设的合理性、建模的创造性、结果的正确性和文字表述的清晰性”为标准。所以，在论文中应努力反映出这些特点。三、

数学建模论文的撰写方法47感谢你的观看2019年8月21对论文的评价主要以“假设的合理性、建模的创造性、结果的正确性论文格式及组成题目摘要，关键词问题重述模型假设分析与建立数学模型数学模型的求解模型检验(总结与评价)模型推广参考文献附录（若有）正文48感谢你的观看2019年8月21论文格式及组成题目模型假设正文7感谢你的观看2019年8月21、题目论文题目是一篇论文给出的涉及论文范围及水平的第一个重要信息。要求简短精练、高度概括、准确得体、恰如其分。既要准确表达论文内容，恰当反映所研究的范围和深度，又要尽可能概括、精练。49感谢你的观看2019年8月211、题目论文题目是一篇论文给出的涉及论文范围及水平的第一个重2、摘要、关键词摘要是论文内容不加注释和评论的简短陈述，其作用是使读者不阅读论文全文既能获得必要的信息。在数学建模论文中，摘要是非常重要的一部分。数学建模论文的摘要应包含以下内容：所研究的实际问题、建立的模型、求解模型的方法、获得的基本结果以及对模型的检验或推广。论文摘要需要概括、简练的语言反映这些内容，尤其要突出论文的优点，如巧妙的建模方法、快速有效的算法、合理的推广等。一般科技论文的摘要要求不列举例证，不出现图、表和数学公式，不自我评价，且字数200以内。前几年，全国大学生数学建模竞赛要求摘要字数应在300字以内。但从2001年开始，为了提高论文评选效率，要求将论文第一页全用作摘要，对字数已无明确限制。在摘要中也可适当出现反映结果数学公式。50感谢你的观看2019年8月212、摘要、关键词摘要是论文内容不加注释和评论的简短陈述，其作3、问题重述数学建模比赛要求解决给定的问题，所以论文中应叙述给定问题。撰写这部分内容时，不要照抄原题，应把握问题的实质，再用较精练的语言叙述问题。51感谢你的观看2019年8月213、问题重述10感谢你的观看2019年8月214、模型假设建模时，要根据问题的特征和建模目的，抓住问题的本质，忽略次要因素，对问题进行必要的简化，做出一些合理的假设。模型假设部分要求用精练、准确的语言列出问题中所给出的假设，以及为了解决问题所做的必要、合理的假设。假设作得不合理或太简单，会导致错误的或无用的模型；假设作得过分详尽，试图把复杂对象的众多因素都考虑进去，会使工作很难或无法继续下去，因此常常需要在合理与简化之间作出恰当的折中。52感谢你的观看2019年8月214、模型假设建模时，要根据问题的特征和建模目的，5、分析与建立模型根据假设，用数学的语言、符号描述对象的内在规律，得到一个数学结构。建模时应尽量采用简单的数学工具，使建立的模型易于被人理解。在撰写这一部分时，对所用的变量、符号、计量单位应作解释，特定的变量和参数应在整篇文章保持一致。为使模型易懂，可借助于适当的图形、表格来描述问题或数据。53感谢你的观看2019年8月215、分析与建立模型12感谢你的观看2019年8月216、模型求解使用各种数学方法或软件包求解数学模型。此部分应包括求解过程的公式推导、算法步骤及计算结果。为求解而编写的计算机程序应放在附录部分。有时需要对求解结果进行数学上的分析，如结果的误差分析、模型对数据的稳定性或灵敏度分析等。54感谢你的观看2019年8月216、模型求解13感谢你的观看2019年8月217、模型检验把求解和分析结果翻译回到实际问题，与实际的现象、数据比较，检验模型的合理性和适用性。如果结果与实际不符，问题常出在模型假设上，应该修改、补充假设，重新建模。这一步对于模型是否真的有用十分关键。55感谢你的观看2019年8月217、模型检验14感谢你的观看2019年8月218、模型推广模型的评价不要流于空泛，需对模型的意义、可信程度、精度等可能问题，需要认真地思考和讨论。将该问题的模型推广到解决更多的类似问题，或讨论给出该模型的更一般情况下的解法，或指出可能的深化、推广及进一步研究的建议。56感谢你的观看2019年8月218、模型推广15感谢你的观看2019年8月219、参考文献在正文中提及或直接引用的材料或原始数据，应注明出处，并将相应的出版物列举在考文献中。需标明出版物名称、页码、著者姓名、出版日期、出版单位等。参考文献按正文中的引用次序列出，其中书籍的表述方式为：[编号]作者，书名[M]，出版地：出版社，出版年。参考文献中期刊杂志论文的表述方式为：[编号]作者，论文名[J]，杂志名，出版年，卷期号：起止页码。参考文献中网上资源的表述方式为：[编号]作者，资源标题，网址，访问时间（年月日）。57感谢你的观看2019年8月219、参考文献在正文中提及或直接引用的材料或原始数据，应注明出10、附录附录是正文的补充，与正文有关而又不便于编入正文的内容都收集在这里。包括：计算机程序、比较重要但数据量较大的中间结果等。为便于阅读，应在源程序中加入足够的注释和说明语句。58感谢你的观看2019年8月2110、附录附录是正文的补充，与正文有关而又不便于编入正文的内59最简单的数学模型之一——“航行问题”用x表示船速，y表示水速，列出方程：答：船速每小时20千米/小时.甲乙两地相距750千米，船从甲到乙顺水航行需30小时，从乙到甲逆水航行需50小时，问船的速度是多少?x=20y=5求解感谢你的观看2019年8月2118最简单的数学模型之一——“航行问题”用x表示船速，y60航行问题建立数学模型的基本步骤

作出简化假设（船速、水速为常数）；

用符号表示有关量（x,y表示船速和水速）；

用物理定律（匀速运动的距离等于速度乘以时间）列出数学式子（二元一次方程）；

求解得到数学解答（x=20,y=5）；

回答原问题（船速每小时20千米/小时）。感谢你的观看2019年8月2119航行问题建立数学模型的基本步骤作出简化假设（船速、水速四、导弹跟踪问题1、实验目的本实验主要涉及常微分方程。通过实验复习微分方程的建模和求解；介绍两种微分方程的数值方法：Euler法和改进的Euler法；还介绍了仿真方法。613、数学模型微分方程建模的方法主要是依据守恒律来建立等量关系式。对于这个问题，寻求等量关系是比较简单的。设坐标系如图3.1所示，取导弹基地为原点O(0,0)，x轴指向正东方,y轴指向正北方。2、实际问题某军的一导弹基地发现正北方向120km处海面上有敌舰一艘以90km/h的速度向正东方向行驶。该基地立即发射导弹跟踪追击敌艇，导弹速度为450km/h，自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌舰？感谢你的观看2019年8月21四、导弹跟踪问题1、实验目的203、数学模型2、实际问题感

当t=0时，导弹位于点O，敌艇位于点A(0,H)，其中H=120（km）。62设导弹在t时刻的位置为P(x(t),y(t))，由题意,感谢你的观看2019年8月21当t=0时，导弹位于点O，敌艇位于点A(0,H方程(3.1),(3.3)连同初值条件构成了一个关于时间变量t的一阶微分方程组的初值问题。其中另外在t时刻,敌艇位置为,其中。由于导弹轨迹的切线方向必须指向敌艇，即直线PM的方向就是导弹轨迹上点P的切线方向,故有或写为63两边对t求导为了寻求x与y的关系,要设法消去变量t,由式(3.2)得感谢你的观看2019年8月21方程(3.1),(3.3)连同初值条件其中22两边对t求导即有把式(3.1)写为代入上式,就得到轨迹方程.这是一个二阶非线性微分方程,加上初值条件,则初值问题64感谢你的观看2019年8月21即有23感谢你的观看2019年8月21就是导弹轨迹的数学模型。值得注意的是，前面导出的一阶微分方程组(3.1)，(3.3)和(3.4)实际上已经是一个数学模型了，不过多一个变量（或说参数）而已。65则（3.5）化为一阶可分离变量方程即易得由初值条件(3.7)即4、解析方法方程（3.5）可以降阶.令感谢你的观看2019年8月21就是导弹轨迹的数学模型。值得注意的是，前面导出的一阶微分方程注意到上式可改写为于是有这样我们又得到一个可分离变量方程积分得利用于是导弹轨迹方程为66感谢你的观看2019年8月21注意到上式可改写为25感谢你的观看2019年8月21设导弹击中敌艇于B(L,H),以Y=H代入(3.9)式,得而导弹击中敌艇的时刻将数据代入(3.10),(3.11)式,得675、数值方法

将初值问题(3.5)~(3.7)化为一阶微分方程组感谢你的观看2019年8月21设导弹击中敌艇于B(L,H),以Y=H代入(3.9)式取自变量y的步长为,于是得节点相应点上的x的值和p的值记为显然,有初值条件我们将介绍两种近似算法来进行数值处理.68设导弹到达处的时刻为那么得到计算的迭代格式ⅠEuler方法

Euler方法十分简单,就是用差商代替微商,即以代之以而代之以这样有感谢你的观看2019年8月21取自变量y的步长为,于于是表3.1是取n=4时的计算结果,读者可以用来检验程序或应用软件的正确性.

表3.1

kykxkpk000013000.052601.50.123905.00.22412011.50.42此时69感谢你的观看2019年8月21于是28感谢你的观看2019年8月21表3.2是对于不同的n值所对应的计算结果.显然,n越大(即h越小),结果就越精确.

表3.2n4812244896120240L11.5215.9617.9720.5522.2523.3323.5824.15T0.1280.1770.2000.2280.2470.2590.2620.268

70注意，由问题（3.1）,(3.3),(3.4)消去t推导出问题（3.5）∼(3.7)是较为巧妙和偶然的.一般而言，一个微分方程组未必能消去一些变元而减少方程的个数。那么，我们能否直接对初值问题(3.1),(3.3),(3.4)进行数值处理呢？答案是肯定的。实际上，只要由方程(3.1),(3.3)解出和的表达式，这样问题变为感谢你的观看2019年8月2129注意，由问题（3.1）,(3.3),(3.4)消去取时间步长对应时导弹轨迹上点的坐标为则Euler格式为当计算到即停止,于是,71感谢你的观看2019年8月21取时间步长对应表3.3和表3.4分别列出了取步长为0.1和0.05时的计算结果:

表3.3ktkxkyk

00.00.000000.0000010.10.0000045.0000020.25.3615489.6794630.322.67495131.21553此时取72感谢你的观看2019年8月21表3.3和表3.4分别列出了取步长为0.

表3.4ktkxkyk

10.050.0000022.5000020.101.0373642.9760730.153.4120567.3504140.207.6461589.4484350.2514.86790110.7579660.3029.19480128.10702此时取表3.5是对应不同的,用Euler法所得相应的步长推进次数n和计算结果.

表3.50.10.050.0050.001n3656278L22.6749529.1948025.6673125.04935T0.251940.324390.285190.2783373感谢你的观看2019年8月21

Euler方法较为简单,但也较为粗糙,对形式较复杂的微分方程更易有较大的误差.人们设计了不少更精确的近似算法,这里我们介绍其中的一种,进一步研究可看参考书.74其中其中而改进的Euler迭代格式则是Ⅱ.改进的Euler方法(预报---校正法)

以一维情况为例,对问题Euler迭代格式是感谢你的观看2019年8月21Euler方法较为简单,但也较为粗糙,对形式由积分表达式的几何意义看,右边为下方的曲边梯形,从图3.2我们可以看出Euler法是用矩形来代替曲边梯形,而改进的Euler法则是用梯形来代替曲边梯形.

对问题(3.18)~(3.20),我们写出相应的改进Euler迭代格式75感谢你的观看2019年8月21由积分表达式34感谢你的观看2019年8月2176感谢你的观看2019年8月2135感谢你的观看2019年8月21表3.6和表3.7分别列出了取步长为0.1和0.05时的计算结果:

表3.6ktkxkyk

00.00.000000.0000010.12.6807744.8397320.212.5752488.2867930.322.07242130.25569此时取表3.7ktkxkyk

10.050.5186822.4880420.102.1059644.9219530.155.0982167.2021340.2010.1609689.06906