数学建模初步_清华大学数学系

数学建模初步_清华大学数学系
数学建模初步_清华大学数学系

清华大学(数学建模获奖论文)

摘 要: 本文以深圳河流域特区为例,模拟了“分流清源”与“混流截排”收集机制的雨水管道与污水管道网络,对两种系统对各种参数的灵敏性和稳定性进行了分析,并对经济效果进行了估算和判断。本文设计了一种建设方案,得到在五年内建造污水管道网络的最佳策略,使得建设完成前排入环境的污水量最少,同时满足政府治污的“一、三、五年目标”,且最节约经济成本。 求解第一问“分流清源”建设方案: 我们根据资料中污水处理厂、污水产生源的位置通过“最小生成树”算法规划出“分流清源”模型中建设的管道主干,通过计算深圳水管网络的分形维数,得到整个系统长度共142公里,并根据资料估算工程费用约7.7亿元。 题目指出“后续管理困难而很难保证不会再出现污水管错接问题”,所以需要对该系统进行了稳定性评估。在已知各污水厂容纳量和污水源产生量和管道连接状况的情况下,通过规划得到了在部分管道错接或损坏的情况下,分别会导致多少污水溢出。结论是在一条管道损毁或错接的情况下,平均溢出增量百分比为6.6%,因此对于该最小生成树管道系统,除了一些关键路径需要额外措施保护之外,其他边稳定性良好。 对于“混流截排”方案: 题目提到“政府已有较大的投入到截排”,我们模拟设定雨水管道(连接污水源和初期雨水池)为三角形网格且已经建好,初期雨水池以及雨污管道已经建好,雨污管道仍然采用“最小生成树”。我们分析了无雨、中小雨、暴雨三种情况下初期雨水池接收雨水和污水的比例。通过“最小溢出固定流”的规划算法,得出各雨水池溢出混水量的最低值。我们通过改变降雨强度、个别雨池容量、排污比例、排污点的参数,来研究这些参数对溢出混水量最低值的影响。 因为环境污染和经济投入量纲不一致,我们用隶属函数和稀释理论两种不同的角度去评估清源和截排的优劣。根据文献中溢出混水的化学需氧量 (COD )就可以估算溢出污水的环境代价,它与污水处理费在五年内共14.1亿元,而清源方案的总代价是18.6亿元。另一方面,我们采取ω(x )=1-e ?(x λ)2的隶属函数来评价污染和花费,把两个拐点分别设置为政府预算20亿和现有污染排放量467万吨这两个位置。实施清源建造期间五年的污染及经济评价分别是0.045、0.35,然而截排则是0.39和0.11,这个评价指标越低越好,通过比较可以看出,清源虽然花费多,但是综合来看还是更占优势的。 求解污水治理的方案。由于建设改造“分流清源”水管系统需要五年的进程,在这五年中,清源与截排共存。我们提出方案:为保证居民秩序不受大规模影响,改造过程必须循序渐进,体现在总体上每年会有完成总长度一定百分比的要求(见第四章),同时对于局部地区也不能把一个节点两条及以上管道在一年内修建。在这种限制条件下,修建进程还必须满足政府治污的“一、三、五年目标”。我们采取搜索剪枝的做法有效解出了每年的修建方案,五年内分别修建水管道长。同时我们算出了在没有修建完成时,每年的污水排放情况并保证五年内它的总数是最低的。由此我们估算出在工程五年期间内环保和经济的总代价为17.1亿元。

什么是数学模型与数学建模

1. 什么是数学模型与数学建模 简单地说:数学模型就是对实际问题的一种数学表述。 具体一点说:数学模型是关于部分现实世界为某种目的的一个抽象的简化的数学结构。 更确切地说:数学模型就是对于一个特定的对象为了一个特定目标,根据特有的内在规律,做出一些必要的简化假设,运用适当的数学工具,得到的一个数学结构。数学结构可以是数学公式,算法、表格、图示等。 数学建模就是建立数学模型,建立数学模型的过程就是数学建模的过程(见数学建模过程流程图)。数学建模是一种数学的思考方法,是运用数学的语言和方法,通过抽象、简化建立能近似刻划并"解决"实际问题的一种强有力的数学手段。 2.美国大学生数学建模竞赛的由来: 1985年在美国出现了一种叫做MCM的一年一度大大学生数学模型(1987年全称为Mathematical Competition in Modeling,1988年改全称为Mathematical Contest in Modeling,其所写均为MCM)。这并不是偶然的。在1985年以前美国只有一种大学生数学竞赛(The william Lowell Putnam mathematial Competition,简称Putman(普特南)数学竞赛),这是由美国数学协会(MAA--即Mathematical Association of America的缩写)主持,于每年12月的第一个星期六分两试进行,每年一次。在国际上产生很大影响,现已成为国际性的大学生的一项著名赛事。该竞赛每年2月或3月进行。 我国自1989年首次参加这一竞赛,历届均取得优异成绩。经过数年参加美国赛表明,中国大学生在数学建模方面是有竞争力和创新联想能力的。为使这一赛事更广泛地展开,1990年先由中国工业与应用数学学会后与国家教委联合主办全国大学生数学建模竞赛(简称CMCM),该项赛事每年9月进行。

清华大学2006数学分析真题参考答案

