计算声学第一章数值计算中的误差分析.ppt

《计算声学第一章数值计算中的误差分析.ppt》由会员分享，可在线阅读，更多相关《计算声学第一章数值计算中的误差分析.ppt（60页珍藏版）》请在三一文库上搜索。

1、前言,课程目的和任务： 通过对一些基本声学和水声学问题的分析和求解，掌握基本声学理论计算与工程研究中常用的数值计算方法，培养综合运用声学专业知识、数学知识和计算机技术解决科学研究中手工所不能解算的问题，具备应用现代计算工具解决工程实际问题的能力。,水声学主要研究声波在水下的辐射、传播与接收，用以解决与水下目标探测和信息传输过程有关的各种声学问题。声波是目前在海洋中唯一能够远距离传播的能量辐射形式。因此作为信息载体的声波，在海洋中所形成的声场时空结构，就成为近代水声学的基本研究内容，而提取海洋中声场信息的结构是我们用来进行水下探测、识别、通信及环境监测等的手段。,前言,前言,海洋环境,波动方程：

2、 波动方程是声学量在声场中满足的基本关系式，反映了波动 特征，也是进行声场计算的基本关系式。在导出波动方程前， 为了使问题简化，需要对介质和声波做一些假设： （1）介质是均匀连续的，即在波长数量级距离内，介质的声学性质保持不变； （2）介质是理想流体介质，声波在其中传播时没有能量损耗，即忽略介质的粘滞性和热传导性； （3）研究小振幅波的传播规律，所谓小振幅波是指各声学量都是一级小量。 波动方程是描述波动运动的数学表达式，它由连续性方程、状 态方程和运动方程推导得到。,前言,波动方程： 理想流体介质中小振幅平面波的波动方程为（沿 轴向传播）： 小振幅声压在三维坐标下的波动方程为 为拉普拉斯算符，

3、在直角坐标系中,前言,海洋声场的数值预报 在建立了能够反映海洋环境因素对声场的制约关系的声 场物理模型（波动方程+定解条件）的基础上，根据可测海 洋环境参数的测定值或预报值，编写程序完成数值计算，给 出相应海洋环境条件下的有关场值。近年来，由于计算机的 快速发展，数值计算声场是一个快速发展的领域。 海洋声场的数值预报方法主要有射线算法、简正波算法、 抛物方程（PE）算法、快速场（FFP）算法等，各自有不同 的适应范围。,前言,前言,函数插值： 已知一组不同深度处的声速值，如何得到任意深度处的 声速值？,前言,前言,数值积分：声线轨迹计算,声线从深度,传播到深度,所经过的水平,距离为,前言,问题

4、：利用射线声学模型对海洋声场进行求解,伪彩色图,前言,前言,三维环境下声传播,前言,三维海洋环境下特征声线求解： 为声线的位置信息，需要求解，其它参数已知。,前言,三维海洋环境下特征声线求解（线性方程组、非线性方程、非线性方程组） 1. 牛顿法迭代法： 泰勒级数展开式的线性部分近似 2. 进化算法： 遗传算法、模拟退化算法、粒子群算法等,前言,曲线拟合：已知目标散射场指向性的实验测量结果如图所示，如何比对其与理论计算结果的误差？,铝球散射声场指向性,频率,kHz,微分方程求解：随机共振系统对微弱信号的检测非线性双稳态随机共振系统 利用四阶龙格库塔算法求解,前言,前言,四阶龙格库塔算法,前言,前

5、言,前言,前言,必要性： 现代科学研究和高技术的发展越来越需要借助计算机进 行数值计算，水声领域也不例外。 讲授的主要内容： 1、数值计算方法：误差分析、方程组求解、非线性方程求解、插值法、最小二乘与曲线拟合、数值微积分、常微分方程求解； 2、进化算法（方程组求解）：量子粒子群算法； 3、虚源法声场建模。,数值计算的对象、任务与特点,对象： 数值计算方法是研究科学与工程技术中数学问题的数值 解及其理论的一个分支，涉及代数、微积分、微分方程等的 数值解问题。 任务： 研究适合在计算机上使用的数值计算方法及相关理论， 如方法的收敛性、稳定性和误差分析等；还要根据计算机的 特点研究如何设计计算方法做

