实验：计算机模拟

实验：计算机模拟

实验：计算机模拟一、实验目的和要求1.掌握计算机模拟的基本原理及方法;2.掌握利用随机数解决问题;二、实验内容1.已知零件C由零件A和零件B连接而成，已知零件A、B的长度为随机变量，数值如下表，试抽取10个样本计算C的平均长度。零件A的长度56789概率0.070.190.380.250.11零件B的长度14151617概率0.230.410.270.09解答：输入xxx.m文件aaa=[];fori=1:10a=rand;ifa<=007a1=5;elseifa>0.07&&a<=0.26a1=6;elseifa>0.26&&a<=0.64a1=7;elseifa>0.64&&a<=0.89a1=8;elseifa>0.89&&a<=1a1=9;endb=rand;ifb<=0.23b1=14;elseifa>0.23&&a<=0.64b1=15;elseifa>0.64&&a<=0.91b1=16;elseifa>0.91&&a<=1b1=17;endc=a1+b1aaa=[aaaa1+b1];end\n运行xxx得到以下结果：xxxc=21c=19c=20c=21c=20c=20c=21c=19c=19c=19>>mean(aaa)ans=19.900010个样本计算C的平均长度为19。2.某服装公司欲预定一批冬装出售，每件冬装的加工费不确切，估计如下：单件成本789101112概率0.050.150.20.30.250.05已知该服装销售量与订价有关预测各销售量概率如下：定价\概率\预计销售500600700800900190.050.150.40.250.15200.30.20.20.20.1210.20.30.350.10.05试用模拟方法决定该公司冬装的订购数与订价，使利润最大？（如订购量大于销量时，每件以5元处理）根据题设编写代码：a.mb=[500600700800900];a=[19,20,21];x1=rand;ifx1<0.05c=7;elseifx1>0.05&x2<0.20c=8;elseifx1>0.20&x1<0.4c=9;elseifx1>0.4&x1<0.7\nc=10;elseifx1>0.7&x1<0.95c=11;elsec=12;endcx2=rand;ifx2<0.05d=500;elseifx2>0.05&x2<0.2d=600;elseifx2>0.2&x2<0.6d=700；elseifx2>0.6&x2<0.85d=800;elsed=900;enddfori=1:5forj=1:3b(i)a(j)ifb(i)<=dz=b(i)*a(j)-b(i)*celsez=d*a(j)-b(i)*c+(b(i)-d)*5endendend运行a.mc=11d=600ans=500ans=19z=4000ans=500ans=20z=4500ans=500ans=21z=5000ans=600ans=19z=4800ans=600ans=20z=5400ans=600ans=21z=6000ans=700ans=19z=4200ans=700ans=20z=4800ans=700ans=21z=5400ans=800ans=19z=3600ans=800ans=20z=4200ans=800ans=21z=4800ans=900ans=19z=3000ans=900ans=20z=3600ans=900ans=21z=4200\n可知：当生产600件定价为21元时，利润最大为6000元。3一个慢跑者在平面上沿椭圆以恒定的速率v=1跑步,设椭圆方程为:x=10+20cost,y=20+5sint．突然有一只狗攻击他．这只狗从原点出发,以恒定速率w跑向慢跑者,狗的运动方向始终指向慢跑者．分别求出w=20,w=5时狗的运动轨迹．（用计算机模拟的方法）建立m3.m文件：w=20;v=1;dt=0.1;d=100;x=0;y=0;fort=0:0.1:100d=sqrt((x1(t)-x)^2+(y1(t)-y)^2);x=w*((x1(t)-x)/d)*dt+x;y=w*((y1(t)-y)/d)*dt+y;plot(x1(t),y1(t),'r+',x,y,'*')holdonend建立x1.mfunctionx=x1(t)x=10+20*cos(t);建立y1.mfunctiony=y1(t)y=20+5*sin(t);运行m3（第二次时将w=20改为w=5）