matlab第二次实验报告重点.doc

《数学建模实验》 上机实验报告 班级：计算机35班 学号： 2130505099 姓名： 田博文 【实验一】 　　一上底面半径2米、下底面半径4米、高4米的圆台形水池内盛满了水，由池底一横截面积为0.001平方米的小孔放水。求在任意时刻的水面高度和将水放空所需的时间。 一、问题分析： 在dt很小一段时间内，水的流速v可以看成是不变的与水面高度h有关的函数， 即v=； 在dt时间内流出的水量，可以近似的为一圆柱， 对应的圆柱高为dh； 由于容器下降的水量与流出的水量相等，令小孔面积为b， 即dV=πdh=bvdt； 可得dh与dt的关系为：dh=bvdt/(π)； 根据相似三角形原理可得， r=4-h/2； 初始条件有； h(0)=4； 为表示方便，引入变量L=4-h 故r=2+L/2 代码： T=0; for L=0:0.001:4-0.001 %选取积分微元为0.001，即dh V=0.001*pi*(2+(1/2*L))^2; %水面处的dV，b=0.001 t=V/(sqrt(2*(4-L)*9.8)*0.001) %dV与dt关系 T=T+t; plot(T,(4-L),red); hold on end 运行结果： 结论： 　　在任意时刻的水面高度和将水放空所需的时间如上图红线所示，横轴为时间，纵轴为水面高度。 【实验二】 　　有A、B、C三个场地，每一个场地都出产一定数量的原料，同时也消耗一定数量的产品，具体数据如下表所示。已知制成每吨产品需要消耗3吨原料，A、B两地，A、C两地和B、C两地之间的距离分别为150千米、100千米和200千米，假设每万吨原料运输1千米的运费为5000元，每万吨产品运输1千米的运费为6000元。由于地区条件的差异，在不同地区设厂的费用不同，由于条件的限制，在B处建厂的规模不能超过5万吨，问：在这三地如何建厂、规模建多大才能使得总费用最小？ 地点年产原料（万吨）年销产品（万吨）生产费用（万元/万吨）A207150B1613120C240100 问题分析： 设nij为i地运往j地的原料量，mij为i地运往j地的产品量，设A地为1地，B为2地，C地为3地。因制成每吨产品需要消耗3吨原料故设A、B、C三地生产规模分别为a、b、c?。 a=(n21+n31)/3, b=(n12+n32)/3, c=(n13+n23)/3, 约束条件为: n12+n13==20, n21+n23==16, n31+n32==24, %年产原料 m21+m31==7, m12+m32==13, m13+m23==0, %年销产品 b=5, %B处建厂规模不能超过5万吨 a==m12+m13, b==m21+m23, c==m31+m32, %保证三地原料无浪费 生产费用： 150a+120b+100c 运输原料费用： 0.5[150(n12+n21)+200(n23+n32)+100(n13+n31)] 运输产品费用： 0.6[150(m12+m21)+200(m23+m32)+100(m13+m31)] 可用矩阵实现 二、代码： c=[75 75 50 50 100 100 150 240 210 120 160 220]; a=[1 -1 1 -1 0 0 3 3 0 0 0 0 -1 1 0 0 1 -1 0 0 3 3 0 0 0 0 -1 1 -1 1 0 0 0 0 3 3 0 0 0 0 0 0 0 0 1 1 0 0]; b=[20;16;24;5]; aeq=[0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 1 0 1 0 1]; beq=[7;13]; vlb=zeros(12,1); [x,final]=linprog(c,a,b,aeq,beq,vlb,[]) 运行结果： 结论： 总方案为：A地生产7万吨，B地生产5万吨，C地生产8万吨， 总费用为：3485万元 【实验三】 　　27个立方形空盒排成3×3×3的三位列阵。如果三个盒子在同一条水平线上，或同一条垂直线上，或同一条对角线上，则认为三盒一线。这样的线共有49条：水平线18条，垂直线9条，水平面对角线6条，垂直面对角线12条，对角面对角线4条。现有白球13个，黑球14个，每个盒子中放入一个球，如何投放，使有单一色球的线数最少？ 问题分析： 每个立方体按顺序编号，成为数组。要想表示单一色球，则盒中装白球置0，黑球置1，如此一来，一条线上均为单一色球则可表示为此线上三个立方体对应值为0或3。 引入函数，功能为算出每种投放方法对应的单一色球的线数