流水渠安装过程_流水渠安装过程视频

喝的水是经过水表里流出来的吗

是的。

       地表水进入水库，水库经由引水渠流向自来水厂，在水厂完成消毒后流入地下自来水管网，最后管网接入各个小区，在水泵水箱的协作下，经过水表后从水龙头流出来。

       毛细血管般的供水管网将连通各个小区。管网末端的水压可以将水送到米左右的高度，而更高的楼层则需要需要小区物业配置水箱，通过水泵二次加压，使之能输送到高层。这些设施需要经过自来水公司的验收和测试，才能办理立户手续，安装计费水表，正式供水。有的家中会安装净水机，净水机出来的水就可以直接饮用了。如果没有净水机，则需要烧开后饮用。

“流水素面”传到中国以后为何没人愿意吃？

        流水素面，顾名思义就是和流水有关的面，大家都知道日本在饮食方面，和我们中国多多少少都会有一些差距。就好比如在日本很火的流水素面，它的结构基本上就是一个类似水渠的东西，有水从上面流下来。然后消费者分别在站在水渠的两侧，用筷子去夹从上游流下来的面。

       一般情况下这种面都是素面，流下来的过程中面也已经凉了，这时候两边的消费者就可以把面夹进自己装满调料的碗里面。然后就可以直接吃了，其实这根面不见得有多好吃，很多日本人也是这么说的。

       他们说这种吃面的方式之所以火，就是因为当地人都喜欢这种氛围，大家围在一起吃面有说有笑的。而且这种吃面的方式算是比较新颖的，日本人都喜欢去尝试一些新的东西，因此流水面出来之后，很快就红遍了日本。

       但是让人感到不能理解的是，火遍日本的流水素面，传入中国后却默默无闻，很多人别说去吃了， 就连听都还没有听说过。而且很多吃货还表示太脏了，这又是因为什么呢？

       很多人可能是这样想的，上游的人他筷子上肯定会沾一些口水，而在夹面的过程中。口水肯定就会沿着水渠流下来，这样下游的人就会接触到在上游人的口水。很多人都担心上面的人，万一他们的身体有什么问题，那下面的人岂不是很危险？

       其实在我看来这种担心完全就是多余的，为什么这么说呢？因为我们平时吃饭的时候，基本上也是一起在碗里面夹。这样对方的口水多多少少也会留在碗里，但即使这样我们吃了这么多年，也没发现有什么问题，所以说这完全是没有必要的担心。

       其次流水素面有水一直在冲下去，这样反而会干净很多。像我们那种传统的吃饭方式，里面的汤汁是一直静置。所以说里面会有更多的口水，不过话又说回来，一般情况下是没有问题的。毕竟我们这么吃了几百年，也不见发生什么事情。

       对此有的网友说：这种面也就是在日本才能火，主要是我想不通这样吃有什么好处，说氛围好像也没有什么氛围大家离的那么远。

