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

计算声学第一章数值计算中的误差分析第1页，共60页，2023年，2月20日，星期四水声学主要研究声波在水下的辐射、传播与接收，用以解决与水下目标探测和信息传输过程有关的各种声学问题。声波是目前在海洋中唯一能够远距离传播的能量辐射形式。因此作为信息载体的声波，在海洋中所形成的声场时空结构，就成为近代水声学的基本研究内容，而提取海洋中声场信息的结构是我们用来进行水下探测、识别、通信及环境监测等的手段。前言第2页，共60页，2023年，2月20日，星期四前言海洋环境第3页，共60页，2023年，2月20日，星期四波动方程：波动方程是声学量在声场中满足的基本关系式，反映了波动特征，也是进行声场计算的基本关系式。在导出波动方程前，为了使问题简化，需要对介质和声波做一些假设：（1）介质是均匀连续的，即在波长数量级距离内，介质的声学性质保持不变；（2）介质是理想流体介质，声波在其中传播时没有能量损耗，即忽略介质的粘滞性和热传导性；（3）研究小振幅波的传播规律，所谓小振幅波是指各声学量都是一级小量。波动方程是描述波动运动的数学表达式，它由连续性方程、状态方程和运动方程推导得到。前言第4页，共60页，2023年，2月20日，星期四波动方程：理想流体介质中小振幅平面波的波动方程为（沿轴向传播）：小振幅声压在三维坐标下的波动方程为为拉普拉斯算符，在直角坐标系中

前言第5页，共60页，2023年，2月20日，星期四海洋声场的数值预报在建立了能够反映海洋环境因素对声场的制约关系的声场物理模型（波动方程+定解条件）的基础上，根据可测海洋环境参数的测定值或预报值，编写程序完成数值计算，给出相应海洋环境条件下的有关场值。近年来，由于计算机的快速发展，数值计算声场是一个快速发展的领域。海洋声场的数值预报方法主要有射线算法、简正波算法、抛物方程（PE）算法、快速场（FFP）算法等，各自有不同的适应范围。前言第6页，共60页，2023年，2月20日，星期四前言函数插值：已知一组不同深度处的声速值，如何得到任意深度处的声速值？

深度（m）0.050.0100.0200.0300.0400.0500.0800.0声速(m/s)1510.51510.41505.81500.81496.01492.01488.11483.21000.01200.02000.03000.04000.01482.61482.41498.01516.61534.8第7页，共60页，2023年，2月20日，星期四前言第8页，共60页，2023年，2月20日，星期四前言数值积分：声线轨迹计算声线从深度传播到深度所经过的水平距离为第9页，共60页，2023年，2月20日，星期四前言问题：利用射线声学模型对海洋声场进行求解

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

利用四阶龙格库塔算法求解前言第16页，共60页，2023年，2月20日，星期四前言四阶龙格库塔算法第17页，共60页，2023年，2月20日，星期四前言第18页，共60页，2023年，2月20日，星期四前言第19页，共60页，2023年，2月20日，星期四前言第20页，共60页，2023年，2月20日，星期四前言必要性：

现代科学研究和高技术的发展越来越需要借助计算机进行数值计算，水声领域也不例外。讲授的主要内容：1、数值计算方法：误差分析、方程组求解、非线性方程求解、插值法、最小二乘与曲线拟合、数值微积分、常微分方程求解；2、进化算法（方程组求解）：量子粒子群算法；3、虚源法声场建模。第21页，共60页，2023年，2月20日，星期四数值计算的对象、任务与特点对象：数值计算方法是研究科学与工程技术中数学问题的数值解及其理论的一个分支，涉及代数、微积分、微分方程等的数值解问题。任务：研究适合在计算机上使用的数值计算方法及相关理论，如方法的收敛性、稳定性和误差分析等；还要根据计算机的特点研究如何设计计算方法做到计算时间短、占用内存小。学习目的：提高应用计算机解决实际问题的能力。

第22页，共60页，2023年，2月20日，星期四数值计算的对象、任务与特点数值计算流程：特点：

既具有数学的抽象性与严格性，又具有应用的广泛性与实际实验的技术性，是一门与计算机紧密结合的实用性很强的有着自身研究方法与理论体系的计算数学课程。第23页，共60页，2023年，2月20日，星期四数值计算中的误差分析内容提要：掌握绝对误差、相对误差、有效数字、数值计算的误差估计以及设计算法的原则。重点内容：绝对误差、相对误差、有效数字的概念，数值计算的误差估计。第24页，共60页，2023年，2月20日，星期四误差与数值计算的误差估计一、误差的来源与分类将一个数的准确值与其近似值之差称为误差。1.分类

过失误差：人为造成，可以避免

非过失误差：

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

设某一量的精确值为，其近似值为，则称为近似值的绝对误差，简称误差。 时称为弱近似值或亏近似值； 时称为强近似值或盈近似值。绝对误差限：如果存在，使得，则称为近似值的绝对误差限，简称误差限或精度（测量时，测量工具最小刻度的一半）。越小，表示近似值的精度越高。在工程技术上常用表示近似值的精度或精确值的范围。§2误差与数值计算的误差估计第28页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计例：用毫米刻度的尺子测得桌子长度近似值为mm，由尺子的精度可以知道，近似值的误差不超过0.5mm，即

