黄河小浪底调水调沙数学建模matlab 程序

黄河小浪底调水调沙问题

摘要：本文利用插值拟合的方法通过Matlab工具模拟出了排沙量与时间、排沙量与水

流量的函数关系，并且求出了总排沙量为1.704亿吨。整个模型简单且方便计算，其中排沙量与水流量的函数关系为分段函数。

关键词：调水调沙 Matlab 插值拟合

一、问题重述

2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：

表1: 试验观测数据单位：水流为立方米/ 秒，含沙量为公斤/ 立方米

(1) 给出估算任意时

篇一：黄河小浪底调水调沙

黄河小浪底调水调沙问题

2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：

表1: 试验观测数据 单位：水流为立方米 / 秒，含沙量为公斤 / 立方米

(1) 给出估算任意时刻的排沙量及总排沙量的方法； (2) 确定排沙量与水流量的变化关系。

篇二：黄河小浪底调水调沙工程数学实验实验报告

《数学实验》实验报告

题目：黄河小浪底调水调沙工程

姓名： 胡 迪学号：201014622 专业：信息与计算

黄河调水调沙利弊之我见(一)

【论文关键词】黄河调水调沙利与弊解析

【论文摘要】黄河泥沙是黄河洪水威胁的症结所在，是黄河治理的难点。黄河的防洪思想、防洪技术、防洪战略决策都着眼于这一点。为确保黄河下游河道不抬高、不淤积、不断流，确保我国经济发展和沿黄地区的生命财产安全。黄河水利委员会创立了举世瞩目用水治沙的历史先河。“调水调沙”的设想经过实践的考验，正一步步走向成功，对防洪、治河、减淤等发挥了巨大作用，也具有较大的效益。但是，笔者认为，在看到有益方面的同时，也存在着这样或那样的不利因素，现结合实际，就调水调沙对黄河下游河道产生利弊问题谈几点粗浅的见解，仅供同仁参考。

1黄河泥沙的危害

黄河水患，根在泥沙。这条世界上含沙量最大的河流素有“斗水七沙”之称，每年携带16亿吨泥沙进入下游，其中有4亿吨淤积在下游河床，使之以每年平均10厘米的速度抬升，目前下游河床平均高出临河地面4－6米，形成举世闻名的“地上悬河”，“悬河”如同悬挂在下游两岸头上的一个大水盆，严重威胁着我国经济发展、社会稳定和沿黄地区的生命财产安全。2什么是调水调沙

所谓调水调沙，就是在充分考虑黄河下游河道输沙能力的前提下，利用水库的调节库容，对水沙进行有效的控制和调节，适时蓄存或泄放，调整天然

小 浪 底

各位游客朋友你们好，一路辛苦了。我代表乐土旅行社对大家的到来表示热烈的欢迎，我是各位本次河南之旅的导游员，大家可以叫我X导。在河南的这段时间将由我向大家提供导游服务。除了我之外，还有一位，就是在我身旁正在专心开车的X师傅。在接下来几天的行程当中，我们将一起给大家提供行程中的各项服务，大家什么有要求或建议，随时可以提出来，我们会尽量满足各位的要求，努力使大家在河南玩的开心。最后，预祝大家在河南度过一个愉快的旅程。

我们今天要参观的景点是黄河小浪底风景区。我们从郑州出发，大概两个半小时的车程就能抵达，在这段时间里，我先简短的向大家介绍一下河南的概况。

河南是中华文明的重要发祥地之一。在中华民族5000年历史长河中，先后有20多个朝代200多位帝王建都或迁都于此，作为全国政治、经济、文化中心长达3000年之久。中国八大古都河南有其四，即九朝古都洛阳、七朝古都开封、殷商古都安阳、商都郑州。在众多的景区中，洛阳的龙门石窟（2000年）、安阳的殷墟（2006）、郑州的“天地之中”登封古建筑群（2010年）被列入世界文化遗产名录。河南被称为“中国历史文化的博物馆”，地下文物和地上文物分别居全国第一位和第二位（陕西的地上文物第一，地下文物第二）。

MATLAB实用程序百例

实例1：三角函数曲线（1） function shili01 h0=figure('toolbar','none',...

'position',[198 56 350 300],... 'name','实例01'); h1=axes('parent',h0,...

'visible','off'); x=-pi:0.05:pi; y=sin(x); plot(x,y);

xlabel('自变量X'); ylabel('函数值Y'); title('SIN( )函数曲线'); grid on

实例2：三角函数曲线（2） function shili02 h0=figure('toolbar','none',...

'position',[200 150 450 350],... 'name','实例02'); x=-pi:0.05:pi; y=sin(x)+cos(x);

plot(x,y,'-*r','linewidth',1); grid on

xlabel('自变量X'); ylabel('函数值Y'); title('三角函数');

file:///E|/Document/发展篇/M&M/竞赛篇/

MATLAB数学建模习题1

一、单项选择题（将选择答案写在答题纸上，每小题2分共20分）

1．在MATLAB命令窗口中键入命令，Vname=prod(7:9)/prod(1:3)，可计算组合数

