1.工程数学研究什么:
利用计算机对各类数学问题进行数值求解的方法,同时对求出的解的精度进行评估,以及如何在计算机上实现求解的过程。
追求稳定的算法和计算复杂度小的算法。
2.误差的来源:
模型误差:数学模型本身存在的,
观察误差:数据来自工程,测量观察存在的误差
方法(截断)误差:n阶泰勒展开,取n项,剩余下称为截断误差。
舍入误差:
3误差
设为精确值,则a为
的一个近似数。
对于近似数a有相对误差和绝对误差
3.1绝对误差
近似数a的绝对误差E(a)为:
![](https://i-blog.csdnimg.cn/blog_migrate/27302009dfebd52590d3bf4458b01283.png)
近似数a的绝对误差等于精确值减去近似值
则近似数a的绝对误差限为:
![](https://i-blog.csdnimg.cn/blog_migrate/e0efa663120c608b12be4f6560fabf4a.png)
近似数a的绝对误差的绝对值小于等于一个值,称该值为近似数a的绝对误差限
绝对误差不能充分说明近似数的精确程度:存在单位,光速和步速比,数值大不代表不精确
3.2相对误差
近似数a的相对误差为:
![](https://i-blog.csdnimg.cn/blog_migrate/d4c5094662b0374eb7b200dc58de8cc1.png)
近似数a的相对误差等于绝对误差与精确值的比
近似数a的相对误差限为:
![](https://i-blog.csdnimg.cn/blog_migrate/3edf4c2a2cc95ed17362b79a8846c663.png)
近似数a的相对误差限
相对误差是属于无量纲的比较:单位上下约去。实际计算精确值取近似值。
3.3有效数字
3.3.1 4舍5入得到的有效数字
四舍五入后得到的近似数,从第一位非零数开始直到最末位,有几位就称该近似数有几位有效数字。
也就是对精确值进行4舍5入,得到的有效数字,有几位是几位。
例1:对进行4舍5入取3位得到3.14,则称3.14近似
有3位有效数字
![](https://i-blog.csdnimg.cn/blog_migrate/47eb304ad4078b09dd7f0288e40582b5.png)
0.002将误差限写成0.5的科学计数法形式
通俗,得到的数据通过4舍5入得到,4舍5入有几位,他的有效数字就有几位。
例2:使3.1416近似,通过4舍5入得到的,有5位有效数字。
![](https://i-blog.csdnimg.cn/blog_migrate/fa7ccaaccf13a71bfb89edbae14987ed.png)
0.000008将误差限写成0.5的科学计数法形式
近似数的有效数字不仅给出了近似数的大小,还指出了他的绝对误差限。
3.4绝对误差与有效数字关系
![](https://i-blog.csdnimg.cn/blog_migrate/c4f1950716a918db021bef0c2d49c448.png)
小数点后第一位不为0,把近似数a写成0.a的形式,在通过绝对误差与有效数字关系求解出n。
3.5例题
3.5.1求绝对误差
已知a和 ,求近似值a的相对误差和绝对误差。
![](https://i-blog.csdnimg.cn/blog_migrate/8d7dc155ddd3004294b33e29510d9b97.png)
![](https://i-blog.csdnimg.cn/blog_migrate/ecf2fe3aa8a2d02bbfc024e9e424d245.png)
a的绝对误差等于精确值减去近似值,同时小于绝对误差限
![](https://i-blog.csdnimg.cn/blog_migrate/535d5bbefa1dbbed1351bfb952b97f41.png)
a的相对误差等与a的绝对误差限除以精确值的近似值
3.5.1求绝对误差
3.5.2求有效数字
![](https://i-blog.csdnimg.cn/blog_migrate/b995a4e1930bb38dbfd3046ffdf98c1e.png)
根据绝对误差与有效数字的关系公式
![](https://i-blog.csdnimg.cn/blog_migrate/fd7167637ae31f4206f1e2feafacb47d.png)
![](https://i-blog.csdnimg.cn/blog_migrate/0e49ba497282ff4e21f3d6784c550b5c.png)
4误差的定性分析和避免误差的危害
为了避免和减少误差,选择稳定的有效的算法
4.1算法的数值稳定性
用一个算法进行计算,初始数据误差(舍入误差),在计算传播过程中,使结果误差增长很快,则称该算法是数值不稳定的,否则称为数值稳定的算法。
4.1.1例题
![](https://i-blog.csdnimg.cn/blog_migrate/ac881761b7ac2908a0cee2e730f400c1.png)
构造En和En-1的式子
![](https://i-blog.csdnimg.cn/blog_migrate/23345543a232c37f69c676dc5283a5ab.jpeg)
![](https://i-blog.csdnimg.cn/blog_migrate/e5bded33703c94c173c6521589fce9fd.png)
![](https://i-blog.csdnimg.cn/blog_migrate/2e7ff8d0dace32f6e846c684feae3201.png)
4.2病态数学问题与条件数
4.2.1病态与条件数介绍
数值稳定性是对算法而言,条件数衡量了函数的相对误差和自变量的变化关系。
病态是数学问题即数学模型本身的性质。
病态数学问题:当输入数据有微小扰动时,会引起解的大扰动。
良态数学问题:当输入数据有微小扰动时,没有引起解的大扰动。
条件数C衡量问题的病态程度,C越大病态可能越严重。
4.2.2条件数定义
y=f(x)的条件数定义为:
![](https://i-blog.csdnimg.cn/blog_migrate/34fcb0cf50299cb4946fc423e4f7451d.png)
![](https://i-blog.csdnimg.cn/blog_migrate/7376ab60ff36671919b37efa402adda6.png)
函数值的相对误差化简成自变量的相对误差与函数值的条件数
例题:
求y=ln(x)的条件数C
![](https://i-blog.csdnimg.cn/blog_migrate/2ecb0f88fddc879056e8ebc888b3316b.jpeg)
![](https://i-blog.csdnimg.cn/blog_migrate/888c09d105ebd12ab2085f2d8b756fff.png)
随着自变量x变化,条件数越接近1越大。
函数值的相对误差与条件数成正比,条件数越大,对数函数值的相对误差就越大。
4.3科学计算避免误差的问题
4.3.1避免两个相近的数相减
![](https://i-blog.csdnimg.cn/blog_migrate/32d6c10500a8d4e823e374559b7c680f.png)
两个相近的数做差导致有效位数的减少
如果无法改变算式,则采用增加有效位数进行计算
![](https://i-blog.csdnimg.cn/blog_migrate/8cbd101edd0792a229a5b7584d2ad390.png)
恒等变换后,有效位数增加了
4.3.2注意简化计算步骤,减少运算次数
4.3.3避免用绝对值很小的数作为除数
较小的误差作分母导致很大的误差。
4.3.4两数相减避免大数吃掉小数
![](https://i-blog.csdnimg.cn/blog_migrate/64d718b33f36c770560964886fdec392.png)
先用数量级高的数做减法运算,降低数量级。在同数量级高的数进行加法运算。
5.参考资料
1.<<数值分析>>北京理工大学出版,感谢曾繁慧 胡行华两位老师
2.辽宁工程技术大学,工程数学考试第一章重点。
1误差与有效数字直接的关系(绝对误差)。
2数值的稳定性(比较哪个更稳定)。
3科学计算需要避免的问题