浙江工业大学第十四届大学生数学建模竞赛
所选赛题： ??A □B □C
我们承诺：
?? 我们仔细阅读了浙江工业大学数学建模竞赛的竞赛规则。我们完全明白，在竞赛
开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与赛题有关的问题。
?? 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的
资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
?? 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛
规则的行为，我们将受到严肃处理。
参赛队员（签名）：
学院年级： 学院年级： 学院年级： 联系方式： 联系方式： 联系方式：
参赛队员（签名）： 参赛队员（签名）：
??我们队伍愿意参加暑期数学建模培训，参加全国大学生数学建模竞赛
题目： 雾霾相关问题探究
摘 要
雾霾天气会造成能见度降低，空气湿度较大，气溶胶污染等环境灾害问题。鉴于雾霾天气已成为一种新的气象及环境灾害性现象，本文将对雾霾引发的一些问题进行探讨，主要包括两个方面，一是建立数学模型描述雾霾通过对植物的直接影响而对生物链造成的危害；二是建立一种算法来处理雾霾天气下摄取的模糊图片。
基于三种群 Volterra 模型，我们通过引入雾霾影响系数，定义了雾霾影响力值，建立了能描述在雾霾环境下，由植物、哺乳动物和爬行动物组成的生物链的常微分方程组。我们主要通过分析常微分方程组的稳定点，以一个在没有雾霾的影响下，三个种群都没有灭绝的情形的数值例子作为研究的参照点，通过增加此情形下的雾霾影响力值，来分析在不同雾霾影响力下生物链的变化。最终，我们能通过数据模拟的方法，得到雾霾影响力值在范围 0 至 1 之间变化时，种群数量的具体变化情况。此外，我们根据雾霾影响力值与种群数量变化之间的关系，建立了评判雾霾影响生物链程度的指标。
本文的三种群 Volterra－雾霾影响模型能在取得实际种群的相关数据的情况下，通过计算得到雾霾影响力值，并以此根据评判指标进行危害性衡量。更进一步地，本模型也适用于评价其他对植物有直接影响的灾害性因素，可见本模型有一定的推广性，实用性以及参考性。
在处理雾霾天气下的降质图片时，我们首先分析了雾霾的成因。通过研究雾霾造成的光线散射的物理过程，提出了霾引起光的衰减模型的 RGB 补偿处理，即在图片中红绿蓝的颜色通道上，增加计算得到的红绿蓝光在霾中的光强衰减值。对于雾霾天气下能见度降低导致图像对比度下降的问题，我们用简单而有效的直方图均衡化方法来处理。此外，在上述方法处理得到的图片上，我们用图像引导滤波方法对细节部分进行加强，以此来处理由雾霾造成的成像轮廓模糊的问题，使图片更加符合人的视觉需求。
本文的雾霾图片处理算法能够较好地提升原图的清晰度，色彩饱和度以及对比度，其去霾效果比较突出，整体效果良好，是一种易于实践，可借鉴且有效的算法。
关键词：三种群 Volterra 模型 常微分方程组 霾引起光的衰减模型
直方图均衡化 图像引导滤波
1
Ⅰ 问题重述
雾天由于湿度大，能见度差，容易造成交通阻塞和引发交通事故；而霾天的细粒子会造成气溶胶污染，粒子被人体吸收将刺激支气管，加重哮喘、鼻炎等呼吸系统病症，可见雾霾天气对人类活动有很大的危害性。而本文主要讨论两类问题，第一类是由于雾霾造成的植物光合作用减弱，间接影响生物链的变化问题；第二类是由于雾霾天气能见度降低，湿度较大造成拍摄的图像降质的问题。所以我们将完成以下 2 项内容： 1、建立数学模型描述雾霾天气对于生物链的影响。
2、建立并评价对雾霾天气下拍摄的图片进行处理的算法。
Ⅱ模型的假设
1、 三种群能够合理地反映生物链的情况。 2、 雾霾的影响时间是持续长久的。 3、 雾霾在大气中是均匀分布的。
Ⅲ雾霾对生物链的影响
3.1 模型的建立
?? 三种群 Volterra－雾霾影响模型
在本问题中，我们建立的雾霾对生物链的影响模型是基于三种群 Volterra 模型[1][2] 的。
我们知道在一个生态系统中，存在哺乳动物吃植物，爬行动物吃哺乳动物的自然现像。当植物、哺乳动物、爬行动物在一个自然环境中生存时, 我们把植物、哺乳动物、爬行动物的数量分别记作 x1(t), x2(t), x3(t)。
对于植物，若不考虑自然资源对植物的限制, 植物独立生存时以指数规律增长, 相 对增长率为 r1 , 即 x?1(t)=r1x1。
考虑到自然资源的有限性，植物，哺乳动物，爬行动物的数量应均遵从 logistic 规律。当植物独立存在时，设 N1，N2，N3 分别为植物，哺乳动物，爬行动物在环境资源
x 容许的最大数量。我们有 x?1(t)=r1x1(1? 1 )。 N 1
而哺乳动物的存在使植物的增长率减小，设减小的程度与哺乳动物的数量和能力成正比，所以我们设单位数量的哺乳动物（相对于 N2）掠取 σ12 倍的单位植物量(相对于 N1)。我们有
x?1(t)=r1x1(1? x1 x2
?? )
12
N1 N 2
在本问题中，植物的生长受到了雾霾的影响。霾会使光线粒子发生散射，雾会使光
线粒子发生折射，致使植物受到的光强减弱，导致其光合作用减弱，其增长率减小（见图 1）。
2
(a)无雾霾时植物接收阳光的情形 (b)有雾霾时植物接收阳光的情形
图 1 植物接收阳光的不同情形
我们在这里引入一个雾霾影响系数 θ，用以表示雾霾对植物增长的抑制作用。我们 有
进一步化简可得：
x?1(t)=r1x1(1x1 ?12 x2 )?θx1
N? N
?
1 2
? x
x? (t)=r1x1(1? ? 1 1 ???x2 )r1 N1
12
N 2
我们定义 k 为雾霾影响力，令 k=θ/r1 ，表示雾霾相对于植物固有增长率的抑制作用。于是我们有
x?1(t)=r1x1(1? x1 ?x2 ?k)
12 ? N1 N 2
(3.1-1)
由于植物不可能会有负增长，所以 k∈[0,1]。
对于哺乳动物，哺乳动物离开植物无法生存，设其独立生存时死亡率为 r2 ，即
?2 (t)= x?r2x2 。
而植物的存在又为哺乳动物提供了食物, 植物的存在相当于使哺乳动物的死亡率降低, 且促使哺乳动物增长, 设单位数量的植物(相对于 N1)供养 σ21 倍的单位哺乳动物量 (相对于 N2)。我们有
?2 (t)= r2x2(?1+??21 x 1 ) xN1
而哺乳动物又为爬行动物提供了食物, 爬行动物的存在使哺乳动物的增长率减小, 设单位数量的爬行动物（相对于 N3）掠取 σ23 倍的单位植物量(相对于 N2) 。我们有
?2 (t)= r2x2(?1+? 21 x1 ???23 x3 ) xN N 3 1 而哺乳动物的增长也应符合 logistic 规律，我们有
3
?2 (t)= r2x2(?1+?21 xx1 x2? ?? 23 x3 ) N 3 N 1 N 2
(3.1-2)
对于爬行动物，爬行动物离开哺乳动物无法生存，设其独立生存时死亡率为 r3 ， 即 x?3 (t) = -r3x3。
而哺乳动物的存在又为爬行动物提供了食物, 哺乳动物的存在相当于使爬行动物的 死亡率降低, 且促使爬行动物增长, 设单位数量的哺乳动物(相对于 N2)供养? 32 倍的单
x2
)。
位爬行动物量(相对于 N3)。我们有 (t)= r3x3(?1+??32 Nx?3
2
而爬行动物的增长也应符合 logistic 规律，我们有
?3 (t) = r3x3(?1+? 32 x2 ? x3 ) xN 2 N 3
(3.1-3)
方程(3.1-1)，(3.1-2)，(3.1-3)构成了在雾霾环境下，植物、哺乳动物、爬行动物三者依存、制约现象的数学模型, 即
??
? x? (t) ? r x (1? x1 ? ? x2 ? k)
1 1 1 ? 12 2 1 ??N N??xx(3.1-4) ? x?2 (t) ? r2 x2 (?1? ? 21 1 ? 2 ? ? 23 x3 ) ??N1 N2 N3 ??x x ? x? (t) ? r x (?1? 2 ? 3 ) ??
3 3 3
???32
N 2 N 3
?? 平衡点及其坐标
虽然非线性微分方程组(3.1-4)无解析解，但是我们可以解出该微分方程组的平衡点及其坐标，结果如下：
P1(0，0，0) P2(0，0，?N3)
P3((1?k)N3，0，0) P4(0，?N2，0)
? N3 (1 ? ? 32 ) ) P5(0， N2 (?1 ? ? 23 ) ， 1 ? ? 23? 32 1 ? ? 23? 32
P6((1-k)N1，0，?N3)
N2[(1 ? k)? 21 ?1]，0) P7( N1(1 ? k ? ?12 ) ， 1 ? ?12? 21 1 ? ?12? 21
N [?1 ? ? 23 ? (1 ? k )? 21 ] P8( N1[1 ? k ? ?12 ? ?12? 23 ? (1 ? k)? 23? 32 ] ， 2，
1 ? ? 12? 21 ? ? 23? 32 1? ?12? 21 ? ? 23? 32
4
雾霾问题论文(终稿)