如果省略了变量名Vname，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名

A）ans； B）pi； C）NaN； D）eps

2．宝石切割问题中，石料左右长度、前后长度、上下高度分别为a1、a2、a3，即a1×a2×a3(cm3)，而精品尺寸为b1×b2×b3(cm3)。操作时，同向切割连续两次再旋转刀具。某一切割方案的切割面积依次为：2a1a2? 2a1b3 ? 2b2b3，则这一切割方案为

A）左右?前后?上下； B）上下?前后?左右； C）前后?上下?左右； D）前后? 左右?上下 3．机场指挥塔位置：北纬30度35.343分，东经104度2.441分，在MATLAB中用变量B=[30 35.343]表达纬度，L=[104 2.441]表达经度。将数据转化为以度为单位的实数，下面正确的语句是

A） P=B(1)+B(2)/60，Q=L(1)+L(2)； B) P =

数统 应数 20121323003 王妍 数统 应数 20121323022 胡可旺

1 生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布； 代码：y=1+sqrt(4)*randn(5) 结果:

2，生成一个20行5列矩阵A，其元素分别以概率0.7、0.1、0.2取自0、1、4这三个整数， 然后将“删除A的全零行”之后的结果赋给矩阵B； for i = 1:20 for j = 1:5

p = rand(); if p<=0.7

A(i,j) = 0;

elseif p>0.7&&p<=0.8 A(i,j) = 1; else

A(i,j) = 4; end end end B = []; for i = 1:20 flag = true; for j = 1:5

if A(i,j)~=0 flag=false; break; end end if ~flag

B = [B;A(i,:)]; end end B

3, （

2、已知速度曲线v(t) 上的四个数据点下表所示

t=[0.15,0.16,0.17,0.18];

v=[3.5,1.5,2.5,2.8];

x=0.15:0.001:0.18

y=i n t e r p1(t,v,x,'s p l i n e')

S=t r a p z(x,y)

p=p o l y f i t(x,y,5);

d p=p o l y d

e r(p);

d p x=p o l y v a l(d p,0.18)

运行结果

S=

0.0687

Dpx=-

3、计算图片文件tu.bmp 给出的两个圆A，B 的圆心，和两个圆的两条外公切线和两条内公切线的切点的坐标。

（1）计算A 圆的圆心坐标

I=imread('tu.bmp');

[m,n]=size(I)

BW=im2bw(I)

BW(:,200:512)=1;

figure, imshow(BW)

ed=edge(BW);

[y,x]=find(ed);

x0=mean(x), y0=mean(y)

r1=max(x)-min(x),r2=max(y)-min(y)

r=(r1+r2)/4

x0 =109.7516

y0 =86.7495

r1 =162

r2 =158

r =80

（2）B圆的圆心坐标和半径

I=imread('tu

matlab 数学建模

1、设A、B两方案的净现金流量（单位：万元）如下表所示： （1）设折现率为10%，计算两个方案的净现值；

（2）计算两个方案的内部收益率。

2、某厂生产的一种电器的销售量y与竞争对手的价格x1和本厂的价格x2有关。下表是该厂在二十个城市的销售记录。

matlab 数学建模

（1） 根据这些数据建立本厂的需求函数模型；

设方程为 y=a + b*x1 + c*x2 Matlab 程序设计如下：

x1=[ 100 110 105 115 118 145 112 165 170 190 156 268 150 200 230 160 270 250 280 286];

x2=[120 140 138 130 136 148 110 150 165 160 130 145 166 145 150 140 180 145 200 220]; y=[102 100 110 115 105 98 95 93 90 89 77 69 92 60 85 82 65 69];

X=[ones(20,1) x1’ x2’];

[b,bint,r,rint,stats]

matlab入门知识

数学建模与数学实验MATLAB入门 入门

后勤工程学院数学教研室

matlab入门知识

●

MATLAB作为线性系统的一种分析和仿真

工具，是理工科大学生应该掌握的技术工具，它 作为一种编程语言和可视化工具，可解决工程 工程、 工程 科学计算和数学学科 科学计算 数学学科中许多问题。 数学学科 向量、数组 矩阵的基 ● MATLAB建立在向量 数组 矩阵 向量 数组和矩阵 础上，使用方便，人机界面直观，输出结果可 视化。 ● ● 矩阵是MATLAB的核心 矩阵 MATLAB的进入与运行方式(两种)

matlab入门知识

MATLAB入门 入门一、变 量 与 函 数 二、数 三、 矩 组 阵

四、 MATLAB编程 编程 五、 实 验 作 业

matlab入门知识

一、变 量 与 函 数1、变量 、 MATLAB中变量的命名规则 命名规则是： 命名规则 (1)变量名必须是不含空格的单个词； (2)变量名区分大小写； (3)变量名最多不超过19个字符； (4)变量名必须以字母打头，之后可以是 任意字母、数字或下划线，变量名中 不允许使用标点符号.

matlab入门知识

特殊变量表特殊变量 ans pi eps flops inf NaN i,j na