辽宁省无资料地区设计暴雨洪水计算方法
的研究
辽宁省无资料地区设~1-
暴雨洪水~1-算75-法的研究
唐继业吴俊秀单丽
（辽宁省水文水资源勘测局）
江秋兰
（辽宁省水文水资源勘测局抚顺分局116000）
【摘要】本文针对辽宁省水工程设计中的实际情况，在认真总结经验的基础上，对流域特大暴雨重现期进行了探
讨;根据不同地区的产流特点，提出了分层扣损的饱卸产漉及非饱和流模型;建立了辽宁中部平厚区的三水”转
亿摸型;提出了综台经验单位线转换为瞬时单位线的流计算方法；在小
流域设计洪永计算上，建立了推理公式辽
宁击和概化过程发法.形成一垂适合辽宁特点的无资料地区设计暴雨
洪水计算方法?
【关键词】重现期模型单位巍
无资料地区暴雨洪水计算问题，一直是国内外水学科专
家学者在不断探索和研究的课题.《辽宁省中小河流（无资料地区）设计暴雨洪水计算方法》一书经过3年的工作编制完成?该书通过对大量水文气象资料分析?全面阐述了辽宁省暴雨，洪水时空变化规律，探人分析了暴雨洪水相关参数，提供出设计洪水计算的新理论，新方法和一系列新图件基础资料详实可靠，计算方法先进，综合成果符合部颁档计洪水计算规范》要求.
l基本资料与系列代表性分析
1.1基本资料
车成果分析暴雨资料的选用时段为最大10rai n,lh,6h,
24h,3d等5个时段.资料系列取自有资料以来截止到1995 年选用站数达306站,年限在25?9O年之间，共有12857 站年?系列最长的站是沈阳，大连，营口，均为91年,起讫时间为1905—1995 年.
1.2亲列代表性分析
首先从定性上开始，绘制各次实测大暴雨等值线图，了解气象成因与天
气系统组合；绘制3d,24h暴雨各站历年实测最高记录图；综合各次大暴雨等值线图，将历次笼罩范围用不同颜色标在一张图上；综合各河洪水调查成果，绘制历史大洪水笼罩范围用不同颜色标在一张图上；通过绘制长系列站正倒累进均值,,刀（一1）等分析图，识别丰枯水的交替变化规律，判别代表性较好的时段和站别?在此基础上,又做了定量分析：长短系列对比分析，典型站改正以及暴雨重现期怙算.
1.3暴雨重现期估算
此方法基本假定：雨洪长系列中特大，次大的个数是相应的;中小流域特大雨洪同次出现的机率基本一致；气候地形一致区的暴雨可以移置.基于以上假定后基本做法为：选取某测站以上流域内有特大暴雨和相应洪水的年份，根据洪水调查和历史文献的考证,先确定洪水的重现期，再作流域内均匀雨量点不考虑重现期的点雨量理论频率计算，求出各雨量点的均值及cv值的流域平均值?又在流域平均附近，假设一组cv 值，在频率格纸上点绘K 一N线,用各点求出的K值?在图上查出不同cv值相应该次暴雨各点值的重现期jv,将各个值相应的重现期相加求算术平均值，取与同攻洪水重现期相等或接近的值，则流域内暴雨中心点雨量的重现期即为所求.
2产流计算方法
在总结经验的基础上，应用填洼，下渗和水动力学原理，分析全省东部
湿润或半湿润，西部干旱或半干旱和中部平原区的产流特点，提出适应不同地区的产流模型.
2.1东部，北部饱和产流模型
辽宁省东部，北部地区属于湿润或半湿润地区，特别是
东部地区年降水量较大，植被好，土壤覆盖层的下渗能力较
大，在流域内的大部分土地上，当土壤缺水量尚未满足时，所
有降雨都被土壤所吸收，不宜超渗产流?这部分地区采用饱和产流模型，即辽宁饱和产流数学模型?蓄水容量曲线采用指数曲线型，产流面积公式为：
F
otz.O=l —e 一(1)
5
式中——产流面积，k;
F――流域面积.k;
P_——降雨量,mln;
?—一流域特征参数.
产流计算公式为：
,1,
R=P—1 一D}(1 —e){2), 口,
式中R ----- 径流量,rain;
o----- 流域韧始蓄水量,mm.
根据已有资料，优选出参数.蒸发采用双层扣损法. 2.2西部非饱和产流模型
辽宁省西部为干旱半干早地区，年降水量偏少，土壤下
43..
渗能力较差，降雨后产流与否与流域包气带缺水量无关.主要取决于雨强与下渗率的关系.根据霍顿下渗公式及其物理意义,结合分析该地区的产流特点，模拟出非饱和产流模型下渗率在流域上分布是不均匀的，可
用指数曲线来模拟下渗率在流域空间分布特征，产流面积公式为：=1 —e(3)
式中b——流域参数，它反映流域下渗率的分配特征
产一下渗率，mm/h.
时段径流深计算公式为：一
R=AP—f(1 一亨)zit(4)
f 一=f—t+f —z=()(,m—P
+
【一()㈤
式中A 一时段降雨量,mm;
, ---- 时段平均下渗率,mm/h;
i----- 时段平均雨强.mm/h;
L----- 流域最大蓄水量,mm;
P. --- 蓄水量,mm.
设计情况下，辽宁经过分析采用P.t=20ram,P.T二
40111/33
经验数据；
, -- 稳定下渗率,mm/h.
根据已知资料.优选出参数
2.3中部平原区”三水”转化产漉模型
辽宁省中部平原区是辽宁省心脏地区?该地区地势平
坦，浅层土壤由亚沙土和亚粘土组成，土壤人渗能力较强，地表径流较小?地下水开采较频繁，地下水运动以垂直交换为主,水平运动微弱?该地区设有台安实验站，主要对天然流域和潜水动态进行实验研究，从平原区”三水转化的物理概念和转化规律出发，针对辽宁省平原区水资源特点，在实验研究的基础上，研制了适应于该区的”三水转化模型?这一模型将包气带土壤计算层分为上下两层，两层各有土壤张力水的蓄水容量（田问持水量）WUM,fiLM作为各自土壤含水量
WUWL的上限.当流域的降水量到达地面，使上层土壤蓄满并满足蒸发后，剩余的降水量RC补给下层，它使下层土壤蓄满井满足土壤蒸发后剩余的水量RT称为总产水量?总产水量RT经过地下水库的调节，超出包气带重力水的蓄水能力的水量称为地表径流量，所剩余的水量RG
数值计算方法课程设计 姓名 学号 成绩
课程实际报告 实验一：秦九韶算法 题目 用选列主元高斯消去法解线性方程组 ???????=+- =-+-=-+-=--02 02 0 21 34343232121x x x x x x x x x x 算法语言： 利用c 语言的知识编写该算法程序 算法步骤叙述： 秦九昭算法的基思路是v[0]=a[0]*x+a[1] v[i]=v[i-1]*x+a[i+1];利用秦九昭算法计算多项式函数。 程序清单： #include
for(i=5;i>=1;i--) {sum=sum*x+a[i-1]; } printf("f(x)=%f/n",sum); } 输出结果计算：
实验总结： 通过运用C 语言，解决了秦九韶算法手写的复杂。为以后的雪地打下基础。 实验二：用选列主元高斯消去法解线性方程组 题目 用选列主元高斯消去法解线性方程组 ???????=+- =-+-=-+-=--02 0 2 0 21 34343232121x x x x x x x x x x 算法步骤叙述 第一步消元——在增广矩阵（A,b ）第一列中找到绝对值最大的元素，将其所在行与第一行交换，再对（A,b ）做初等行变换使原方程组的第一列元素除了第一行的全变为0； 第二步消元——在增广矩阵（A,b ）中第二列中（从第二行开始）找到绝对值最大的元素，将其所在行与第二行交换，再对（A,b ）做初等行变换使原方程组的第二列元素除了第一和第二行的全变为0； 第三步消元——在增广矩阵（A,b ）中第三列中（从第三行开始）找到绝对值最大的元素，将其所在行与第三行交换，再对（A,b ）做初等行变换使原方程组的第三列第四行元素为0； 第四，按x4-x3-x2-x1的顺序回代求解出方程组的解,x[n]=b[n]/a[n][n],x[i]=(b[i]-Σa[i][j]x[j])/a[i][i],i=n-1,…,2,1 程序清单： #include
3.用梯度法求下列无约束优化问题：MinF(X)=x12+4x22，设初始点取为X(0)={2,2}T，以梯度模为终止迭代准则，其收敛精度为5。 1)求初始点梯度▽F(X) ▽F(X)={2x1,8x2}T▽F(X(0))={4,16}T (2)第一次搜索 |▽F(X(0))|=16.5,S(0)=- ▽F(X(0))/16.5=-{0.243,0.97}T α(0)=2.157 X(1)=X(0)+α(0)S(0)={1.476,-0.923}T ▽F(x(1))={2.952,-0.738}T |▽F(x(1))|=3.043<5.0 故满足要求，停止迭代。 最优点X*={1.476,-0.0923}T 最优值F(X*)=2.21 4.
5.
6. 用外点法求解约束优化问题： ()()12211221min ..0()0 f X x x s t g X x x g X x =+=-≤=-≤ ， 收敛准则：(1) ()0.10.01k k X X εδ+-≤=，约束容限＝ 解：（1）利用外点法惩罚法构造无约束优化问题 () ( ) 12()22()212121(min ,()() k k k x x X r x x r x x r x +??Φ=?++-+-??可行域内）（可行域外） （2）此例只是为了说明外点法的思路，用微分法求解上述无约束优化问题。 用极值条件求解： 在可行域内：偏导数不可能等于0，即可行域内无极值 在可行域外，令： ()2()11211 ()2122 14()2012()0k k k r x x x r x x r x x x ?Φ =+-+=??Φ =--=?
第一章通风除尘与气力输送系统的设计 第一节概述 在食品加工厂中，车间的通风换气、设备和物料的冷却、粉尘的清除等都需要通风除尘系统来完成。粉状、颗粒状的物料(如奶粉、谷物等)的输送都可借助气力输送系统实现。通风除尘和气力输送系统是食品加工厂的常用装置。 食品加工厂中粉尘使空气污染，影响人的身体健康。灰尘还会加速设备的磨损，影响其寿命。灰尘在车间或排至厂房外，会污染周围的大气，影响环境卫生。由于粉尘的这些危害性，国家规定工厂中车间部空气的灰尘含量不得超过10mg/m3，排至室外的空气的灰尘含量不得超过150mg/m3，为了达到这个标准，必须装置有效的通风除尘设备。 图1是食品加工厂常见的通风除尘装置。主要由通风机、吸风罩、风管和除尘器等部分组成。当通风机工作时，由于负压的作用，外界空气从设备外壳的缝隙或专门的风管引入工作室，把设备工作时产生的粉尘、热量和水汽带走，经吸风罩沿风管送入除尘器净化，净化后的空气排出室外。 气力输送系统的形式与通风除尘系统相似，但其目的是输送物料，主要由接料器（供料器）、管道、卸料器、除尘器、风机等部分组成。气力输送系统除了起到输送作用外，还可以在输送过程中对物料进行清理、冷却、分级和对作业机完成除尘、降温等。小型面粉厂气力输送工艺流程如图2。
风机 气力输送具有设备简单、一次性投资低、可以一风多用等特点，与机械输送相比，气力输送的缺点主要是能耗较大，对颗粒物料易造成破碎。 通风除尘和气力输送都是利用空气的流动性能来进行空气的净化或物料的搬运的，因此，流体力学是本章的基础知识。有关流体力学的知识可参阅相关书籍资料，在此不再敷述。本章主要讨论食品加工厂通风除尘和气力输送系统的设计。 第二节通风除尘系统的设计与计算 1 通风除尘系统的设计原则和计算容 通风除尘系统也叫除尘网路或风网。通风除尘网路有单独风网和集中风网两种形式。在确定风网形式时，当： 1)吸出的含尘空气必须作单独处理； 2)吸风量要求准确且需经常调节； 3)需要风量较大；或设备本身自带通风机；
重庆邮电大学 数学与应用数学 专业 《数值计算方法》课程设计 姓名： 李金徽 王莹 刘姝楠 班级： 1131001 1131002 1131002 学号： 2010213542 2010213570 2010213571 设计时间： 2012-6-4 指导教师： 朱伟
一、课程设计目的 在科学计算与工程设计中，我们常会遇到求解线性方程组的问题，对于系数矩阵为低阶稠密矩阵的线性方程组，可以用直接法进行消元，而对于系数矩阵为大型稀疏矩阵的情况，直接法就显得比较繁琐，而迭代法比较适用。比较常用的迭代法有Jacobi 迭代与Gauss － seidel 迭代。本文基于两种方法设计算法，并比较他们的优劣。 二、课程设计内容 给出Jacobi 迭代法和Gauss-Seidel 迭代法求解线性方程组的算法思想和MATLAB 程序实现，并对比分析这两种算法的优劣。 三、问题的分析（含涉及的理论知识、算法等） Jacobi 迭代法 方程组迭代法的基本思想和求根的迭代法思想类似，即对于线性 方程组Ax = b( 其中n n n R b R R A ∈?∈,），即方程组 )1(2211222221211 1212111?? ???? ?=+?++??=+?++=+?++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 将系数矩阵A 写为 )2(000000 21122 12122 11U L D a a a a a a a a a A n n n n nn --≡??? ?? ? ? ??---- ??????? ??----??????? ??= 若选取D M =，则U L A M N +=-=，方程组)1(转化为等价方程组 b x U L Dx ++=)(
现代设计方法综合训练题 第二部分 非选择题 三、填空题 请在每小题的空格中填上正确答案。错填、不填均无分。1．计算机辅助设计（CAD）是指人们在计算机的辅助下，对产品或工程进行设计、绘图、分析计算或编写技术文件以及显示、输出的一种设计方法。 2．CAD系统的软件根据其用途可分为三类，它们是：软件、支撑软件专用软件。 3．在特征模型中，形状特征是其它特征的载体，非几何特征信息一般作为附加在形状特征上。 4．在单峰搜索区间[a,b]内，任取两个试算点a1,a2，若两点的函数值 F(a1)>F(a2)，则缩小后的区间为。 5．当有两个设计变量时，目标函数与设计变量之间的关系是中的一个曲面。 6．在有限元方法中，求总刚度矩阵的方法主要有两种，其中一种方法是利用刚度系数集成的方法获得总刚度矩阵的，该方法应用了 原理。 7．单元刚度矩阵具有对称性、性和奇异性。 8．2/3表决系统中各子系统的可靠度为R，则该系统的可靠度为 。 9．可靠度是对产品可靠性的度量。 10．某串联机电系统由N个子系统组成，各子系统的可靠度服从指数分布，且第i个子系统的失效率为λ，则该系统的平均寿命为。11．优化设计亦称最优化设计，它是以为基础，以电子计算机为辅助工具的一种设计方法。 12．优化设计的数学模型一般由设计变量、和约束条件三个基本要素组成。 13．国标GB3187-82将可靠性定义为：“产品在和规定的时间内，完成规定功能的能力。” 14．所谓特征指的是反映、可按一定原则分类的、具有相对独立意义的典型结构形状。 15．梯度是函数对各个设计变量的所组成的列矢量，并以符 号“”或grad表示。 16．二次插值法的基本思想是：在选定的单峰区间内取一点，连同两端点，利用这三点的函数值构成一个，作为原函数的近似，求出近似二次多项式的极小点作为原函数的近似最优点。
淮海工学院计算机工程学院课程设计报告书 课程名：《数值分析》 题目：数值分析课程设计 班级： 学号： 姓名：
数值分析课程设计 课程设计要求 1、研究第一导丝盘速度y与电流周波x的关系。 2、数据拟合问题运用样条差值方法求出温度变化的拟合曲线。 课程设计目的 1、通过编程加深对三次样条插值及曲线拟合的最小二乘法的理解； 2、学习用计算机解决工程问题，主要包括数据处理与分析。 课程设计环境 visual C++ 6.0 课程设计内容 课程设计题目1： 合成纤维抽丝工段中第一导丝盘的速度对丝的质量有很大的影响，第一丝盘的速度和电流周波有重要关系。下面是一组实例数据： 其中x代表电流周波，y代表第一导丝盘的速度 课程设计题目3： 在天气预报网站上获得你家乡所在城市当天24小时温度变化的数据，认真观察分析其变化趋势，在此基础上运用样条差值方法求出温度变化的拟合曲线。然后将该函数曲线打印出来并与原来的温度变化数据形成的曲线进行比较，给出结论。写出你研究的心得体会。 课程设计步骤 1、利用最小二乘法写出题1的公式和算法； 2、利用excel表格画出数据拟合后题1的图像； 3、在Visual C++ 6.0中编写出相应的代码； 4、搜索11月12日南通当地一天的温度变化数据； 5、在Visual C++ 6.0中编写出相应的代码； 6、利用excel表格画出数据拟合后题3的图像 课程设计结果 课程设计题目1 数值拟合
解：根据所给数据，在excel窗口运行： x=[49.2 50.0 49.3 49.0 49.0 49.5 49.8 49.9 50.2 50.2] y=[16.7 17.0 16.8 16.6 16.7 16.8 16.9 17.0 17.0 17.1] 课程设计题目3 数据为：X=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]; Y=[12,12,11,12,12,12,12,12,13,15,16,17,17,18,17,17,17,16,15,15,15,15,14,14]; 源代码为： 第一题： #include
《现代设计方法》课程习题集 西南科技大学成人、网络教育学院 版权所有 习题 【说明】：本课程《现代设计方法》（编号为09021）共有单选题,计算题,简答题, 填空题等多种试题类型，其中，本习题集中有[ 填空题,单选题]等试题类型未进入。 一、计算题 1. 用黄金分割法求解以下问题（缩小区间三次）。 342)(m in 2+-=x x x f ，给定初始区间[][]3,0,=b a ，取1.0=ε。 2. 用黄金分割法求解以下问题（缩小区间三次） 32)(m in 2+=x x f ，给定[][],1,2a b =-，取1.0=ε 3. 用黄金分割法求解以下问题（缩小区间三次） 432+=x )x (f min ，给定[][]40,b ,a =，取10.=ε。 4. 用黄金分割法求解以下问题（缩小区间三次）。 12)(m in 3+-=x x x f ，给定初始区间[][]3,0,=b a ，取5.0=ε 5. 用黄金分割法求解以下问题（缩小区间三次）。 107)(m in 2+-=x x x f ，给定初始区间[][]3,0,=b a ，取1.0=ε 6. 用梯度法求解无约束优化问题： 168)(m in 22221+-+=x x x X f ，取初始点[]T X 1,1)0(= ，计算精度1.0=ε。 7. 用梯度法求解96)(m in 12221+-+=x x x X f ，[]T X 1,1)0(= ，1.0=ε。 8. 用梯度法求解44)(m in 22221+-+=x x x X f ，[]T X 1,1)0(=，1.0=ε 。
9. 用梯度法求解无约束优化问题：1364)(m in 222 121+-+-=x x x x X f ，取初始点[]T X 1,1)0(=，计算精度1.0=ε。 10. 用梯度法求解1212221422)(m in x x x x x X f --+=,[]T X 1,1)0(=，1.0=ε 。（请迭代两次） 11. 有三个可靠度均为0.9的子系统组成的并联系统，试比较纯并联及2/3[G]表决系统的可靠度。 12. 一个由2个子系统组成的系统，其可靠度指标为0.85，试按等同分配法分配子系统的可靠度：（1）组成串联系统，（2）组成并联系统。 13. 已知某零件的应力和强度均呈正态分布，零件强度：MPa 516=δμ（均值），MPa S 2.24=δ（标准差），应力：MPa 378=σμ（均值），Mpa S 5.41=σ（标准差），试计算零件的可靠度与失效概率。 14. 由应力分析表明，某零件所承受的应力是拉应力，可用正态分布来描述，MPa T 3500=μ，标准差MPa S T 400=。该零件在制造过程中所引起的残余应力也可用正态分布来描述，其均值MPa C 1000=μ，标准差MPa S C 150=。由强度分析表明，该零件的强度也服从正态分布，其均值MPa 5000=δμ。现要求出当保证该零件的可靠度不低0.999时，零件强度的标准差的最低值应为多少? 15. 由应力分析表明，某零件所承受的应力是拉应力，可用正态分布来描述，MPa T 3500=μ，标准差MPa S T 400=。该零件在制造过程中所引起的残余应力也可用正态分布来描述，其均值MPa C 1000=μ，标准差MPa S C 150=。由强度分析表明，该零件的强度也服从正态分布，其均值MPa 5000=δμ。现要求出当保证该零件的可靠度不低0.999时，零件强度的标准差的最低值应为多少?
5 低压吸运气力输送系统设计计算示例 （1）单管气力输送系统设计计算示例 例7.3 如图7.78所示，由压榨车间将破碎饼粕送至浸出车间的气力输送系统。浸出车间日处理25 T/d （1）设计输送量G 计的确定 根据浸出车间要求处理饼25T/d ，按24h 计，则 G =25/24=1000（kg/h ） 由公式7-25，得： G 计=α×G =1.1×1000=1100（kg/h ） （2）输送风速V 的选择 由表7.56，取V 为21m/s 。 （3）输送浓度μ的选择 取μ=0.4。 （4）输送风量Q a 的确定 由公式7-27，得： 29924 .02.11100 =?= = μ ρa a G Q 计 （m 3/h ） （5）确定管径D 的确定 由公式7-28，得： 195.021 14.336002992 4.36004=???= = V Q D a π（m ） 取200mm 。则实际输送浓度为： 39.02378 2.11100=?==a a Q G ρμ计 （6）压力损失计算 输料输送压力损失H 物 ①空气通过作业机的压力损失H 机 由表7.1，H 机=0 ②接料器压力损失H 接 采用诱导式接料器，由表7.57，阻力系数为0.7。由公式7-31，得： g V H a j 22 ρζ=接 9.1881.92212.17.02 =???= （mmH 2O ） ③加速物料压力损失H 加 查表7.60得，i 谷粗=17mmH 2O/t ，由公式7-， H 加= i 谷粗G 算=17×1.1=18.7 （mmH 2O ） ④摩擦压力损失H 摩 查表7.65，R =2.21mmH 2O/m ，K 粗=0.669；由公式7-35，得： 236)39.0669.01(70.8421.2)1(=?+?=+=μm K RL H 摩（mmH 2O ） ⑤弯头压力损失H 弯 采用弯头90°，曲率半径为6D ，ζw 为0.083，查表7.60，K w =1.6，由公式7-45，得： 6.3)39.06.11(81 .92212.1083.0)1(22 2=?+???=+=μρζw a w K g V H 弯（mmH 2O ） ⑥恢复压力损失H 复 查表7.61和表7.62，△=0.35，β=1.5，由公式7-47，得： H 复=βΔΗ加=1.5×0.35×18.7=9.8 （mmH 2O ）
JAVA实现计算器课程设计 计算机科学系 计应0701班 指导老师:刘其昌 设计人员:陈秀桃 设计日期:2009年11月10日——2009年12月20日 计算器的设计 目录 第一章绪 论 ..................................................................... .. (2)
1.1 开发环 境 ..................................................................... . (2) 1.2 基本功能介 绍 ..................................................................... ......... 2 第二章系统设 计 ..................................................................... (3) 2.1 系统流程 图 ..................................................................... . (3) 2.2 系统功能框 图 ..................................................................... (3) 2.3 需求分 析 ..................................................................... ................ 4 第三章软件设 计 ..................................................................... (5) 3.1 界面设 计 ..................................................................... . (5) 3.2 代码设 计 .....................................................................
浙江省1月自学考试现代设计方法试题及答案
浙江省 1月自学考试现代设计方法试题 课程代码：09323 一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在题干的括号内。每小题1分，共30分) 1.对第Ⅱ象限中的一个点P 实施-?????? ? ???100010001坐标变换，则变换后P 点位于( ) A.Ⅰ象限 B.Ⅱ象限 C.Ⅲ象限 D.Ⅳ象限 2.滚筒式绘图仪上来自x 方向的脉冲信号使得( ) A.抬笔或落笔 B.笔架左右移动 C.滚筒带动图纸移动 D.停机和开机 3.图形变换矩阵T= 200020001??????? ???,则变换后的图形是原来的( ) A.2倍 B.1倍 C.3倍 D.4倍 4.扫描仪与主机接口，常见接口为( ) A.并行 B.RS-232 C.SCSI D.通信 5.用来定义用户在二维平面或三维世界中的物体，并符合右手定则的直角坐标系是( ) A.设备坐标系 B.世界坐标系
C.物理坐标系 D.规格化坐标系 6.CAD 支撑软件是在CAD 系统中，支撑( )进行CAD 工作的实用性功能软件。 A.用户 B.硬件 C.绘图机 D.打印机 7.对单凸物体进行面可见性测试，当外法线矢量N 和视线矢量S 的夹角( )时，面为可见的。 A.小于60° B.小于90° C.大于90° D.大于60° 8.数据库不但能为多个用户服务，而且能同时被多个用户使用，这种共享是( ) A.自动化的 B.并发的 C.智能化的 D.较大冗余的 9.标准件的图形构成分为4个层次，其中最基本的通用几何元素( ) A.G 类构件 B.K 类构件 C.B 类构件 D.A 类构件 10.在透视投影变换中，斜透视有( )个主灭点。 A.1 B.2 C.3 D.4 11.判断矩阵2014-???? ? ?，它应是( ) A.负定矩阵 B.正定矩阵 C.不定矩阵 D.对称矩阵
气力输送系统的设计要点 【摘要】本文简要介绍了气力输送系统的分类和组成，并对气力输送系统设计中存在的一些重要问题进行归纳总结，为以后的工程设计提供参考。 【关键词】气力输送；分类；组成；设计要点 0.前言 气力输送是借助负压或正压气流通过管道输送粉料的技术。与其他机械输送方式如斗提、皮带等相比，具有设备简单、布置灵活、占地面积小、操作及维修方便等特点，在钢铁、煤炭、电力、化工、粮食等行业得到广泛应用[1]。气力输送系统设计的合理与否，对输送效率、运行成本和使用寿命都有重要影响，因此本文对气力输送系统设计中着重考虑的问题进行归纳总结，希望引起工程设计同行的重视，为将来的工程设计提供参考。 1.气力输送系统 1.1气力输送的分类 根据输送管中物料的密集程度，气力输送可分为稀相输送和密相输送。稀相输送的混合比一般为0.1～25，输送气速为18～30m/s，高于浓相输送[2]。 根据输送管中气体的压力大小，气力输送可分为吸送式和压送式。吸送式的输送管内压力低于大气压，能自吸进料，缺点是必须负压卸料，而且物料输送距离较短；压送式的输送管内压力高于大气压，卸料方便，物料输送距离较长，其缺点是须用给料器将物料送入带压的管道中[3]。 1.2气力输送系统的组成 气力输送系统主要包括给料系统、输料系统、集料系统、动力系统和控制系统五大部分。 给料系统的作用是保证粉尘能够连续、均匀地进入输送管中，主要包括粉料缓冲斗、插板阀、旋转给料阀、给料器等。由于吸送式气力输送的输送管内存在一定负压，能够自吸进料，故其给料器通常采用L型或V型给料器，压送式的给料器较复杂，一般采用船型给料器或仓泵。 输料系统是粉料输送的关键环节，由输送直管、弯管、吸气口、吹扫口等组成，输送管的布置对气力输送系统的压力损失、连续稳定运行有至关重要的影响。 集料系统的作用是使料气分离，并将粉料收集后集中处理，主要包括集料器、卸料阀、粉料储罐等。集料器即除尘器，烟尘粒径小、混合比大时，应采用二级
数理学院2014级信息与计算科学 课程设计 姓名：刘金玉 学号： 3141301240 班级： 1402 成绩：
实验要求 1．应用自己熟悉的算法语言编写程序，使之尽可能具有通用性。2．上机前充分准备，复习有关算法，写出计算步骤，反复检查，调试程序。（注：在练习本上写，不上交） 3．完成计算后写出实验报告，内容包括:算法步骤叙述，变量说明，程序清单，输出计算结果，结构分析和小结等。（注：具体题目 具体分析，并不是所有的题目的实验报告都包含上述内容！）4．独立完成，如有雷同，一律判为零分！ 5．上机期间不允许做其他任何与课程设计无关的事情，否则被发现一次扣10分，被发现三次判为不及格！非特殊情况，不能请 假。旷课3个半天及以上者，直接判为不及格。
目录 一、基本技能训练 (4) 1、误差分析 (4) 2、求解非线性方程 (6) 3、插值 (12) 4、数值积分 (12) 二、提高技能训练 (16) 1、 (16) 2、 (18) 三、本课程设计的心得体会（500字左右） (21)
一、基本技能训练 1、误差分析 实验1.3 求一元二次方程的根 实验目的： 研究误差传播的原因与解决对策。 问题提出：求解一元二次方程20ax bx c ++= 实验内容： 一元二次方程的求根公式为 1,22b x a -+= 用求根公式求解下面两个方程： 2210(1)320(2)1010 x x x x +-=-+= 实验要求： （1） 考察单精度计算结果（与真解对比）； （2） 若计算结果与真解相差很大，分析其原因，提出新的算法（如先求1x 再 根据根与系数关系求2x ）以改进计算结果。 实验步骤： 方程（1）： 根据求根公式，写出程序： format long a=1;b=3;c=-2; x1=((-1)*b+sqrt(b^2-4*a*c))/2*a x2=((-1)*b-sqrt(b^2-4*a*c))/2*a
《数值分析》课程设计实验报告 龙格—库塔法分析Lorenz 方程 200820302033 胡涛 一、问题叙述 考虑著名的Lorenz 方程 () dx s y x dt dy rx y xz dt dz xy bz dt ?=-???=--???=-?? 其中s ，r ，b 为变化区域内有一定限制的实参数，该方程形式简单，表面上看并无惊人之处，但由该方程揭示出的许多现象，促使“混沌”成为数学研究的崭新领域，在实际应用中也产生了巨大的影响。 二、问题分析 Lorenz 方程实际上是一个四元一阶常微分方程，用解析法精确求解是不可能的，只能用数值计算，最主要的有欧拉法、亚当法和龙格- 库塔法等。为了得到较高精度的，我们采用经典四阶龙格—库塔方法求解该问题。 三、实验程序及注释 （1）算法程序 function [T]=Runge_Kutta(f,x0,y0,h,n) ％定义算法，其中f 为待解方程组， x0是初始自变量，y0是初始函数 值，h 是步长，n 为步数 if nargin<5 n=100; ％如果输入参数个数小于5，则步数 n=100 end r=size(y0);r=r(1); ％返回初始输出矩阵的行列数，并将 值赋给r(1) s=size(x0);s=s(1); ％返回初始输入矩阵的行列数，并 将值赋给s(1) r=r+s; T=zeros(r,n+1); T(:,1)=[y0;x0]; for t=2:n+1 ％以下是具体的求解过程 k1=feval(f,T(1:r-1,t-1)); k2=feval(f,[k1*(h/2)+T(1:r-1,t-1);x0+h/2]); k3=feval(f,[k2*(h/2)+T(1:r-1,t-1);x0+h/2]); k4=feval(f,[k3*h+T(1:r-1,t-1);x0+h]); x0=x0+h; T(:,t)=[T(1:r-1,t-1)+(k1+k2*2+k3*2+k4)*(h/6);x0]; end
全国2011年7月高等教育自学考试考前练习题 现代设计方法 （课程代码：02200） 一、单项选择题（本大题共20小题，每小题1分，共20分） 1．数字化仪在早期CAD 作业中经常使用，它是一种（ ） A ．输入设备 B ．输出设备 C ．存储设备 D ．绘图设备 2．滚筒式绘图仪一般采用插补原理进行绘图，但在绘图过程中有时不需要插补。绘制以下方向的线段时需要插补的是（ D ） A ．X 方向线段 B ．Y 方向线段 C ．±45o 方向线段 D ．一般方向线段 3．在采用Cohen —Sutherland 编码裁剪算法剪裁线段时，如果一条线段跨越两个区，则当其两端点编码的逻辑乘为零时，则该线段( ) A ．全部可见 B ．全部不可见 C ．部分可见 D ．不确定 4．在消隐处理中，当进行平面可见性测试时，若平面的法线矢量和视线矢量的夹角大于90o ，则该面（ ） A ．可见 B ．不可见 C ．部分可见 D ．不确定 5．一逐行扫描光栅显示器的水平扫描频率为15600H Z ，垂直扫描频率为50H Z ，则该显示器（ ） A ．水平方向有15600个像素点 B ．垂直方向有312条水平扫描线 C ．一条水平扫描线上有312个像素点 D ．垂直方向上有50条扫描线 6．F （x ）为单值、连续、可微且无约束的一元函数，则在点x=x * 处有极大值的充分条件是（ ） A ．0)(='* X F B ．0)(0)(>''='** X F X F C ．0)(=''*X F D ．0)(0 )(<''='** X F X F 7．外点罚函数法的特点是（ ） A ．适合处理等式约束的最优化问题 B ．探索过程在可行域内进行 C ．适合处理不等式约束的最优化问题 D ．需要先求出初始内点 8．在单峰搜索区间[a ，b]内任取两点a 1、a 2，若函数值F(a 1)=F(a 2)，则缩小后的区间为（ ） A ．[a ，a 1] B ．[a 1，b] C ．[a 1，a 2] D ．[a ，a 2] 9．以下关于函数的梯度的说法不正确的是（ ） A ．函数的梯度是标量 B ．函数值沿梯度方向变化最剧烈 C ．函数的梯度是矢量 D ．求函数的极小值时常沿负梯度方向搜索 10．多元函数F(x)在点X * 附近一阶偏导数连续，则该点为极大值点的充分条件为（ ） A ．F （x * ）=0 B ． F （x *）=0，H （x *）正定 C ．F （x * ）=0 D ． F （x * ）=0，H （x * ）负定 11．平面桁架结构中，某单元局部编码依次对应的总体编码为8，6，则单元则度矩阵中的元素 k 34应放入总体刚度矩阵[K]中的（ ） A ．16行15列 B ．11行12列 C ．12行16列 D ．11行15列 12．平面三解形单元的刚度矩阵阶数为（ ） A ．2×2 B ．3×3 C ．4×4 D ．6×6 13．平面刚架结构中的杆单元，单元中的位移分量个数为（ ） A ．3 B ．6 C ．9 D ．12 14．平面问题的弹性矩阵与材料的（ ） A ．弹性模量和硬度有关 B ．弹性模量和泊松比有关 C ．弹性模量和强度有关 D ．强度和泊松比有关 15．若把平面应力问题的弹性矩阵改成平面应变问题的弹性矩阵，只需将（ ）E /（12 μ-） μ- A ．E 换成E /（1μ-），μ换成μ/（1μ-） B ．E 换成E /（12 μ-），μ换成μ（12 μ-） C ．E 换成E /（12μ-），换成μ/（1μ-）
气力输送系统的设计原则与程序 在设计压送式气力输送装置时，首先必须要对被输送物料的性质和料粒形状，输送条件，现场状况等进行了解和研究，在此基础上充分发挥气力输送的优点，正确选择气力输送的类型，以利于提高生产效率。 一、设计原则 1、输送物料的性质和料粒形状物料的粒度常取平均粒度作为物料的计算粒度，并要了解物料粒度的分布情况。物料的流动性一般用堆积角和摩擦角的大小来间接表示。同一种物料由于含水量不同，流动性有很大的差别，对物料的含水量需考虑是内部水分还是表面水分，要考虑物料的粘附作用。 ●物料的密度和堆密度是直接影响气力输送装置的外形尺寸、结构形式及功率 消耗的大小。 ●物料破碎率决定气力输送的布置路线、输送距离和选定合适的气流速度。 ●物料的腐蚀性对输送管道的材质提出特殊的要求。 ●物料有静电效应时，要安装必要的地线和防止带电装置，防止产生静电。
●对爆炸性物料，除防止静电外，必须采取防爆安全措施。 ●对输送有害物料，必须考虑采取密闭的搬运安全措施，防止管道和设备磨损 或损坏而外泄。 2、输送量在压送式气力输送装置设计时，要根据单位时间的输送量来确定装置的容量及规格。气力输送装置往往是成套设备中的一部分，必须与其他主机及辅机匹配，如果在输送量的大小上发生矛盾，可以采取中间料斗贮存缓冲的办法予以解决。输送量还与工艺有关，根据工艺要求决定采用间歇式还是连续式的装置，在选用压送式气力输送形式还应考虑装置的可靠性，要估计气力输送一旦发生故障对生产的影响。 3、输送起点和终点的状况在保证工艺的前提下尽可能缩短输送距离，充分发挥压送式气力输送的优势。装置的安装高度和给料方式要允分考虑周围的环境，必须不阻碍交通，便于检修，并减少设备维护费用。 4、降噪及环保气源机械的噪声影响环境，在气源进口及出口处，必须采取降低噪声措施。如风机或空气压缩机安装在单独的房间内，采用消声器等。气力输送装置必须考虑排气的除尘效果，采用各种类型适合于气力输送特点的除尘器，防止对大气的污染，若采用湿法除尘器时，要考虑污水处理。 5、自动化水平程度气力输送装置可实现集中自动控制，由中央控制室进行远程控制。这不仅减少操作人员，而且实现自动连锁，防止事故发生。 6、安装要点气力输送装置安装在室外时要考虑防雨防冻措施。岔道、增压器、气动或电气控制元件、阀、限位开关等必须要有箱体，防止雨淋而失灵。 7、特殊条件的要求输送高温物料需考虑冷却因素，输送管道要考虑保温和加热。气源机械(如空压机)要考虑水冷条件及排水措施。
课 程 设 计 我再也回不到大二了， 大学是那么短暂 设计题目 数值分析 学生姓名 李飞吾 学 号 x x x x x x x x 专业班级 信息计x x x x x 班 指导教师 设 计 题 目 共15题如下 成绩
数值分析课程设计 1．1 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？（15621） 试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题 解：算法分析：解该问题主要使用递推算法，关于椰子数目的变化规律可以设起初的椰子数为0p ，第一至五次猴子在夜里藏椰子后，椰子的数目分别为01234,,,,p p p p p 再设最后每个人分得x 个椰子,由题: 14 (1)5 k k p p +=- （k=0,1,2,3,4）51(1)5 x p =- 所以551p x =+,11k k p p +=+利用逆向递推方法求解 15 1,4 k k p p +=+ （k=0,1,2,3,4） MATLAB 代码： n=input('n= '); n= 15621 for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p==fix(p), break end end disp([x,p]) 1．2 设，1 5n n x I dx x =+? （1）从0I 尽可能精确的近似值出发，利用递推公式： 11 5(1,2,20)n n I I n n -=-+= 计算机从1I 到20I 的近似值； （2）从30I 较粗糙的估计值出发，用递推公式：
1.试用进退法确定函数f（x）=的一维优化初始区间[a,b],给定初始 点=-1，初始步长h=1。 解：==-1，=f（）=9 =+h=-1+1=0，=f（）=5 比较，，由于>，作前进计算： =+2h=-1+2=1，=f（）=3 比较，，由于>，再作前进计算： =，==5 ==1，==3 =+4h=-1+4=3，=f（）=5 由于，可知初始区间已经找到，即[a,b]=[0,3]。 2.设某种单元的可靠度(t)=，其中=0.001/h，试求出： （1）由这种单元组成的二单元串联系统，二单元并联系统及2/3（G）表决系统的平均寿命； （2）当t=100h、500h、1000h时，一单元、二单元串联、二单元并联及2/3(G)表决系统的可靠度，并加以比较。 解：（1）一个单元与系统的平均寿命分别为： 单=1/ =1000h 串 =1/2=500h 并=3/2=1500h （） =5/6=833.3h （2）当t=100h时，一个单元与系统的可靠度分别为： 单==0.905 串 = 单 ==0.819 并=（ 单 ）=1-（）=0.991 = 单 - 单 =0.975 当t=500h时，一个单元与系统的可靠度分别为： 单==0.6065 串 = 单 ==0.3678 并=1-（ 单 ）=1-（）=0.8452 （）= 单 - 单 =0.6575 当t=1000h时，一个单元与系统的可靠度分别为： 单==0.368 串 = 单 ==0.135 并=1-（ 单 ）=1-（）=0.600= 单单 =0.306 从计算结果可以看出：
（1）一个单元的可靠度高于二单元串联系统的可靠度，但低于二单元并联系统的可靠度； （2）2/3（G）系统的平均寿命为一个单元的平均寿命的5/6 倍，明显低于一个单元的平均寿命。 3.已知约束优化问题：minf（x）=+ S．t．-1≥0 试写出内点罚函数与外点罚函数的表示式。 解：内点罚函数： Φ（）=++ 外点罚函数： 当-+1≤0，Φ（x，）=+ 当--+1＞0，Φ（x，）=++（） 4.现在要用钢板制作一个有盖的长方本储水箱，要求各边长均不超过20厘米， 且长度为宽度的2倍，试确定三边长度值，使该储水箱的容积最大，要求其表面积不超过400平方厘米。 解：（1）建立数学模型用复合形法迭代3次。 取储水箱长和高为设计变量，，则其宽0.5，数学模型为 maxF（X）=0.5 s．t．+3≤400 0≤≤20 0≤≤20 （2）用复合形法求解 求得的近似结果为=F(X*)=509 1已知右上图所示等腰直角三角形的单元刚度矩阵为： （）= 对 称 右图所示薄板结构中节点2处所受载荷以 及材料的弹性模量和板厚分别为： 100KN，E=2×N/cm，t=0.1cm 求节点2处的各位移分量。 解得.7×N,=-70.7×N {}=
数值计算方法课程设计 学号 班级
实验要求 1．应用自己熟悉的算法语言编写程序，使之尽可能具有通用性。2．上机前充分准备，复习有关算法，写出计算步骤，反复检查，调试程序。（注：在练习本上写，不上交） 3．完成计算后写出实验报告，容包括：所用的算法语言，CPU时间，算法步骤叙述，变量说明，程序清单，输出计算结果，结构分析和小结等。（注：具体题目具体分析，并不是所有的题目的实验报告都包含上述容！） 4．至少需要选择5道必做题目。其余的也可以选择，如果多选，可酌情加分！ 5．独立完成，如有雷同，一律判为零分！ 6．上机期间不允许做其他任何与课程设计无关的事情，否则被发现一次扣10分，被发现三次判为零分！
上机实习题目 1． 编写九韶算法程序，并用该程序计算多项式623)(3 5+-+=x x x x f 在1.3 1.2, ,1.1=x 的值。 public class Qinjiushao { public double result; public double x=-2;//定义一个未知数x public double[] b; public double[] c; public int i; public String abc; public void calculate(){ // abc="x^5+3*x^3-2*x+6";多项式，可以截取字符串获取系数 double[] a={2,0,-3,3,-4};//多项式的系数 double[] b=new double[a.length]; double[] c=new double[b.length]; for(i=0;i
相关主题
文本预览