是想写一篇数学建模论文。题目是“水渠流量的设计” 1某地计划修一条水渠，使水渠流量最大.

       本题主要目的是建立相关模型解决在修建水渠过程中的诸多问题，从而实现

       工程量最优化。

       针对问题一，为了求得开掘水渠的土石方量，本文通过对比分段三次

       Hermite 插值与三次样条插值，最终采用分段三次 Hermite 插值的方法对已知数 据点进行插值拟合，得到关于水渠的曲线方程 y f x ，对水渠曲线方程积分即

       得到水渠长度 dx yL ,利用 MATLAB 求解得到水渠长度为： m 5. 。

       因此最终解得开掘水渠的总土石方量为： 3 mLSV 。

       针对问题二，在问题一的基础上，本文建立积分上限函数模型：令 a ，

       1x 满足 1 2 1

       6

       x

       a

       V

       S y dt ，求出 i x 后， 1 ix 满足 1 2 1

       6

       i

       i

       x

       x

       V

       S y dt ，从而将

       总土石方量的六等分，得到 7 .x 1 ， 2 .x 2 ， x3 ， x4 ，

       x 5 ，进而确定了六等分点的坐标 y,x 。

       针对问题三，设在沿水渠的公路上有三个变量，分别为 k ji x ,x,x ，为使得

       运输工作量最小，本文建立了无约束规划模型，利用 MATLAB 求解得到最小运输

       量为 4 .7 m 。并给出了修建两条公路时水渠上的位置坐标 7.,B 和 4.C ， 。

       关键词：Hermite 插值 MATLAB 积分上限函数 无约束规划

       一、问题重述

       在某地区开掘水渠，已知该水渠经过的若干点。

       问题一，求解水渠施工的总石方量；

       问题二，如果将水渠的分成 6 截，每截土石方量相同，分段点应该取在何位

       置；

       问题三，设平行于水渠修一条路。河道中挖出的土石方要运往 A(,)

       处为了方便运输，计划在沿水渠的公路上选择两点修建通往 A 处的临时公路，使

       得总的土石方运输工作量最小。

       二、问题的分析

       针对问题一，本题要求开掘水渠的总土石方量，已知水渠截面积，则主要目

       的在于求得水渠长度。已知水渠经过的若干点的位置，要得到水渠的长度，本文

       想到用插值拟合可以得到水渠曲线，对曲线积分则得到水渠长度。插值与拟合的

       方法有多种，样条插值会较光滑，但不一定能保持原有形状，考虑到要更好的保

       持水渠的形状，于是，本文选用 Hermite 方法进行插值拟合。

       针对问题二， 要将水渠六等分且每段的土石方量相同，此问题为函数的反

       解问题，因此，在已知水渠曲线函数的情况下，本文可以考虑到用积分上限函数

       求解，从而确定 x 点，进而得到 y 点。

       针对问题三，要修建公路以运输土石方，从而使运输量工作量最小。此问题

       为规划问题，在问题二中，本文已知 x 与土石方量 V 存在关系，又因为运输工作

       量等于土石方量与距离的乘积，因此，本文使用无约束规划模型，求工作量最小

       值即可。

       三、模型假设

       1、修建的两条临时公路为直线。

       2、沿水渠的公路函数曲线近似与水渠的曲线函数相同。

       四、符号说明 xf 水渠曲线方程

       V 土石方量

       S 水渠截面积 L 水渠长度

       ix 水渠上点的横坐标

       iy 水渠上点的纵坐标

       iW 土石方运输工作量

       1L 临时公路 2L 临时公路

       五、模型的建立与求解

       5.1 问题一

       5.1.1 插值与拟合

       由已知水渠经过的点，做出散点图（图 1）

       0.8 0.9 1 1.1 1.2 1.3 1.4 x 4

       X/m

       Y/m

       水渠散点图

       图 1.水渠散点图

       方法 1、利用 Hermite 方法对已知数据点进行插值。 3 设 已 知 函 数 xfy 在 1 n 个互异节点 n x ,L,x,x 上 的 函 数 值 ii xfy n,L,1,0i 和导数值 i ' i ' x fy ，要求一个至多 2 n +1 次的多项 式 xH ，使得 i i yxH i ' i ' y xH n,1,0i Hermite 插值多项式为： 2 ' i i i i i i H x h x x a y y y

       其中，

       2

       n

       ij 0j j i

       j i x x xx h ， n ij 0j j i i x x 1 a 。

       利用 MATLAB 进行插值，得到插值曲线（图 2）。

       0.8 0.9 1 1.1 1.2 1.3 1.4 x 4

       Hermite插值曲线与原始数据点

       X/m

       Y/m

       Hermite插值曲线 原始数据点

       图 2.Hermite 插值曲线与原始数据点

       方法 2、利用样条差值对已知数据点进行插值。 3

       定义样条函数：

       数学上将具有一定光滑性的分段多项式称为样条函数。具体的说，给定区间 a,b 的一个划分

       0 1 1nn :a x x x x b

       如果函数 () sx满足： 1. 在每个小区间 1 , ( 0,1, , 1) ii x x i n 上 () sx是k 次多项式；

       2. () sx在 a,b 上具有 1 k 阶连续导数。

       则称 () sx为关于划分的k 次样条函数，其图形称为k 次样条曲线。 , , , n x x x 称

       为样条节点， 1 2 1 , , , n x x x 称为内节点， 0, n xx称为边界点，这样样条函数的全体

       记作 ( , ) p Sk ，称为k 次样条函数空间。

       显然，折线是一次样条曲线。

       若 ( ) ( , ) p s x S k ，则 () sx是关于分划的k 次多项式样条函数。k 次多项式

       样条函数的一般形式为

       1

        ( ) ( ) !! i kn j ki kj ij x s x x x ik

       其中 ( 0,1, , ) i ik 和 ( 1,2, , 1) j jn 均为任意常数，而

       ( ) , ( ) , ( 1,2, , 1) 0, k jjk j j x x x x x x j n xx 本文使用 3 k 的情况：即为三次样条函数。 三次样条函数：对于 a,b 上的划分 0 1 1nn :a x x x x b ，则

       1 2 3 3

       3 0 1

       1 ( ) ( ) ( ,3) 2! 3! 3! n j jp j aa s x x x x x x S

       其中

       3 3 ( ) , ( ) , ( 1,2, , 1) 0, jj j j x x x x x x j n xx

       三次样条函数差值：

       由于 3( ) ( ,3) ps x S中含有 3 n 个待定系数，故应需要 3 n 个插值条件，已

       知插值节点 i x 和相应的函数值 ( ) ( 0,1,2, , ) ii f x y i n ，这里提供了 1 n 个条件，

       还需要 2 个边界条件。

       常用的三次样条函数的边界条件有 3 中类型：

       （1） 3 0 3 ( ) , ( ) n s a y s b y 。由这中边界条件建立的样条插值函数称为 () fx的

       完备三次样条插值函数。

       特别的， 0'0 n yy 时，样条曲线在端点处呈水平状态。

       如果 () fx 不知道，可以要求 3() sx 与 () fx 在端点处近似相等。这时以

       0 1 2 3 , , , x x x x 为节点作一个三次 Newton 插值多项式 () a Nx，以 1 2 3 , , , n n n n x x x x 作一

       个三次 Newton 插值多项式 () b Nx，要求

       ( ) ( ), ( ) ( ) ab s a N a s b N b

       由这种边界条件建立的三次样条称为 () fx的 Lagrange 三次样条插值函数。

       （2） 3 0 3 3 ( ) , ( ) s a y s b y 。特别的 0 nn yy 时，称为自然边界条件。

       （3） 3 3 3 3 ( 0) ( 0), ( 0) ( 0) s a s b s a s b ，（这里要求 ( 0) ( 0) s a s b ）

       此条件称为周期条件。

       利用 MATLAB 进行三次样条插值，得到插值曲线（图 3）。

       0.8 0.9 1 1.1 1.2 1.3 1.4 x 4

       X/m

       Y/m

       三次样条插值曲线与原始数据点

       三次样条插值曲线 原始数据点

       图 3.三次样条插值曲线与原始数据点

       Hermite 插值与三次样条插值的对比

       5

       :

       SPLINE 提供的函数 s(x)的构建方法和 PCHIP 里面的函数 p(x)完全相同，只

       O x

       y

       0 AM

       1M

       2M

       1nM

       n BM

       图 4

       不过在 X(j)处的斜率的选择方法不一样，

       SPLINE 函数的 s(x)在 X(j)的二阶导数 D^2s(x)也是连续的，这导致了如下

       结果：

       （1） SPLINE 更加光滑，即，D^2s(x)是连续的。

       （2） 如果数据是一个光滑函数的值，则 SPLINE 更加精确。

       （3） 如果数据不是光滑的，则 PCHIP 没有 overshoots，也不太震荡。

       （4） PCHIP 建立的难度较小。

       （5） 这两种函数估计的难度是一样的。

       三次样条比 Hermite 插值光滑，样条的两阶导数连续，而 Hermite 插值一阶

       导数连续。不连续的两阶导数隐含着不连续的曲率。人的眼睛可以检测出图形上

       曲率的不连续。另一方面，Hermite 插值是保形状的，而样条插值不一定保形状。

       通过对比 Hermite 插值与三次样条插值，针对本题并无明显差异。为了更好

       的保证图形形状，减小误差，本文采用 Hermite 插值。

       5.1.2 求解水渠长度

       圆的周长可以利用圆的内接正多边形的周长当边数无限增多时的极限确定。

       类似的方法，可以用来建立平面连续曲线的弧长，应用定积分来计算弧长。

       设 AB 、 是曲线弧的两个端点。在弧 AB 上以此取分点：

       0 1 2 1 1 , , , , , , , , i i n n A M M M M M M M B ，并以此连接相邻分点得一折线（图

       4）。

       当分点的数目无限增加且每小段 1ii MM 都缩向一点时，如果此折线的长

       1

       1

       n

       ii

       i

       MM 的极限存在，则称此极限为曲线弧 AB 的弧长，并称此曲线弧 AB 是

       可求长的。

       由于光滑曲线弧是可求长的，故可应用定积分来计算弧长。

       设曲线弧由参数方程：

        () , () xt t yt

       给出，其中 ( ), ( ) tt 在 ,

       上具有连续导数，且 ( ) ( ) tt 、 不同时为零，现计

       算该曲线弧的长度。 取参数t为积分变量，它的变化区间为 ,

       。相应于 ,

       上任一小区间

        , t t dt 的小弧段的长度 s 近似等于对应的弦的长度 ( ) ( ) xy ，因为

       ( ) ( ) ( ) x t dt t dx t dt

       ( ) ( ) ( ) y t dt t dy t dt

       所以， s 的近似值（弧微分）即弧长元素为

       2 2 2 2 2 2 2 2 ( ) ( ) ( )( ) ( )( ) ( ) ( ) ds dx dy t dt t dt t t dt

       于是所求弧长为

        ( ) ( ) s t t dt

       当曲线弧由直角坐标方程

       ( ) ( ) y f x a x b

       给出，其中 () fx在 , ab上具有一阶连续导数，这时曲线弧由参数方程

       ()

       () xx

       a x b

       y f x

       从而所求的弧长为

       b a s y dx

       利用插值后得到的水渠的曲线函数，对其进行积分，则为水渠长度。 dx yL 用 MATLAB 求解得到 m 5.L

       5.1.3 求解土石方量

       已知，水渠长度，水渠截面积。

       则：

        2 mS

       3mLSV

       5.2 问题二

       设函数 () fx在区间 , ab上连续，并且设x为 , ab上的一点。观察 () fx在部 分区间 , ax上的定积分

       ()

       x

       a f x dx

       首先，由于 () fx在 , ax上依旧连续，因此该定积分存在。这里，x即表示

       定积分的上限，又表示积分变量。因为定积分与积分变量的记号无关，所以，为

       了明确起见，可以吧积分变量改用其他符号，例如用t表示，则上面的定积分可

       以写成

       ()

       x

       a f t dt

       如果上限x在区间 , ab上任意变动，则对于每一个取定的x值，定积分有一 个对应值，所以它在 , ab上定义了一个函数，记作 () x ：

       ( ) ( ) ( ) x a x f t dt a x b

       () x 便为积分上限函数。

       本文针对问题二建立积分上限函数模型：

       6

       x

       a

       V

       S y dt 通过起点a作为积分下限，求得第一个积分上限，即第一个等分点，第一个

       等分点为积分下限，求得第二个积分上限，即第二个等分点，以此类推。改变积

       分上下限，确定5个等分点，将水渠六等分，且每段土石方量相同。

       利用 MATLAB 求解（见附录 8.2），得到等分点坐标为： 4.,7. ， 2.,2. ， 6. ， ， 2., ， ,

       且每段的土石方量为： 3 m8.

       5.3 问题三 由问题二知土石方量 V 与水渠曲线函数存在关系。首先建立 xFV 模型。

       设在沿水渠的公路上有三个变量为 k ji x ,x,x ，修建的临时公路需要保证运输工

       作量最小，因此，在 D 点左边开掘水渠的土石方均运到 B 处，在 D 右边开掘水渠

       的土石方都运往 C 处。最终将土石方由 B、C 两处运往 A 处（示意图见图4）

       y = f(x)

       （xk）

       （xj）D

       （xi）

       L2

       L1

       A

       C

       B

       图4.水渠临时公路修建示意图

       运输工作量等于土石方量乘以距离，因此对于水渠曲线上的运输工作量本文

       建立的模型为：

       以 i 0 x ~x 段为例，设：该段水渠长度为

       n L L,L i 0 i0

       ，该段土石方量为

       n SL

       n V V,V i 0i0 i0

       ，

       则：

       )LnL(V)L2L(V)LL(VW i0i0i0i0

       )n(LVVLn i0

       n2 1n

       LV

       当 i 0 i 0 x x

       2

       x

       x

       2 i0 dx y1dxy1S 2 1 LV 2 1 W,n

       同理可得 1 kjkij W,W,W

       则：

       ijikjk1 WWWWW

        2 2 2 2 1 1 1 1 kk j j k k xx x x x x S y dx y dx S y dx y dx j i j i ii x x 2x x 2x 2x 2 dx y1dxy1S 2 1 dxy1dxy1S 2 1 对于由 B、C 两点运往 A 处的运输工作量本文建立的模型为：

        2 0k 2 0k x 0i 2 0i

       x

       2 2 y yxxdxy1Syyxxdxy1SW

       j

       j

        要使运输工作量最小，即 1 W 、 2 W 之和达到最小，因此，本文建立无约束规

       划模型： kji2kji1 x,x,xWx,x,xWMin

       即：

        2 0i 2 0i x 2 2x x 2 2x 2 y yxxdxy1dxy1 2 1 dxy1 2 1 Min j j i i

        2 0k 2 0k x 2 x 2 2 x x 2 y yxxdxy1dxy1 2 1 dxy1 2 1 jk k j

       运用 MATLAB 求解得：（程序代码见附录 8.3）

       最小运输工作量为： 4 .7 m B、C 两点的坐标为： 7.,B 和 4.C ， 。

       六、模型的优缺点

       优点：

       1、 通过对比 Hermite 插值与三次样条插值，发现求得的水渠长度分

       别为 .5m 和 .1m，对本题无明显差异。

       2、 对于运输量的规划问题，准确的反应了最优解。

       缺点：

       1、 对于插值函数的曲线积分，近似了曲线的导数，存在一定误差。

       2、 规划问题的运算量较大。利用 MATLAB 算法优势不明显。