数学建模计算机模拟实验

1、实验6：计算机模拟班级：2015级电科班，学号：222015333210187，姓名：吴京宣，第1组一、实验目的：学习计算机模拟的基本过程与方法实验器材：配有网卡的PC机二台一组，安装有WindowsXP以上版本的操作系统二、实验内容1模拟的概念2产生随机数的计算机命令3计算机模拟实例三、实验步骤对以下问题,编写M文件:1编一个福利彩票电脑选号的程序2.某报童以每份0.03元的价格买进报纸,以0.05元的价格出售.根据长期统计,报纸每天的销售量及百分率为销售量200210220230240250百分率0.100.200.400.150.100.05已知当天销售不出去的报纸，将以每份0.02元的

2、价格退还报社.试用模拟方法确定报童每天买进多少份报纸,能使平均总收入最大?某设备上安装有4只型号规格完全相同的电子管，已知电子管寿命服从10002000h之间的均匀分布.电子管损坏时有两种维修方案，一是每次更换损坏的那只；二是当其中1只损坏时4只同时更换已知更换时间为换1只时需1h,4只同时换为2h.更换时机器因停止运转每小时的损失为20元，又每只电子管价格10元，试用模拟方法确定哪一个方案经济合理？导弹追踪问题：设位于坐标原点的甲舰向位于x轴上点A(l,0)处的乙舰发射导弹，导弹头始终对准乙舰如果乙舰以最大的速度(常数)沿平行于y轴的直线行驶，导弹的速度为5.模拟导弹运行的轨迹.乙舰行驶多远

3、时，导弹将它击中？三、实验结果要求模拟彩票选号，则由最常见的双色球摇号为例。双色球由红球和蓝球两部份组成，从33个红球号码(0133)中选择6个，再从16个蓝球号码(0116)中选择1个。开奖时，在红色球中随机摇出六个红号，在蓝色球中随机摇出一个蓝号。双色球的投注方法2）2）可选择单式或复式。一注彩票只选取六个红号和一个蓝号来组成为单式投注；复式：一注彩票选取多个红号或多个蓝号，可提高中奖率。因此选号的方式也契合于单式或者复式的摇号，第一段程序是单式摇号（6红1蓝），第二段程序是复试摇号（15红3蓝）Red=O00000：fori=l:6k=rand;forj=0:32if(kj/33&k.i

4、/16&k:/i+Qy1&)Blue=i+l;endendBlueendRed=30Blue=g531-21410Red=000000000000000;fori=l:15k=rand;forj=0:32ifBISkfj+i0劈)ReafiJj+l；endendRedendBlue=.000;for1=1：芬k=rand;forj=0:15if(k険jj/16止k-6)Bluej/16也Kb(购进价)c(退回价)售出一份赚a-b；退回一份赔b-c每天购进多少份可使收入最大？购进太多-卖不完退回-赔钱购进太少不够销售-赚钱少应根据需求确定购进量每天需求量是随机的-每天收入是随机的优化问题的目标函

5、数应是长期的日平均收入，等于每天收入的期望BUYMIN：每天的最小购买量；BUYMAX:每天的最大购买量；SIMUDAY:模拟时间;sell_amount:报童销售量；buy_amount：报童购买量percentage：销售百分率；ave_profit:总平均利润；loop_buy:当天购买量；oop_day:当天时间；调查需求量的随机规律每天需求量为r的概率f(r),r=0,1,2设每天购进n份，日平均收入为G(n)已知售出一份赚a-b；退回一份赔b-c，可得函数与程序如下：EraiN=200;BUYMAX250;.SIMUDAY=1.Oe+5.e11_aniount=200:10:250

6、:percentage=0.10.30.70.850.951:buy_ajnourLt=0;ave_profit=0:forloop_buy=BUYKIN:BUYJUXuiR_profit=0;forloop_day=l:SIMUDAYindeK.：=find(percentage=rand);sujr_profit=si.un_profit+GetProfiop_buy?sell_amount(indes)；endbuy_amountLjloop_buy;怨ave_profiitf.SlMUDAY;sndbuy_amouiit.(O=;aCve_profit(1;val,id=maK(ave

7、_profbuy=buy_amount(idXxlabel=?每天的购买重；ylabel=平均利润；plot(buy_amountnave_profit?)fijnctiondy=eql(k?y)dy=ZBros(2n1);dy(l)=y(2);dy(25=l/5*sqrt(1+y(1)2)/1-k)s=750money2=结果如下：说明每天进货在220左右可以达到最大，最大在4.2793左右。val=4.2793id=21buy=2203由于要求使用时间随即分布，所以，使用了unifrnd语句实现使用时间的估计,然后使用while循环计算不同的方案下会花费的成本。其中s为方案一的花money

8、2为方案二的花费。由以下结果可知方案一的花费低于方案二。moneyl=0,0,0,0;fori=l:4r=unifrnd(1000?2000);t=r;whilet=10000moneyl(i)=moneyl(i)+30r=urnfrnd(1000?2000);t=t+r;endends=sujn(moneyl)money2=0;r=unifrnd(1000?2000);tr;whilet=10000money2=money2+200;rifrndflOOO,2000);t=七+匸；endinoney24.导弹飞行始终对准乙舰，则在任意时间导弹的位置与斜率之间的关系为dy/dx=(v*t-y)

9、/(l-x)，编写函数表示此表达式，然后使用odel5s语句拟合出导弹运动的路径，同时作出乙船的运动路径，由此可得一个交点坐标，就是相遇点。在设置参数时xf=1无法正常运行，因此将其改为0.9999得以正常运算。fijnctiondy=eqly)dy=zeros(2,1);dy(l)=y(2);dy(2)=l/5*sqrt(1+y(1)2)/(1-k)k0=0;sf=0.9999;k,y=odel5seql?k0sf?00);ploty?1)/b.);holdony=0:0.01:2;plot(1?y?b*?)结果如下：导弹大致在(1，0.2)处击中乙舰实验总结1本次实验旨在利用Matlab的操作界面和基本的数据