基本数学运算1.pptVIP

* * 第一章　基本的数值运算 §1.1 数值微分 §1.2 数值积分 §1.3 函数方程求根 §1.4 双原子分子振动能级的半经典量子化 §1.5 中心位势的散射 在数值运算的大多数计算中, 有三种最基本的数值运算, 那就是数值微分, 数值积分, 数值求根. §1.0 插值方法 § 1.0 插值方法 当我们需要从一些分离的数据推断他的局部信息时, 就需要插值. 拉格朗日(Lagrange)插值 已知分离点 要求 y(x) 的最简单方法是线性插值 如果用三点 来插值 如果是 n +1 点插值 拉格朗日 n 次插值公式 牛顿插值法 由拉格朗日法, 如果增减节点数, 要重新计算基函数. 用牛顿法可以解决这个问题. 系数 a 由以下插值条件决定 一阶差商 同理可得 二阶差商 若定义 k 阶差商为 那么有 这样我们就得到牛顿插值公式 牛顿法的优点是: 当增加一个节点时, 只要增加一项就行了, 且有递推关系 如果不仅给出各点的函数值, 而且给出各点导数值.这时可使用 2n+1 阶 Hermite 插值. §1.1 数值微分 利用泰勒级数来求数值微分公式 差分公式 ①对(1.2)式取一级近似，得 ②如果对(1.2)式取二级近似(假设f”i存在)，所得两式相减，得 ③如果对(1.2)、(1.3)式取四级近似，容易证明下述公式 ④对(1.2)式取三级近似，容易证明，二阶导数的三点中心差分公式为 ⑤若函数的各阶导数都存在，利用(1.2)、(1.3)式的适当组合，就可得到各阶导数的4点、5点差分公式，其表达式见表1.2。高阶差分公式不一定比低阶差分公式好。如果函数具有连续高阶导数选取高阶差分公式较好，如果不具高阶导数选取低阶差分公式较好。 表1.2 ⑥如果不是等间隔格点，或者所求导数不在格点上，则常用拉格朗日插值公式来计算微分。三点拉格朗日插值公式为 §1.2 数值积分 ①梯形积分：对于积分 ②定步长辛卜生积分：若利用泰勒展开式(1.2)和差分公式(1.6)(1.8)式，得 这就是定步长辛卜生积分公式，其误差关于h是五阶的。它要求被积函数必须存在1阶和2阶导数。 辛卜生积分方法：定步长方法，变步长方法，自适应积分方法。更高阶的积分方法：龙贝格积分方法，高斯型积分方法，… 。 积分方法选取原则： (1). 若被积函数具有任意阶导数，且无很大的尖峰存在，通常使用高阶的龙贝格积分方法。若想使积分花最少的时间，最好使用固定的N点高斯积分(见第三章)方法。因为它的积分精度最高。但条件是能预知N取多大，就能满足精度要求。 (2). 若被积函数具有少数几个较高较窄的峰或谷，往往采用自适应辛卜生积分方法。因为这种方法是专门为这种情况设计的。 (3). 若被积函数不具有高阶导数，但仍具有1阶导数，通常采用变步长辛卜生积分方法。 (4). 若被积函数不光滑，甚至不连续，一般采用梯形积分方法较为可靠。 另外，对于高振荡函数，还有专门的积分方法。 §1.3 函数方程求根 简单搜索法、牛顿法、弦割法和连分式方法。其它十几种方法。 一、简单搜索法(即对分法或二分法)。 SQR1.FOR CALL SQR1(2.0,1.E-5) END FUNCTION F(X) F=X*X-5. N=N+1 WRITE(7,*)N,X,F END SUBROUTINE SQR1(X,EPS) DX=.5 F1=F(X) 10 IF(DX.LE.EPS)GOTO 20 X=X+DX I=I+1 F2=F(X) IF(F1*F2.GT.0.)GOTO 10 X=X-DX DX=.5*DX GOTO 10 20 CONTINUE END N x f(x) 1 2.000000 -1.000000 2 2.500000 1.250000 3 2.250000 6.250000E-02 4 2.125000 -4.843750E-01 5 2.250000 6.250000E-02 6 2.187500 -2.148438E-01 7 2.250000 6.250000E-02 8 2.218750 -7.714844E-02 9 2.250000 6.250000E-02 10 2.234375 -7.568359E-03 11 2.2