6、到计算时间短、占用内存小。 学习目的： 提高应用计算机解决实际问题的能力。,数值计算的对象、任务与特点,数值计算流程： 特点： 既具有数学的抽象性与严格性，又具有应用的广泛性与 实际实验的技术性，是一门与计算机紧密结合的实用性很强 的有着自身研究方法与理论体系的计算数学课程。,数值计算中的误差分析,内容提要： 掌握绝对误差、相对误差、有效数字、数值计算的误差 估计以及设计算法的原则。 重点内容： 绝对误差、相对误差、有效数字的概念，数值计算的误 差估计。,误差与数值计算的误差估计,一、误差的来源与分类 将一个数的准确值与其近似值之差称为误差。 1. 分类 过失误差： 人为造成，可以避免 非过失

7、误差： 无法避免，分析产生原因，限制在许可范围之内,2 误差与数值计算的误差估计,2. 误差来源（非过失误差） 模型误差： 数学模型是通过对实际问题进行抽象和简化建立的，是一种近似描述。 观测误差： 测量工具精度与测量手段的限制。 舍入误差： 计算机位数的限制，由于计算机的字长是有限的，对参与计算的数据和最后得到的计算结果，都必然用有限位小数代替无穷位小数。,2 误差与数值计算的误差估计,截断误差： 由数值方法求得的数学问题的近似解与数学模型的精确解之间的误差，是数值计算方法固有的。 取部分和作近似 截断误差：,绝对误差与绝对误差限 绝对误差： 设某一量的精确值为 ，其近似值为 ，则称 为近似

8、值 的绝对误差，简称误差。 时称 为弱近似值或亏近似值； 时称 为强近似值或盈近似值。 绝对误差限：如果存在 ，使得 ，则称 为近似值 的绝对误差限，简称误差限或精度（测量时， 测量工具最小刻度的一半）。 越小，表示近似值 的精度 越高。在工程技术上常用 表示近似值的精度或精确 值的范围。,2 误差与数值计算的误差估计,2 误差与数值计算的误差估计,例：用毫米刻度的尺子测得桌子长度近似值为 mm，由尺子的精度可以知道，近似值的误差不超过0.5mm，即 表明精确值 在区间 内，可以写成 绝对误差限 mm，即绝对误差限是末位的半个单位。,相对误差和相对误差限 相对误差：绝对误差与精确值之比，即 称

9、为近似值的相对误差。 实际中，由于精确值 一般无法知道，所以常取 作为近似值 的相对误差。 相对误差限：若存在 ，使得 ，则称 为近似值 的相对误差限。 注意：绝对误差和绝对误差限与 有相同的量纲，相对误差 和相对误差限是无量纲的，工程中常以百分数来表示。,2 误差与数值计算的误差估计,例1.1 国际大地测量学会建议光速采用 其含义是绝对误差限为多少？而其相对误差限为多少？,2 误差与数值计算的误差估计,2 误差与数值计算的误差估计,绝对误差限： 近似值： 相对误差限：,有效数字 如果近似值 的绝对误差限是某一位的半个单位，就称其 “准确”到这一位，且从该位开始直到 的第一位非零数字共 有n位

10、，则称近似数 有n位有效数字。 有效数字既能表示近似值的大小，又能表示其精确程度 （绝对误差限）。 例1.2 设 ，其近似值 ，问 有 几位有效数字？如果 ， 有几位有效数字？,2 误差与数值计算的误差估计,练习题,1.指出如下有效数的绝对误差限、相对误差限和有效数字位数。 2.将22/7作为 的近似值，它有几位有效数字？绝对误差限和相对误差限各为多少？,2 误差与数值计算的误差估计,有效数字 设 的近似值 可以写成如下的标准形式 所以当其绝对误差限满足 时，称近似值 具 有 位有效数字。,2 误差与数值计算的误差估计,结论：如果 ，有 位有效数字，则其相对误差限为 反之，如果 的相对误差限满

11、足 则 至少有 位有效数字。,2 误差与数值计算的误差估计,例1.3 要使 的近似值的相对误差小于 ，至少需取几位有效数字？,误差的传播与估计 实际的数值计算中，参与运算的数据往往都是近似值，带有误差。而在进一步运算中都会产生舍入误差或截断误差，这些误差在运算过程中会进行传播，影响计算结果。,2 误差与数值计算的误差估计,一元函数的泰勒（Taylor）中值定理： 如果函数 在区间 内有直到 阶导数， ，则有 其中，拉格朗日型余项 （ 介于 之间）。,2 误差与数值计算的误差估计,2 误差与数值计算的误差估计,泰勒（Taylor）公式估计误差的方法： 以二元函数 为例，设 和 分别是 和 的近

