§ 避免绝对值小的数作除数 § 避免两个相近的数据相减 § 要防止大数“吃掉”小数 § 尽量减少计算工作量 § 选用数值稳定性好的算法

Romberg求积方法是在积分区间逐次分 半的过程中利用外推法产生的一种数值积 分方法,当被积函数的光滑性条件满足时, 可以得到较精确的积分近似法

本文简要回顾了计算数学和科学计算的发展历史，从数值代数、数值逼近、最优化方法和微分方程计算方法等4个方面简要介绍了计算数学的主要研究内容和最新进展，最后提出进一步发展我国科学计算的建议

一、求积公式的代数精度 二、插值型求积公式 三、复化求积法

提出了点接触热弹流润滑的简化数值解。在等温弹流润滑完全数值解的计算结果基础上,通过联立求解能量方程和热界面方程等,计算出重载条件下点接触弹流润滑状况下的三维温度分布。研究了热弹流问题的温升规律,并提出了油膜和接触体表面的最大温升公式。该公式得到了弹流测温实验的初步验证

由微积分学基本定理,当f(x)在[a,b]上连续时,存在原函数F(x) 由 NewtonLeibnitsI-式if(x)df()-F(a) 有时用上面的方法计算定积分有困难

以某电厂的脱硝工程为背景,对烟道、反应器和支架这种大型复杂结构进行了整体分析,并以烟道出口和反应器作为研究对象,分析了支架变形对烟道出口和反应器的影响.结果表明,由于支架变形的影响,烟道和反应器的变形分布由对称分布变为非对称,平均变形值要比不考虑时的大,但是支架的变形对烟道和反应器的应力值及分布基本无影响.因此在以强度为控制因素的设计中,烟道和反应器进行单独计算的设计方法是可行的

根据中锰钢热轧组织结构确立两相区奥氏体化的几何模型和初始条件,利用DICTRA动力学分析软件对中锰钢马氏体基体奥氏体化过程进行计算分析.在奥氏体化初期的形核过程中,马氏体中过饱和的碳锰元素从铁素体迅速转移到奥氏体并在相界面奥氏体一侧聚集.后续的相变过程中,碳在奥氏体中快速均化,但锰在相界面奥氏体一侧的聚集加剧.相变初期奥氏体界面推移速度比中后期高出若干个数量级,但随时间推移迅速衰减.相变初期相界面推移是碳扩散主导,相变后期界面推移受到锰在奥氏体中扩散速度制约.温度升高可显著提高相界面推移速度.达到相同数量奥氏体的情况下,低温长时退火有利于锰从铁素体向奥氏体转移并提高其在奥氏体中的富集度,从而提高奥氏体的稳定性

可以说,误差分析是计算方法研究的无底洞,所以大家既 要提得起,也要放得下。提得起,就是基本概念,思想,方法 要领会得了;放得下,就是不主张钻牛角尖。 本章的基本概念有:截断误差,舍入误差,绝对误差,相 对误差,误差限,有效数字等。 四则运算的误差估计仅理解几个重要结论即可;利用微分 估算误差既有理论意义,又有实用价值,应当牢记;有效数字 与相对误差的关系在一般的教科书中都很重要,我并不这样认 为:既没有理论价值,程序设计也几乎不会用到它们;本次礼 包给了几个典型例题,都不难,应当没有问题

§9.1引言 §9.2 Euler方法 §9.3 Runge-Kutta公式 §9.4单步法的进一步讨论 §9.5线性多步法