数学建模预测模型
数学建模讲座--预测模型
年份
1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973
时序 （ t） 12 13 14 15 16 17 18 19 20 21 22
总额 （ yt ） 604.5 638.2 670.3 732.8 770.5 737.3 801.5 858.0 929.2 1023.3 1106.7
k
(一) 直线趋势外推法
适用条件：时间序列数据（观察值）呈直线 上升或下降的情形。 该预测变量的长期趋势可以用关于时间 的直线描述，通过该直线趋势的向外延伸 （外推），估计其预测值。 两种处理方式：拟合直线方程与加权拟合直线 方程
例 3.1 某家用电器厂 1993～2003 年利润额数据资料如表 3.1 所示。试预测 2004、2005年该企业的利润。
二 、趋势外推法经常选用的数学模型
根据预测变量变动趋势是否为线性，又分为线性趋势外推法 和曲线趋势外推法。
ˆt b0 b （一）线性模型y 1t （二）曲线模型 1.多项式曲线模型 2.简单指数曲线模型 3.修正指数曲线模型 4.生长曲线模型 （龚珀资曲线模型）
2
ˆt b0 b1t b2t bk t y 多项式模型一般形式：
预测模型简介
数学模型按功能大致分三种： 评价、优化、预测 最近几年，在大学生数学建模竞赛常常出 现预测模型或是与预测有关的题目：
1.疾病的传播； 2.雨量的预报； 3.人口的预测。
统计预测的概念和作用
（一）统计预测的概念
概念: 预测就是根据过去和现在估计未来，预测未来。 统计预测属于预测方法研究范畴，即如何利用科学的统计 方法对事物的未来发展进行定量推测.
数学建模地震预测模型
数学建模竞赛论文题目：地震预测数学建模：志鹏学号：12291233 学院：电气工程学院：鑫学号：10291033 学院：电气工程学院：书铭学号：12291232 学院：电气工程学院目录摘要 (3)一、问题重述 (4)二、问题的分析 (4)三、建模过程 (5)问题1：地震时间预测 (5)1、问题假设 (5)2、参数定义 (6)3、求解 (6)问题2：地震地点预测 (7)1、问题假设：72、参数定义83、求解过程：8四、模型的评价与改进 (12)参考文献 (13)摘要大地振动是地震最直观、最普遍的表现。
在海底或滨海地区发生的强烈地震，能引起巨大的波浪，称为海啸。
在大陆地区发生的强烈地震，会引发滑坡、崩塌、地裂缝等次生灾害。
对人们的生产生活成巨大影响，严重威胁人们的生命和财产安全，所以，对地震的预测是十分必要的。
本文根据从1900年以来中国发生的八级以上地震的时间和地点分析，利用合理的数学建模方法，对下一次中国可能发生的八级以上地震的和时间和地点进行合理的预测。
建模方法分为对于时间的预测和地点的预测两个方面。
问题1：对于时间的预测采用的方法为指数平滑法，它是通过计算指数平滑值，配合一定的时间序列预测模型对现象的未来进行预测。
其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。
问题2：对于地点的预测根据长久的数据表明，八级以上地震主要发生在东经70°——110°，北纬20°——50°这个围，据此将整个地震带划分为100个区域，按顺序进行编号。
建立时间与地震区域编号的数学模型，利用线性回归的方法对下次地震地点预测。
关键词：地震，预测，数学建模，指数平滑法，线性回归一、问题重述地震预报问题，大地震的破坏性是众所周知的，为了减少大地震带来的灾难，人们提出了各种预报地震的方法，以求减少大地震产生的破坏。
本赛题请大家用数学建模的方式预报下一次大地震发生的时间和地点。
数学建模预测类模型
数学建模预测类模型
数学建模预测类模型是一种采用数学模型的预测技术，可以通过对现有的数据进行建模、分析和推断，来预测未来的潜在发展情况。
这种技术的主要应用领域有金融、经济、市场预测、风险评估、贸易分析和营销策略等。
数学建模预测类模型的基本思想是利用现有的数据建立一个模型，来表示一种潜在的发展情况，以便预测未来发展趋势。
这种技术比较灵活，可用于很多不同的问题。
首先，我们需要确定所要求的模型建模形式。
有些模型可以采用统计方法，例如回归分析、相关分析等，有些模型可以采用数学优化方法，如对偶理论、回归分析、支配集理论等，还有些模型可以采用模糊逻辑、分类规则、神经网络等，具体选择哪种形式要根据实际情况而定。
其次，建立数学建模预测类模型是需要考虑的一个重要因素是数据的质量。
通常来说，模型的准确性取决于所使用的数据是否足够准确和完整。
因此，在建立模型之前，必须首先确定数据的收集与准备方式，确保数据的准确性和完整性。
最后，模型建立完毕后，还需要进行模型检验，以评估模型的预测准确性和可靠性。
一般来说，可以采取两种方法来进行模型检验：使用已有数据进行模型训练和验证，或使用未来的数据来测试模型的准确性。
如果模型具有较高的准确性，则可以放心使用该模型进行未来的预测。
综上所述，数学建模预测类模型是一种采用数学模型分析已有数据，预测未来发展趋势的技术。
为了有效地使用现有数据建立准确的数学模型，在建立模型之前，需要确定模型的建模形式，同时保证数据的准确性和完整性；建立完毕模型后，还要进行模型检验，以评估模型的准确性和可靠性。
希望通过本文的介绍，能够帮助读者更好地理解数学建模预测类模型，并熟练地应用它们。
数学建模中的预测方法：时间序列分析模型
自相关函数
k 满足 ( B) k 0
它们呈指数或者正弦波衰减，具有拖尾性
3）ARMA（ p, q）序列的自相关与偏自相关函数均是拖尾的
（2）模型的识别
自相关函数与偏自相关函数是识别ARMA模型的最主 要工具，B-J方法主要利用相关分析法确定模型的阶数. 若样本自协方差函数 k 在 q 步截尾，则 X t 是MA（ q ）序列
注：实参数 1 ,2 ,
,q 为移动平均系数，是待估参数
引入滞后算子，并令 (B) 1 1B 2 B2 q Bq 则模型【3】可简写为
X t ( B)ut
【4】
注1：移动平均过程无条件平稳 注2：滞后多项式的根都在单位圆外时，AR过程与MA过程 能相互表出，即过程可逆，
2
N 为样本大小，则定义AIC准则函数
用AIC准则定阶是指在
得 AIC (S )
p, q
最小的点
ˆ,q ˆ) (p
作为
( p, q)
的估计。
2p N 2( p q ) 2 ( p , q ) ˆ ARMA 模型 ： AIC ln N
AR（ p ）模型 ：
ˆ2 AIC ln
应用案例：
（1）CUMCM2004-A:奥运临时超市网点设计；
（2）CUMCM2004-B:电力市场的输电阻塞管理；
（3）CUMCM2005-A:长江水质的评价与预测；
（4）CUMCM2006-B:艾滋病疗法的评价与预测； （5）CUMCM2008-B:高校学费标准探讨问题。
3.灰预测GM(1,1)：小样本的未来预测 应用案例
k 在
2） kk 的截尾性判断 作如下假设检验：M N
H0 : pk , pk 0, k 1, , M H1 : 存在某个 k ，使kk
数学建模——预测模型简介
数学建模——预测模型简介在数学建模中，常常会涉及⼀些预测类问题。
预测⽅法种类繁多，从经典的单耗法、弹性系数法、统计分析法，到现在的灰⾊预测法、专家系统法和模糊数学法、甚⾄刚刚兴起的神经元⽹络法、优选组合法和⼩波分析法等200余种算法。
下⾯将简要介绍⼏类预测⽅法：微分⽅程模型、灰⾊预测模型、差分⽅程预测、马尔可夫预测、插值与拟合、神经元⽹络。
⼀、下⾯是这⼏种类型的使⽤场景对⽐：模型⽅法适⽤场景优点缺点微分⽅程模型因果预测模型，⼤多为物理、⼏何⽅⾯的典型问题，其基本规律随着时间的增长呈指数增长，根据变量个数确定微分⽅程模型。
适⽤于短、中、长期的预测，既能反映内部规律以及事物的内在关系，也嫩能够分析两个因素之间的相关关系，精度⾼便与改进。
由于反映的内部规律，⽅程建⽴与局部规律的独⽴性为假定基础，长期预测的偏差性较⼤。
灰⾊预测模型该模型不是使⽤原始数据，⽽是通过求累加、累减、均值中的两种或者全部⽅法⽣成的序列进⾏建模的⽅法。
不需要⼤量数据，⼀般四个数据即可，能够解决历史数据少、序列完整性及可靠性低的问题。
只适⽤于指数增长的中短期预测。
差分⽅程预测常根据统计数据选⽤最⼩⼆乘法拟合出差分⽅程的系数，其稳定性依赖于代数⽅程的求根。
差分⽅程代替微分⽅程描述，在⽅程中避免了导函数，可以⽤迭代的⽅式求解。
精度较低（⽤割线代替切线。
）马尔可夫预测某⼀系统在已知情况下，系统未来时刻的情况只与现在时刻有关，与历史数据⽆关的情况。
对过程的状态预测效果良好，可考虑⽤于⽣产现场危险状态的预测。
不适宜于中长期预测。
插值与拟合适⽤于物体轨迹图像的模型。
例如，导弹的运动轨迹测量的预测分析。
分为曲线拟合和曲⾯拟合，通过找到⼀个函数使得拟合原来的曲线，这个拟合程度可以⽤⼀个指标来进⾏判断。
神经元⽹络在控制与优化、预测与管理、模式识别与图像处理、通信等⽅⾯有⼗分⼴泛的应⽤。
多层前向BP⽹络适⽤于求解内部机制复杂的问题，有⼀定的推⼴、概括能⼒。
数学建模之预测模型总结
数学建模之预测模型总结数学建模是一种通过数学方法解决实际问题的过程，它可以帮助我们理解和预测各种现实世界中的现象。
在数学建模中，预测模型是一个非常重要的部分，它可以帮助我们预测未来的趋势和结果，为决策提供重要的参考依据。
本文将从数学建模的角度出发，总结预测模型的基本原理和常见方法。
预测模型的基本原理。
预测模型的基本原理是通过已知的数据来建立一个数学模型，然后利用这个模型来预测未来的结果。
在建立模型的过程中，我们需要首先确定预测的目标，然后收集相关的数据，进行数据分析和处理，最后选择合适的数学方法建立模型。
预测模型的建立过程需要考虑到多种因素，如数据的可靠性、模型的可解释性和预测的准确性等。
常见的预测模型方法。
在数学建模中，有许多常见的预测模型方法，其中最常见的包括线性回归模型、时间序列分析、神经网络模型和机器学习模型等。
下面将对这些方法进行简要介绍。
线性回归模型是一种基本的预测模型方法，它假设自变量和因变量之间存在线性关系，并通过最小二乘法来估计模型参数。
线性回归模型简单易懂，但对数据的要求较高，需要满足一些前提条件才能得到可靠的结果。
时间序列分析是一种专门用于处理时间序列数据的预测模型方法，它包括自回归模型、移动平均模型和ARIMA模型等。
时间序列分析适用于具有一定规律性和周期性的数据，可以很好地捕捉数据的趋势和季节性变化。
神经网络模型是一种基于人工神经网络的预测模型方法，它通过模拟人脑神经元之间的连接来实现对复杂非线性关系的建模。
神经网络模型适用于大规模数据和复杂问题，但需要大量的数据和计算资源来训练模型。
机器学习模型是一种基于数据驱动的预测模型方法，它包括决策树、随机森林、支持向量机和深度学习等。
机器学习模型适用于大规模数据和复杂问题，可以自动学习数据的特征和规律，但对数据的质量和标注要求较高。
预测模型的应用领域。
预测模型在各个领域都有着广泛的应用，如经济学、金融学、管理学、环境科学、医学和工程等。
数学建模的预测模型
11周的大米销售数量。
表1 4
§1 时间序列预测法
•分析： –大米是日常生活必需品，不受季节、周期的影响；
–数据记录的时间单位为周，时间间隔短很少受趋势的长时期 因素的影响。
t 1 t t
F t 1为 第 t+ 1 时 期 的 时 间 序 列 预 测 值 yt Ft 为第t 为第t 时间的时间序列的实际值 时间的时间序列的预测值
（2）
为平滑系数
0
1
求解例1
分析: 为了预测第11周的大米销售量，除了要知道前10周的实际销量 外，还要知道第10周的预测值。而要知道第10周的预测值，必须知道第9
在此用移动平均法预测。
步骤：
1.选定n的取值，取n为3； 2.选取距离第11周最近的3周数据。第8、9、10周的数据分别为54， 63，73； 3.按公式计算，得 第11周销售量预测值为：
54 63 73 3
管 理 运 筹 学
5
6 3.3 3
§1 时间序列预测法
（续例1 ）
4.用同样方法，获得第4—10周各周的预测值。
147.93
158.76 170.78
13
§1 时间序列预测法
三、用时间序列趋势进行预测
假定时间序列趋势为线性。
例2 某种品牌的冰箱最近十年的销售数量,如表5所示：
表5
年（t） 销量（万台）（yt） 年（t） 销量（万台）（yt）
1
2
40.3
44.2
6
7
54.8
64.1
评价模型预测模型优化模型数理统计模型
评价模型预测模型优化模型数理统计模型1.引言1.1 概述概述本文旨在评价模型预测模型优化模型数理统计模型，并探讨这些模型在实际应用中的价值和局限性。
模型在科学研究和实践中扮演着重要的角色，它们被广泛运用于各个领域，包括金融、医学、工程等。
通过对模型的评价、预测、优化和数理统计的研究，我们可以更好地理解和预测系统的行为，提高系统的性能和效率。
在本文中，我们将分别介绍评价模型、预测模型、优化模型和数理统计模型的概念、方法和应用。
评价模型主要关注模型的准确性、鲁棒性和可解释性，通过评估模型的性能，可以判断模型在实际应用中的可行性和可靠性。
预测模型则旨在预测未来的趋势和结果，它可以通过历史数据和统计方法来建立，并对未来的情况进行预测和分析。
优化模型则致力于寻找最优解或最优策略，通过优化模型，我们可以在给定的约束条件下达到最佳的效果。
数理统计模型是一种基于数学和统计学原理的理论模型，它能够以概率和统计的方式分析和描述数据的规律和特征。
在本文的结论部分，我们将对评价模型预测模型优化模型数理统计模型进行总结和回顾。
通过对这些模型的研究，我们可以看到它们在实际应用中的重要性和优势。
同时，我们也需要认识到这些模型存在的局限性和挑战，例如数据的质量问题、模型假设的合理性等。
在未来的研究中，我们需要继续优化和改进这些模型，以更好地应对实际问题和需求。
总之，本文将对评价模型预测模型优化模型数理统计模型进行深入研究和探讨，并总结它们在实际应用中的价值和局限性。
通过对这些模型的理解和应用，我们可以推动科学研究和实践的发展，并提高系统的性能和效率。
文章结构部分的内容可以如下编写：1.2 文章结构本文分为引言、正文和结论三个部分。
具体结构如下：引言部分首先对文章的主题进行了概述，介绍了评价模型、预测模型、优化模型和数理统计模型这四个主要内容，并指出了本文的目的。
正文部分主要分为四个部分，分别是评价模型、预测模型、优化模型和数理统计模型。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
§4 名词解释与符号说明
一、名词解释 1．土壤背2．土壤重金属污染：由于人类活动将重金属带入到土壤中，致使土壤中重金属含 量明显高于背景含量，并可能造成现存或潜在的土壤质量退化、生态与环境恶化的现象
[2]。
3．污染源：造成环境污染的污染物发生源，通常指向环境排放有害物质或对环境 产生有害影响的场所、设备、装置或人体。 4．地质环境：自然环境的一种，指由岩石圈、水圈和大气圈组成的环境系统。 二、符号说明 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 符号 符号说明 Pij 第 i 区域土壤中第 j 种元素的污染指数（ i = 1, L 5 ， j = 1,2 L8 ）
Pi C ij S ij
第 i 区域综合污染指数（ i = 1, L 5 ） 第 i 区域土壤中第 j 种重金属元素的平均实测浓度 （ i = 1, L 5 ，j = 1,2 L8 ） 第 i 区域土壤中第 j 种重金属元素的平均背景值（ i = 1, L 5 ， j = 1,2 L8 ） 参与评价的重金属种类总数 综合污染指数值 重金属浓度的特征向量矩阵 分别代表母质、废弃物、废水和沉降因子（ i = 1,L 4 ） 不同年份的地质环境综合评价值（ i = 1, 2,L , 30 ） 地质环境评价指标( i = 1, 2,..., 7 ) 非负约束因子的矩阵迭代运算终止的判定条件参数 非负因子载荷矩阵 第 i 种评价因子的相关性比值( i = 1, 2,...,8 ) 修正内梅罗指数法的第 i 个指标权重值( i = 1, 2,...,8 )
§2 问题的分析
城市表层土壤重金属污染状况涉及多方面问题，比如如何确定污染源、污染原因以 及评价城区污染程度。为了解决此类问题，更好研究城市地质环境的演变模式，我们分
2
析题目及附录中所给的相关信息，从不同角度建立多个模型分别加以讨论。 一、对问题一的分析 问题一涉及到两个小问题，首先要给出 8 种主要重金属元素在该城区的空间分布， 以直观了解不同区域重金属的污染程度。 通过附表数据， 利用 MATLAB 画出各种重金属 污染物在五个不同区域环境的分布图。对该城区土壤 319 个采样点的 8 种重金属浓度实 测均值与背景均值（见图 1）进行比较，可以发现 8 种重金属浓度平均值都高于背景值 上限，说明该城区存在一定程度的污染。其次要定量解决该城区内不同区域重金属的污 染程度问题，需建立一个指标进行综合评价。采用单因子污染指数法和内梅罗综合污染 指数法，构建城市土壤污染的评价指标体系，利用其计算式得到对应的指数值，定量评 价不同区域土壤重金属污染的程度。
§3 模型的假设
1．由于该城区海拔变化较小，在建模过程中不考虑海拔对重金属浓度含量的影响； 2．在因子分析中，假设公共因子彼此不相关且具有单位方差；
3
3. 假设通过 MATLAB 运用插值法所得数据能较好反映各种污染物的实际分布； 4．假设地质环境综合评价模型中各因素所需数据均能合理的渠道获取； 5．我们在确定污染源位置时，假设污染物由污染源向四周连续扩散； 6．在研究地质环境演变模式时，假设所考虑地区地势没有较大的变化。
1
§1 问题的提出
一、背景知识 1．概述 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响 日显突出，导致城市土壤污染日益严重。对城市土壤地质环境异常的查证，以及如何应 用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境 的演变模式，日益成为人们关注的焦点。因此，对城市土壤污染的研究分析工作显得尤 为重要。 2．现状与对策 城市土壤作为城市环境中重金属的主要蓄积库，很大程度上反映了城市环境受重金 属污染状况。 土壤重金属污染主要来源于人为污染源的输入， 主要是工业生产活动中“三 废”物质的排放、交通运输过程中产生的废物、居民生活中丢弃的废弃物质等。考虑到 不同的区域环境受人类活动影响的程度不同，将城区按功能划分为生活区、工业区、山 区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、……、5 类区。污染源在 城市区域分布的差异性导致了不同功能区土壤中重金属的含量有所不同。土壤重金属含 量影响植物、动物的基本元素的组成 ，直接影响人类和其它生物的健康。 我们针对某城市城区土壤重金属污染状况，选取 319 个土壤采样点作为样本，将所 考察的城区划分为间距 1 公里左右的网格子区域，按照每平方公里 1 个采样点对表层土 （0~10 厘米深度）进行取样、编号，并用 GPS 记录采样点的位置。应用专门仪器测试 分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照 2 公里的间距 在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 利用这些数据开展该城市环境质量评价，研究该城市地质环境的演变模式，可为城市生 态环境调控及其可持续发展提供基础科学决策依据。 二、相关试验数据 针对某城市城区土壤重金属污染状况，得到该城市城区表层土壤的相关信息： 1．表层土壤 319 个采样点位置、海拔高度及其所属功能区的信息数据（见附件的 表 1） ； 2．表层土壤 319 个采样点 8 种主要重金属元素浓度的信息数据（见附件的表 2） ； 3．表层土壤 8 种主要重金属元素背景值的信息数据（见附件的表 3） 。 三、要解决的问题 1．问题一：根据附表的相关数据，给出 8 种主要重金属元素在该城区的空间分布， 并分析该城区内不同区域重金属的污染程度； 2．问题二：通过对数据分析，说明重金属污染的主要原因； 3．问题三：分析重金属污染物的传播特征，由此建立模型，确定污染源的位置； 4．问题四：通过对前三个问题讨论，分析所建立模型的优缺点，为更好地研究城 市地质环境的演变模式，搜集相关信息，并利用这些信息建立模型解决问题。
n
P* V Xi
Fi Ai
α
C Ri Wi
§5 模型的建立与求解
从所要解决的问题和对问题所做的假设出发，分别对四个问题进行详细的分析与求 解。 一、模型Ⅰ 重金属污染综合评价模型 重金属污染综合评价模型 1．模型的分析 问题一涉及到两个小问题， 对于第一小问， 需借助 MATLAB 软件分别画出 8 种重金 [1] 属浓度的空间分布图 。由于该城区的海拔为 0~308m，变化幅度较小。为了研究方便，
Cd 的空间分布
Cr 的空间分布
Cu 的空间分布
Zn 的空间分布
Ni 的空间分布
5
Pb 的空间分布 图2
Zn 的空间分布 8 种重金属元素的分布
由图 2 可以看出，8 重金属元素没有一个固定的变化趋势，在各区域均有不同特点 的分布，无法看出各采样点的综合污染水平，整体看污染源呈零星分布。 对于第二小问，评价一个地区污染程度的方法有很多，目前常用方法有模糊贴近度 法、综合评价法、主成分分析法等，但这些方法都存在一定的缺陷性，不能明确的分出 不同区域环境的污染程度。这里我们采用单因子污染指数法和内梅罗综合指数法[2][3]， 其中单因子污染指数法是计算出 8 种重金属污染指数，可以表示出某个区域环境的不同 重金属污染程度，评价较为具体；内梅罗综合指数法是将 8 种重金属污染性进行综合分 析，表示某个区域各种重金属的总体污染程度。为了具体全面的评价，模型要以各种重 金属背景值为评价参照标准，计算出五个功能区域 8 种重金属的单因子污染指数和内梅 罗综合指数。 2．模型的准备 模型的准备 ⑴首先计算出五个区域环境的八种重金属单因子污染指数和内梅罗综合指数。 其计算公式分别为： Pij = Cij / Sij
城市表层土壤重金属污染的 城市表层土壤重金属污染的探析
摘 要
本文针对城市表层土壤重金属污染问题，以表层土壤重金属浓度为切入点，构建 城市土壤污染的评价指标体系，建立重金属污染综合评价模型、污染因子模型、污染源 主成分模型和地质环境综合评价模型。依据附件数据，运用数值分析、多元统计分析和 模糊综合评价等相关知识与算法，得出有关城市环境质量评价的相关结论。最后经过讨 论分析，对模型做出改进和推广。 模型Ⅰ 模型Ⅰ 建立重金属污染综合评价模型。借助 MATLAB 软件，采用四维样条插值法 画出该城区 8 种重金属浓度的空间分布图。同时采用单因子污染指数法，对不同区域的 8 种重金属污染程度进行定量分析。考虑到不同区域重金属污染程度不同，我们定义了 土壤污染综合评价指标，运用内梅罗综合指数法得出生活区、工业区、山区、主干道路 区及公园绿地区的综合评价值为 3.170，13.533，1.253，9.426，2.734，表明山区重金属 污染程度最轻，工业区污染最严重，评价结果合理。 模型Ⅱ 模型Ⅱ 建立重金属污染因子模型。依据因子分析原理，分别用几个公共因子来解 释该城区不同区域土壤 8 种重金属污染指标，借助 SPSS 软件得出不同区域污染因子载 荷值， 并结合不同功能区域污染的具体情况， 进行综合分析得出整个城区重金属污染的 主要原因：工矿企业污染、燃煤污染、交通污染和居民生活污染。 模型Ⅲ 模型Ⅲ 建立污染源主成分模型。利用主成分分析法，得到 8 种重金属污染指标的 四个主成分指标，即母质因子、废物因子、废水因子和沉降因子。借助 MATLAB 软件 编程计算出土壤 319 个采样点不同因子得分， 利用排名结果判断出各属性因子的污染源 位置。其中废物因子污染源在以主干道路区点（3299，6018，4）为中心的一平方公里 、工业区点（1647， 范围内；废水因子污染源在以主干道路区点（18134，10046，41） 2728， 6） 为中心的一平方公里范围内； 沉降因子污染源在以主干道路区点 （18134， 10046， 41）为中心的一平方公里范围内。 模型Ⅳ 模型Ⅳ 建立地质环境综合评价模型。结合前三个模型，针对城市地质环境演变模 式，总结区域地质、水文环境、土壤质量及地质灾害等主要因素，最终选择 7 个地质环 境评价指标，构建地质环境综合评价体系。通过综合评价值的变化趋势，可以完成对城 市地质环境演变模式的分析。 在模型改进中，建立模型Ⅴ（非负因子约束模型） ，消除了模型Ⅱ中因子载荷矩阵 出现负值的情况，使模型结果贴近实际，具有较强的解释力。 最后，本文对模型进行了误差分析、优缺点评价与推广。 关键词：重金属污染；地质环境演变；内梅罗综合指数法；因子分析； MATLAB 关键词：