科学计算与数学建模
中南大学 郑洲顺教授
课程来源https://www.xuetangx.com/course/CSU07011000630/5882763?channel=learn_title
数学思维、科学计算
案例和课程导入
案例
贷款金额的公式推导,利率R,第k个月剩余多少金额
+思维,工程问题,生活问题,科学、快速、准确—— 转换成数学考虑
利用计算机解决面临的实验研究或理论研究,难以验证的、复杂的问题
+区分,解析解,与数值解
+区分,精确解,近似解(含误差)
数学建模步骤
- 分析问题,所有相关因素,所有因素的关系——转换为数学模型 例如买房—— 总价、贷款金额、首付、利率、还款额、贷多少时间
- 求解后——对解进行分析,可靠性,数值还是解析,数值是否是稳定的、收敛的 检验后再应用;
- 检验时不相符,再进行假设
- 重复上述步骤,直到与实际问题相符
计算和建模形成中心课题
多项式加法
更好的算法?——秦九韶,改写计算的顺序,减少乘法的次数
+问题,用什么标准评判最好准则?
- 算法不同,对初始误差的传播不同
误差
误差定义和分类
认识误差,有些误差不会扩大和缩小,比如模型误差
模型误差,取决于数学方法描述的能力,实际生活的了解能力(专业知识的理解能力)
观察误差(测量误差),初始误差,
不同误差的传播,算法的稳定性定义
初始误差小——传播后误差误差也小;因此算法稳定性与观测误差有关
用有限逼近无限——截断误差,比如有限的泰勒展开、对无理数的分式展开吗。要求,算法具备收敛性
构造算法—— 计算次数少,计算机舍入误差小
认识总结
模型误差,不会扩大缩小;观察误差,算法稳定性关联;截断误差,算法收敛性关联;舍入误差,受计算次数影响
如何度量误差?
绝对误差,误差
绝对误差的上限
相对误差,5公斤差0.2公斤,和100米差0.2米
相对误差,刻画误差的程度,但计算困难
相对误差近似值,绝对误差/x近似值;证明出用相对误差近似值去估计相对误差是可靠的
误差对于数值计算的重要性
误差的传播,算法的稳定性
加法相对误差传播公式—— 推导n个数的和的函数相对误差不超过这个max(n个数的函数相对误差)—— 加法是稳定的;而减法会扩大传播
+重点,减法传播的误差取决于|x1 - x2|,如果两个相近的近似数相减,其相对误差可能会很大
+重点,设计算法中,应尽量避免相近近似数相减
乘法误差传播公式
+重点,绝对值大,误差大
除法误差传播公式——利用泰勒级数展开近似—— 传播系数为函数在近似点的偏导数
+重点,影响出发绝对误差传播的关键因素—— 分母,除数,除数在(0,1)之间时,误差传播产生放大效应,避免用一个极小数作分母;绝对误差与分母大小有重要关系;而相对误差和乘法传播一样是代数和
乘方和开方,一个扩大一个缩小
+重点,开方可以加强精度
+知识,大数吃小数
有些算法正算误差大,反着算误差小。+思考,传播的规律
算法不稳定,失之毫厘谬之千里
总结
尽量减少运算次数,多项式换成秦九韶;
乘方幂次低,减少乘法和加法次数
合理安排运算顺序
级差悬殊,大数小数,通常小数先处理,积累到一定量级
减法,避免相近的近似数相减,用恒等变换解决
乘法,避免绝对值很大的数值作相乘
除法,避免绝对值很小的数作分母
延伸,减法误差的传播解决:
利用恒等变换解决相近近似数相减后,带来的误差传播增大的问题
参考来源:百度文库 > 误差的传播与估计 > https://wenku.baidu.com/view/f2a401895b0102020740be1e650e52ea5418ce43.html