数学建模短学期7.docVIP

数学建模短学期作业7 1. 某报童以每份0.03元的价格买进报纸,以0.05元的价格出售. 根据 长期统计,报纸每天的销售量及百分率为 销售量200210220230240250百分率0.100.200.400.150.100.05已知当天销售不出去的报纸，将以每份0.02元的价格退还报社.试确定报童每天买进报纸数量,使报童的平均总收入为最大? 解：1、问题假设： （1） 模拟时间充分大； （2） 报童购买报纸量介于销售量最小值与最大值之间； （3）不考虑有重大事件发生时卖报的高峰期，也不考虑风雨天气时卖报的低谷期。 2、符号假设： BUYMIN：每天的最小购买量 BUYMAX：每天的最大购买量 SIMUDAY：模拟时间 sell_amount：报童销售量 buy_amount： 报童购买量 percentage：销售百分率 ave_profit：总平均利润 loop_buy ：当天购买量 loop_day ：当天时间 利用函数GetProfit.m： Code：nction re=GetProfit(a,b)if ab % 供不应求：报童购买量小于销售量??? re=a*(0.05-0.03); else % 供过于求：报童购买量大于销售量??? re=b*(0.05-0.03)+(a-b)*(0.02-0.03); end 在命令窗口输入： BUYMIN=200;?? %? 每天的最小购买量BUYMAX=250;?? %? 每天的最大购买量SIMUDAY=1.0e+5;? % 模拟时间 sell_amount=200:10:250;?????? % 销售量percentage=[0.1 0.3 0.7 0.85 0.95 1];? %? 百分率buy_amount=0;ave_profit=0;for loop_buy=BUYMIN:BUYMAX??? sum_profit=0;??? for loop_day=1:SIMUDAY??????index=find(percentage=rand); % 产生随机数，用于决定当天的销售量??????sum_profit=sum_profit+GetProfit(loop_buy,sell_amount(index(1)));??? end??? buy_amount=[buy_amount,loop_buy]; % 循环嵌套??? ave_profit=[ave_profit,sum_profit/SIMUDAY]; % 循环嵌套endbuy_amount(1)=[]; % 第一个元素置空ave_profit(1)=[];[val,id]=max(ave_profit)??? % 显示最大平均收入valbuy=buy_amount(id)??? % 显示在平均收入最大情况下的每天的购买量buyxlabel=每天的购买量;ylabel=平均利润;plot(buy_amount,ave_profit,*:); 得： 该结果说明当报童每天买进报纸数量为220，报童的平均总收入为最大，且最大为4.2801。 毕 2、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据： 养护时间x234579121417212856抗压强度y354247535965687376828699试求型回归方程. 解：利用matlab拟合工具箱，可找出x与y的函数关系： x=0.2995*exp（0.05278*y） （1） 曲线和散点图基本拟合，转化（1）式： y=22.843+18.947*lnx 毕 3、在研究化学动力学反应过程中，建立了一个反应速度和反应物含量的数学模型，形式为 其中是未知参数，是三种反应物（氢，n戊烷，异构戊烷）的含量，y是反应速度.今测得一组数据如下表，试由此确定参数，并给出置信区间. 其中的参考值为（1，0.05, 0.02, 0.1, 2）. 序号反应速度y氢x1n戊烷x2异构戊烷x318.554703001023.79285801034.8247030012040.02470