《计算方法(B).ppt》由会员分享,可在线阅读,更多相关《计算方法(B).ppt(27页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS计算方法(B)主讲:张瑞E-Mail:http:/ 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS第0章 绪论l计算方法的作用l计算方法的内容l误差l一些例子数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS现实中,具体的科学、工程问题的解决:实际问
2、题实际问题物理模型物理模型数学模型数学模型数值方法数值方法计算机求结果计算机求结果计算方法是一种研究并解决数学问题的数值近似解近似解方法数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS数值数值分析分析输入复杂问题或运算输入复杂问题或运算 计算机计算机近似解近似解数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS计算方法的特性l理论性:数学基础l实践性计算方法连接了模型到结果的重要环节 随着计算机
3、的飞速发展,数值分析方法已深入到计算物理、计算力学、计算化学、计算生物学、计算经济学等各个领域。本课仅限介绍最常用的数学模型的最基本的数值分析方法。数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS学习的目的、要求l会套用、修改、创建公式l编制程序完成计算课程评分方法课程评分方法(GradingPolicies)总分总分(100)=平时作业平时作业(20)+上机作业上机作业(15)+期末期末(65)数 学 系University of Science and Technology of Chin
4、aDEPARTMENT OF MATHEMATICS3、到网站 http:/ 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS内容2、数值代数线性代数的数值求解,如解线性方程组、数值代数线性代数的数值求解,如解线性方程组、逆矩阵、特征值、特征向量逆矩阵、特征值、特征向量3、微分方程常微分,、微分方程常微分,Runge-Kutta法、积分法法、积分法1、数值逼近数学分析中的数值求解,如微分、积分、数值逼近数学分析中的数值求解,如微分、积分、100亿/秒,算3,000年,而Gauss消元法2660次数
5、 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差l绝对误差设为精确值,为近似值,为误差或绝对误差例如:作Taylor展开,舍弃,即为误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSl相对误差称为相对误差150分满考139,100分满考90,两者的绝对误差分别 为11和10,优劣如何?前者相对误差(150139)/150=0.073,后者相对误差(100-90)/100=0.100数 学
6、系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS有效位数l当x的误差限为某一位的半个单位,则这一位到第一个非零位的位数称位x的有效位数。有效位的多少直接影响到近似值的绝对误差和相对误差 的近似值3.141具有几位有效位数?的近似值3.142具有几位有效位数?数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差来源l原始误差模型误差(忽略次要因素,如空气阻力)物理模型,数学模型l方法误差截断误差(算法本
7、身引起)l计算误差舍入误差(计算机表示数据引起)数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差的运算1、两相近数相减,相对误差增大2、数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS例子例子求根数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS3、小数作除数,绝对误差增大误差的运算数 学
8、系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS一些例子1、则,我们有构造方法如下:1.2.数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSn00.1820.1820.18210.0880.0900.08820.0580.0500.05830.04310.0830.043140.0343-0.1650.034350.02841.0250.028460.024-4.9580.02470.02124.93
9、30.02180.019-124.5400.019数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS原因:对格式1,如果前一步有误差,则被放大5倍加到这一步称为不稳定格式稳定格式,对舍入误差有抑制作用在我们今后的讨论中,在我们今后的讨论中,误差误差将不可回避,将不可回避,算法的算法的稳定性稳定性会是一个非常重要的话题。会是一个非常重要的话题。数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS2、有
10、时候,模型本身就是病态(系数引入小变化,解产生大变化)数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS例:例:蝴蝶效应蝴蝶效应 纽约的一只蝴蝶翅膀一拍,风和日丽的北纽约的一只蝴蝶翅膀一拍,风和日丽的北京就刮起台风来了?!京就刮起台风来了?!NYBJ以上是一个以上是一个病态问题病态问题 /*ill-posed problem*/*ill-posed problem*/关于本身是病态的问题,我们还是留给数学家去头痛吧!关于本身是病态的问题,我们还是留给数学家去头痛吧!数 学 系University
11、 of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSLab01.级数计算级数计算Hamming(1962)x=0.0,0.5,1.0;,10.0,100.0,300.00.绝对误绝对误差小于差小于1.0e-6.输出输出两列两列输输出:出:x和和(x)如如Cfprintf:fprintf(outfile,“%6.2f,%16.12en,x,psix);/*hererepresentsaspace*/数 学 系University of Science and Technology of ChinaDEPARTMENT OF MA
12、THEMATICSSampleOutput(representsaspace)0.00,0.644934066848e+0010.10,0.534607244904e+001.1.00,0.99999900000e+00110.00,0.000000000000e+001.300.00,0.20942212934e-002数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSH.W.给出计算如下式子的方法,以达到相当的精度其中,()、()中x接近,()中xa数 学 系University of Sc
13、ience and Technology of ChinaDEPARTMENT OF MATHEMATICS一些基本数学定理介值定理若f(x)在a,b上连续,则任意C在f(a)与f(b)之间,都存在使 f(c)=C若f(x)在a,b上连续,x1,xn为a,b内的点,g1,gn为同号的实数,则存在 使数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS积分均值定理数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSLagrange中值定理 若f(x)在a,b上连续,在(a,b)内有有限导数f(x),则存在a与b之间的点c,满足数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSTaylor展开
限制150内