数学建模计算机模拟.ppt

在Matlab软件包中编程，共需三个Ｍ－文件:randlp.m,mylp.m,lpconst.m.主程序为randlp.m.%mylp.mfunctionz=mylp(x)%目标函数z=2*x(1)^2+x(2)^2-x(1)*x(2)-8*x(1)-3*x(2);%转化为求最小值问题第30页,共37页，星期六，2024年，5月%randlp.mfunction[sol,r1,r2]=randlp(a,b,n)%随机模拟解非线性规划debug=1;a=0;%试验点下界b=10;%试验点上界n=1000;%试验点个数r1=unifrnd(a,b,n,1);%n?1阶的［a,b］均匀分布随机数矩阵r2=unifrnd(a,b,n,1);sol=[r1(1)r2(1)];z0=inf;fori=1:nx1=r1(i);x2=r2(i);lpc=lpconst([x1x2])；iflpc==1z=mylp([x1x2])；ifzz0z0=z；sol=[x1x2]；endendendToMatlab(randlp)返回第31页,共37页，星期六，2024年，5月实验作业1、编一个福利彩票电脑选号的程序。第32页,共37页，星期六，2024年，5月4.某设备上安装有四只型号规格完全相同的电子管，已知电子管寿命为1000--2000小时之间的均匀分布。当电子管损坏时有两种维修方案，一是每次更换损坏的那一只；二是当其中一只损坏时四只同时更换。已知更换时间为换一只时需1小时，4只同时换为2小时。更换时机器因停止运转每小时的损失为20元，又每只电子管价格10元，试用模拟方法决定哪一个方案经济合理？5.导弹追踪问题：设位于坐标原点的甲舰向位于x轴上点A(1,0)处的乙舰发射导弹，导弹头始终对准乙舰.如果乙舰以最大的速度(是常数)沿平行于y轴的直线行驶，导弹的速度是5，模拟导弹运行的轨迹.又乙舰行驶多远时，导弹将它击中？第33页,共37页，星期六，2024年，5月初始化:i=0,k1=0,k2=0,k3=0i=i+1骰子点数?k1=k1+1k2=k2+1k3=k3+1k1=k1+1i＜20?E=(k2+k3)/20E1=0*k1/20+1*k2/20+2*k3/20停止硬币正面?YNNY1，2，34，56第34页,共37页，星期六，2024年，5月投掷硬币的计算机模拟1、产生服从U（0，1）的随机数R12、将区间[0，1]两等分：若，则对应硬币正面若，则对应硬币反面第35页,共37页，星期六，2024年，5月掷骰子的计算机模拟1、产生服从U（0，1）的随机数R22、将区间[0，1]六等份：若，则对应骰子点数为1若，则对应骰子点数为2若，则对应骰子点数为3若，则对应骰子点数为4若，则对应骰子点数为5若，则对应骰子点数为6第36页,共37页，星期六，2024年，5月初始化:i=0,k1=0,k2=0,k3=0i=i+1R2=?k1=k1+1k2=k2+1k3=k3+1k1=k1+1i＜20?E=(k2+k3)/20E1=0*k1/20+1*k2/20+2*k3/20停止R1=0.5YNNYR23/6其它R25/6ToMatlab(liti1)第37页,共37页，星期六，2024年，5月实验目的实验内容学习计算机模拟的基本过程与方法。1、模拟的概念。4、实验作业。3、计算机模拟实例。2、产生随机数的计算机命令。第2页,共37页，星期六，2024年，5月连续系统模拟实例:追逐问题离散系统模拟实例:排队问题用蒙特卡洛法解非线性规划问题返回计算机模拟实例第3页,共37页，星期六，2024年，5月模拟的概念模拟就是利用物理的、数学的模型来类比、模仿现实系统及其演变过程，以寻求过程规律的一