7.1 引言 7.3 欧拉法 7.5 单步法的收敛性与稳定性 7.2 离散变量法 7.6 线性多步法 7.4 龙格-库塔法 7.7 方程组与高阶方程初值问题的数值解法 7.8 数值实验

6.1 初值问题的Euler方法 6.2 Runge-Kutta方法 6.3 线形多步法 6.4 一阶常微分方程组数值解法 6.5 常微分方程边值问题的数值解法

1. 函数f(x)以离散点列给出时，而要求我们给出导数值， 2. 函数f(x)过于复杂 这两种情况都要求我们用数值的方法求函数的导数值 微积分中，关于导数的定义如下：

Tavlor公式 多项式是一类很重要的函数,其明显特点是结构 简单,因此无论是数值计算还是理论分析都比较方便 从计算的角度看,只须加、减、乘三种运算,连除法 都不需要,这是其它函数所不具备的优点。 用多项式近似地表示给定函数的问题不仅具有实 用价值,而且更具有理论价值。一般的函数不好处理 先用较好处理的多项式近似替代,然后通过某种极限 手续再过渡到一般的函数

1引言 问题的提出 在实际问题中常遇到这样的函数 J=(x),其在某个区间[a,b上 是存在的。但是,通过观察或测量或 试验只能得到在[,b区间上有限个 离散点x0x1xn上的函数值 y;=f(x;), (=0,…,n)或者∫(x)的函数表达 式是已知的,但却很复杂而 不便于计算,希望用一个简单的函数 来描述它

工程实践中有多种振动问题,如桥梁或建筑物的振 动,机械机件、飞机机翼的振动,及一些稳定性分析和 相关分析可转化为求矩阵特征值与特征向量的问题 1已知A=(an)mn,求代数方程q(1)=det(I-A)=0 的根。φ(λ)称为舶特征多项式,一般有n个零点,称 为的特征值 2设为伯的特征值,求相应的齐次方程(4/-A)x=0 的非零解(即求Ax=λx的非零解),x称为矩阵A对应 于孔的特征向量

工程实际计算中，线性方程组的系数矩阵常常具有对 称正定性，即其各阶顺序主子式及全部特征值均大于零。 矩阵的这一特性使它的三角分解也有更简单的形式，从而 导出一些特殊的解法，如平方根法与改进的平方根法

周期性非均质复合材料具有微观结构特征，需要均匀化理论进行宏观和微观的多尺度分析来研究其性能表现。针对其耐久强度性能，应用塑性极限安定下限定理，特别分析了其在长期交变载荷下的安定状态。结合工程应用目标，提出一种全新的代表性单元边界条件，结合圆锥二次优化算法进行数值计算，可以从材料微结构和组分性能出发，经过弹性应力场求解确定位移边界载荷数值，最终由优化求解得到复合材料板材的面内塑性性能容许域。所求得的应力域以单向应力为基，可根据结构宏观的单向应力状态变化幅值直接进行安定状态与否的判定。通过文中的多个算例，验证了所编写的软件及计算流程的可行性及数值准确性，展示了该方法在工程模型中的应用场合和工程实践意义

工程实际计算中,线性方程组的系数矩阵常常具有对 称正定性,即其各阶顺序主子式及全部特征值均大于零 矩阵的这一特性使它的三角分解也有更简单的形式,从而 导出一些特殊的解法,如平方根法与改进的平方根法 定理:设A是对称正定矩阵,则存在唯一的非 奇异下三角阵L,使得 A=LL 且L的对角元素皆为正

一、基本QR方法 60年代出现的QR算法是目前计算中小型矩阵的全部特征值与 特征向量的最有效方法。实矩阵、非奇异。 理论依据:任一非奇异实矩阵都可分解成一个正交矩阵Q和 一个上三角矩阵R的乘积,而且当R的对角元符号取定时,分解是 唯一的