数学模型与实验报告姓名：王珂班级：121111学号：442指导老师：沈远彤数学模型与实验一、数学规划模型某企业将铝加工成A,B两种铝型材，每5吨铝原料就能在甲设备上用12小时加工成3 吨A 型材，每吨A 获利2400 元，或者在乙设备上用8 小时加工成4 吨B 型材，每吨B 获利1600 元。
现在加工厂每天最多能得到250 吨铝原料，每天工人的总工作时间不能超过为480 小时，并且甲种设备每天至多能加工100 吨A,乙设备的加工能力没有限制。
（1 ）请为该企业制定一个生产计划，使每天获利最大。
（2）若用1000 元可买到1 吨铝原料，是否应该做这项投资若投资，每天最多购买多少吨铝原料（3）如果可以聘用临时工人以增加劳动时间，付给工人的工资最多是每小时几元（4）如果每吨A 型材的获利增加到3000 元，应否改变生产计划题目分析：每5 吨原料可以有如下两种选择：在甲机器上用12 小时加工成3 吨A 每吨盈利2400 元在乙机器上用8 小时加工成4 吨B 每吨盈利1600 元限制条件：原料最多不可超过250 吨，产品A 不可超过100 吨。
工作时间不可超过480 小时线性规划模型：设在甲设备上加工的材料为x1 吨，在乙设备上加工的原材料为x2 吨，获利为z，由题意易得约束条件有：Max z = 7200x1/5 +6400x2/5x1 + x2 三25012x1/5 + 8x2/5 三4800 三3x1/5 三100, x2 三0用LINGO 求解得：VARIABLE VALUE REDUCED COSTX1X2ROW SLACK OR SURPLUS DUAI PRICE1234做敏感性分析为：VARIABLE CURRENT ALLOWABLE ALLOWABLECOFF INCREASE DECREASEX1X2ROW CURRENT ALLOWABLE ALLOWABLERHS INCREASE DECREASE234INFINITY可见最优解为x1= 100, x2=150, MAXz=33600Q因此最优解为在甲设备上用100吨原料生产A产品，在乙设备上用150吨原料生产B产品。
最大盈利为336000. 由运算结果看约束条件1 (原料)的影子价格是960,即每增加1吨原料可收入960,小于1000元，因此不购入。
同理可得，每小时的影子价格是40元，因此聘用员工的工资不可超过每小时40 丿元。
4、由敏感性分析可得，在最优解不变的前提下，x1予许的变化范围上限是1920, 下限是1280。
若每吨A获利增加到3000,价值系数变为1800,在允许范围内，所以保持原计划不变。
二、微分方程模型在鱼塘中投放n0尾鱼苗，随着时间的增长，尾数将减少而每尾的重量将增加。
设尾数n(t)的(相对)减少率为常数；由于喂养引起的每尾鱼重量的增加率与鱼的表面积成正比，由于消耗引起的每尾鱼重量的减少率与重量本省成正比。
分别建立尾数和每尾鱼重的微分方程，并求解。
用控制网眼的办法不捕小鱼，到时刻T才开始捕捞，捕捞能力用尾数的相对减少量|n/n|表示，记作E,即单位时间捕获量是En(t)。
问如何选择T和E,使从T 开始的捕获量最大。
基本假设：1. 鱼塘里的鱼无繁殖，且不会自然死亡。
2. 鱼苗尾数相对减少率为常数。
3. 由于喂养引起的每尾鱼重量的增加率与其表面积成正比；由于消耗引起的每尾鱼重量的减少率与本身重量成正比。
4. 将鱼简化为椭球体，且其密度分布均匀，初始状态相同。
2由基本假设：鱼苗尾数n t 相对减少率为常数，则可得以下微分方程:由基本假设:由于喂养引起的每尾鱼重量的增加率与其表面积成正比; 重量的减少率与重量本身成正比。
可得以下微分方程：k .s t k 2m t 2dt又因为要通过设定渔网网格面积来确定最大捕获量， 而渔网网格面积由每尾鱼的 最小横截面相关，又每尾鱼的横截面面积与鱼的表面积相关。
由基本假设中鱼群 的表面积服从正态分布，即：2s t u2eh其中U 为st 的均值，2为st 的方差。
dnt dtkn t由于消耗引起的每尾鱼则在此条件下：NTPstsT nt 4又由N t Ent 5得：E P s t s T 6模型的求解：关于鱼尾数随时间变化的微分方程组:dn tkn tdtn 0 n0可直接求解得：n t n0e kt7又椭球体的体积为:V 4 abc38表面积近似为：2s 4 abc 39又m V10则可得:413不妨设W b T c T此时s T 4aTbTcT则s T 4Wa T由基本假设s t 服从正态分布，则11则将 11 式代入式可得:dmtdt 4 k 13m t 4k 2m t12m b所以求解可得4 k 1k 213 m°4 k 1k 2则:141516P s t s T1 P s t s Tst u s T u1 Ps T u14Wa T u1其中t为标准正态分布函数则由此将渔网网眼面积和单位时间最大捕获量联系起来，此时仅需将通过调查将a t函数进行研究，进而使得P s t sT 取得最大值，则此时N t P s t s T nt取得最大值又E P s t s T则可通过查找标准正态分布表求得结论。
三、统计回归模型下表列出了某城市18位35岁一44岁经理的年平均收入X1千元，风险偏好度X2和人寿保险额y千元，其中风险偏好度是根据发给每个经理的问卷调查表综合评估得到的，它的数值越大，就越偏爱高风险。
研究人员想研究此年龄段中的经理所投保的人寿保险额与年平均收入及风险偏好度之间的关系。
研究者预计，经理的年均收入和人寿保险额之间存在着二次关系，并有把握地认为风险偏好度对人寿保险额有线性效应，但对风险偏好度对人寿保险额是否有二次效应以及两个自变量是否对人寿保险额有交互效应，心中没底。
请你通过表中的数据来建立一个合适的回归模型，验证上面看法，并给出进一步的分析。
数学模型解：为大致分析y与x1和x2关系，首先利用表1的数据分别作岀y对于x1和x2的散点图(见图1和图2中的圆点)x1=[ ]；>> y仁[196 63 252 84 126 14 49 49 266 49 105 98 77 14 56 245 133 133];>> p=polyfit(x1,y1,2)P =+000 +001>> x2=0::85;y2=polyval(p,x2); plot(x1,y1,'o',x2,y2)y对X1的散点图从图中可以发现，随着治的增加，y的值有明显向上弯曲的二次增长趋势，图中的曲线是用二次函数模型2y 0i x i 2 x i(1)拟合的。
(其中是随机误差)>> x3=[7 5 10 6 4 5 4 6 9 5 2 7 4 3 5 1 8 6]; >> q=polyfit(x3,y1,1)+001 +001>> x4=0::15;y3=polyval(q,x4); plot(x3,y1,'o',x4,y3)y 对X 2的一次的散点图从图中可以发现，随着X 2的增加，y 的值比较明显的线性增长趋势, 是用线性函数模型y 01X 2（2）拟合的。
（其中是随机误差）综合上面的分析，结合模型（1）和（2）建立如下的回归模型2y 01X 1 2x 2 3X 11015图中的曲线 2C 0150WO50(3)(3)式右端的x i和X2称为回归变量，° i X i 2X2 3X2是给定年平均收入X i、风险偏好度X2时，人寿保险额y的值，其中的参数0, 1, 2, 3称为回归系数。
还有影响y的其它因素作用都包含在随机误差中。
模型求解：使用MATLAB统计工具箱的命令regress求解，求解过程如下>> x1=[ ];x2=[7 5 10 6 4 5 4 6 9 5 2 7 4 3 5 1 8 6]; x3=x1.*x1; x0=o nes(18,1); x=[x0 xi' x2'x3'];y=[196 63 252 84 126 14 49 49 266 49 105 98 77 14 56 245 133 133];>> [b,b in t,r,ri nt,stats]=regress(y',x,b =+001+000 bint =+001 +001 +000+000 +000 stats =+004 +000由此得到模型（3）的回归系数估计值及其置信区间（置信水平0.05）、检验统计量R2,F,p的结果见下表结果分析：R2=指因变量y （保险额）的％可由模型确定，F的值远远超过F的检验的临界值，p远小于，因此模型（3）从整体来看是可用的。