黄河小浪底调水调沙工程数学实验实验报告
小浪底实习报告
随着我国水电事业的蓬勃发展，小浪底水利枢纽工程作为世界上最大的水利枢纽工程之一，其建设和管理对于我国乃至全球的水电事业都具有重要的示范意义。
为了更好地了解水利枢纽工程的建设和管理，提升自身的专业素养，我于2021年7月至9月在小浪底水利枢纽工程进行了为期两个月的实习。
二、实习单位简介小浪底水利枢纽工程位于河南省洛阳市孟津县，是黄河流域治理的关键工程，也是我国水利工程的典范。
工程于1997年开工建设，2003年全面建成，总装机容量为1,820万千瓦，是我国目前最大的水利枢纽工程。
三、实习内容在实习期间，我主要参与了以下工作：1. 施工现场观摩通过实地观摩施工现场，我了解了水利枢纽工程的施工过程和工艺，包括大坝建设、电站安装、导流洞施工等。
这让我对水利枢纽工程的建造有了直观的认识。
2. 工程管理学习在工程管理方面，我学习了工程项目的组织架构、项目管理流程、质量安全管理等知识。
通过参与工程例会，我了解了工程项目的进度、质量、安全等方面的控制情况。
3. 设备操作培训在设备操作方面，我接受了水轮发电机组、水工建筑物、监测系统等设备的操作培训。
通过实际操作，我掌握了设备的操作技能和安全注意事项。
4. 数据分析与处理在数据分析与处理方面，我参与了水库调度、发电量统计、设备运行状态分析等工作。
通过这些工作，我熟悉了水利枢纽工程的数据收集、处理和分析方法。
5. 现场调研与记录在实习期间，我还参与了现场调研活动，对工程周边的环境、水资源、社会经济状况进行了调查，并进行了详细的记录和分析。
1. 专业知识提升通过实习，我对水利枢纽工程的建设和管理有了更加深入的了解，掌握了相关的专业知识和技能。
2. 实践能力增强实习过程中，我参与了实际工程项目的操作和管理，锻炼了我的实践能力，提高了我的问题解决能力。
3. 团队协作意识增强在实习过程中，我与团队成员共同完成了各项任务，增强了团队协作意识。
4. 社会责任感提升通过实习，我深刻认识到水利枢纽工程对于国家和社会的重要性，增强了社会责任感。
黄河小浪底调水调沙问题数学建模
黄河小浪底调水调沙问题数学建模
黄河小浪底调水调沙问题是指通过调整黄河水流的水量和输沙量来解决黄河小浪底河道的淤积问题。
数学建模可以帮助我们分析和预测黄河小浪底的水流和沙传输规律，从而提出合理的调水和调沙方案。
以下是数学建模中可能涉及的一些步骤和方法：
1. 数据收集和处理：收集黄河小浪底相关的水文数据、地质资料和历史数据，对数据进行整理和处理，建立合适的数据模型。
2. 建立水流模型：通过流体力学理论和水流实验数据，建立黄河小浪底水流的数学模型，包括水流速度、水动力和水力调控方面的参数。
3. 建立沙传输模型：根据黄河小浪底河道的地质特征和沙传输规律，建立沙传输的数学模型，包括输沙通道的沙动力和沙质输运规律方面的参数。
4. 模型验证和参数拟合：利用已有的观测数据和实验数据验证建立的水流和沙传输模型，并通过参数拟合来优化模型的准确性和适用性。
5. 模拟预测和优化调控：利用建立的数学模型，进行水流和沙传输的模拟预测，通过调整输水和输沙量来优化黄河小浪底的调水和调沙方案，以达到降低淤积和维护航道的目的。
数学建模可以辅助相关专业的研究人员和决策者做出科学的决策，使调水和调沙方案更加合理和有效，减少淤积和保护黄河流域的生态环境。
黄河小浪底调水调沙问题数学建模
黄河小浪底调水调沙问题数学建模黄河是中国第二长河流，也是中国北方主要的水源之一。
然而，由于年际变化和人类活动的影响，黄河水沙特性的变化对地区社会经济和生态环境产生了巨大影响。
黄河小浪底是黄河下游的一个关键水文站点，对黄河的水沙调控起着重要作用。
因此，对于黄河小浪底的调水调沙问题进行数学建模具有重要意义。
数学建模是通过数学方法分析和解决实际问题的过程。
对于黄河小浪底的调水调沙问题，我们可以从以下几个方面进行数学建模：1. 水量平衡模型：黄河小浪底是一个重要的水源供给站点，掌握黄河的水量情况对于调水调沙至关重要。
因此，我们可以建立一个水量平衡模型，根据入库、出库等因素来估计黄河在小浪底的流量。
这个模型可以包括如下因素：入流量（降雨、地表径流、地下径流等）、出流量（供水、排水等）以及河道水量的变化。
通过这个模型，可以对黄河小浪底的水量进行预测和调控。
2. 水沙关系模型：黄河的水沙关系对于调水调沙具有重要影响。
水沙关系模型可以通过分析黄河不同断面的水位和水沙含量之间的关系，来估计黄河的河床输沙量。
这个模型可以包括如下因素：断面形态特征、流量、水沙含量等。
通过这个模型，可以了解到黄河的水沙变化规律，并对黄河小浪底的调沙情况进行预测和控制。
3. 沉积模型：黄河的床面沉积是一个长期过程，对于调水调沙有着重要影响。
沉积模型可以通过分析黄河不同断面的沉积速率、沉积厚度等变化，来估计黄河的床面沉积情况。
这个模型可以包括如下因素：流率、输沙率、流态等。
通过这个模型，可以对黄河小浪底的沉积情况进行预测和控制。
4. 排沙方案优化模型：为了减少黄河小浪底的沙泥淤积问题，需要设计科学合理的排沙方案。
排沙方案优化模型可以通过考虑沙泥淤积的成因、河道特征、水流特性等因素，来确定最佳的排沙方案。
这个模型可以包括如下因素：流态、输沙率、河道形态等。
通过这个模型，可以设计出最优的排沙方案，从而实现黄河小浪底的水沙调控。
综上所述，黄河小浪底的调水调沙问题可以通过数学建模的方式来研究和解决。
小浪底水利工程实训报告
一、实训背景随着我国水利事业的发展，水利工程在国民经济和社会发展中扮演着越来越重要的角色。
为了更好地了解水利工程的实际情况，提高自己的专业素养，我于近日参加了小浪底水利工程实训。
小浪底水利枢纽工程作为我国黄河流域治理的重要工程，其建设过程、运行管理以及社会经济效益等方面都具有典型意义。
二、实训目的1. 了解小浪底水利枢纽工程的基本情况，包括工程规模、建设背景、主要功能等。
2. 学习水利工程的基本原理、施工技术和管理方法。
3. 培养实际操作能力，提高解决实际问题的能力。
4. 增强团队合作意识和沟通能力。
三、实训内容1. 工程概况小浪底水利枢纽工程位于河南省洛阳市以北40千米，是黄河干流三门峡以下唯一能取得较大库容的控制性工程。
该工程以防洪、防凌、减淤为主要功能，兼顾供水、灌溉和发电，是治理开发黄河的关键性工程。
2. 工程建设小浪底工程自1991年4月开工，历时11年，共完成土石方挖填9478万m³，混凝土348万m³，钢结构3万T，安置移民20万人。
在工程建设过程中，面对塌方、设计变更、施工管理等各种困难，建设者以高度的主人翁责任感，强烈的爱国主义情怀，沉着应对，奋勇拼搏，最终取得了工期提前、投资节约、质量优良的好成绩。
3. 工程运行与管理小浪底水利枢纽工程自2001年底主体工程全面完工以来，已正常运行多年。
在运行管理方面，工程严格执行国家法律法规和行业标准，确保工程安全、可靠、高效运行。
4. 调水冲沙实验2002年，我国在小浪底水利枢纽工程成功进行了调水冲沙实验，以减轻黄河下游河床的淤沙程度。
实验结果表明，小浪底工程在减轻河床淤积、提高黄河水资源利用效率等方面取得了显著成效。
四、实训体会1. 工程规模宏大，技术复杂小浪底水利枢纽工程是我国水利史上最具挑战性的项目之一，技术复杂，施工难度大。
在工程建设过程中，建设者充分发挥了聪明才智，克服了重重困难，为我国水利事业树立了典范。
2. 以人为本，关注民生小浪底工程建设过程中，高度重视移民安置工作，确保移民利益得到充分保障。
小浪底调水调沙
1.摘要本文综合利用试验观测的数据，根据曲线拟合的方法合理的估算任意时刻黄河小浪底水库的排沙最和总排沙量，同时根据排沙帚:和水流最之间的关系图像大致拟合出两者之间的关系。
小浪底观测站从6月29到7月10 口，每天分别在早晨8点和夜晚8点对水流量”和水的含沙量Y,做出检测，针对问题一，我们知道，任意时刻的排沙量Y等于该时刻的水的含沙最与水流最之积，利用表一中各数据点，把时间点依次记作t（l<=t<=24）利用曲线拟合，得到排沙量和时间之间的关系，再利用辛普森公式求得某一时间段的总排沙量。
针对问题二，从排沙量和水流量的关系图像上我们可以大致看出一•者成一•次关系，同样利用曲线拟合就可以得出排沙最和水流量的关系了。
关键词：含沙星水流星排沙量曲线拟合辛普森积分公式2 •问题重述在小浪底水库蓄水后，黄河水利委员会进行了多次试验，特别是2004年6月到7月进行的黄河第三次调水调沙试验具有典型的意义。
这次试验首次由小浪底、三门峡和万家寨三大水库联合调度，进行接力式防洪预泄放水，形成人造洪峰进行调水调沙试验成功。
这次试验的一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底的库区的沉积泥沙。
在小浪底开闸泄洪以后，从6月2711开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29口先后到达小浪底，7月3口达到最大流星为2720m'3/s,使小浪底水库的排沙量也不断的增加。
表1:试验观测数据单位：水流为立方米/秒，含沙最为公斤/立方米（1）根据表一中的观测数据，合理的建立模型來估算任意时刻的排沙量和总的排沙量。
（2）确定排沙皐和水流量之间的关系。
3•问题分析含沙量即是单位体积的水中沙的含量，排沙量是指单位体积的水中带走的沙的星,因此可知:排沙量二水的存沙量X水流量静水时，水的含沙量是很少的，由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底的库底的沉积泥沙，使得水屮的含沙量增多，随着人造洪峰的逐渐平息，使水流量的减少，也使得含沙星也呈逐渐减少的趋势，所以排沙量也是减少的，由丁•洪峰在正常情况下是随着时间逐渐平息的，所以, 水的含沙量和排沙最都和时间存在着一定的关系，根据表一，水的含沙量和水流量的积求出该时刻水的排沙量，将时间按记录次库依次记作1 2 3 ........................................................... 24,然后按二次函数拟合即可得到排沙暈和时间的关系式，然后插值得到任意时刻的排沙量，对排沙最和时间的关系式在任意时间段上积分求得值就是这段时间总的排沙量。
黄河小浪底调水调沙
•
图4
所以拟合后的函数为V= 95*t^3-5.5e+003*t^2+7.7e+004*t3.2e+004，通过图像可以看出排沙量与时间近似服从正态分 布，进行拟合。
计算总含沙量
通过Matlab编程可以计算出定积分，结果如下
%jisuan.m syms t; S=0.014*t^3-1.3*t^2+21*t+16; v=0.13*t^3-14*t^2+2.4e+002*t+1.5e+003; V=v*S; simple(V); syms t; V=95*t^3-5.5e+003*t^2+7.7e+004*t-3.2e+004; int(12*60*60*V,t,0,24) ans =170366976000 即总含沙量为1.704亿吨
表1: 试验观测数据
日期 时间 水流量 8:00 1800
单位：水流为立方米 / 秒，含沙量为公斤 / 立方米
6.29 20:00 1900 8:00
6.30 20:00 2200 8:00 2300
7.1 20:00 2400 8:00 2500
7.2 20:00 2600 8:00 2650
7.3 20:00 2700 8:00 2720
因为某一时刻的排沙量V=v(t)S(t)，所以将所拟合出来的多项 式带入上式，通过Matlab进行计算可以得到下面答案 即排沙量与时间的关系为： V=0.0018*t^6-0.365*t^5+24.29*t^4-582.92*t^3 +2866*t^2+35340*t+24000
由于这里的多项式次数过高，对其再进行一次拟合， 有下面结果：
最新-2019小浪底调水调沙原理 精品
2019小浪底调水调沙原理篇一：黄河小浪底调水调沙黄河小浪底调水调沙问题2019年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为755亿立方米，在这之前，小浪底共积泥沙达1415亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：表1试验观测数据单位：水流为立方米秒，含沙量为公斤立方米(1)给出估算任意时刻的排沙量及总排沙量的方法；(2)确定排沙量与水流量的变化关系。
篇二：黄河小浪底调水调沙工程数学实验实验报告《数学实验》实验报告题目：黄河小浪底调水调沙工程姓名：胡迪学号：201914622专业：信息与计算科学黄河小浪底调水调沙问题2019年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功。
整个试验期为20多天，小浪底从6月19日开始预泄放水，至到7月13日恢复正常供水结束。
小浪底水利工程按设计拦沙量为755亿3,在这之前，小浪底共积泥沙达1415亿。
这次调水调沙试验一个重要的目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700，使小浪底水库的排沙量也不断地增加。
黄河小浪底调水调沙灰色数学研究
第１ ９卷 第 ４期 ２０ ０ ７年 １ ２月
河南纺织 高等专科 学校 学报
Ｊ ＯＵＲＮＡ Ｌ ＯＦ ＨＥＮ ＥＸ ’ Ｅ ＣＯ Ｌ ＧＥ ＡＮ Ｔ １ Ｌ Ｌ Ｅ Ｉ
Ｖ １１ Ｎ ． ０． ９． ｏ ４ Ｄ ｅ ，０ ７ ｅ ．２ ０
束 。这次调 水调沙试 验 的一个 重要 目的就是 冲刷小 浪底库 区沉 积 的泥 沙 。在 小 浪 底 水 库 开 闸 泄 洪 以
收 稿 日期 ：０ ７— ９—１ ２０ ０ ７
根据上 述试 验数据 ， 我们 利用 灰色系统 理论 ， 研
作者简介： 饶明贵（９ ５一） 男， 南光 山人 ， １６ 、 河 讲师 ， 主要研 究应用数学。
不断增 加 。表 １是 由小 浪 底 观 测站 从 ２０ ０４年 ６月 ２ 日到 ７月 １ ９ Ｏ日检测 到 的试验数据 。
表 １ 试 验 观 测 数 据
信息不完备的系统 ， 通过 已知信 息来研 究和预测 未 知 领 域 从 而 达 到 了 解 整 个 系 统 的 目 的 。 灰 色 数学与概率论 、 模糊数 学一起并 称为研 究不确定 性 系统 的三种 常用 方 法 ， 有 能 够 利 用 “ 数 据 ” 具 少 建模 寻 求 现 实规 律 的 良好 特 性 ， 服 了 数 据 不 足 克
或 系统 周 期 短 的 问题 。
目前 ， 色 数 学 理 论 … 得 到 了 极 为 广 泛 的应 灰 用 ， 仅 成功 地 应 用 于 工 程 控 制 、 济 管 理 、 会 不 经 社 系统 、 态 系 统 等 领 域 ， 且 在 复 杂 多 变 的 农 业 生 而
系统 ， 如在水利、 气象 、 生物 防治 、 机决策 、 农 农业 规划 、 农业经济等方 面 的应用 也取得 了可喜 的成 果 。笔 者 利用 灰 色 系统 理 论 ， 以黄 河 小 浪 底 第 三 次调水调沙试验数 据 为依 据 ， 研究排 沙量 与水 流 量的函数关系 ， 以期使 灰色数 学能够 在制黄工程 调 水调 沙 中发 挥作 用 。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《数学实验》实验报告题目：黄河小浪底调水调沙工程：胡迪学号： 201014622专业：信息与计算科学黄河小浪底调水调沙问题2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功。
整个试验期为20多天，小浪底从6月19日开始预泄放水，至到7月13日恢复正常供水结束。
小浪底水利工程按设计拦沙量为75.5亿m 3,在这之前，小浪底共积泥沙达14.15亿t 。
这次调水调沙试验一个重要的目的就是由小浪底上游的和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6月27日开始水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700 ，使小浪底水库的排沙量也不断地增加。
表1是由小浪底观测站从6月29日到7月10日检测到的试验数据。
表1 试验观测数据 （ 单位：水流为s m/3，含沙量为3kg/m ）现在，根据试验数据建立数学模型研究下面的问题： （1）给出估算任意时刻的排沙量及总排沙量的方法； （2）确定排沙量与水流量的变化关系。
关键词：拟合，SAS ，Matlab ，线性回归，调水调沙实验问题分析：1、对于问题一，所给数据中水流量x 和含沙量h 的乘积即为该时刻的排沙量y 即：y=hx 。
2、对于问题二，研究排沙量与排水量的关系，从实验数据中可以看出，开始排沙量随水量增加而增加，而后随水流量的增加而减少，显然变化关系并非线性的关系，为此，把问题分为两部分，从水流量增加到最大值为第一阶段，从水流量最大值到结束为第二阶段，分别来研究水流量与排沙量之间的函数关系。
模型假设：1、水流量和排沙量都是连续的，不考虑上游泄洪所带来的含沙量和外界带来的含沙量。
2、时间是连续变化的，所取时间点依次为1,2,3，…,24,单位时间为12h 。
模型的建立与求解：<一>对于问题一，因为排沙量与时间的散点图基本符合正态曲线，如图二所示。
所以，排沙量的对数与时间的函数关系就应该符合二次函数关系，因而排沙量取对数后，再与时间t 进行二次回归，排沙量取自然后的数据见表2. 假设排沙量与时间函数关系的数学模型是两边取对数得 Lny=at^2+bt+c先由表二做出排沙量的自然对数lny 与时间t 的散点图见图一，并利用SAS 软件进行拟合，得到排沙量的自然对数与时间的回归方程为： Lny=-0.0209t^2+0.4298t+10.6321由回归拟合参数表可知回归方程是显著的，因为相关系数人R^2=0.9629,误差均方S^2=0.0543，说明回归曲线拟合效果很好。
所以排沙量与时间之间的函数关系式为ec bt at y ++=2^et t y 6312.104289.02^0209.0++-=图二：排沙量对时间的曲线图最后对所求出的函数关系在区间[0,24]之间进行积分结果为总排沙量1.93962亿吨，此与媒体报道的排沙量几乎一样。
dtet t ⎰++-246312.104289.02^0209.0*60*60*12<二>对于第二个问题，两个阶段的数据如表三、表四所示序号 1 2 3 4 5 6 7 8 1800 1900 2100 2200 2300 2400 2500 2600 水流量x32 60 75 85 90 98 100 102 含沙量h序号 1 2 3 4 5 6 7 8 2650 2600 2500 2300 2200 2000 1850 1820 水流量x116 118 120 118 105 80 60 50 含沙量h对于第一阶段，有表四用MATLAB作图（如图三）可以看出其变化趋势，我们用多项式做最小二乘拟合。
设三次拟合函数关系h=a0+a1x+a2x^2+a3x^3其中a0,a1,a2,a3,为待定系数。
四次拟合函数关系h= a0+a1x+a2x^2+a3x^3+a4x^4其中a0,a1,a2,a3,a4为待定系数。
图三：第一阶段水流量与排沙量之间的关系图三次多项式拟合由MATLAB拟合函数求解出a0=a1=0,a2=0.0032,a3=-2.4929.则拟合函数h=0.0032x^2-2.4929x^3,拟合效果如图四所示图四：三次多项式拟合效果，红线为拟合曲线类似的四次多项式拟合由MATLAB拟合函数求解出a0=a1=a2=0,a3=0.0121,a4= -7.4347则拟合函数h=0.0121x^3-7.4347x^4,拟合效果如图五所示图五：四次多项式拟合效果，蓝线线为拟合曲线对于第二阶段，有表五用MATLAB作图可以看出其变化趋势，我们用多项式做最小二乘拟合。
设三次拟合函数关系h=a0+a1x+a2x^2+a3x^3其中a0,a1,a2,a3,为待定系数。
四次拟合函数关系h= a0+a1x+a2x^2+a3x^3+a4x^4其中a0,a1,a2,a3,a4为待定系数。
三次多项式拟合由MATLAB拟合函数求解出a0=a1=0,a2=-0.9475,a3= 464.9601.则拟合函数h=-0.9475x^2+464.9601x^3,拟合效果如图图六所示图六：三次拟合函数拟合效果类似的四次多项式拟合由MATLAB拟合函数求解出a0=a1=0，a2=-0.0013，a3= 1.1219a4=-354.5952则拟合函数h=-0.0013x^2+1.1219x^3-354.5952x^4,拟合效果如图七所示图七：四次拟合函数拟合效果结论以及分析检验：<一>用SAS 软件做线性回归得到排沙量与时间的函数关系式为:再利用所求函数在区间[0,24]上进行积分得到总排沙量1.93962亿吨，这与现实情况基本相符。
<二>对于第一阶段三次多项式拟合由MATLAB 拟合函数求解a0=a1=0,a2=0.0032,a3=-2.4929则拟合函数h=0.0032x^2-2.4929x^3对于第一阶段四次多项式拟合由MATLAB 拟合函数求解出a0=a1=a2=0,a3=0.0121,a4= -7.4347则拟合函数h=0.0121x^3-7.4347x^4对于第二阶段三次多项式拟合由MATLAB 拟合函数求解出a0=a1=0,a2=-0.9475,a3= 464.9601则拟合函数h=-0.9475x^2+464.9601x^3对于第二阶段四次多项式拟合由MATLAB 拟合函数求解出a0=a1=0，a2=-0.0013，a3= 1.1219 a4=-354.5952则拟合函数h=-0.0013x^2+1.1219x^3-354.5952x^4讨论与推广：1、对于第一个问题排沙量与时间不是严格的正态函数关系可能与实际有些偏差，此外还可以用SAS 软件进行高次的多向式回归2、对于第二个问题，由于MATLAB 软件的计算可能有些偏差导致拟合的函数关系可能与实际有稍微偏差，此外，还可以进行高次的拟合。
附录：1、排沙量与时间的关系图像的MATLAB 程序： >>t=1:1:24;>>y=[57600,114000,157500,187000,207000,235200,250000,265200,2862000,2400,312800,307400,306800,300000,271400,231000,160000,111000,91000,54000,45500,30000,8000,4500]; >> plot(t,y,'r')2、对排沙量求自然对数的MATLAB 程序与结果： y3=log(y) y3 =Columns 1 t hrough 1710.9613 11.6440 11.9672 12.1389 12.2405 12.3682 12.4292 12.4882 12.5644 12.6195 12.6533 12.6359 12.6340 12.6115 12.5113 12.3502 11.9829Columns 18 through 2411.6173 11.4186 10.8967 10.7255 10.3090 8.9872 8.4118 3、第一阶段的排沙量与水流量之间的关系MATLAB 程序：>>x=[1800,1900,2100,2200,2300,2400,2500,2600,2650,2700,2720]; >> h=[32,60,75,85,90,98,100,102,108,112,115];>> x1=[2650,2600,2500,2300,2200,2000,1850,1820,1800,1750,1500,1000,900];et t y 6312.104289.02^0209.0++-=>> h1=[116,118,120,118,105,80,60,50,40,32,20,8,5];>> plot(x,h,'r:')4、第一阶段三次多项式拟合函数以及拟合效果程序与结果：>> A1=polyfit(x,h,3)>> in polyfit at 80A1 =0.0000 -0.0000 0.0032 -2.4929>> z1=polyval(A1,x);plot(x,h,'k+',x,z1,'r')5、第一阶段四次多项式拟合函数以及拟合效果程序与结果：>>A2=polyfit(x,h,4)>> In polyfit at 80A2 =-0.0000 0.0000 -0.0000 0.0121 -7.4347>> z2=polyval(A2,x);>> plot(x,h,'*',x,z2,'r')6、第二阶段三次多项式拟合函数以及拟合效果程序与结果：>>A3=polyfit(x1,h1,3)>> In polyfit at 80A3 =-0.0000 0.0006 -0.9475 464.9601>> z3=polyval(A3,x1);>> plot(x,h,'*',x1,z3,'b')7、第二阶段四次多项式拟合函数以及拟合效果程序与结果：>>A4=polyfit(x1,h1,4)>> In polyfit at 80A4 =-0.0000 0.0000 -0.0013 1.1219 -354.5952>> z4=polyval(A4,x1);>> plot(x1,h1,'k*',x1,z4,'r:')参考文献：【1】启源，数学模型（第三版），高等教育【2】楼顺天,MATLAB程序设计语言（第二版），电子科技大学【3】娜，在SAS中拟合ARCH/GARCH模型。