一、微分方程的解析解方法 二、微分方程问题的数值解法 1 微分方程问题算法概述 2 四阶定步长 Runge-Kutta-算法及 MATLAB实现

为进一步揭示远场条件下金属矿山崩落矿岩运移演化机理，综合利用物理试验、数值模拟和理论分析等手段，构建单口放矿模型开展近?远场崩落矿岩流动特性研究。首次基于离散元软件PFC3D和刚性块体模型构建放矿数值模型，并通过近场放矿物理试验与模拟结果的对比分析，证明了刚性块体模型在崩落矿岩流动特性研究中的可靠性与优越性。在此基础上，对远场条件下松动体形态变化规律、矿岩流动体系内的应力演化规律及其力学机理进行了量化研究。研究结果表明：1）近?远场条件下的松动体形态变化均符合倒置水滴理论。在放矿初始阶段，松动体最大宽度随高度增大呈幂函数形式快速增加；随后，松动体最大宽度随高度增大而近似线性增加。2）崩落矿岩流动过程中存在明显的应力拱效应。随着矿岩散体松动范围不断扩大，松动体外围一定范围内的垂直应力均呈明显下降趋势，水平应力逐渐增大并在松动区域到达前出现激增现象；而松动体内的水平应力与垂直应力则急剧下降至较低水平

在冷连轧无取向硅钢薄带过程中，为了实现锥形工作辊窜动自动控制边降，需要合理的确定功效系数与策略。这种系数的获得，不只需要研究本道次的轧辊弹性变形、薄带横向流动、机架间变形对窜辊效率的影响，更重要的是需研究上游机架窜辊对下游机架的影响。这就需要高效的仿真模型来完成以上计算。基于边降区域的金属横向流动理论，建立了将横向流动视为纯剪切增量的数值模型，避免了沿带宽方向建立刚度矩阵，从而提高了计算效率。同时考虑了薄带在机架间发生的轧后屈服流动，由于锥形工作辊窜动，打破了带钢断面的等比例遗传关系，使得轧后带钢在边部区域需要缩宽并减薄来补偿边部延伸率差。所建立的数值模型通过工业现场实验验证，相比于原有模型具有更高的精度。完成了两个机架连续计算，研究了上游机架窜辊对下游机架出口边降的影响。研究发现，第一机架的边降控制范围最宽，第二、三机架控制范围逐渐变窄。根据该规律设计了根据三点边降偏差的配合调控策略，相比单点策略在工业应用中取得了更好效果

1、引言 一阶常微分方程的初值问题:

武汉大学：《数值分析》课程教学资源（PPT课件讲稿）第二章 求解线性方程组的数值解法 2.2 解线性方程组的迭代法

只要近似地算出右边的积分Ikfx,y(x)dx,则可通 过yi1=y+k近似(x).而选用不同近似式Ik,可得到不同的计算公式

定义在假设y=y(x),即第i步计算是精确的前提下,考 虑的截断误差R=y(x)-y称为局部截断误差鬥 local truncation error

高次插值有 Runge现象,故采用分段低次插值→分段低次合成的 Newton- Cotes复合求积公式

利用插值多项式P(x)则积分易算在[a,b]上取a≤x