环境系统分析数学模型的参数估计及灵敏度分析-PPT课件.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,环 境 系 统 分 析,第,8,讲,1,第四章 数学模型的参数估计及灵敏度分析,前章所述的一些解析模型常用于环境质量的模拟预测和控制规划,一维解析模型广泛地用于各种河流的水质模拟和预测中,三维解析模型在大气质量的预测中普通采用,在流动均匀稳定的条件下，二维解析模型可用来模拟河流的水质,在模型具体应用时，必须首先对模型中的参数进行估值和进行灵敏度的分析。,2,一、模型参数的估值方法,有经验公式，图解法，最小二乘法和最优化方法等估值方法,除经验公式外，其余方法均是利用系统输入输出数据和数学模型本身确定合理的参数数值。,1,、,图解法,对经适当处理后以转换为直线的公式，均可用图解法估计参数，其误差取决于点位的精度和绘制直线的精度。,3,2,、一元线性回归分析法,亦称最小二乘法,该法有两个假定：,所有自变量的值均不存在误差，因变量的值则含有测量误差；,与各测量点拟合最好的直线为能使各点到直线的竖向偏差（因变量偏差）的平方和最小的直线。,4,偏差的平方和最小意味着各个点的偏差均很小。,最佳的,b,和,m,的估计值：（,y=mx+b,）,由,5,3,、多元线性回归分析（原理相同）,以二元为例,6,4,、最优化估值方法,函数一般式,：,建立目标函数：,使其最小（,Z min,）。,对一个连续可微的目标函数可采用最速下降法（一阶梯度法）。,7,梯度法的步骤如下：,第一步,:,设,1,2,m,的初值为,1,2,m,允许迭代误差为,.,第二步：计算目标函数的初值,第三步：计算目标函数对参数的梯度。,8,大家学习辛苦了，还是要坚持,继续保持安静,9,在函数的形式比较复杂，不易求得梯度的解析式时，可以计算其数值梯度,.,第四步：计算参数修正步长,10,二阶梯度矩阵,H,（,）,亦称海森矩阵,。,11,对于复杂的数学表达式，海森矩阵的解析值很难计算，可以数值梯度来近似的解析值。,对于海森矩阵的对角元素：,对于非,对角元素：,12,第五步：计算参数,i,的修正值,i,1,第六步，计算新的目标函数值,Z,1,第七步,，比较,Z,1,和,Z,0,若 ，,则停止运算，并输出参数的估计值,i,1,（,i=1,2n,）,否则令 返回第三步,。,13,若以相对误差表示则可取,|(Z1-Z0)/Z1|,否则计算的允许选代误差（也称截断误差）要视目标函数的绝对值大小而定。用最优化方法估值时，要由经验给定参数的初值。,14,例：已知河流沿程的溶解氧（,DO,）的测定数据如下：,若起点的,BOD,（,L,0,）为,20mg/l,，饱和溶解氧（,Cs,）为,10.0mg/l,，河流平均流速为,Ux=4.0km/h,由,S-P,模型可知河流溶解氧的变化规律符合下述方程,:,X(km),0 8 28 36 56,DO(mg/l),10.0 8.5 7.0 6.1 7.2,15,试确定其中的耗氧速度常数,Kd,和得氧速度,常数,Ka,。,解：首先，建立目标函数,16,用一阶梯度法，据前述的七步，编制计算机程序，给定初值，,K,0,d,=1.0d,-1,=0.042h,-1,，,K,0,a,=2.0d,-1,=0.083h,-1,当目标,Z=0.4681,时，得到参数的最优估计值,:,K,d,=0.053 h,-1,=1.27 d,-1,K,a,=0.19 h,-1,=4.67 d,-1,。（,取的是,0.0001,）。,17,5,、网格法,假定有,n,个等定参数，且已知各参数的取值范围，把各搜索区间（取值范围）分成若干个等分，则参数空间,=(,1,2,n,),T,就被划分成若干网格，计算所有网格顶点上的目标函数值，并取其中最小的值所对应的参数值作为最优估计值。,若精度还不够，则可再分细些。,18,6,、经验公式计算法,如：河流的复氧速度常数，大气扩散方程中的方差等。除经验公式计算法外，其余方法均应有自度量和因变量的实测输入输出数据，注意使用条件，范围。,二、模型的验证与误差分析,在模型建立且参数估值之后，还应对模型进行验证和误差分析方可投入应用。,19,验证所用的数据应与参数估值时所用数据独立，以模型的计算结果和实测数据之间的吻合程度来判断。,常用方法：,1,、图形表示法,观测值为横坐标，计算值为纵坐标，据各自变量可得上面相应的两值。,由于环境系统问题的复杂性，对于大系统，有的文献认为，对于观测值和计算值在,2,倍误差范围内都认为满意。,20,21,2,、相关系数法,统计学上衡量曲线拟合程度的量。,y,和,y,分别为观测值和计算值的平均值。,r,越大相关关系越好（,0,r,1,）。,22,当对,y=+y+,作回归分析证明,=0,和,=1,时用相关系数验证才有实际意义。,表示计算值,y,和实测值,y,之间的误差。,3,、相对误差法,e,i,=y,i,-y,i,/y,i,n,组观测值与相应计算值数据可得,n,个误差值，将这,n,个误差值从小到大排列，可以求得小于某一误差值的误差的出现频率以及累积频率为,10%,、,50%,和,90%,的误差。,23,通常采用中值误差（累积频率为,50%,）作为衡量模型精确度的度量。,中值误差与统计学上的概率误差是一致的。,中值误差可从误差分布的累积曲线上求出，也可按下式计算：,24,常用,e,0.5,的,10%,作为水质模型验证标准，还有用绝对中值误差的。,(,公式分母中,y,i,去掉,),利用相关系数、相对中值误差和绝对中值误差等验证方法还可验证所用参数估值方法哪种效果更好些。,25,三、数学模型的灵敏度分析,由于环境系统是一个开放性系统，各种影响非常复杂，很难精确定量，各种数学模型存在着不确定性（有许多假设），模型中的参数也有误差，因此，利用模型进行的模拟和规划的真实性，可靠性究竟如何，如何对此做出估计，换言之，状态变量对参数的灵敏度如何，目标函数对参数的灵敏度如何以及目标函数对状态变量的灵敏度如何，需进行分析。,26,灵敏度分析可以估计模型计算结果的偏差，且还有助于建立低灵敏度系统，（这种系统在运行上比较可靠），有助于确定合理的设计裕量，这比盲目给定安全系数要合理得多。,（希望是低灵敏度高预测精度的模型）,27,误差分析是直接验证模型计算结果与实测值的差异，针对一些零散值而作的，而灵敏度分析是从另一角度考虑该模型参数的误差大小对状态变量所引起的计算误差和对目标函数所引起的误差的一种敏感程度。,28,下面仅介绍一下状态变量和参数的数目都是,1,时的灵敏度分析。,若决策变量（污染物排放量等）保持不变，则状态变量,x,和目标,Z,均可表示为参数,的函数：,x*=f(,0,),Z*=f(,0,),x*,和,Z*,分别表示参数,取,0,值的状态变量值和目标函数值。,29,灵敏度的定义为：,在,=,0,附近，状态变量（或目标）相对于原值的变化率和参数,相对于,0,的变化率的比值称为状态变量（或目标函数）对参数的灵敏度，即：,状态变量对参数的灵敏度为：,30,目标函数对参数的灵敏度为：,式中,x=x,x*,z=z,z*,=,0,当,0,时，忽略高阶微分项得：,31,例：已知某河段的,BOD,降解规律可用下式表示：,L=L,0,e,-K,d,t,若已知河段初始的,BOD,浓度,L,0,=15mg/l,BOD,衰减速度常数,K,d,=0.1 d,-1,假定,K,d,的变化幅度在,10%,试求,t=2d,时的,BOD,值及其变化幅度。,解：,K,d0,=0.1d,-1,时,t=2d,的,BOD,浓度为：,L*=L,0,e,-K,d0,t,=15,e,-0.1,2,=12.28mg/l,32,BOD,对,K,d,的一阶灵敏度系数为：,BOD,对,k,d,的灵敏度为：,33,已知：,K,d,/,K,d0,=10%,，所以,BOD,的变化幅度为：,变化与,K,d,的变化方向相反。,因为,2%,10%,所以属低灵敏度模型。,34,1,双曲线函数,2,指数函数,y=A e,bx,式中,A,0,ln y=ln A+b x,令,Y=ln y ,a=ln A ,X=x,则,Y=a+b X,35,y=A e,b/x,式中,A,0,ln y=ln A+b/x,令,Y=ln y ,a=ln A ,X=1/x,则,Y=a+b X,3,对数函数,y=a+b lnx,令,Y=y ,X=ln x,则,Y=a+b X,36,4.,幂函数,y=AX,b,(A,0),lny=lnA+blnx,令,Y=ln y ,a=ln A ,X=ln x,则,Y=a+b X,5.S,曲线,y=1/(a+be,-x,),1/y=a+be,-x,令,Y=1/y ,X=e,-x,则,Y=a+b X,37,值得注意的是，广义线性函数的剩余平方和、剩余标准差和相关系数应以原,y,模式求。,还有一些,广义线性函数，不再列举。可按此思路处理。,38,