05灵敏度分析

《05灵敏度分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、环境系统分析第四章数学模型的参数估计及灵敏度分析前章所述的一些解析模型常用于环境质量的模拟预测和控制规划一维解析模型广泛地用于各种河流的水质模拟和预测中三维解析模型在大气质量的预测中普通采用在流动均匀稳定的条件下，二维解析模型可用来模拟河流的水质在模型具体应用时，必须首先对模型中的参数进行估值和进行灵敏度的分析。一、模型参数的估值方法有经验公式，图解法，最小二乘法和最优化方法等估值方法除经验公式外，其余方法均是利用系统输入输出数据和数学模型本身确定合理的参数数值。1、图解法对经适当处理后以转换为直线的公式，均可用图解法估计参数，其误差取

2、决于点位的精度和绘制直线的精度。2、一元线性回归分析法亦称最小二乘法该法有两个假定：①所有自变量的值均不存在误差，因变量的值则含有测量误差；②与各测量点拟合最好的直线为能使各点到直线的竖向偏差（因变量偏差）的平方和最小的直线。偏差的平方和最小意味着各个点的偏差均很小。最佳的b和m的估计值：（y=mx+b）由3、多元线性回归分析（原理相同）以二元为例4、最优化估值方法函数一般式：建立目标函数：使其最小（Zmin）。对一个连续可微的目标函数可采用最速下降法（一阶梯度法）。梯度法的步骤如下：第一步:设θ1,θ2,…,θm的初值为θ1°,θ2°

3、,…θm°,允许迭代误差为ℇ.第二步：计算目标函数的初值第三步：计算目标函数对参数的梯度。在函数的形式比较复杂，不易求得梯度的解析式时，可以计算其数值梯度.第四步：计算参数修正步长λ二阶梯度矩阵H（θ°）亦称海森矩阵。对于复杂的数学表达式，海森矩阵的解析值很难计算，可以数值梯度来近似的解析值。对于海森矩阵的对角元素：对于非对角元素：第五步：计算参数θi的修正值θi1第六步，计算新的目标函数值Z1第七步，比较Z1和Z0若，则停止运算，并输出参数的估计值θi1（i=1,2……n）否则令返回第三步。若以相对误差表示则可取

4、(Z1-Z0)/Z1

5、

6、≤ε否则计算的允许选代误差（也称截断误差）要视目标函数的绝对值大小而定。用最优化方法估值时，要由经验给定参数的初值。例：已知河流沿程的溶解氧（DO）的测定数据如下：若起点的BOD（L0）为20mg/l，饱和溶解氧（Cs）为10.0mg/l，河流平均流速为Ux=4.0km/h,由S-P模型可知河流溶解氧的变化规律符合下述方程:X(km)08283656DO(mg/l)10.08.57.06.17.2试确定其中的耗氧速度常数Kd和得氧速度常数Ka。解：首先，建立目标函数用一阶梯度法，据前述的七步，编制计算机程序，给定初值，K0d=1.0d

7、-1=0.042h-1，K0a=2.0d-1=0.083h-1当目标Z=0.4681时，得到参数的最优估计值:Kd=0.053h-1=1.27d-1Ka=0.19h-1=4.67d-1。（ℇ取的是0.0001）。5、网格法假定有n个等定参数，且已知各参数的取值范围，把各搜索区间（取值范围）分成若干个等分，则参数空间θ=(θ1,θ2,…,θn)T就被划分成若干网格，计算所有网格顶点上的目标函数值，并取其中最小的值所对应的参数值作为最优估计值。若精度还不够，则可再分细些。6、经验公式计算法如：河流的复氧速度常数，大气扩散方程中的方差等。除经

8、验公式计算法外，其余方法均应有自度量和因变量的实测输入输出数据，注意使用条件，范围。二、模型的验证与误差分析在模型建立且参数估值之后，还应对模型进行验证和误差分析方可投入应用。验证所用的数据应与参数估值时所用数据独立，以模型的计算结果和实测数据之间的吻合程度来判断。常用方法：1、图形表示法观测值为横坐标，计算值为纵坐标，据各自变量可得上面相应的两值。由于环境系统问题的复杂性，对于大系统，有的文献认为，对于观测值和计算值在2倍误差范围内都认为满意。2、相关系数法统计学上衡量曲线拟合程度的量。y和y'分别为观测值和计算值的平均值。r越大相关

9、关系越好（0≤r≤1）。当对y=α+βy‘+ε作回归分析证明α=0和β=1时用相关系数验证才有实际意义。ε表示计算值y和实测值y’之间的误差。3、相对误差法ei=∣yi-yi∣/yin组观测值与相应计算值数据可得n个误差值，将这n个误差值从小到大排列，可以求得小于某一误差值的误差的出现频率以及累积频率为10%、50%和90%的误差。通常采用中值误差（累积频率为50%）作为衡量模型精确度的度量。中值误差与统计学上的概率误差是一致的。中值误差可从误差分布的累积曲线上求出，也可按下式计算：常用e0..5的10%作为水质模型验证标准，还有用绝对

10、中值误差的。(公式分母中yi去掉)利用相关系数、相对中值误差和绝对中值误差等验证方法还可验证所用参数估值方法哪种效果更好些。三、数学模型的灵敏度分析由于环境系统是一个开放性系统，各种影响非常复杂，很难精确定