计算方法公式总结_计算公式总结
计算方法公式总结由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“计算公式总结”。
计算方法公式总结
绪论
exx,x为准确值,x为近似值。绝对误差绝对误差限
r|e||xx|,ε为正数,称为绝对误差限
xxe表示相对误差 通常用exxrxxe相对误差e*xxr相对误差限|er|r或|e|r 有效数字
一元函数y=f(x)
'e(y)f(x)e(x)绝对误差e(y)f(x)'e(x)xf'(x)e(y)er(x)相对误差ryyf(x)二元函数y=f(x1,x2)绝对误差 f(x1,x2)f(x1,x2)e(y)dx1dx2
x1x2f(x1,x2)x1f(x1,x2)x2e(y)er(x1)er(x2)相对误差rx1yx2y
机器数系
注:1.β≥2,且通常取2、4、6、8 2.n为计算机字长
3.指数p称为阶码(指数),有固定上下限L、U 4.尾数部 s0.a1a2an,定位部p
n112(1)(UL1)5.机器数个数机器数误差限
1np舍入绝对 |xfl(x)|截断绝对|x2fl(x)|np
|xfl(x)||xfl(x)|11n1n舍入相对截断相对
|x||x|2
秦九韶算法
方程求根
f(x)(xx)mg(x),g(x)0,x*为f(x)=0的m重根。
二分法
迭代法
f(x)0xk1(xk)
k=0、1、2……
**lim{x}x(x){xk}为迭代序列,(x)为迭代函数,kk
局部收敛
注:如果知道近似值,可以用近似值代替根应用定理3判断是否局部收敛
牛顿迭代法
f(x)f(xk)f(xk)(xxk)0
f(xk)xk1xk'(k0,1,2,)f(xk)注:牛顿迭代对单根重根均局部收敛,只要初值足够靠近真值。
'
牛顿迭代法对初值要求很高,要保证初值在较大范围内也收敛,加如下四个条件
注:证明牛顿迭代法大范围收敛性,要构造一个区间[ε,M(ε)],其中f()M()',在这个区间内验证这四个条件。
f()
如果知道根的位置,构造[ε,M(ε)]时应该包括根,即ε+常数
线性方程组求解
有两种方法:消去法和迭代法
高斯消去法 利用线性代数中初等行变换将增广矩阵转化为等价上三角矩阵。
注意:第一行第一列为0,将第一列不为0的某一行与第一行交换位置,继续初等行变换。对角占优矩阵
a11aA21an1na12a22an2a1na2n ann则称A为按行严格对角占优矩阵 |aii||aij|(i1,2,,n)j1jin|ajj||aij|(j1,2,,n)i1ij则称A为按列严格对角占优矩阵
aijaji(i1,jn)xR,x0,(x,Ax)0
则称A是对称正定的。
当A是上面三种情况时,用高斯消去法消元时追赶法是高斯消元法的一种特例
nakk0,不用换行。
列主元高斯消元法
|aik|,即第k次消元把k~n行第k列绝对值当|ask|maxkin最大的行(s行)调到第k行,再进行高斯消元。(k)(k)
迭代序列构造
AxbxBxfx第三个等式为迭代序列,B为迭代矩阵。迭代收敛判别
1.充分条件:迭代矩阵范数小于1,B1
结论:Ax=b有唯一解x*
(k1)Bx(k)f
2.充要条件:迭代矩阵谱半径小于1,(B)1 Jacobi迭代法
ALDU其中L(low)为下三角,U为上三角,D为对角线元素
迭代格式:x(k1)D(LU)x(k)D1b
1
迭代矩阵JD(LU)
1收敛性判据:
|IJ|0|D||LDU|0|LDU|0
求出最大值小于1(J的谱半径小于1)即迭代格式收敛.1Gau-Seidel迭代法
迭代格式
x(k1)D(Lx1(k1)Ux(k)b)
(k)x(k1)(DL)Ux11(DL)1b
迭代矩阵:G(DL)U
常数矩阵:g(DL)1b
收敛性判据:
|IG|0|(DL)||(DL)U|0|(DL)U|0
求出最大值小于1(G的谱半径小于1)即迭代格式收敛.结论:当A是严格对角占优的,则Jacobi和Gau-Seidal迭代法均是收敛的1插值法
用插值多项式p(x)代替被插函数f(x)
nP(x)aaxax插值多项式:,01nn+1个点P(xi)yi(i0n)
插值区间:[a,b],插值点满足
ax0x1xnb
求插值多项式P(x),即求多项式系数的过程为插值法
带入可知求系数的插值点行列式为范德蒙行列式,不为0,有唯一解。即n+1插值条件对应的不超过n次的插值函数P(x)只有一个。一次线性插值nxx0xx1Py0y1y0l0(x)y1l1(x)1(x)x0x1x1x0(xxi)lk(x)i0(xx)(xkxi)ikki
ni0iki0ikn(xxi)Lagrange插值多项式
Ln(x)yklk(x)k0k0 nnxxi()yki0xxiikkn插值余项
非插值节点上Lagrange插值多项式为被插函数f(x)的近似值
f(n1)()nRn(x)f(x)Ln(x)(xxi)(n1)!i0(a,b)
带导数插值条件的余项估计
注:推导过程用罗尔中值定理构造辅助函数
(t)Rn(t)K(x)Wn1(t)
第二条性质用于可以证明阶数不大于n的f(x)的插值余项为0.差商和Newton插值法
记忆方法:先记分母,最后一个减去第一个,对应的分子第一项是最后一个临近k元素的差商,第二项是第一个临近k个元素的差商。
牛顿插值多项式
通常记作Nn(x)分段样条插值
分段二次样条插值
讨论n为奇偶情况时的三个点 余项估计式
三次样条插值函数
第一类边界条件(端点一阶导数已知)
D0等于第一个式子,dn等于第二个式子
自然边界条件(端点二阶导数已知二阶导数和M0,Mn=0)
曲线拟合最小二乘原理
函数关于n个点线性无关
23n1,x,x,x,,x注:线性无关的函数为才是最小二乘多项式
注:记住公式即可。
数值积分和数值微分
xk为求积节点,Ak为求积系数。
插值求积公式
梯形公式
Simpson公式
Cotes公式
截断误差
代数精度
当f(x)为不超过m次多项式时上式成立,f(x)为m+1多项式时上式不成立。则称为求积公式有m次代数精度。
梯形公式代数精度为1,Simpson公式代数精度为3,Cotes公式代数精度为5
截断误差 梯形公式
Simpson公式
Cotes公式
Gau求积公式
求积公式代数精度为2n+1 [-1,1]上的两点Gau公式(3次代数精度)
111f(x)dxf(3)f(3)1[-1,1]上的三点Gau公式(5次代数精度)
538531f(x)dx9f(5)9f(0)9f(5)1
记住 xktk,AkAk的关系,tkAk查表即可
复化梯形公式2阶,复化Simpson公式4阶,复化Cote公式6阶
计算机通过不断把区间二分,所得前后两次积分差值满足精度条件即可
1|I2n(f)In(f)|时 给定精度ε,p211|I(f)I2n(f)|p|I2n(f)In(f)|21因而可以取I2n(f)为I(f)的近似值。
梯形
Simpson数值微分
数值微分截断误差
中点公式:
f(x0h)f(x0h)D(h) 2h常微分方程数值解法
Euler方法
欧拉公式(单步显式公式)求出的近似解
局部截断误差
Euler公式的局部截断误差(一阶精度)
后退Euler公式
梯形公式(二阶精度)
改进Euler公式(二阶精度)
截断误差(推导要求掌握,利用梯形和Euler公式的截断误差)
材料费用的分配1.按原材料定额消耗比例分配原材料费用某产品原材料定额消耗量=该产品实际产量X单位产品原材料定额消耗量 原材料消耗量分配率=原材料实际消耗总量/各种产品......
培训讲义 —— 灭火器材的配置与计算灭火器 的配置 与计算 §1.火灾的分类和灭火器配置场所的危险等级:一、火灾的分类(即根据火灾中燃烧物质的特征进行分类。按国家标准(GBl......
第一章:基本概念x1x2...xm.xm1xm2...xmn 1.xx1x2...xm.xm1xm2...xmnxmn1x若xx1mn及其以前的非零数字称为准确数字。 准确到n位小数,x10n,称x2各位数字都准确的近似数称为有效数......
1.何为有根区间给定一个方程f(x)=0,如果f(x)在[a,b]上连续,又f(a).f(b)3.作图法寻找有根区间适用于哪种情况函数f(x)比较简单时适用4.对于已知方程,如何利用逐步搜索法在区间......
高中化学计算方法总结高中化学教师,在开展计算教学时,应该引导学生掌握常见的解题方法与解题技巧,以促进教学效果的提升。下面为大家总结了高中化学几种计算方法,希望帮助到......
