污染物扩散模型

《污染物扩散模型》由会员分享，可在线阅读，更多相关《污染物扩散模型（6页珍藏版）》请在装配图网上搜索。

1、污染物扩散模型污染物扩散模型 一、问题分析 题目要求利用马氏链模型来解决该问题.由题目条件知,要让各城市污染物浓度在无论时间有多大都要小于某一个特定值,可将各城市下一刻点污染物浓度与目前的污染物浓度表示出来,得到一个关于污染物浓度变化的递推公式,对该公式进行利用递推法可得到污染物浓度的表达式,令其小于题目中给出的特定即可实现对问题的求解. 二、模型假设 1各城市污染物浓度仅与浓度扩散的转移概率有关. 2扩散到给出城市之外的污染物不会再回来. 三、符号约定 vi(i=1,2,L,k) pij(i,j=1,2,L,k) 不同的城市 污染物从vi扩散到vj的概率 ci(t)(i=1,2,L,k) d

2、i(i=1,2,L,k) t时间点城市vi的污染物浓度 城市vi的污染源排出的污染物数量 各城市污染物浓度最大限度 ci(i=1,2,3) *四、模型建立与求解 根据题目条件可知,各城市下一刻的污染物浓度是在目前污染物浓度在各个城市之间转移后的浓度再加上这一时刻该城市污染源排出的污染物量,即 c(t+1)=c(t)Q+d 其中c(t)=(c1(t),c2(t)Lck(t)为由各地区污染物浓度组成的k维向量,d=(d1,d2,Ldk)为由排除污染物组成的k维向量. 下面对式进行递推： 由式可得到 c(t)=c(t-1)Q+d c(t-1)=c(t-2)Q+d 将式带入到式中有 1 c(t)=c(

3、t-2)Q+dQ+d 同理可得 2c(t)=c(t-3)Q+dQ+dQ+d 32M 依次类推,可得个城市污染物浓度的表达式为 t-1 c(t)=c(0)Q+dQ s=0ts将这k个城市以及城市中的污染物看做一个系统,如果k个城市的污染物浓度视为该系统的k个状态,并增加一个状态0表示污染物扩散到k个城市之外将不再回来,污染物扩散的无后效性表明可用马氏链模型描述其变化过程,那么污染物在k+1个状态间的转移矩阵可表示为 10P= RQ其中第一行对应状态0,由污染物一旦离开这k个城市将不会再回来可知状态0是一个吸收状态,现假设各地区均对应于非吸收状态,并且由这些状态出发最终可到达0状态,从而形成一个吸

4、收链,由于(I-Q)可逆,并且有 (I-Q)t-1=s=0Q s因此可得到当时间t时,有Q0.这样在式中令t可得 c()=d(I-Q) 题目中给出当时间t充分大时必有 -1ci(t)ci 式可以表示为 ci()ci 结合式与式有 d由 *(I-Q)-1*c 131Q=30可以得出 01323131 313 2 (I*-Q)-13=3336634.5 6将上式以及c=(25,25,25)代入到式中即 3 d33式可表示为 36634.5(25,25,25) 63d1+3d2+3d3253d1+6d2+6d325 3d+4.5d+6d25231由上面的不等式组可以看出: 对于d1,d2,d30,只要3d1+6d2+6d325就可满足题目要求. 综上知,当污染物排出量d1,d2,d3满足3d1+6d2+6d325时可以时整个系统内的污染物浓度控制在给定范围之内. 3