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