清华大学2006数学分析真题参考答案 1.若数列{}n x 满足条件11221n n n n x x x x x x M ----+-++-≤g g g 则称{}n x 为有界变差数列,证:令10y =,11221n n n n n y x x x x x x ---=-+-++-g g g (n=2,3,….) 那么{}n y 单调递增,由条件知{}n y 有界, {}n y ∴收敛 ,从而0,0N ε?>?>,使当n m N >>时,有 n m y y ε-<,此即:11211n n n n m m x x x x x x ε---+--+-++-,考虑1()f x 和 3()f x 。 (i)若()132()()()f x f x f x <<,由于()f x 在12[,]x x 上连续,由介值定理,必存在 412[,]x x x ∈,使43()()f x f x =,定与一一映射矛盾。 (ii) ()312()()()f x f x f x <<,这时考虑23[,]x x ,必存在523[,]x x x ∈使得 51()()f x f x =,也得到矛盾。 (2)若存在123,,x x x I ∈且123x x x <<,123()()()f x f x f x ><。由介值定理,存在 412[,]x x x ∈,523[,]x x x ∈,使得42()()f x f x =,也与一一映射矛盾。 ∴f(x)在I 必严格单调。 3.证:设()f x 在(,)a b 内两个不同实根为12x x <,即12()()0f x f x ==。 由罗尔定理,存在12(,)c x x ∈,使()0f c '= (1) 因为()0f x ≥,从而为()f x 极小值点,由费马定理 12()()0f x f x ''∴== (2) 由(1),(2)对()f x '在1[,]x c 和2[,]c x 用罗尔定理,则存在3144(,),(,),x x c x c x ∈∈ 使34()()0f x f x ''''==。再一次对()f x ''在34[,]x x 上应用罗尔定理, 34[,](,)x x a b ξ?∈?,使(3)()0f ξ=。 4.证:令t=a+b-x,则 ()()()b b b a a a f x dx f a b t dt f a b x dx =+-=+-? ??。对6 a π = ,

清华大学数学建模竞赛题目

2011年清华大学“华罗庚杯”数学建模竞赛论文格式规范●论文(答卷)用白色A4纸,上下左右各留出2.5厘米的页边距。 ●论文第一页为承诺书,具体内容和格式见本规范第二页。 ●论文题目和摘要写在论文第二页上,从第三页开始是论文正文。 ●论文从第二页开始编写页码,页码必须位于每页页脚中部,用阿拉伯数字从 “1”开始连续编号。 ●论文不能有页眉,论文中不能有任何可能显示答题人身份的标志。 ●论文题目用3号黑体字、一级标题用4号黑体字,并居中。论文中其他汉字 一律采用小4号黑色宋体字,行距用单倍行距。 ●提请大家注意:摘要在整篇论文评阅中占有重要权重,请认真书写摘要(注 意篇幅不能超过一页)。 ●引用别人的成果或其他公开的资料(包括网上查到的资料) 必须按照规定的 参考文献的表述方式在正文引用处和参考文献中均明确列出。正文引用处用方括号标示参考文献的编号,如[1][3]等;引用书籍还必须指出页码。参考文献按正文中的引用次序列出,其中书籍的表述方式为: [编号] 作者,书名,出版地:出版社,出版年。 参考文献中期刊杂志论文的表述方式为: [编号] 作者,论文名,杂志名,卷期号:起止页码,出版年。 参考文献中网上资源的表述方式为: [编号] 作者,资源标题,网址,访问时间(年月日)。

2011年清华大学“华罗庚杯”数学建模竞赛 承诺书 我们仔细阅读了清华大学数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛的题目是:商务电梯运载效率的分析与建模 参赛队员(打印并签名) :1. 2. 3. 日期: 2011 年 5 月 16 日

数学建模常用方法

数学建模常用方法 建模常用算法,仅供参考: 1、蒙特卡罗算法(该算法又称随机性模拟算法,是通过计算机仿真来解决问题的算法,同时可以通过模拟可以来检验自己模型的正确性,是比赛时必 用的方法) 2、数据拟合、参数估计、插值等数据处理算法(比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用M a t l a b作为工具) 3、线性规划、整数规划、多元规划、二次规划等规划类问题(建模竞赛大多数问题属于最优化问题,很多时候这些问题可以用数学规划算法来描述,通 常使用L i n d o、L i n g o软件实现) 4、图论算法(这类算法可以分为很多种,包括最短路、网络流、二分图等算法,涉及到图论的问题可以用这些方法解决,需要认真准备) 5、动态规划、回溯搜索、分治算法、分支定界等计算机算法(这些算法是算法设计中比较常用的方法,很多场合可以用到竞赛中) 6、最优化理论的三大非经典算法:模拟退火法、神经网络、遗传算法(这些问题是用来解决一些较困难的最优化问题的算法,对于有些问题非常有帮助,但是算法的实现比较困难,需慎重使用) 7、网格算法和穷举法(网格算法和穷举法都是暴力搜索最优点的算法,在很多竞赛题中有应用,当重点讨论模型本身而轻视算法的时候,可以使用这种 暴力方案,最好使用一些高级语言作为编程工具) 8、一些连续离散化方法(很多问题都是实际来的,数据可以是连续的,而计 算机只认的是离散的数据,因此将其离散化后进行差分代替微分、求和代替积分等思想是非常重要的) 9、数值分析算法(如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比如方程组求解、矩阵运算、函数积分等算法就需要额外编写库函数进行调用) 10、图象处理算法(赛题中有一类问题与图形有关,即使与图形无关,论文 中也应该要不乏图片的,这些图形如何展示以及如何处理就是需要解决的问题,通常使用M a t l a b进行处理) 一、在数学建模中常用的方法: 1.类比法 2.二分法 3.量纲分析法 4.差分法 5.变分法 6.图论法 7.层次分析法 8.数据拟合法 9.回归分析法 10.数学规划(线性规划、非线性规划、整数规划、动态规划、目标规划) 11.机理分析 12.排队方法

清华大学数值分析A第一次作业

7、设y0=28,按递推公式 y n=y n?1? 1 100 783,n=1,2,… 计算y100,若取≈27.982,试问计算y100将有多大误差? 答:y100=y99?1 100783=y98?2 100 783=?=y0?100 100 783=28?783 若取783≈27.982,则y100≈28?27.982=0.018,只有2位有效数字,y100的最大误差位0.001 10、设f x=ln?(x? x2?1),它等价于f x=?ln?(x+ x2?1)。分别计算f30,开方和对数取6位有效数字。试问哪一个公式计算结果可靠?为什么? 答: x2?1≈29.9833 则对于f x=ln x?2?1,f30≈?4.09235 对于f x=?ln x+2?1,f30≈?4.09407 而f30= ln?(30?2?1) ,约为?4.09407,则f x=?ln?(x+ x2?1)计算结果更可靠。这是因为在公式f x=ln?(x? x2?1)中,存在两相近数相减(x? x2?1)的情况,导致算法数值不稳定。 11、求方程x2+62x+1=0的两个根,使它们具有四位有效数字。 答:x12=?62±622?4 2 =?31±312?1 则 x1=?31?312?1≈?31?30.98=?61.98 x2=?31+312?1= 1 31+312?1 ≈? 1 ≈?0.01613

12.(1)、计算101.1?101,要求具有4位有效数字 答:101.1?101= 101.1+101≈0.1 10.05+10.05 ≈0.004975 14、试导出计算积分I n=x n 4x+1dx 1 的一个递推公式,并讨论所得公式是否计算稳定。 答:I n=x n 4x+1dx 1 0= 1 4 4x+1x n?1?1 4 x n?1 4x+1 dx= 1 1 4 x n?1 1 dx?1 4 x n?1 4x+1 dx 1 = 1 4n ? 1 4 I n?1,n=1,2… I0= 1 dx= ln5 1 记εn为I n的误差,则由递推公式可得 εn=?1 εn?1=?=(? 1 )nε0 当n增大时,εn是减小的,故递推公式是计算稳定的。

系统的描述与数学建模

系统的描述与数学建模 [摘要]数学建模就是利用数学方法将系统的文字语言描述转化成数学方式表达。由于影响系统的因素多种多样,当用数学表达系统时,我们要求尽可能要使得数学建模既能从本质上反映系统,又能使得系统的数学模型具有简单性。 [关键词]系统的建模数学建模 数学建模就是利用数学方法将系统的文字语言描述转化成数学方式表达。由于影响系统的因素多种多样,当用数学表达系统时,我们要求尽可能要使得数学建模既能从本质上反映系统,又能使得系统的数学模型具有简单性。一个极其复杂的数学模型对于分析系统毫无帮助。 为了说明这种数学建模的方法,我们举一个简单的例子。比如我们研究某一地区人口的健康状况。假定在我们的研究时段内没有人口的自然死亡,按照自然规律人口总是以一定的概率,变成亚健康、或者患上某种轻疾病、或者患上重疾病。在一定的环境和医疗条件下,部分亚健康者和患者会得以康复,这是一种简单运算的系统描述,并没有具体地给出定量表达。为了能用数学的方法表达这个描述,我们按照以下方式将人口分类:(1)健康人。(2)亚健康人。(3)患轻病人。(4)患重病人。 根据上面的关系和一些假定条件,我们可以得到相应的微分方程,至于方程的详细导出我们以后再讨论。这里我们需要指出,前面我们只是一种说明性的举例,在实际建模过程中,要依赖于系统所在的环境,按照前面方法得到的应是确定性模型,在随机环境中,上面所述的量应当对应成相应状态的概率。 再比如排队系统,是最常见的一种系统,这类系统主要描述顾客到达,接受服务然后离开这一过程。系统由顾客与服务员两个单元组成。这类问题主要由以下四个因素决定:(1)顾客来到窗口的频率。(2)窗口的个数。(3)排队规则。(4)服务时间分布;所以我们必须对它们作适当的假定。 在单个服务台的排队系统模型M/M/1,即系统只设一个服务台床的情况。假定顾客是相互独立地到达系统,而且顾客到达系统的间隔时间服从负指数分布 F(t)=1-e -λt (输入过程),又服务窗为每一位顾客的服务时间也同时服从负指 数分布H(t)=1-e -μt (运行方式)。对这种最简单的排队模型,我们将依照不同的系统规则确定排队系统所满足的微分方程。 M/M/1损失制排队模型是指系统内只设一个服务窗,系统容量为1(即有一个排队位置而无排队等待位置),顾客到达和窗口服务时间均为负指数分布,且

清华大学数学课介绍

数学科学系 00420033数学模型3学分48学时 Mathematical Modelling 建立数学模型是用数学方法解决实际问题的关键步骤。本课程从日常生活的有趣问题入手,介绍数学模型的一般概念、方法和步骤,通过实例研究介绍一些用机理分析方法建立的非物理领域的模型及常用的建模数学方法,培养同学用建模方法分析和解决实际问题的意识和能力。 00420152数学建模引论2学分32学时 Introduction of Mathematical Modelling 本课程以案例分析的方式组织教学,主要面向低年级的学生,各个学期根据对学生数学基础的不同要求,选择案例。我们这里所选择的都是实际应用价值非常突出的案例。 00420163数理科学与人文3学分48学时 Mathematical and Physical Sciences and Humanities 本课程旨在加强学生以通识教育为目标的思维和训练,提高学生的科学素质。该课程虽然以知识为载体,却并不以传授理论知识为主要目的,而是以启迪思想,养成思考的习惯,以提升学生的创新意识。 00420183博弈论3学分48学时 Game Theory This is an introductory course on the basic concepts of Game Theory. Topics to be covered are:Combinatorial Game Theory, Games in Extensive Form, 2-person 0-sum games, Bimatrix games, Nash Equilibrium, Correlated equilibrium, Evolutionary Game Theory, Repeated Prisoner’s Dilemma, Bargaining Problems, Games in Coalition form, Shapley value, Nucleolus, 2-side matching problem. 10420095微积分(1)5学分80学时 Calculus(1) 内容包括:实数,函数,极限论,连续函数,导数与微分,微分中值定理,L'Hospital法则,极值与凸性,Taylor公式,不定积分与定积分,广义积分,积分应用,数项级数,函数级数,幂级数,Fourier级数。 10420115微积分(2)5学分80学时 Calculus(2) n维空间中的距离、邻域、开集与闭集,多元函数的极限与连续,多元函数微分学,空间曲线与曲面,重

清华大学电力系统 同步发电机的数学模型21

长江三峡水电枢纽

同步汽轮发电机的转子同步水轮发电机的转子气隙 定子 同步发电机的FLASH.SWF 11

定子上3个等效绕组 a 相绕组 b 相绕组 c 相绕组 转子上3个等效绕组 同步发电机简化为:定子3个绕组、转子3个绕组、气隙、定子铁心、转子铁心组成的6绕组电磁系统励磁绕组 d 轴等效的阻尼绕组轴等效的阻尼绕组Q 15d 轴 q 轴120度 120度 120度 定子、转子铁心同轴(忽略定、转θ sin )M F =磁动势零点 θ 的,无饱和,无磁滞和涡流损耗,

19 磁链与电流、电压的参考正方向 1、设转子逆时针旋转为旋转正方向; 3、定子三相绕组端电压的极性与相电流正方向按发电机惯例来定义,即 正值电流i a 从端电压u a 的正极流出发电机,b 、c 相类似。 定子绕组的正电流产生负的磁链!! 2、定子三相绕组磁链ψa ,ψb ,ψc 的正方向与a 、b 、c 三轴正方向一致; + -21 5、d轴上的励磁绕组f、阻磁链正方向与d轴磁链正方向与q轴的正方向一致;正电流由端电压,因此绕组电阻: a 相绕组 b 相绕组 c 相绕组 +

26 励磁绕组d 轴阻尼绕组 轴阻尼绕组 绕组、 28 绕组的磁链方程-6个 定子绕组的磁链a 相绕组的磁链= a 相绕组电流产生的自磁链+ b 相绕组电流产生的互磁链+ c 相绕组电流产生的互磁链+励磁绕组电流产生的互磁链+D 绕组电流产生的互磁链 + Q 绕组电流产生的互磁链

31 转子绕组的磁链励磁绕组的磁链= a 相绕组电流产生的互磁链+ b 相绕组电流产生的互磁链+ c 相绕组电流产生的互磁链+励磁绕组电流产生的自磁链+D 绕组电流产生的互磁链+ Q 绕组电流产生的互磁链 36 a 相绕组磁路磁阻(磁导)的变化与转子d 轴与a 相绕组轴线的夹角θa (=ωt )有关 磁路的磁导λaa ,自感L aa 为θa 的周 期函数,周期为π。 θa θa =±π/2 磁路磁导最小,自感最小 a θa =0,π磁路磁导最大,自感最大 a

中国大学生数学建模竞赛历年试题

中国大学生数学建模竞赛(CUMCM)历年赛题一览! CUMCM历年赛题一览!! CUMCM从1992年到2007年的16年中共出了45个题目,供大家浏览 1992年A)施肥效果分析问题(北京理工大学:叶其孝) (B)实验数据分解问题(复旦大学:谭永基) 1993年A)非线性交调的频率设计问题(北京大学:谢衷洁) (B)足球排名次问题(清华大学:蔡大用) 1994年A)逢山开路问题(西安电子科技大学:何大可) (B)锁具装箱问题(复旦大学:谭永基,华东理工大学:俞文此) 1995年:(A)飞行管理问题(复旦大学:谭永基,华东理工大学:俞文此) (B)天车与冶炼炉的作业调度问题(浙江大学:刘祥官,李吉鸾) 1996年:(A)最优捕鱼策略问题(北京师范大学:刘来福) (B)节水洗衣机问题(重庆大学:付鹂) 1997年:(A)零件参数设计问题(清华大学:姜启源) (B)截断切割问题(复旦大学:谭永基,华东理工大学:俞文此) 1998年:(A)投资的收益和风险问题(浙江大学:陈淑平) (B)灾情巡视路线问题(上海海运学院:丁颂康) 1999年:(A)自动化车床管理问题(北京大学:孙山泽) (B)钻井布局问题(郑州大学:林诒勋) (C)煤矸石堆积问题(太原理工大学:贾晓峰) (D)钻井布局问题(郑州大学:林诒勋) 2000年:(A)DNA序列分类问题(北京工业大学:孟大志) (B)钢管订购和运输问题(武汉大学:费甫生) (C)飞越北极问题(复旦大学:谭永基) (D)空洞探测问题(东北电力学院:关信) 2001年:(A)血管的三维重建问题(浙江大学:汪国昭) (B)公交车调度问题(清华大学:谭泽光) (C)基金使用计划问题(东南大学:陈恩水) (D)公交车调度问题(清华大学:谭泽光) 2002年:(A)车灯线光源的优化设计问题(复旦大学:谭永基,华东理工大学:俞文此) (B)彩票中的数学问题(解放军信息工程大学:韩中庚) (C)车灯线光源的优化设计问题(复旦大学:谭永基,华东理工大学:俞文此))

数学建模 自习室管理系统

一.问题重述: 近年来,大学用电浪费比较严重,集中体现在学生上晚自习上,一种情况是去某个教室上自习的人比较少,但是教室的灯却全部打开,第二种情况是晚上上自习的总人数比较少,但是开放的教室比较多,这要求提供一种最节约、最合理的管理方法。根据题目所给出的数据,有以下问题。数据见表。 1.假如学校有8000名同学,每个同学是否上自习相互独立,上自习的可能性为0.7. 要使需要上自习的同学满足程度不低于95%,开放的教室满座率不低于4/5,同时尽量不超过90%。问该安排哪些教室开放,能达到节约用电的目的。 2.在第一问基础上,假设这8000名同学分别住在10个宿舍区,现有的45个教室分为9个自习区,按顺序5个教室为1个区,即1,2,3,4,5为第1区,…, 41,42,43,44,45为第9区。这10个宿舍区到9个自习区的距离见表2。学生到各教室上自习的满意程度与到该教室的距离有关系,距离近则满意程度高,距离远则满意程度降低。假设学生从宿舍区到一个自习区的距离与到自习区任何教室的距离相同。请给出合理的满意程度的度量,并重新考虑如何安排教室,既达到节约用电目的,又能提高学生的满意程度。另外尽量安排开放同区的教室。3.假设临近期末,上自习的人数突然增多,每个同学上自习的可能性增大为0.85,要使需要上自习的同学满足程度不低于99%,开放的教室满座率不低于4/5,同时尽量不超过95%。这时可能出现教室不能满足需要,需要临时搭建几个教室。 假设现有的45个教室仍按问题2中要求分为9个区。搭建的教室紧靠在某区,每个区只能搭建一个教室,搭建的教室与该区某教室的规格相同(所有参数相同),学生到该教室的距离与到该区任何教室的距离假设相同。问至少要搭建几个教室,并搭建在什么位置,既达到节约用电目的,又能提高学生的满意程度。

2014年清华大学数学建模竞赛

2014年清华大学“华罗庚杯”数学建模竞赛论文格式规范●论文(答卷)用白色A4纸,上下左右各留出2.5厘米的页边距。 ●论文第一页为承诺书,具体内容和格式见本规范第二页。 ●论文题目和摘要写在论文第二页上,从第三页开始是论文正文。 ●论文从第二页开始编写页码,页码必须位于每页页脚中部,用阿拉伯数字从 “1”开始连续编号。 ●论文不能有页眉,论文中不能有任何可能显示答题人身份的标志。 ●论文题目用3号黑体字、一级标题用4号黑体字,并居中。论文中其他汉字 一律采用小4号黑色宋体字,行距用单倍行距。 ●提请大家注意:摘要在整篇论文评阅中占有重要权重,请认真书写摘要(注 意篇幅不能超过一页)。 ●引用别人的成果或其他公开的资料(包括网上查到的资料) 必须按照规定的 参考文献的表述方式在正文引用处和参考文献中均明确列出。正文引用处用方括号标示参考文献的编号,如[1][3]等;引用书籍还必须指出页码。参考文献按正文中的引用次序列出,其中书籍的表述方式为: [编号] 作者,书名,出版地:出版社,出版年。 参考文献中期刊杂志论文的表述方式为: [编号] 作者,论文名,杂志名,卷期号:起止页码,出版年。 参考文献中网上资源的表述方式为: [编号] 作者,资源标题,网址,访问时间(年月日)。

2014年清华大学“华罗庚杯”数学建模竞赛 承诺书 我们仔细阅读了清华大学数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛的题目是: 参赛队员(打印并签名) :1. 2. 3. 日期:年月日

清华大学杨顶辉数值分析第6次作业

清华大学杨顶辉数值分析第6次作业

9.令*()(21),[0,1]n n T x T x x =-∈,试证*{()}n T x 是在[0,1]上带权 2 ()x x x ρ= -****0123(),(),(),()T x T x T x T x . 证明: 1 1 **2 1 1 * *20 12 2 1**20 ()()()(21)(21)211()()()()()211()22 ()()1()1()()()()()1n m n m n m n m n m n n m n m x T x T x dx x T x dx x x t x x T x T x dx t T t dt t t t T t dt t T x x x T x T x dx t T t t ρρρ---=---=-=++-= --= -???? ?令,则 由切比雪夫多项式1 01=02 m n dt m n m n ππ ≠??? =≠??==??? 所以*{()}n T x 是在[0,1]上带权2 ()x x x ρ= - *00*11* 2 2 2 2*33233()(21)1()(21)21 ()(21)2(21)188()(21)4(21)3(21)3248181 T x T x T x T x x T x T x x x x T x T x x x x x x =-==-=-=-=--=-=-=---=-+- 14.已知实验数据如下: i x 19 25 31 38 44 i y 19.0 32.3 49.0 73.3 97.8 用最小二乘法求形如2y a bx =+的经验公式,并求均方误差 解: 法方程为

办公室电话系统模拟(数学建模)

排队论在电话问题中的应用 摘要 本文建立一个模拟办公室电话系统模型,解决由三个电话机占线而可能打不进电话的问题。根据该办公室的电话系统状况得知其服从排队论模型规律,则应用排队论知识建立模型。 用)(t Pn 表示在时刻t ,服务系统的状态为n (系统占线条数为n )的概率。通过输入过程(顾客打进电话),排队规则,和服务机构的具体情况建立关于)(t Pn 的微分差分方程求解。令0)('=t P n 把微分方程变成差分方程,而不再含微分了, 把)(t Pn 转化为与t 无关的稳态解。关于标准的M/M/s 排队模型各种特征的规定于标准的M/M/1模型的规定相同。另外规定各服务器工作是相互独立(不搞协作)且平均服务率相同 .==...==s 21μμμμ于是整个服务机构的平均服务率为μs 。令ρ=λ/su 只有当时λ/su<1时才不会排成无限的队列,成这个系统为服务强度,各顾客服务时间服从相同的负指数分布 ' 通过模型我们可以得到:无占线、一条占线、两条占线、三条占线的概率分别 是%,%,%,%。 · 关键词:泊松分布,指数分布,概率,期望,Little 公式

… 一、问题重述 一个办公室有三条电话线可打进,也就是说在任意时刻最多能接待三个顾客,顾客打电话是随机的,其时间服从上午9点至下午5点的均匀分布,每次电话持续时间是均值为6分钟的随机变量。 经理关心由于三个电话机占线而可能打不进电话的顾客数。他们当中部分人稍后可能重拨电话,而其他人则可能放弃通话,一天中接通的电话平均数是70。 请你建立一个模型模拟办公室电话系统,帮助经理在休息时思考这个问题,用你的模型做下述估计: (1)} (2)无电话占线、有一条、两条占线和三条都占线的时间百分比; (3)未打进电话的顾客所占百分比。 二、问题的分析 这是一个多服务台混合制模型M/M/s/K,顾客的相继到达时间服从参数为的负指数分布(即顾客的到达过程为Poisson流),服务台的个数为s,每个服务台的服务时间相互独立,且服从参数为的负指数分布,系统的空间为K。求平稳分布,考虑系统处的任一状态n。假设记录了一段时间内系统进入状态n和离开状态n的次数,则因为“进入”和“离开”是交替发生的,所以这两个数要么相等要么相差1。但就这两件事件平均发生率来说,可以认为是相等的。 三、基本假设 ①顾客的相继到达时间服从参数为λ的负指数分布; ②服务时间服从参数μ的负指数分布; ③顾客选择打进哪一条线是随机的而且是等可能的; ④, ⑤某条线接通时,其他顾客不能接通,则称为占线 四、符号定义及变量说明 ①:顾客的相继到达时间服从参数为λ的负指数分布,服务时间服从参数μ的负指 数分布; ②:) Pn表示在时刻t服务系统的状态为n(系统中顾客数为n)的概率,(t

数学建模A题系泊系统设计完整版

数学建模A题系泊系统 设计 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

系泊系统的设计 摘要 本题要求观测近海观测网的组成,建立模型对其中系泊系统进行设计,在不同风速和水流的情况下确定锚链,重物球,钢管及浮标等的状态,从而使通讯设备的工作效果最佳。求解的具体流程如下: 针对问题一,分别对系统中的受力物体在水平方向和竖直方向上的力进行分析,找出锚链对锚无拉力时的临界风速,运用力矩平衡求出钢管与钢桶的倾斜角度。对于锚链,将其等效为悬链线模型,根据风速不同判断锚链的状态,从而求出结果。 ?时能够正常工作。为针对问题二,需要调节重物球的质量,使通讯设备在36m m 了确定重物球的质量,首先将实际风速与临界风速进行比较,判断此时系统中各物体的状态,与题目中已知数据进行比较。在钢桶倾斜角度达到临界角度时,计算锚链与海床的夹角并于题中数据进行比较,计算重物球的质量。在浮标完全没入海面时,计算相应条件下重物球的质量,从而确定满足条件的重物球的质量范围。 针对问题三,要求在不同条件下,求出系泊系统中各物体的状态。以型号I锚链为例,当水流方向与风速方向相同时,系统条件最差,分析在不同水深条件下的系泊系统设计。由题中已知条件确定系统设计的限制条件,对系统各物体进行受力分析,以使整体结果最小,即可得出最优的系泊系统设计。 关键词:悬链线多目标非线性规划 一、问题重述 近浅海观测网的传输节点由浮标系统、系泊系统和水声通讯系统组成(如图1所示)。某型传输节点的浮标系统可简化为底面直径2m、高2m的圆柱体,浮标的质量为1000kg。系泊系统由钢管、钢桶、重物球、电焊锚链和特制的抗拖移锚组成。锚的质量为600kg,锚链选用无档普通链环,近浅海观测网的常用型号及其参数在附表中列出。钢管共4节,每节长度1m,直径为50mm,每节钢管的质量为10kg。要求锚链末端与锚的链接处的切线方向与海床的夹角不超过16度,否则锚会被拖行,致使节点移位丢失。水声通讯系统安装在一个长1m、外径30cm的密封圆柱形钢桶内,设备和钢桶总质量为100kg。钢桶上接第4节钢管,下接电焊锚链。钢桶竖直时,水声通讯设备的工作效果最佳。若钢桶倾斜,则影响设备的工作效果。钢桶的倾斜角度(钢桶与竖直线的夹角)超过5度时,设备的工作效果较差。为了控制钢桶的倾斜角度,钢桶与电焊锚链链接处可悬挂重物球。 系泊系统的设计问题就是确定锚链的型号、长度和重物球的质量,使得浮标的吃水深度和游动区域及钢桶的倾斜角度尽可能小。

清华大学杨顶辉数值分析第6次作业

9.令*()(21),[0,1]n n T x T x x =-∈,试证*{()}n T x 是在[0,1] 上带权()x ρ=的正交多项式,并求****0123(),(),(),()T x T x T x T x . 证明: 1 1 * *0 1 1 * *011**0 ()()()(21)(21)211()()()()()2()()()()()()()()n m n m n m n m n m n n m n m x T x T x dx x T x dx t x x T x T x dx t T t dt t T t dt T x x T x T x dx t T t ρρρ---=--=-== = ???? ?令,则 由切比雪夫多项式1 01=02 m n dt m n m n ππ ≠??? =≠??==??? 所以*{()}n T x 是在[0,1] 上带权()x ρ= *00*11* 22 2 2*33233()(21)1()(21)21 ()(21)2(21)188()(21)4(21)3(21)3248181 T x T x T x T x x T x T x x x x T x T x x x x x x =-==-=-=-=--=-=-=---=-+- 14.已知实验数据如下: 用最小二乘法求形如2y a bx =+的经验公式,并求均方误差 解: 法方程为

22222(1,)(1,1)(1,)(,)(,1)(,)a y x b x y x x x ?????? =???? ?????? ?? 即 5 5327271.453277277699369321.5a b ??????=???????????? 解得 0.972579 0.050035a b =?? =? 拟合公式为20.9725790.050035y x =+ 均方误差 2 4 2 2 0[]0.015023i i i y a bx σ==--=∑ 21.给出()ln f x x =的函数表如下: 用拉格朗日插值求ln 0.54的近似值并估计误差(计算取1n =及2n =) 解:1n =时,取010.5,0.6x x == 由拉格朗日插值定理有 1 100.60.5 0.693147 0.510826 0.50.(60.60.51.82321)0 1.()6047()52 j j j x x x L x f x l x ==------=-=∑ 所以1ln0.54(0.54)0.620219L ≈=- 误差为ln 0.54(0.620219)= 0.004032ε=-- 2n =时,取0120.4,0.5,0.6x x x === 由拉格朗日插值定理有

人口模型预测 数学建模作业

上传是为了分析数学的乐趣,请粘贴复制的时候也多思考哈。为了更多的 学子们。 2014年数学建模论文 第二套 题目:人口增长模型的确定 专业、姓名:土木135 提交日期: 2015/7/2晚上

题目:人口增长模型的确定 摘要 对美国人口数据的变化进行拟合,并进行未来人口预测,在第一个模型中,考虑到人口连续变化的规律,用微分方程的方法解出其数量随时间变化的方程,用matlab里的cftool工具箱求出参数,即人口净增长率r=0.02222,对该模型与实际数据进行对比,并计算了从1980年后每隔10年的人口数据,与实际对比,有很大出入。因此又改进出更为符合实际的阻滞增长模型,应用微分方程里的分离变量法和积分法解出其数量随时间变化的方程,求出参数人口增长率r=0.02858和人口所能容纳最大值m x=258.9,与实际数据对比,拟合得很好,并预测出1980年后每隔10年的人口数据,与实际对比,比较符合。为了便于比较两个模型与实际数据的描述情况作对比,又做出了两个模型与实际数据的对比图,以及两个模型的误差图。 关键词:人口预测微分方程马尔萨斯人口增长模型阻滞增长模型 一、问题重述 1790-1980年间美国每隔10年的人口记录如下表所示。 表1 人口记录表 试用以上数据建立马尔萨斯(Malthus)人口指数增长模型,并对接下来的每隔十年 预测五次人口数量,并查阅实际数据进行比对分析。 如果数据不相符,再对以上模型进行改进,寻找更为合适的模型进行预测。 二、问题分析 由于题目已经说明首先用马尔萨斯人口增长模型来刻划,列出人口增长指数增长方 程并求解,并进行未来50年内人口数据预测,但发现与实际数据有较大出入。考虑到 实际的人口增长率是受实际情况制约的,因此,使人口增长率为一变化的线性递减函数, 列出人口增长微分方程,求出其方程解,并预测未来五十年内人口实际数据。 三、问题假设

数学建模答案 (5)

一、解释下列词语,并举例说明(每小题满分5分,共15分) 1.模型 模型是系统知识的抽象表示。我们不能仅仅通过语言来描述一个系统,也不能仅仅通过记忆来记录关于系统的知识。知识是通过某种媒介来表达的,这种媒介所表达的内容就是模型。而知识形成媒介的过程就是建模,或者称为模型化。通常模型可以使用多种不同的媒介来表达,比如纸质或电子文档、缩微模型/原型、音像制品等等。而表达模型的体现方式也是多种多样的,常见的有图表、公式、原型、文字描述等等。 2.数学模型 由数字、字母、或其他数学符号组成的,描述现实对象(原型)数量规律的数学结构。具体地说,数学模型也可以描述为:对于现实世界的一个特定对象,为了一个特定的目的,根据特有的内在规律,做出一些简化假设后,运用适当的数学工具,得到的一个数学结构称之为数学模型.如概率的功利化定义。 3.抽象模型 通过人们对原型的反复认识,将获取的知识以经验的形式直接存储在大脑中的模型称之谓思维模型。从实际的人、物、事和概念中抽取所关心的共同特性,忽略非本质的细节把这些特性用各种概念精确地加以描述。 二、简答题(每小题满分8分,共24分) 1.模型的分类 按照模型替代原型的方式,模型可以简单分为形象模型和抽象模型两类.形象模型:直观模型、物理模型、分子结构模型等;抽象模型:思维模型、符号模型、数学模型等。 2.数学建模的基本步骤 1)建模准备:确立建模课题的过程; 2)建模假设:根据建模的目的对原型进行抽象、简化。有目的性原则、简明性原则、真实性原则和全面性原则; 3)构造模型:在建模假设的基础上,进一步分析建模假设的各条款,选择恰当的数学工具和构造模型的方法对其进行表征,构造出根据已知条件和数据,分析模型的特征和模型的结构特点,设计或选择求解模型的数学刻划实际问题的数学模型.; 4)模型求解:构造数学模型之后,方法和算法,并借助计算机完成对模型的求解; 5)模型分析:根据建模的目的要求,对模型求解的数字结果,或进行稳定性分析,或进行系统参数的灵敏度分析,或进行误差分析等。;

清华大学高等数值分析实验设计及答案

高等数值分析实验一 工物研13 成彬彬2004310559 一.用CG,Lanczos和MINRES方法求解大型稀疏对称正定矩阵Ax=b 作实验中,A是利用A= sprandsym(S,[],rc,3)随机生成的一个对称正定阵,S是1043阶的一个稀疏阵 A= sprandsym(S,[],0.01,3); 检验所生成的矩阵A的特征如下: rank(A-A')=0 %即A=A’,A是对称的; rank(A)=1043 %A满秩 cond(A)= 28.5908 %A是一个“好”阵 1.CG方法 利用CG方法解上面的线性方程组 [x,flag,relres,iter,resvec] = pcg(A,b,1e-6,1043); 结果如下: Iter=35,表示在35步时已经收敛到接近真实x relres= norm(b-A*x)/norm(b)= 5.8907e-007为最终相对残差 绘出A的特征值分布图和收敛曲线: S=svd(A); %绘制特征值分布 subplot(211) plot(S); title('Distribution of A''s singular values');; xlabel('n') ylabel('singular values') subplot(212); %绘制收敛曲线 semilogy(0:iter,resvec/norm(b),'-o'); title('Convergence curve'); xlabel('iteration number'); ylabel('relative residual'); 得到如下图象:

为了观察CG方法的收敛速度和A的特征值分布的关系,需要改变A的特征值: (1).研究A的最大最小特征值的变化对收敛速度的影响 在A的构造过程中,通过改变A= sprandsym(S,[],rc,3)中的参数rc(1/rc为A的条件数),可以达到改变A的特征值分布的目的: 通过改变rc=0.1,0.0001得到如下两幅图 以上三种情况下,由收敛定理2.2.2计算得到的至多叠代次数分别为:48,14和486,由于上实验结果可以看出实际叠代次数都比上限值要小较多。 由以上三图比较可以看出,A的条件数越大,即A的最大最小特征值的差别越大,叠代所需要的步骤就越多,收敛越慢。 (2)研究A的中间特征值的分布对于收敛特性的影响: 为了研究A的中间特征值的分布对收敛速度的影响,进行了如下实验: 固定A的条件数,即给定A的最大最小特征值,改变中间特征值得分布,再来生成A,具体的实现方法是,先将原来的生成A进行特征值分解: [U,S]=svd(A);

相关文档
最新文档