表明精确值在区间内，可以写成绝对误差限mm，即绝对误差限是末位的半个单位。第29页，共60页，2023年，2月20日，星期四相对误差和相对误差限相对误差：绝对误差与精确值之比，即称为近似值的相对误差。实际中，由于精确值一般无法知道，所以常取作为近似值的相对误差。相对误差限：若存在，使得，则称为近似值的相对误差限。注意：绝对误差和绝对误差限与有相同的量纲，相对误差和相对误差限是无量纲的，工程中常以百分数来表示。§2误差与数值计算的误差估计第30页，共60页，2023年，2月20日，星期四例1.1

国际大地测量学会建议光速采用其含义是绝对误差限为多少？而其相对误差限为多少？§2误差与数值计算的误差估计第31页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计绝对误差限：近似值：相对误差限：第32页，共60页，2023年，2月20日，星期四有效数字如果近似值的绝对误差限是某一位的半个单位，就称其“准确”到这一位，且从该位开始直到的第一位非零数字共有n位，则称近似数有n位有效数字。有效数字既能表示近似值的大小，又能表示其精确程度（绝对误差限）。例1.2

设，其近似值，问有几位有效数字？如果，有几位有效数字？§2误差与数值计算的误差估计第33页，共60页，2023年，2月20日，星期四练习题1.指出如下有效数的绝对误差限、相对误差限和有效数字位数。

2.将22/7作为的近似值，它有几位有效数字？绝对误差限和相对误差限各为多少？第34页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计有效数字

设的近似值可以写成如下的标准形式

所以当其绝对误差限满足时，称近似值具有位有效数字。第35页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计结论：如果

，有位有效数字，则其相对误差限为反之，如果的相对误差限满足则至少有位有效数字。第36页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计例1.3

要使的近似值的相对误差小于，至少需取几位有效数字？第37页，共60页，2023年，2月20日，星期四误差的传播与估计实际的数值计算中，参与运算的数据往往都是近似值，带有误差。而在进一步运算中都会产生舍入误差或截断误差，这些误差在运算过程中会进行传播，影响计算结果。§2误差与数值计算的误差估计第38页，共60页，2023年，2月20日，星期四一元函数的泰勒（Taylor）中值定理：如果函数在区间内有直到阶导数，，则有其中，拉格朗日型余项（介于之间）。

§2误差与数值计算的误差估计第39页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计泰勒（Taylor）公式估计误差的方法：以二元函数为例，设和分别是和的近似值，是函数值的近似值。函数在点处的Taylor展开式为第40页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计式中一般都是小量值，如果忽略它们的高阶无穷小量，则上式简化为因此的绝对误差为第41页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计系数分别是一阶偏导数在处的值，称为对的绝对误差的增长因子，分别表示绝对误差经过传播后增大或缩小的倍数。第42页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计的相对误差

分别是对的相对误差的增长因子，表示相对误差经过传播后增大或减小的倍数。

第43页，共60页，2023年，2月20日，星期四

由此可以得到两近似数的和、差、积、商的误差估计（绝对误差）为§2误差与数值计算的误差估计第44页，共60页，2023年，2月20日，星期四例1.4

经过四舍五入得到，，问他们分别具有几位有效数字？，，，的绝对误差限分别是多少？§2误差与数值计算的误差估计第45页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计解：记和的精确值分别是和，则分别具有5位有效数字第46页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计例1.5

测得某电阻两端的电压和流过的电流分别为伏、安，求电阻的阻值，并求及。

第47页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计解：有，已知伏，安，得欧的绝对误差：由于，，所以从而的相对误差第48页，共60页，2023年，2月20日，星期四例1.6设，，都精确到2位小数，估计的相对误差。§2误差与数值计算的误差估计第49页，共60页，2023年，2月20日，星期四§2误差与数值计算的误差估计解：所以第50页，共60页，2023年，2月20日，星期四在由误差估计式得出绝对误差限和相对误差限的估计时，由于取了绝对值并用三角不等式放大，是按照最坏情形得出的，所以结果是保守的。 一般来说，为了保证计算结果的精确度，在计算过程中，比结果中所要求的有效数字位数多取1位或2位就可以了。§2误差与数值计算的误差估计第51页，共60页，2023年，2月20日，星期四计算机只能对有限位数进行计算，从而在运算中产生误差是不可避免的。许多实际问题的求解往往需要进行成千上万次的数值计算，为了保证计算结果的可靠性，必须防止误差的产生、传播与扩大。一个好的算法应该是计算量小、精度高，算法稳定，在计算过程中占用计算机的存储单元和工作单元少。§3选用和设计算法时应遵循的原则第52页，共60页，2023年，2月20日，星期四§3选用和设计算法时应遵循的原则选择算法应遵循的原则：1.算法是否稳定；2.算法的逻辑结构是否简单；3.算法的运算次数和算法的存储量是否尽量少。第53页，共60页，2023年，2月20日，星期四减少运算误差的几项措施：1.选用数值稳定的计算公式，控制舍入误差的传播；在数值计算中，对于某一问题选用不同的算法，所得到的结果往往不同，有时甚至大不相同。这主要是由于初始数据的误差或计算时的舍入误差在计算过程中的传播因算法的不同而异。对某一算法，如果初始数据的误差或舍入误差对计算结果的影响较小，则称该算法是数值稳定的；否则，称为数值不稳定算法。

§3选用和设计算法时应遵循的原则第54页，共60页，2023年，2月20日，星期四例1.7

计算积分解：由算法一

将代入递推公式分别计算，由于§3选用和设计算法时应遵循的原则第55页，共60页，2023年，2月20日，星期四所以，因此算法一误差过大，原因是的近似值所产生的误差每计算一次放大10倍，误差随n的增大迅速递增算法二

由后项递推前项，中的误差传递到时下降为原来的1/1