计算方法公式总结_第1页
计算方法公式总结_第2页
计算方法公式总结_第3页
计算方法公式总结_第4页
计算方法公式总结_第5页
已阅读5页,还剩26页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、计算方法公式总结绪论绝对误差 ,为准确值,为近似值。绝对误差限 ,为正数,称为绝对误差限相对误差 通常用表示相对误差相对误差限或有效数字一元函数y=f(x)绝对误差 相对误差 二元函数y=f(x1,x2)绝对误差 相对误差机器数系注:1. 2,且通常取2、4、6、82. n为计算机字长3. 指数p称为阶码(指数),有固定上下限L、U4. 尾数部 ,定位部5. 机器数个数机器数误差限舍入绝对 截断绝对舍入相对 截断相对秦九韶算法方程求根 ,为f(x)=0的m重根。二分法迭代法 k=0、1、2为迭代序列,为迭代函数, 局部收敛注:如果知道近似值,可以用近似值代替根应用定理3判断是否局部收敛牛顿迭代

2、法注:牛顿迭代对单根重根均局部收敛,只要初值足够靠近真值。牛顿迭代法对初值要求很高,要保证初值在较大范围内也收敛,加如下四个条件注:证明牛顿迭代法大范围收敛性,要构造一个区间,M(),其中,在这个区间内验证这四个条件。 如果知道根的位置,构造,M()时应该包括根,即+常数线性方程组求解有两种方法:消去法和迭代法高斯消去法利用线性代数中初等行变换将增广矩阵转化为等价上三角矩阵。注意:第一行第一列为0,将第一列不为0的某一行与第一行交换位置,继续初等行变换。对角占优矩阵则称A为按行严格对角占优矩阵 则称A为按列严格对角占优矩阵则称A是对称正定的。当A是上面三种情况时,用高斯消去法消元时,不用换行。

3、追赶法是高斯消元法的一种特例列主元高斯消元法当,即第k次消元把kn行第k列绝对值最大的行(s行)调到第k行,再进行高斯消元。迭代序列构造第三个等式为迭代序列,B为迭代矩阵。迭代收敛判别1. 充分条件:迭代矩阵范数小于1,结论:Ax=b有唯一解x* 2. 充要条件:迭代矩阵谱半径小于1,Jacobi迭代法其中(low)为下三角,为上三角,为对角线元素迭代格式:迭代矩阵收敛性判据:求出最大值小于1(J的谱半径小于1)即迭代格式收敛.Gauss-Seidel迭代法 迭代格式 迭代矩阵:常数矩阵:收敛性判据:求出最大值小于1(G的谱半径小于1)即迭代格式收敛.结论:当A是严格对角占优的,则Jacobi

4、和Gauss-Seidal迭代法均是收敛的插值法用插值多项式p(x)代替被插函数f(x)插值多项式:,n+1个点插值区间:,插值点满足求插值多项式P(x),即求多项式系数的过程为插值法带入可知求系数的插值点行列式为范德蒙行列式,不为0,有唯一解。即n+1插值条件对应的不超过n次的插值函数P(x)只有一个。一次线性插值Lagrange插值多项式插值余项非插值节点上Lagrange插值多项式为被插函数f(x)的近似值带导数插值条件的余项估计注:推导过程用罗尔中值定理构造辅助函数第二条性质用于可以证明阶数不大于n的f(x)的插值余项为0.差商和Newton插值法记忆方法:先记分母,最后一个减去第一个

5、,对应的分子第一项是最后一个临近 k元素的差商,第二项是第一个临近k个元素的差商。牛顿插值多项式通常记作Nn(x)分段样条插值分段二次样条插值讨论n为奇偶情况时的三个点余项估计式三次样条插值函数第一类边界条件(端点一阶导数已知) D0等于第一个式子,dn等于第二个式子自然边界条件(端点二阶导数已知 二阶导数和M0,Mn=0)曲线拟合最小二乘原理函数关于n个点线性无关注:线性无关的函数为才是最小二乘多项式注:记住公式即可。数值积分和数值微分为求积节点,为求积系数。插值求积公式梯形公式Simpson公式Cotes公式截断误差代数精度当f(x)为不超过m次多项式时上式成立,f(x)为m+1多项式时上

6、式不成立。则称为求积公式有m次代数精度。梯形公式代数精度为1,Simpson公式代数精度为3,Cotes公式代数精度为5截断误差梯形公式Simpson公式Cotes公式Gauss求积公式求积公式代数精度为2n+1-1,1上的两点Gauss公式(3次代数精度)-1,1上的三点Gauss公式(5次代数精度)记住,的关系,查表即可复化梯形公式2阶,复化Simpson公式4阶,复化Cote公式6阶计算机通过不断把区间二分,所得前后两次积分差值满足精度条件即可给定精度,时因而可以取为的近似值。梯形Simpson数值微分数值微分截断误差中点公式:常微分方程数值解法Euler方法欧拉公式(单步显式公式)求出的

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论