12、似值， 是函数值 的近似值。函数 在点 处 的Taylor展开式为,2 误差与数值计算的误差估计,式中 一般都是小量值，如 果忽略它们的高阶无穷小量，则上式简化为 因此 的绝对误差为,2 误差与数值计算的误差估计,系数 分别是一阶偏导数在 处的 值，称为 对 的绝对误差的增长因子，分别 表示绝对误差 经过传播后增大或缩小的 倍数。,2 误差与数值计算的误差估计,的相对误差 分别是 对 的相对误差 的增长因子，表示相对误差 经过传 播后增大或减小的倍数。,由此可以得到两近似数 的和、差、积、商的误 差估计（绝对误差）为,2 误差与数值计算的误差估计,例1.4 经过四舍五入得到 ， ，问他们 分别

13、具有几位有效数字？ ， ， ， 的 绝对误差限分别是多少？,2 误差与数值计算的误差估计,2 误差与数值计算的误差估计,解：记 和 的精确值分别是 和 ，则 分别具有5位有效数字,2 误差与数值计算的误差估计,例1.5 测得某电阻两端的电压和流过的电流分别 为 伏、 安，求电阻的阻值 ， 并求 及 。,2 误差与数值计算的误差估计,解：有 ，已知 伏， 安，得 欧 的绝对误差： 由于 ， ，所以 从而 的相对误差,例1.6 设 ， ， 都精确到2位小数， 估计 的相对误差。,2 误差与数值计算的误差估计,2 误差与数值计算的误差估计,解： 所以,在由误差估计式得出绝对误差限和相对误差限 的估计

14、时，由于取了绝对值并用三角不等式放大， 是按照最坏情形得出的，所以结果是保守的。 一般来说，为了保证计算结果的精确度，在计 算过程中，比结果中所要求的有效数字位数多取1位 或2位就可以了。,2 误差与数值计算的误差估计,计算机只能对有限位数进行计算，从而在运算 中产生误差是不可避免的。许多实际问题的求解往 往需要进行成千上万次的数值计算，为了保证计算 结果的可靠性，必须防止误差的产生、传播与扩大 。一个好的算法应该是计算量小、精度高，算法稳 定，在计算过程中占用计算机的存储单元和工作单 元少。,3 选用和设计算法时应遵循的原则,3 选用和设计算法时应遵循的原则,选择算法应遵循的原则： 1.算法

15、是否稳定； 2.算法的逻辑结构是否简单； 3.算法的运算次数和算法的存储量是否尽量少。,减少运算误差的几项措施： 1.选用数值稳定的计算公式，控制舍入误差的传播； 在数值计算中，对于某一问题选用不同的算法，所得到 的结果往往不同，有时甚至大不相同。这主要是由于初始数 据的误差或计算时的舍入误差在计算过程中的传播因算法的 不同而异。对某一算法，如果初始数据的误差或舍入误差对 计算结果的影响较小，则称该算法是数值稳定的；否则，称 为数值不稳定算法。,3 选用和设计算法时应遵循的原则,例1.7 计算积分 解：由 算法一 将 代入递推公式分别计算，由于,3 选用和设计算法时应遵循的原则,所以 ，因此算

16、法一误差过大，原因是 的近似值所 产生的误差每计算一次放大10倍，误差随n的增大迅速递增 算法二 由后项递推前项， 中的误差传递到 时下降为原来的 1/10。由于 令,3 选用和设计算法时应遵循的原则,代入递推公式分别计算，计算结果如表所示。,3 选用和设计算法时应遵循的原则,3 选用和设计算法时应遵循的原则,2.尽量简化计算步骤以便减少运算次数； 3.避免两个相近数相减； 4.绝对值太小的数不宜作除数； 5.合理安排运算数序，防止大数淹没小数；多个数相加时，绝对值小者先加；多个数相乘时，有效位数多者先乘。,课后题,1.要使 的相对误差不超过 ，至少需要保留多少位 有效数字？ 2.设 为 的近似数，证明 的相对误差约为 的相对误差 的 倍。 3.求 的近似有效数： （1）使绝对误差不超过0.01； （2）是相对误差不超过0.01； 4.正方形的边长约为10cm，问测量边长的误差界多大时才能 保证面积的误差不超过0.1cm2？,课后题,5.设 的相对误差为1%，求 的相对误差。 6.为使积分 的近似值的相对误差不超过1%，问至 少要取记为有效数字？ 7.设 ，假定 是准确的，而对 的测量有0.1秒的 误差，试证当 增加时， 的绝对误差增加，而相对误差减少。,