哈工大数学建模课作业第一题和第三题
数学建模第一次作业
14-15(2)数学建模第一次作业注意事项：提交时间截至3月27日课前，请将电子文档发送至邮箱sxjm@。
两个题目做到一个word文档里，文档和邮件标题均以“学号+姓名”命名。
请注意提交时间(顺序会影响给分结果)。
一、(必做题)ppt的思考题(1)~(4),由学号的后两位除以4的余数来确定；二、(必做题)本文档里的题目1~5，由学号的后两位除以5的余数来确定；三、(选做题)对于“生猪价格下降1%”理解的，0.65(11%)tp=-请根据ppt课件上的过程给出相应的结果(包括图形和灵敏性分析等)。
1油污清理问题一处石油泄漏污染了200英里的太平洋海岸线，所属石油公司被责令在14天内将其清除，预期则要被处以10000美元/天的罚款。
当地的清洁队每周可以清理5英里的海岸线，耗资500美元/天，额外雇佣清洁队则要付每支清洁队18000美元的费用和500美元/天的清洁费用.(1). 为使公司的总支出最低，应该额外雇佣多少支清洁队？采用5步方法，并求出清洁费用。
(2). 讨论清洁队每周清洁海岸线长度的灵敏性。
分别考虑最优的额外雇佣清洁队的数目和公司的总支出。
(3). 讨论罚金数额的灵敏性。
分别考虑公司用来清理漏油的总天数和公司的总支出。
(4). 石油公司认为罚金过高而提出上诉。
假设处以罚金的唯一目的是为了促使石油公司及时清理泄漏的石油，那么罚金的数额是否过高？*(5). (选做题)即使一开始采取围堵措施，海浪仍导致油污以每天0.5英里的速度沿海岸线扩散，这将导致最终清理的海岸线超过200海里，请分析扩散速度对公司总支出的影响。
2报刊价格问题一家有80000订户的地方日报计划提高其订阅价格。
现在的价格为每周1.5美元，据估计如果每提高定价10美分，就会损失5000订户。
(1)采用五步法，求使利润最大的订阅价格(2)对(1)中所得结论讨论损失5000订户这一参数的灵敏性。
分别假设这个参数值为3000，4000，5000，6000或7000，计算最优订阅价格(3)设n=5000为提高定价10美分而损失的订户数，求最优订阅价格p作为n的函数关系。
哈工大图论习题
1．画出具有4个顶点的所有无向图(同构的只算一个)。
2．画出具有3个顶点的所有有向图(同构的只算一个)。
3．画出具有4个、6个、8个顶点的三次图。
4．某次宴会上，许多人互相握手。
证明：握过奇数次手的人数为偶数(注意，0是偶数)。
5．证明：哥尼斯堡七桥问题无解。
6．设u与v是图G的两个不同顶点。
若u与v间有两条不同的通道(迹)，则G中是否有回路？7．证明：一个连通的(p，q)图中q ≥p-1。
8．设G是一个(p，q)图，δ(G)≥[p/2]，试证G是连通的。
9．证明：在一个连通图中，两条最长的路有一个公共的顶点。
10．在一个有n个人的宴会上，每个人至少有m个朋友(2≤m≤n)。
试证：有不少于m+1个人，使得他们按某种方法坐在一张圆桌旁，每人的左、右均是他的朋友。
11．一个图G是连通的，当且仅当将V划分成两个非空子集V1和V2时，G总有一条联结V1的一个顶点与V2的一个顶点的边。
12．设G是图。
证明：若δ(G)≥ 2，则G包含长至少是δ(G)+1的回路。
13．设G是一个(p，q)图，证明：(a)q≥p，则G中有回路；(b)若q≥p+4，则G包含两个边不重的回路。
14．证明：若图G不是连通图，则G c 是连通图。
15．设G是个(p，q)图，试证：(a)δ(G)·δ(G C)≤[(p-1)/2]([(p+1)/2]+1)，若p≡0，1，2(mod 4)(b) δ(G)·δ(G C)≤[(p-3)/2]·[(p+1)/2]，若p≡3(mod 4)16．证明：每一个自补图有4n或4n+1个顶点。
17．构造一个有2n个顶点而没有三角形的三次图，其中n≥3。
18.给出一个10个顶点的非哈密顿图的例子，使得每一对不邻接的顶点u和v，均有degu+degv≥919．试求Kp中不同的哈密顿回路的个数。
20．试证：图四中的图不是哈密顿图。
21．完全偶图Km，n为哈密顿图的充分必要条件是什么？22．菱形12面体的表面上有无哈密顿回路？23．设G是一个p(p≥3)个顶点的图。
建模习题答案
田佳王伊陈鹏《数学建模入门》练习题练习题1：发现新大陆！发现新大陆！人人都能做到，可是最终哥伦布做到了。
为什么哥伦布能做到呢？(参考答案：有兴趣、能想到、去做了、坚持到底。
)答： 1）从其主观条件分析：他具有一个优秀水手的素质：对大海的热爱，具有宝贵的航海经验，接触过航海所必不可少的宇宙学和数学，并且学会了绘制地图和使用各种航海工具。
更为重要的事，在航海强国葡萄牙，哥伦布在思想上为远航做好了准备。
他阅读了《马可·波罗游记》，对东方的富饶遐想无限，使他产生了到东方区的想法；他接触了学者托斯勘内里，接受了“地圆学说”，坚定了从海上到达东方的信念。
2）从客观条件分析：出于共同的对黄金的追求，哥伦布与西班牙王室达成了一致（签订《圣塔菲协定》），西班牙为其提供了自己的船队、自己的船员。
当时中国的指南针也已传到航海界，这一发明对其也有及其重要的作用练习题2：棋盘问题有一种棋盘有64个方格，去掉对角的两个格后剩下62个格（如下图），给你31块骨牌，每块是两个格的大小。
问能否用这些骨牌盖住这62个方格？答：这个问题涉及到数学上的一个典型排列：完美覆盖31张不重叠的多米诺牌则盖住31个白方格和31个黑方格。
因此，这副被剪过的棋盘没有完美覆盖，上述推理可总结为：31黑白 32黑＋30白更一般地．可以将棋盘上的方格交替徐成黑色和自色，切除一些方格，得到一块切过的棋盘什么时候能有一个完美覆盖？为使完美覆盖存在，这块被切过的棋盘必须又有相等的黑方格数和白方格数但是，这个条件却不是充分的，最后是不能够用这些骨牌盖住其余方格的练习题3：硬币游戏如果你和你的对手准备依次轮流地将硬币放在一个长方形桌子上，使得这些硬币不重叠。
最后放上硬币的人为胜者，在开始时你有权决定先放还是后放。
为了能赢得这场比赛，你决定先放还是后放呢？答：决定先放。
首先将硬币放在长方形桌子的中心，然后根据对手所放的硬币，找一桌子中心为对称中心的位置，直至对手没有地方放硬币为止，由于长方形的对称性，只有中心不存在对称位置，故先放者赢。
《数学建模》课程第一章自测练习及解答提示
《数学建模》课程第一章自测练习及解答提示一、填空题：1．设年利率为0.05，则10年后20万元的现值按照复利计算应为 . 解：根据现值计算公式：10)05.01(20)1(+=+=n R S Q 2783.1221201011≈=（万元） 应该填写：12.2783万元．2．设年利率为0.05，则20万元10年后的终值按照复利计算应为 . 解：根据终值计算公式：10)05.01(20)1(+=+=n R P S =5779.322021910=（万元） 应该填写：32.57793．所谓数学建模的五步建模法是指下列五个基本步骤，按一般顺序可以写出为 .解：应该填写：问题分析，模型假设，模型建立，模型求解，模型分析．4．设某种商品的需求量函数是,1200)(25)(+-=t p t Q 而供给量函数是3600)1(35)(--=t p t G ，其中)(t p 为该商品的价格函数，那麽该商品的均衡价格是 .解： 由商品的均衡价格公式：80352536001200)(=++=++=c a d b t p 应该填写：80.5．一家服装店经营的某种服装平均每天卖出110件，进货一次的批发手续费为200元，存储费用为每件0.01元/天，店主不希望出现缺货现象，则最优进货周期与最优进货量分别为 .解：根据经济订购批量公式：1911001.020022*≈⨯⨯==R c c T s b 209701.011020022*≈⨯⨯==s b c R c Q 应该填写：.2097,19**=≈Q T二、分析判断题1. 从下面不太明确的叙述中确定要研究的问题，需要哪些数据资料（至少列举3个），要做些甚麽建模的具体的前期工作（至少列举3个） ，建立何种数学模型：一座高层办公楼有四部电梯，早晨上班时间非常拥挤，该如何解决．解：（1）要研究的问题：如何设置四部电梯的停靠方式，使之发挥最大效益．（2）所需资料为：每天早晨乘电梯的总人数、各层上、下电梯的人数、电梯的速度、楼层的高度、层数等．（3）要做的具体建模前期工作：观察和统计所需资料，一般讲，需要统计一周内每天的相关资料．（4）可以建立概率统计模型，亦可在适当的假设下建立确定性模型．2.一条公路交通不太拥挤，以至人们养成“冲过”马路的习惯，不愿意走临近的“斑马线”．交管部门不允许任意横穿马路，为方便行人，准备在一些特殊地点增设“斑马线”，以便让行人可以穿越马路．那末“选择设置斑马线的地点”这一问题应该考虑哪些因素？试至少列出3种．解：（1）车流的密度（2）车的行驶速度（3）道路的宽度（4）行人穿越马路的速度（5）设置斑马线地点的两侧视野等．3．怎样解决下面的实际问题．包括需要哪些数据资料，要作些什么观察、试验以及建立什么样的数学模型等．（1）估计一个人体内血液的总量．（2）为保险公司制定人寿保险计划（不同年龄的人应缴纳的金额和公司赔偿的金额）．（3）估计一批日光灯管的寿命．（4）确定火箭发射至最高点所需的时间．（5）决定十字路口黄灯亮的时间长度．（6）为汽车租赁公司制订车辆维修、更新和出租计划．（7）一高层办公楼有4部电梯，早晨上班时间非常拥挤，试制订合理的运行计划解：（1）注射一定量的葡萄糖，采集一定容量的血样，测量注射前后葡萄糖含量的变化，即可估计人体的血液总量．注意采集和测量的时间要选择恰当，使血液中的葡萄糖含量充分均匀，又基本上未被人体吸收．（2）调查不同年龄的人的死亡率，并估计其在未来一定时期的变化，还应考虑银行存款利率和物价指数，保险金与赔偿金之比大体上应略高于死亡率．（3）从一批灯管中取一定容量的样本，测得其平均寿命，可作为该批灯管寿命的估计值．为衡量估计的精度，需要从样本寿命确定该批灯管寿命的概率分布，即可得到估计值的置信区间．还可试验用提高电压的办法加速寿命测试，以缩短测量时间．（4）根据牛顿第二定律建立火箭向上发射后的运动方程，初速已知，若不考虑空气阻力，很容易算出到达最高点（即速度为零）时间；若考虑空气阻力，不妨设其与火箭速度（或速度的平方）成正比，并有试验及拟合方法确定阻力系数，再解方程得到结果．（5）司机看到黄灯后停车要有一定的刹车距离S1，设通过十字路口的距离为S2，汽车行驶速度为v，则黄灯的时间长度t应使距停车线S1之内的汽车能通过路口，即t （S1+S2）/v．S1可由试验得到，或按照牛顿第二定律解运动方程，进一步可考察不同车重、不同路面及司机反应灵敏程度等因素的影响．（6）根据资料和经验确定维修费用随着车龄和行驶里程的增加而增加的关系，再考虑维修和更新费用，可以以一年为一个时段，结合租金决定应该维修或更新．（7）统计在各层上班的人数，通过数据或计算确定电梯运行时间，以等待的人数与时间乘积为目标，建立优化模型，确定每部电梯运行的楼层（有的从大厅直接运行到高层）．4．为了培养想象力、洞察力，考察对象时除了从正面分析外，还常常需要从侧面或反面思考，试尽可能迅速地回答下列的问题：（1）某甲早8：00从山下旅馆出发，沿一条路径上山，下午5：00到达山顶并留宿．次日早8：00沿同一路径下山，下午5：00回到旅馆．某乙说，甲必在2天中的同一时刻经过路径中的同一地点．为什么？（2）甲乙两站之间有电车相通，每隔10分钟甲乙两站相互发一趟车，但发车时刻不一定相同，甲乙之间有一中间站丙，某人每天在随机的时刻到达丙站，并搭乘最先经过丙站的那趟车，结果发现100天中约有90天到达甲站，约有10天到达乙站．问开往甲乙两站的电车经过丙站的时刻表是如何安排的？（3）某人住T 市在他乡工作，每天下班后乘火车于6：00抵达T 市车站，他的妻子驾车准时到车站接他回家．一日他提前下班搭乘早一班火车于5：30抵T 市车站，随即步行回家，他的妻子像往常一样驾车前往，在半路上遇到他，即接他回家，此时发现比往常提前10分钟．问他步行了多长时间．解：（1）设想有两个人一人上山，一人下山，同一天同时出发，沿同一路径，必定相遇．（2）不妨设从甲站到乙站经过丙站的时刻表是：8：00，8：10，8：20，…，那么从乙站到甲站经过丙站的时刻表应该是：8：09，8：19，8：29，…．（3）步行了25分钟．设想他的妻子驾车遇到他后，先带他去车站，再回家，汽车多行驶了10分钟，于是带他去车站这段路程汽车跑了5分钟，而到车站的时间是6：00，所以妻子驾车遇到他的时刻是5：55．三、计算题1．下面是众所周知的智力游戏：人带猫、鸡、米过河，船除需要人划之外，至多能载猫、鸡、米三者之一，而当人不在场时猫要吃鸡、鸡要吃米．试设计一个安全过河方案，并使渡河次数尽量地少．解：人、猫、鸡、米分别记为i =1, 2, 3, 4，当i 在此岸时记x i =1，否则记x i =0，则此岸的状态可用s =（x 1, x 2, x 3, x 4）表示．记s 的反状态为s '=（1-x 1, 1-x 2, 1-x 3, 1-x 4），允许状态集合为S ={（1, 1, 1, 1），（1, 1, 1, 0），（1, 1, 0, 1），（1, 0, 1, 1）（1, 0, 1, 0）及它们的5个反状态}．决策为乘船方案，记作d =（u 1, u 2, u 3, u 4），当i 在船上时记u i =1，否则记u i =0，允许决策集合为D ={（1, 1, 0, 0），（1, 0, 1, 0），（1, 0, 0, 1），（1, 0, 0, 0）}．记第k 次渡河前的状态为s k ，第k 次渡河的决策为d k ，则状态转移律为s k +1=s k +（-1）k d k ，设计安全过河方案归结为求决策序列d 1, d 2, …, d n ∈D ，使状态s n ∈S 按状态转移律由初始状态s 1=（1, 1, 1, 1）经n 步到达s n +1=（0, 0, 0, 0）．一个可行方案如下：2．假定人口的增长服从这样的规律：时间t 的人口为x (t )，t 到t +∆t 时间内人口的增长与x m - x (t )成正比 （其中x m 为最大容量）．试建立模型并求解．作出解的图形并与指数增长模型、阻滞增长模型的结果进行比较．解 )(d d x x r t x m -=，r 为比例系数,0)0(x x =， 解为rt m m x x x t x ---=e )()(0，如图1中粗实线所示．当t 充分大时，它与Logistic 模型相近． 图13．在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗？比如洁银牙膏50g 装的每支1.50元，120g 装的每支3.00元，二者单位重量的价格比是1.2：1．试用比例x x方法构造模型解释这个现象．（1）分析商品价格c 与商品重量w 的关系．价格由生产成本、包装成本和其它成本决定，这些成本中有的与重量w 成正比，有的与表面积成正比，还有与w 无关的因素．（2）给出单位重量价格c 与w 的关系，画出它的简图，说明w 越大c 越小，但是随着w 的增加c 减小的程度变小．解释实际意义是什么？ 解：（1）生产成本主要与重量w 成正比，包装成本主要与表面积s 成正比，其它成本也包含与w 和s 成正比的部分，上述三种成本中都含有与w 和s 无关的成分．又因为形状一定时一般有s ∝w 2/3，故商品的价格可表为C = αw +β w 2/3+γ（α，β，γ为大于0的常数）．（2）单位重量价格131--++==w w w C c γβα，其 图2 简图如图2所示．显然c 是w 的减函数，说明大包装商品比小包装商品便宜；曲线是下凸的,4．用宽w 的布条缠绕直径d的圆形管道，要求布条 不重叠，问布条与管道轴线的夹角α应多大（如图3）． 若知道管道长度，需用多长布条（可考虑两端的影响）． 如果管道是其它形状呢？解：将管道展开如图4，可得απcos d w =，若d 一 图3 定，0→w ，2πα→；d w π→，0→α．若管道长度为l ，不考虑两端的影响时布条长度显然为w dl π，若考虑两 端的影响，则应加上απsin dw ．对于其它形状管道，只需将d π 改为相应的周长即可．5．建立不允许缺货的生产销售存贮模型．设生产速率为常数k ，销售速率为常数r ，k >r ．在每一生产周期T 内，开 图4 始的一段时间（0<t <T 0）一边生产一边销售，后来的一段时间（T 0<t <T ）只销售不生产，画出贮存量)(t q 的图形．设每次生产准备费为1c ，单位时间每件产品贮存费为2c ，以总费用最小为目标确定最优生产周期．讨论k 》r 和k ≈ r 的情况． 解： 贮存量)(t q 的图形如图5．单位时间总费用 KT r k r c T c T c 2)()(21-+=， 使)(T c 达到最小值的最优周期)(221r k r c k c T -=*． 图5 当k 》r 时，rc c T 212=*，相当 于不考虑生产的情况．当k ≈ r 时，∞→*T ，因为0产量被销量抵消，无法形成贮存量．四、综合应用题1．试建立方桌问题在四条腿脚呈长方形情形时的数学模型，以说明方桌能否在地面上放稳的问题．（ 提示：要求按照五步建模法进行建模工作，本题至少应给出前四个步骤．） 解：问题分析所谓方桌可否在地面上放稳，可视为其四个桌脚可否同时着地，从而可将问题归结为桌脚与地面的距离是否同时为零，故构造这个距离函数是建模的关键，而证明四个距离函数同时为零这个命题是建模的最终目的．模型假设(1) 四条桌腿同长,视四个桌脚为四个几何点,四脚的连线呈长方形；(2) 地面的高度是连续变化的，即将地面看作数学上的连续曲面；(3)模型建立 如图6，以长方形的两条对角线的交点为原点建立平面直角坐标系，且不妨设A ，C 两桌脚开始时位于横轴上，则问题与旋转角度θ有关．注意到假设3，设A ，B 两个桌脚与地面距离之和为0)(≥θf ，另外两个桌脚与地面距离之 和为,0)(≥θg 则)(,θθf ∀与)(θg 中至少有一个为零，当 图6 0=θ时不妨假设0)(,0)(>=θθg f ．又由假设2，以上两个函数均为旋转角度的连续函数，于是有命题：已知,0)0(,0)0(,0)()()(),(>==∀g f g f g f 且，的连续函数，对是θθθθθθ则0θ∃，使得.0)()(00==θθg f上述命题即为所建立的数学模型．模型求解将桌子旋转0180)(π，则A 、B 两点与D 、C 两点恰好交换位置．由假设便有，)(,0)(ππg f >.0=又由前述假设，.0)0(,0)0(>=g f令),()()(θθθg f h -=则有.0)(,0)0(><πh h 由于)(),(θθg f 的连续性知)(θh 也是连续函数．依据连续函数的基本性质（零点定理），必至少存在一个角度0θ，,00πθ<<使得0)(0=θh ，即).()(00θθg f =又根据θθθ∀=,0)()(g f 成立，故有.0)()(00==θθg f 模型分析由于本问题结论简单，符合实际，故分析过程从略．2．试建立确定情形下允许缺货的存储问题的数学模型．提示： 所谓的确定情形下的存储模型是指文字教材第一章提到过的不允许缺货的存储模型；所谓允许缺货是在不允许缺货模型假设条件下，再考虑因缺货造成的损失建立相应的模型．（要求按照五步建模法进行建模工作，本题应给出五个步骤．）解： 问题分析由题设，只须在不允许缺货模型条件下，考虑因缺货造成的损失即可．而缺货损失按天计算与下列因素有关：货物总需求量、缺货量、缺货时刻、每单位的缺货费用等． 模型假设 （1）每次定货费为C 1，每天每单位货物的存储费为C 2 （2）每天货物的需求量为r 单位.（3） 每T 天定货Q 单位，所定货物可在瞬间到达．（4）允许缺货，每天每单位货的缺货费为C 3缺货时，存储量q 视为负值，则)(t q 的图形变为,Q rt q +-=如图7所示．模型建立 图7 货物在1T t =时售完，则必有一段时间缺货．又在T t =时下一次定货量Q 到达，于是有1rT Q = （1）在一个定货周期内的总费用包括定货费1C 、存储费Q T C dt t q C T 102221)(1⎰=和缺货费.)(13dt t q C T T ⎰其中21)(2)()(11T T r dt Q rt dt t q TT T T -=-=⎰⎰ 其中用到了（1）式．于是总费用应为2/)(2/213121T T r C QT C C C -++= （2） 则由（1）式解出r Q T /1=并代入（2）式可得r Q rT C r Q C C C 2/)(2/23221-++= （3）每天的平均总费用便是rT Q rT C rT Q C T C T C Q T C 2/)(2///),(23221-++== （4）（4）式即为所求的数学模型．模型求解对（4）式分别求总费用对定货周期和定货量的偏导数，并令其为零解得0)()(22322322221=-+----=∂∂Q rT T C Q rT rT C rT Q C T C T C0)(32=--=∂∂Q rT rTC rT Q C Q C 由3230C C rT C Q Q C +=⇒=∂∂，代入0=∂∂TC 便可解出 32321*33221*2;2C C C C r C Q C C C rC C T +=+=. （5） （5）式就是在允许缺货情形下，最佳定货周期与最佳定货量公式．模型分析当3C 远远超过2C 时，（5）式就转化为不允许缺货模型中的相应结论，这也说明所建模型是合理的，结论也是正确的．。
哈工大图论习题
1．画出具有4个顶点的所有无向图(同构的只算一个)。
2．画出具有3个顶点的所有有向图(同构的只算一个)。
3．画出具有4个、6个、8个顶点的三次图。
4．某次宴会上，许多人互相握手。
证明：握过奇数次手的人数为偶数(注意，0是偶数)。
5．证明：哥尼斯堡七桥问题无解。
6．设u与v是图G的两个不同顶点。
若u与v间有两条不同的通道(迹)，则G中是否有回路？7．证明：一个连通的(p，q)图中q ≥p-1。
8．设G是一个(p，q)图，δ(G)≥[p/2]，试证G是连通的。
9．证明：在一个连通图中，两条最长的路有一个公共的顶点。
10．在一个有n个人的宴会上，每个人至少有m个朋友(2≤m≤n)。
试证：有不少于m+1个人，使得他们按某种方法坐在一张圆桌旁，每人的左、右均是他的朋友。
11．一个图G是连通的，当且仅当将V划分成两个非空子集V1和V2时，G总有一条联结V1的一个顶点与V2的一个顶点的边。
12．设G是图。
证明：若δ(G)≥ 2，则G包含长至少是δ(G)+1的回路。
13．设G是一个(p，q)图，证明：(a)q≥p，则G中有回路；(b)若q≥p+4，则G包含两个边不重的回路。
14．证明：若图G不是连通图，则G c 是连通图。
15．设G是个(p，q)图，试证：(a)δ(G)·δ(G C)≤[(p-1)/2]([(p+1)/2]+1)，若p≡0，1，2(mod 4)(b) δ(G)·δ(G C)≤[(p-3)/2]·[(p+1)/2]，若p≡3(mod 4)16．证明：每一个自补图有4n或4n+1个顶点。
17．构造一个有2n个顶点而没有三角形的三次图，其中n≥3。
18.给出一个10个顶点的非哈密顿图的例子，使得每一对不邻接的顶点u和v，均有degu+degv≥919．试求Kp中不同的哈密顿回路的个数。
20．试证：图四中的图不是哈密顿图。
21．完全偶图Km，n为哈密顿图的充分必要条件是什么？22．菱形12面体的表面上有无哈密顿回路？23．设G是一个p(p≥3)个顶点的图。
数学建模第一次作业作业
(i)取定 x0 3.9, t0 1790, ，拟合待定参数 r .
t=1790:10:2000; c=[3.9,5.3,7.2,9.6,12.9,17.1,23.2,31.4,38.6,50.2,62.9,76.0,92.0,106.5,123.2,131.7,150.7,179.3,20 4.0,226.5,251.4,281.4]; f1=@(r,t)3.9.*exp(r.*(t-1790)); r0=0.02; r=nlinfit(t,c,f1,r0), se1=sum((c-f1(r,t)).^2), plot(t,c,'k+',1780:1:2010,f1(r,1780:1:2010),'k') (ii)取定 t0 1790 ，拟定待定参数 t0 、 x0 、 r
数学建模第一次作业 1、绘制图形 （1）程序及图形如下： n=500; t=linspace(0,2*pi,n); x1=cos(t); y1=sin(t); x2=2*cos(t); y2=2*sin(t); x3=2*cos(t); y3=sin(t); plot(x1,y1,'k',x2,y2,'k',x3,y3,'k') axis equal;title('参数方程画 x^2+y^2=1, x^2+y^2=4, x^2/4+y^2=1 的图像'); gtext('x^2+y^2=1') gtext('x^2+y^2=4') gtext('x^2/4+y^2=1')
2
黎曼函数的图像 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 分 母 P的 最 大 值 n =36时
数学建模作业题+答案
数学建模MATLAB 语言及应用上机作业11. 在matlab 中建立一个矩阵135792468101234501234A ⎡⎤⎢⎥⎢⎥=⎢⎥-----⎢⎥⎣⎦答案：A = [1,3,5,7,9;2,4,6,8,10;-1,-2,-3,-4,-5;0,1,2,3,4]2. 试着利用matlab 求解出下列方程的解（线性代数22页例14）123412423412342583692254760x x x x x x x x x x x x x x +-+=⎧⎪--=⎪⎨-+=-⎪⎪+-+=⎩ 答案：A=[2 ,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6]; B=[8;9;-5;0]; X=A\B 或A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6] b=[8,9,-5,0]' X=inv(A)*b3. 生成一个5阶服从标准正态分布的随机方阵，并计算出其行列式的值，逆矩阵以及转置矩阵。
答案:A=randn(5) det(A) inv(A) A'4. 利用matlab 求解出110430002A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦的特征值和特征向量。
答案：A=[-1,1,0;-4,3,0;0,0,2] [V,D]=eig(A)5.画出衰减振荡曲线3sin3t y et -=在[0,4]π上的图像。
要求，画线颜色调整为黑色，画布底面为白色。
（在实际中，很多打印机时黑白的，因此大多数作图要考虑黑白打印机的效果。
） 给出恰当的x ，y 坐标轴标题，图像x 轴的最大值为4π。
6. 生成一个0-1分布的具有10个元素的随机向量，试着编写程序挑选出向量中大于0.5的元素。
数学建模和Matlab 上机作业2（2016-9-20）跟老师做（不用整合进作业中）：上机演示讲解：函数，递归的两个例子的写法。
附：1. Fibonacci Sequence （斐波那契数列）在数学上，费波那西数列是以递归的方法来定义： F1= 1；F2= 1；F （n ）=F （n-1）+F （n-2） 2. 阶乘举例：数学描述：n!=1×2×……×n ；计算机描述：n!=n*(n-1)!自己做（需要整合进作业中，提交到系统中）：1. 写一个m 文件完成分值百分制到5分制的转换（即输入一个百分制，转换后输出一个5级对应的得分，联系条件控制语句）。
数学模型-哈尔滨工业大学
系统建模的抽象过程
（可观测） 输入变量 （可观测） 输出变量
真实系统
观测屏障
抽象
真实系统
ω（t） ρ（t）
黑箱 灰箱 白箱
数学描述 ω（t）、ρ（t）---输入输出变量对
17
数学建模信息源
数据
2．先验知识
3．试验数据
先验 知识
建模
模型应用
目的 达到否？
目标
1．目标和目的
18
图2.1 数学建模的信息源
在定义一个系统时，首先要确定系统的边界。尽管世界上的事物是 相互联系的，但当我们研究某一对象时，总是要将该对象与其环境区别开 来。边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系 统对边界以外的环境的作用称为系统的输出。
边界 输出
环境
输入
系统
8
研究上述系统时，可分为三类问题： （1）系统的分析问题 （2）系统的控制问题 （3）系统的辨识问题
K 0 s W0 ( s ) e (Ts 1) n
15
1.1.3 数学模型
本质上讲，系统数学模型是从系统概念出发的关 于现实世界的一小部分或几个方面的抽象的“映像”。
系统数学模型的建立需要建立如下抽象：输入、 输出、状态变量及其间的函数关系。这种抽象过程称 为模型构造。
系统数学建模就是将真实系统抽象成相应的数学 表达式（一些规则、指令的集合）。
2
主要内容安排
第一章 绪论 第二章 系统辨识常用输入信号 第三章 系统数学描述及经典辨识法 第四章 最小二乘法辨识
第五章 极大似然法辨识
第六章 系统结构辨识
3
实验1 白噪声和M序列的产生 实验2 相关分析法辨识脉冲响应
2023年全国数学建模a题第一题
2023年全国数学建模A题第一题一、概述在2023年全国数学建模比赛中，A题的第一题是一个涉及到矩阵运算和线性代数知识的题目。
该题目具有一定的挑战性和实用性，要求参赛选手能够灵活运用数学知识解决实际问题，展现出扎实的数学建模能力。
二、题目要求本题要求选手根据给定的矩阵运算规则，进行定性分析和定量分析，并给出相应的理论模型。
1. 给定一个m×n的实矩阵A=[本人j]，其中1≤i≤m，1≤j≤n。
对矩阵A中的每个元素进行如下的操作：若本人j是正数，则将本人j的值减半；若本人j是负数，则将本人j取绝对值。
2. 对经过上述操作后的矩阵A进行转置得到矩阵B，即B=A^T。
3. 给定一个m×n的实矩阵C=[cij]，其中1≤i≤m，1≤j≤n。
对矩阵C 中的每个元素进行如下的操作：若cij是正数，则将cij的值加倍；若cij是负数，则将cij取相反数。
4. 对经过上述操作后的矩阵C进行转置得到矩阵D，即D=C^T。
5. 求证：矩阵B与矩阵D的乘积得出的矩阵E满足以下性质：矩阵E 是一个对称矩阵。
三、解题思路针对题目要求，参赛选手需要按照以下步骤进行解题：1. 对给定矩阵A进行定性分析，分析经过操作后的矩阵B的特点和性质。
2. 对给定矩阵C进行定量分析，计算经过操作后的矩阵D的数值。
3. 推导并证明矩阵B与矩阵D的乘积得出的矩阵E是一个对称矩阵。
四、解题方法针对题目要求，可以采用以下解题方法：1. 对矩阵A进行定性分析，可以根据给定的操作规则，通过数学推导和分析，得出矩阵B的特点和性质，包括矩阵B的元素取值规律、矩阵B的转置性质等。
2. 对矩阵C进行定量分析，可以利用线性代数知识，计算经过操作后的矩阵D的数值，并对矩阵D的元素取值规律进行分析。
3. 推导并证明矩阵B与矩阵D的乘积得出的矩阵E是一个对称矩阵，可以通过矩阵乘法运算和对称矩阵的性质进行推导和证明。
五、解题过程针对题目要求，参赛选手可以按照以下步骤进行解题过程的详细描述：1. 对矩阵A进行定性分析，包括对矩阵A中每个元素进行操作后的取值规律的推导和分析，以及矩阵B的转置性质的分析。
数学建模习题及答案课后习题
第一部分课后习题1. 学校共1000名学生，235人住在A宿舍，333人住在B宿舍，432人住在C宿舍。
学生们要组织一个10人的委员会，试用下列办法分配各宿舍的委员数：(1)按比例分配取整数的名额后，剩下的名额按惯例分给小数部分较大者。
(2) 2.1节中的Q值方法。
(3)d' Hondt方法：将A , B, C各宿舍的人数用正整数n=1 , 2, 3,…相除，其商数如下表：将所得商数从大到小取前10个(10为席位数)，在数字下标以横线，表中A, B, C行有横线的数分别为2, 3, 5,这就是3个宿舍分配的席位。
你能解释这种方法的道理吗。
如果委员会从10人增至15人，用以上3种方法再分配名额。
将3种方法两次分配的结果列表比较。
(4)你能提出其他的方法吗。
用你的方法分配上面的名额。
2. 在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗。
比如洁银牙膏50g装的每支1.50元，120g装的3.00元，二者单位重量的价格比是 1.2: 1。
试用比例方法构造模型解释这个现象。
(1)分析商品价格C与商品重量w的关系。
价格由生产成本、包装成本和其他成本等决定，这些成本中有的与重量w成正比，有的与表面积成正比，还有与w无关的因素。
(2)给出单位重量价格c与w的关系，画出它的简图，说明w越大c越小，但是随着w 的增加c减少的程度变小。
解释实际意义是什么。
3. 一垂钓俱乐部鼓励垂钓者将调上的鱼放生，打算按照放生的鱼的重量给予奖励，俱乐部只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。
假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据(胸围指鱼身的最大周长) ：4. 用宽w的布条缠绕直径d的圆形管道，要求布条不重叠，问布条与管道轴线的夹角应多大(如图)。
若知道管道长度，需用多长布条(可考虑两端的影响) 。
如果管道是其他形状呢。
5. 用已知尺寸的矩形板材加工半径一定的圆盘，给出几种简便、有效的排列方法，使加工 出尽可能多的圆盘。
数学建模作业及答案
数学建模作业姓名：叶勃学号：班级：024121一：层次分析法1、 分别用和法、根法、特征根法编程求判断矩阵1261/2141/61/41A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦11/2433217551/41/711/21/31/31/52111/31/5311A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦的特征根和特征向量（1）冪法求该矩阵的特征根和特征向量 程序为：#include<iostream> #include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20 #define err 0.0001 //幂法求特征值特征向量 void main(){cout<<"**********幂法求矩阵最大特征值及特征向量***********"<<endl; int i,j,k;double A[n][n],X[n],u,y[n],max;cout<<"请输入矩阵:\n"; for(i=0;i<n;i++) for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 cout<<"请输入初始向量:\n"; for(i=0;i<n;i++)cin>>X[i]; //输入初始向量 k=1; u=0;while(1){ max=X[0]; for(i=0;i<n;i++) {if(max<X[i]) max=X[i]; //选择最大值 }for(i=0;i<n;i++)y[i]=X[i]/max; for(i=0;i<n;i++)X[i]=0;for(j=0;j<n;j++)X[i]+=A[i][j]*y[j]; //矩阵相乘}if(fabs(max-u)<err){cout<<"A的特征值是 :"<<endl; cout<<max<<endl; cout<<"A的特征向量为:"<<endl; for(i=0;i<n;i++) cout<<X[i]/(X[0]+X[1]+X[2])<<" ";cout<<endl;break;}else{if(k<N) {k=k+1;u=max;} else {cout<<"运行错误\n";break;}}} }程序结果为：（2）和法求矩阵最大特征值及特征向量程序为：#include<stdio.h>#include<iostream>#include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j,k;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********和法求矩阵的特征根及特征向量*******"<<endl;cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 //计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;} //求特征向量w[0]=0;w[1]=0;w[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){w[i]+=W[i][j];}cout<<"特征向量为:"<<endl; for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征根为："<<endl;cout<<max/n<<endl; }运行结果为：（3）根法求矩阵最大特征值及特征向量：程序为：#include<stdio.h>#include<iostream>#include<math.h>using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********根法求矩阵的特征根及特征向量*******"<<endl; cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵//计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;}//求特征向量//w[0]=A[0][0];w[1]=A[0][1];w[2]=A[0][2];w[0]=1;w[1]=1;w[2]=1;for(i=0;i<n;i++){for(j=0;j<n;j++){w[i]=w[i]*W[i][j];}w[i]=pow(w[i], 1.0/3);}cout<<"特征向量为:"<<endl;for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征值为:"<<endl; cout<<max/n;}运行结果为：2、编程验证n阶随机性一致性指标RI：运行结果：3、考虑景色、费用、居住、饮食、旅途五项准则，从桂林、黄山、北戴河三个旅游景点选择最佳的旅游地。
数学建模课后习题作业
选修课——数学建模部分习题详细解答【陈文滨】1、在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？【模型假设】（1）椅子四条腿一样长，椅脚与地面接触处视为一点，四脚的连线呈长方形．（2）地面高度是连续变化的，沿任何方向都不会出现间断 (没有像台阶那样的情况)，即从数学的角度看，地面是连续曲面．这个假设相当于给出了椅子能放稳的必要条件．（3）椅子在任何位置至少有三只脚同时着地．为保证这一点，要求对于椅脚的间距和椅腿的长度而言，地面是相对平坦的．因为在地面上与椅脚间距和椅腿长度的尺寸大小相当的范围内，如果出现深沟或凸峰(即使是连续变化的)，此时三只脚是无法同时着地的。
【模型建立】在上述假设下，解决问题的关键在于选择合适的变量，把椅子四只脚同时着地表示出来．首先，引入合适的变量来表示椅子位置的挪动．生活经验告诉我们，要把椅子通过挪动放稳，通常有拖动或转动椅子两种办法，也就是数学上所说的平移与旋转变换．然而，平移椅子后问题的条件没有发生本质变化，所以用平移的办法是不能解决问题的．于是可尝试将椅子就地旋转，并试图在旋转过程中找到一种椅子能放稳的情形．注意到椅脚连线呈长方形，长方形是中心对称图形，绕它的对称中心旋转180度后，椅子仍在原地．把长方形绕它的对称中心O旋转，这可以表示椅子位置的改变。
于是，旋转角度θ这一变量就表示了椅子的位置．为此，在平面上建立直角坐标系来解决问题．如下图所示，设椅脚连线为长方形ABCD，以对角线AC所在的直线为x轴，对称中心O为原点，建立平面直角坐标系．椅子绕O点沿逆时针方向旋转角度θ后，长方形ABCD转至A1B1C1D1 的位置，这样就可以用旋转角θ（0≤θ≤π）表示出椅子绕点O旋转θ后的位置．其次，把椅脚是否着地用数学形式表示出来．我们知道，当椅脚与地面的竖直距离为零时，椅脚就着地了，而当这个距离大于零时，椅脚不着地．由于椅子在不同的位置是θ的函数，因此，椅脚与地面的竖直距离也是θ的函数．由于椅子有四只脚，因而椅脚与地面的竖直距离有四个，它们都是θ的函数．而由假设(3)可知，椅子在任何位置至少有三只脚同时着地，即这四个函数对于任意的θ，其函数值至少有三个同时为0．因此，只需引入两个距离函数即可．考虑到长方形ABCD是中心对称图形，绕其对称中心 O沿逆时针方向旋转180°后，长方形位置不变，但A,C和B,D对换了．因此，记A、B两脚与地面竖直距离之和为f（θ），C、D两脚与地面竖直距离之和为g（θ），其中θ∈[0，π]，从而将原问题数学化。
数学建模资料 2004哈尔滨工业大学数学建模竞赛试题
2004哈尔滨工业大学数学建模竞赛试题A：西大直街的交通线联动信号控制问题城市的交通管理是众所周知的难题。
当你乘车时令你烦心的是刚通过了一个十字路口就被不远处另一十字路口的红灯止住。
能否通过交通信号的联动控制使得在主干线上最大可能的解决这一问题，请你对西大直街从护军街到通达街一线的考察研究，用数学建模的方法给出一个最优联动红绿灯控制方案。
要求：1.对西大直街从护军街到通达街一线17时至18时的交通现状进行现场调查，采集建模所需的数据；2.对现行交通信号系统的合理性进行分析；3.给出你的最优联动红绿灯控制原理和方案；4.给当地报纸写一篇短文（不超过两页），阐述你的方案的可行性。
B：股市全流通方案的设想全国人大常委会副委员长成思危日前在"第八届(2004)资本市场论坛"上指出,股权的流动性分裂给资本市场的发展带来了很多弊病,因此股市要在规范的同时,重视发展,在发展中实现全流通。
股市的全流通问题应该考虑到：一、有利于资本市场的改革开放和稳定发展；二、有利于保护中小投资者的合法权益。
推进全流通,要注意可行性,同时一定要有周密的策划（附件一）。
张卫星提出了中国股市的全流通改造三条原则（附件二）。
依据以上原则对股市全流通方案通过数学建模提出你的设想，要求：1.对上证50（1.浦发银行、2.白云机场、……、50.长江电力）的最近20日均价与总股本、流通股（A、B股）占总股本的比例、03年每股收益、净资产、概念（国企大盘、民企、全流通、其他类）做出相关性分析。
2.分析哪些因素对股市全流通方案是至关重要的，提出你的全流通方案设想。
要求兼顾国家、企业、中小投资者的利益，为稳定市场，设定一个方案实施后的股价最大振幅（比如3%）。
3.按照你的方案给出以下股票的具体实施办法：01浦发银行、08民生银行、10宝钢股份、18中国联通、20清华同方、25安阳钢铁、32申能股份、36哈药集团、37上海石化、43东方集团、46四川长虹、49张江高科。
《数学建模》习题及参考答案 第一章 建立数学模型
第一章部分习题3(5). 决定十字路口黄灯亮的时间长度.4. 在1.3节“椅子能在不平的地面上放稳吗”的假设条件中，将四角的连线呈正方形改为长方形，其余不变，试构造模型并求解.5. 模仿1.4节商人过河问题中的状态转移模型，作下面这个众所周知的智力游戏：人带着猫、鸡、米过河，船除希望要人计划之外，至多能载猫、鸡、米三者之一，而当人不在场时猫要吃鸡、鸡要吃米，设计一个安全过河方案，并使渡河次数尽量地少.6. 利用1.5节表1和表3给出的1790-2000年的美国实际人口资料建立下列模型： (1) 分段的指数增长模型. 将时间分为若干段，分别确定增长率r. (2) 阻滞增长模型. 换一种方法确定固有增长率r 和最大容量x m .7. 说明1.5节中Logistic 模型（9）可以表示为()()01t t r mex t x --+=，其中t 0是人口增长出现拐点的时刻，并说明t 0与r ，x m 的关系.8. 假定人口的增长服从这样的规律：时刻t 的人口为x (t)，t 到t +△t 时间内人口的增量与x m -x (t)成正比（其中为x m 最大容量）. 试建立模型并求解. 作出解的图形并与指数增长模型、阻滞增长模型的结果进行比较.9(3). 甲乙两站之间有电车相通，每隔10分钟甲乙两站相互发一趟车，但发车时刻不一定相同。
甲乙之间一中间站丙，某人每天在随机的时刻到达丙站，并搭乘最先经过丙站的那趟车，结果发现100天中约有90天到达甲站，约有10天到达乙站。
问开往甲乙两站的电车经过丙站的时刻表是如何安排的。
参考答案3(5). 司机看到黄灯后停车要有一定的刹车距离1s ，设通过十字路口的距离为2s ，汽车行驶速度为v ，则黄灯的时间长度t 应使距停车线1s 之内的汽车能通过路口，即()vs s t 21+≈其中s 1可由试验得到，或按照牛顿第二定律解运动方程，进一步可考察不同车重、不同路面及司机反应灵敏程度等因素的影响.4. 相邻两椅脚与地面距离之和分别定义为()()θθg f 和，将椅子旋转ο180，其余作法与1.3节相同.5. 人、猫、鸡、米分别记为4,3,2,1=i ，当i 在此岸时记1=i x ，否则记0=i x ，则此岸的状态可用()4321,,,x x x x s =表示。
数学建模练习与思考题
数学建模练习与思考题第⼀部分练习与思考题第1章建⽴数学模型1.1 在稳定的椅⼦问题中，如设椅⼦的四脚连线呈长⽅形，结论如何？（稳定的椅⼦问题见姜启源《数学模型》第6页）1.2 在商⼈们安全过河问题中，若商⼈和随从各四⼈，怎样才能安全过河呢？⼀般地，有n 名商⼈带n 名随从过河，船每次能渡k ⼈过河，试讨论商⼈们能安全过河时，n 与k 应满⾜什么关系。
（商⼈们安全过河问题见姜启源《数学模型》第7页）1.3 ⼈、狗、鸡、⽶均要过河，船需要⼈划，另外⾄多还能载⼀物，⽽当⼈不在时，狗要吃鸡，鸡要吃⽶。
问⼈、狗、鸡、⽶怎样过河？1.4 有3对夫妻过河，船⾄多载两⼈，条件是任⼀⼥⼦不能在其丈夫不在的情况下与其他的男⼦在⼀起。
问怎样过河？1.5 如果银⾏存款年利率为5.5%，问如果要求到2010年本利积累为100000元，那么在1990年应在银⾏存⼊多少元？⽽到2000年的本利积累为多少元？1.6 某城市的Logistic 模型为2610251251N N dt dN ?-=，如果不考虑该市的流动⼈⼝的影响以及⾮正常死亡。
设该市1990年⼈⼝总数为8000000⼈，试求该市在未来的⼈⼝总数。
当∞→t 时发⽣什么情况。
1.7 假设⼈⼝增长服从这样规律：时刻t 的⼈⼝为)(t x ，最⼤允许⼈⼝为m x ，t 到t t ?+时间内⼈⼝数量与)(t x x m -成正⽐。
试建⽴模型并求解，作出解的图形并与指数增长模型和阻滞增长模型的结果进⾏⽐较。
1.8 ⼀昼夜有多少时刻互换长短针后仍表⽰⼀个时间？如何求出这些时间？1.9 你在⼗层楼上欲乘电梯下楼，如果你想知道需要等待的时间，请问你需要有哪些信息？如果你不愿久等，则需要爬上或爬下⼏个楼层？1.10 居民的⽤⽔来⾃⼀个由远处⽔库供⽔的⽔塔，⽔库的⽔来⾃降⾬和流⼊的河流。
⽔库的⽔可以通过河床的渗透和⽔⾯的蒸发流失。
如果要你建⽴⼀个数学模型来预测任何时刻⽔塔的⽔位，你需要哪些信息？第2章初等模型2.1 学校共1000名学⽣，235⼈住在A 宿舍，333⼈住在B 宿舍，432⼈住在C 宿舍。
数学建模基础练习一及参考答案
数学建模基础练习一及参考答案练习1 matlab练习一、矩阵及数组操作1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）,然后将正态分布矩阵中大于1的元素变为1，将小于1的元素变为0。
2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。
3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。
4．随机生成10阶的矩阵，要求元素值介于0~1000之间，并统计元素中奇数的个数、素数的个数。
二、绘图5．在同一图形窗口画出下列两条曲线图像，要求改变线型和标记y1=2x+5；y2=x^2-3x+1，并且用legend标注。
6．画出下列函数的曲面及等高线z=sinxcosyexp(-sqrt(x^2+y^2))． 7．在同一个图形中绘制一行三列的子图，分别画出向量x=[1 5 8 10 12 5 3]的三维饼图、柱状图、条形图。
三、程序设计8．编写程序计算（x在[-8，8]，间隔0.5）先新建的，在那上输好，保存，在命令窗口代数；9．用两种方法求数列前15项的和。
10．编写程序产生20个两位随机整数，输出其中小于平均数的偶数。
11．试找出100以内的所有素数。
12．当时，四、数据处理与拟合初步1随机产生由10个两位随机数的行向量A,将A中元素按降序排列为B,再将B重排为A。
14．通过测量得到一组数据t 1 2 3 4 5 6 7 8 9 10 y 842 362 754 368 169 038 034 016 012 005 分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合，并画出散点及两条拟合曲线对比拟合效果。
15．计算下列定积分16．（1）微分方程组当t=0时，x1(0)=1，x2(0)=-0.5，求微分方程t在[0，25]上的解，并画出相空间轨道图像。
数学建模课后参考答案
数学建模课后参考答案数学建模课后参考答案数学建模是一门应用数学的学科，旨在解决实际问题。
在学习数学建模的过程中，课后作业是巩固知识、提高能力的重要环节。
然而，由于数学建模问题的多样性和复杂性，有时候我们可能会遇到一些难以解决的问题，或者对于某些题目的答案不够确定。
因此，提供一份数学建模课后参考答案是很有必要的。
1. 问题描述假设有一座小岛，岛上有一座高度为h的灯塔，灯塔的光照范围是一个圆形区域。
现在有一只船在岛外的海上，船上的人想知道距离灯塔多远的位置才能看到灯塔。
请问，船上的人应该停留在哪个位置才能看到灯塔？2. 建模过程首先，我们可以根据几何知识得出，船上的人能够看到灯塔的条件是船在灯塔的光照范围内。
因此，我们需要确定灯塔的光照范围。
灯塔的光照范围是一个圆形区域，半径为r。
根据几何知识，我们可以得出光照范围的半径与灯塔的高度之间的关系：r = √(2hR)其中，R为地球半径，h为灯塔的高度。
接下来，我们需要确定船在哪个位置才能看到灯塔。
我们可以假设船位于距离灯塔为d的位置，且船与灯塔连线与地球表面垂直。
此时，船与灯塔连线与地球表面的夹角为θ。
根据三角函数的定义，我们可以得出：tan(θ) = h/d解出θ后，我们可以得到船位于距离灯塔为d的位置时，船与灯塔连线与地球表面的夹角。
3. 答案求解根据上述建模过程，我们可以得到船位于距离灯塔为d的位置时，船与灯塔连线与地球表面的夹角θ。
如果我们已知灯塔的高度h和地球半径R，我们可以使用数学软件或计算器来计算出θ的近似值。
例如，假设灯塔的高度h为100米，地球半径R为6371千米。
我们可以使用计算器来计算出船位于距离灯塔为d的位置时，船与灯塔连线与地球表面的夹角θ的近似值。
根据公式tan(θ) = h/d，我们可以解出θ的近似值为θ ≈ 0.0157 弧度。
4. 结论根据上述计算结果，船位于距离灯塔为d的位置时，船与灯塔连线与地球表面的夹角θ的近似值为θ ≈ 0.0157 弧度。
建模仿真作业
题1,使用Simulink 创建系统,求解非线性微分方程 .其初始值为 ,绘制函数的波形.题2 力-质量系统，要拉动一个箱子（拉力f=1N)，箱子质量为M(1kg)，箱子与地面存在摩擦力[(b=0.4N(/m/s)],其大小与车子的速度成正比。
拉力作用时间为2s. 用Simulink 仿真箱体位移以及箱体速度的变化情况。
题3，力-弹簧-阻尼系统，假设箱子与地面无摩擦存在，箱子质量为M(1kg)，箱子与墙壁间有线性弹簧(k=1N/m)与阻尼器(b=0.3N/ms-1)。
阻尼器主要用来吸收系统的能量，吸收系统的能量转变成热能而消耗掉。
现将箱子拉离静止状态2cm 后放开，用Simulink 仿真试求箱子的运动轨迹。
题4，下图所示简单的单摆系统,假设杆的长度为L,且质量不计,钢球的质量为m.单摆的运动可以以线性的微分方程式来近似,但事实上系xx x x x 44)23(2=--2)0(,0)0(==x x统的行为是非线性的,而且存在粘滞阻尼,假设粘滞阻尼系数为bkg/ms-1.选取b=0.03,g=9.8,L=0.8,m=0.3, 试用Simulink仿真求取单摆角度以及振幅变化情况。
题5:蹦极跳系统：当你系着弹力绳从桥上跳下来时。
自由下落的物体满足牛顿运动定律：F=ma.假设绳子的弹性系数为k，它的拉伸影响系统的动力响应，如果定义人站在桥上时绳索下端的初始位置为0位置，x为拉伸位置，用b(x)表示绳子的张力。
m为物体的质量，g 是重力加速度，a1,a2是空气阻尼系数且均为1。
未伸长时绳索的顶端部距地面为50m。
设绳索长度30m，起始速度为0；物体质量为90kg，g为9.8m/s2,弹性系数k为20，a1和a2均为1。
试用Simulink 仿真仿真判断绳索长度是否安全？题6: 上述蹦极跳系统中，如果绳索长度不变，要保证安全，在其他条件不变的情况下，绳索弹性系数至少应该为多少？题7: 上述蹦极跳系统中，如果绳索弹性系数不变，要保证安全，绳索长度应该满足何条件？题8，重置弹力球。
哈工大画法几何与机械制图答案
第一章 制图的基本知识
1-8 按规定方法标注尺寸（数值从图中量取，取整数）。
第一章 制图的基本知识
1-9 已知椭圆长轴为70，短轴为40，作椭圆。
（1） 同心圆法
（2） 四心圆弧法
第一章 制图的基本知识
1-10 按已知斜度和锥度，完成下列图形，并标注斜度与锥度值。
′
′
′
n ′ c
n
八解
第二章 点、直线和平面
2-14 判断两直线的相对位置。
′ ′
′ ′
′ ′
′ ′
′ ′
′ ′
(1)平行
′′
′′
(2) 交叉
′′
′ ′
(3)
交
叉
′a
c
′
′
′
b
d
(4)
相交
(5)
交叉
(6)
交叉
第二章 点、直线和平面
2-15 AB与CD相交，补出线段所缺的投影。
（1）
′
′ ′
′
′
′
′ ′
′
′
′
′
′ ′
d
k l
e
′ ′
′
k′
′
k
′
′
k′
′
′
k
第三章 直线与平面、平面与平面的相对位置
3-15 过点A作平面垂直于已知直线AB。
3-16 求点M到平面 ABCD的距离。
c′
′
d′
′
c
d ′ ′
′
′
′
k′
′
′
k
第三章 直线与平面、平面与平面的相对位置
数学建模第一次实验
春季开学典礼校长致辞范文尊敬的各位领导、亲爱的老师、亲爱的同学们：大家好！在这美丽的春天里，我们迎来了新学期的开始，为了庆祝我们的开学，今天我们举行了这庄重而隆重的开学典礼。
我作为学校的校长，非常荣幸能够在这个特殊的时刻与大家相聚，并向大家致以热烈的欢迎！首先，我要向新进学校的同学表示最真诚的欢迎！希望你们能够尽快适应新的学习环境，在这里度过愉快而充实的时间。
作为学校的成员，你们将肩负着更多的责任，希望你们能够珍惜这个机会，努力学习，做一个有梦想、有担当的人。
我要感谢所有的教职员工，特别是辛勤工作了整个暑假的老师们。
正是你们的辛勤付出，为学校的发展和学生的成长提供了坚实的保障。
你们用心去教育每一个学生，你们不仅是他们的老师，更是他们的朋友和榜样。
我们要共同努力，为学生创造更好的学习环境，让每一个孩子都能够充分展示自己的才华和潜力。
在新的学期里，我们将迎来新的挑战和机遇。
作为一所优秀的学校，我们要注重培养学生的品格素养和创新能力。
只有在良好的品德教育和严格的学风氛围下，学生才能够健康成长，并发挥出自己的才华。
同时，我们要关注学生的全面发展，注重培养他们的实践能力和团队合作精神。
我们将积极开展各类实践活动和社团活动，让学生在实践中获得成长和收获。
在新的学期里，我们还将继续加强学校的教育教学工作。
教育是无形的火焰，而教师是点燃火焰的火种。
我们要努力提高教师的教学能力和专业素养，使他们能够更好地教育学生，引导他们发展潜能。
我们要注重教学方法的创新和教育资源的共享，提高教育教学质量，为学生的发展创造更好的条件。
最后，我要向全体同学提出要求：首先，要树立远大的人生目标，并为之努力奋斗。
每个人都有自己的梦想，只有明确了目标，才能够更好地前进。
其次，要把握机会，勇往直前。
人生充满了无数的选择和机遇，关键在于我们是否有勇气去抓住和利用。
再次，要注重个人素质的培养。
学习知识只是人生的一部分，更重要的是培养自己的品德和能力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、某银行经理计划用一笔资金进行证券投资业务，可供购进的证券及其相应信息如下表所示，且有如下规定和限制：（1）市政证券的收益可以免税，其它证券的收益需要按50%的税率纳税；（2）政府及代办机构的证券总共至少购进400万元；（3）所购证券的平均信用等级不超过1.4（信用等级越小，信用程度越高）；请回答下列问题：（1）若该经理有1000万资金，应如何投资？（2）如果能够以2.75%的利率借到不超过100万元，该经理应该如何操作？（3）在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？注：为简化问题起见，题中的税前收益率和利率都与年限无关，即都为固定值。
问题分析这个优化问题的目标是使投资获得最大利润，要做的决策就是投资计划，即在每个证劵下投资多少金额，决策受到信用等级、到期年限、到期税前收益率以及每种证劵的纳税税率的约束。
基本模型决策变量：设每种证劵分别投资A、B、C、D、E（万元），平均信用等级为X，平均到期年限为Y。
目标函数：设投资总金额为Q，投资的利润为W（万元），根据条件有W=A×4.3％+B×5.4％×50％+C×5.0％×50％+D×4.4％×50％+E×4.5％=0.043×A+0.027×B+0.025×C+0.022×D+0.045×E约束条件：平均信用等级X=(2×A+2×B+C+D+5×E)/ Q≤1.4平均到期年限Y=(9×A+15×B+4×C+3×D+2×E)/Q≤5非负约束所有的证劵投资均为非负值附加约束B+C+D≥400模型分析与假设每种证劵投资资金均为连续变量取值，税前收益率和利率都与年限无关；每种证劵投资资金符合比例性、可加性、连续性。
模型求解根据题设的条件，针对问题一有如下函数关系及约束条件W=0.043×A+0.027×B+0.025×C+0.022×D+0.045×EA+B+C+D+E=1000=QB+C+D≥4002×A+2×B+C+D+5×E≤1.4×Q=14009×A+15×B+4×C+3×D+2×E≤5×Q=50000≤A≤10000≤B≤10000≤C≤10000≤D≤10000≤E≤1000模型求解用LINGO软件求解，程序如下从程序运行的结果中可以得到最优解为：A=218.1818，B=0，C=736.3636，D=0，E=45.45455,W=29.83636。
即对A中证劵投资218.1818万元，对C中证劵投资736.3636万元，对E中证劵投资45.45455万元，对B和D中证劵不进行投资，此种投资方案可使投资利润最大为29.83636万元。
针对问题二在问题一的条件下，从程序运行的结果中可以看到，程序的第二条所对应的影子价格为2.983636％，大于2.75％，所以应该借足100万元进行投资。
用LINGO软件求解，程序如下从程序运行的结果中可以得到最优解为：A=240，B=0，C=810，D=0，E=50,W=32.82。
即对A中证劵投资240万元，对C中证劵投资810万元，对E中证劵投资50万元，对B和D中证劵不进行投资，此种投资方案可使投资利润最大为32.82万元。
针对问题三在问题一的条件，从程序的运行结果中，从灵敏度分析的结果中可以看到A的系数在0.0417～0.0465之间变化时，投资的最优解仍然不变，因此证劵A的税前收益增加为 4.5％时，投资不应改变；C的系数在0.02444～0.0423333之间变化时，投资的最优解仍然不变，因此C的税前收益减少为4.8％时，即0.024已经小于0.02444，投资方案应改变。
3、某公司有三个加工厂，某月公司接到4份分别来自于A，B，C，D客户的订单依次为15,17,22,12件，三个加工厂开工费用（固定值，与生产量无关）和生产量及各加工厂的至各个客户的每件运输费用如下面两个表格所示，请制定生产和运输费用，使得总费用最小。
单件运输费用（元）问题分析：这个优化问题的目标是使生产与运输费用总和最少，即要求每个客户投给每个工厂的订单量做出一个最优分配决策，即每个客户给每个加工厂分配加工量做出的决策要考虑每个加工厂的开工费用、每个工厂生产量、每个工厂到每个客户单件运输费用和每个客户的订单量这四个因素。
基本模型决策变量：设A 、B 、C 、D 四个客户分别给每个工厂的分配变量为A1、A2、A3、B1、B2、B3、C1、C2、C3、D1、D2、D3。
目标函数：设生产与运输总费用为W （元），开工费用为W1（元），运输费用为W2（元），W=W1+W2。
约束条件：开工费用：根据每个工厂生产量及A 、B 、C 、D 四个客户的订单总量66件，可判断至少有两个工厂要开工，即W1至少是两个工厂开工费用的总和。
运输费用：根据已有的费用关系可有以下关系，W2=6×A1+4×A2+8×A3+2×B1+9×B2+8×B3+6×C1+5×C2+C3+7×D1+3×D2+5×D3订单量：根据要求有，A1+A2+A3=15，B1+B2+B3=17，C1+C2+C3=22，D1+D2+D3=12 工厂生产量：根据要求有，A1+B1+C1+D 1≤39，A2+B2+C2+D2≤35，A3+B3+C3+D3≤31 非负约束：每个客户分配给每个工厂的分配额不能取负值且为整数。
模型求解：上述决策可分为两部分实现：一、只开两个工厂时的最优决策；二、三个工厂一起开工的最优决策。
一、 只开两个工厂1、 只开工厂一、工厂二根据已知条件有：A3=B3=C3=D3=0，A1+A2=15 B1+B2=17 C1+C2=22 D1+D2=12A1+B1+C1+D 1≤39 A2+B2+C2+D 2≤35 W1=91+70=161W2=6×A1+4×A2+2×B1+9×B2+6×C1+5×C2+7×D1+3×D2将上述的W2变为如下式子W2=6×A1+4×(15-A1）+2×B1+9×(17-B1)+6×C1+5×(22-C1)+7×D1+3×(12-D1)=2×A1-7×B1+C1+4×D1+359即变为如下条件W=W1+W2=2×A1-7×B1+C1+4×D1+520A1+B1+C1+ D1≤39 A1+B1+C1+ D1≥31 0≤A1≤150≤B1≤170≤C1≤220≤D1≤12用LONGO软件求解，程序如下此种情况下为A1=0，B1=17，C1=14，D1=0，A2=15，B2=0，C2=8，D2=12，A3=B3=C3=D3=0，W=4152、只开工厂二、三根据已知条件有：A1=B1=C1=D1=0，A2+A3=15B2+B3=17C2+C3=22D2+D3=12A2+B2+C2+D2≤35A3+B3+C3+D3≤31W1=70+24=94W2=4×A2+8×A3+9×B2+8×B3+5×C2+C3+3×D2+5×D3 将上述的W2变为如下式子W2=4×A2+8×（15-A2）+9×B2+8×(17-B2)+5×C2+22-C2+3×D2+5×(12-D2) =-4×A2+ B2+4×C2-2×D2+338即为W=W1+W2=-4×A2+ B2+4×C2-2×D2+432A2+B2+C2+ D2=350≤A2≤150≤B2≤170≤C2≤220≤D2≤12用LINGO软件求解，程序如下此种情况下为：A2=15，B2=8，C2=0，D2=12，A3=0，B3=9，C3=22，D3=0，A1=B1=C1=D1=0，W=3563、只开工厂一、三根据已知条件有：A2=B2=C2=D2=0，A1+A3=15B1+B3=17C1+C3=22D1+D3=12A1+B1+C1+D1≤39A3+B3+C3+D3≤31W1=91+24=115W2=6×A1+8×A3+2×B1+8×B3+6×C1+C3+7×D1+5×D3 将上述的W2变为如下式子W2=6×A1+8×（15-A1）+2×B1+8×(17-B1)+6×C1+22-C1+7×D1+5×(12-D1) =-2×A1-6×B1+5×C1+2×D1+338即为W=W1+W2=-2×A1-6×B1+5×C1+2×D1+453A1+B1+C1+ D1≤39A1+B1+C1+ D1≥350≤A1≤150≤B1≤170≤C1≤220≤D1≤12用LINGO软件求解，程序如下此种情况下为：A1=15，B1=17，C1=0，D1=3，A3=0，B3=0，C3=22，D3=9，A2=B2=C2=D2=0，W=327二、三个工厂同时开根据已知条件有：A1+A2+A3=15B1+B2+B3=17C1+C2+C3=22D1+D2+D3=121≤A1+B1+C1+D1≤391≤A2+B2+C2+D2≤351≤A3+B3+C3+D3≤31W1=91+70+24=185W2=6×A1+4×A2+8×A3+2×B1+9×B2+8×B3+6×C1+5×C2+C3+7×D1+3×D2+5×D3将上述的W2变为如下式子W2=6×A1+4×A2+8×（15-A1-A2）+2×B1+9×B2+8×(17-B1-B2)+6×C1+5×C2+22-C1-C2+7×D1+3×D2+5×(12-D1-D2)=-2×A1-4×A2-6×B1+ B2+5×C1+4×C2+2×D1-2×D2+338 即为W=W1+W2=-2×A1-4×A2-6×B1+ B2+5×C1+4×C2+2×D1-2×D2+523 1≤A1+B1+C1+ D1≤391≤A2+B2+C2+ D2≤3531≤A1+B1+C1+ D1+ A2+B2+C2+ D2≤650≤A1+ A2≤150≤B1+B2≤170≤C1+C2≤220≤D1+ D2≤12用LINGO软件求解，程序如下此种情况下为：A2=15，B1=17，C3=22，D2=12，其余均为0，W=337结果分析从上面的求解中可以得出最优解为只开工厂一、三，此时的总费用最少为327，最终的决策为：工厂一加工客户A和B的所有订单以及D客户的3件，工厂三加工D客户的9件。