08-数学实验之八—计算机模拟课件.pptVIP

数学实验之八—计算机模拟;3、培养和树立采用计算机模拟解决实际问题尤其是复杂问题的能力。;概述; 状态随时间连续变化的系统称为连续系统。通常该系统的模型一般可以用微分方程的形式表达，通过一些物理机理推导出来。模拟结果往往是近似的。 例如，追逐问题。;1、追逐问题;模拟方法： 1. 建立平面直角坐标系； 2. 以时间间隔△t 进行采样，并且计算各个 时刻下的状态： A: (x1(t), y1(t)) → (x1(t+ △t ), y1(t+ △t )) B: (x2(t), y2(t)) → (x2(t+ △t ), y2(t+ △t )) (x1(t+ △t ), y1(t+ △t )) （由几何原理） ≈(x1(t)+v △t cos(θ), y1(t) +v △t sin(θ);选取足够小的△t ，模拟到任意两人的距离小于v △t 为止。;v=1; dt=0.05; d=20; % x=[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8)] x=[0 0 0 10 10 10 10 0]; x(9)=x(1); x(10)=x(2); hold on axis(equal); axis([0 10 0 10]); while(d0.1) for i=1:2:7 d=sqrt((x(i)-x(i+2))^2+(x(i+1)-x(i+3))^2); x(i)=x(i)+v*dt*(x(i+2)-x(i))/d; x(i+1)=x(i+1)+v*dt*(x(i+3)-x(i+1))/d; plot(x(i),x(i+1),.) end x(9)=x(1);x(10)=x(2); end hold on; 某军的一导弹基地发现正北方向120km处海面上有敌艇一艘以90km/h的速度向正东方向行驶。该基地立即发射导弹跟踪追击敌艇，导弹速度为450km/h，自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌艇？;P(x(t), y(t));消去t：;如何决定导弹位置 P2(x2, y2)?;function [num,y_j,L,T]=simu1(x,y,t,eps) k=0; while k1000 p=90*k*t-x; q=120-y; d1=p/(p^2+q^2)^0.5; d2=q/(p^2+q^2)^0.5; x=x+450*t*d1; y=y+450*t*d2; if (abs(q)eps) break; end k=k+1; end num=k; L=x; y_j=y; T=L/90;;（2）如果当基地发射导弹的同时，敌艇立即由仪器发觉。假定敌艇为高速快艇，它即刻以135千米/小时的速度与导弹方向成垂直的方向逃逸，问导弹何时何地击中敌艇？ 不停地改变逃跑策略，运动轨迹如何？ 修正敌艇速度：90 → 135千米/小时;离散系统;排队问题;输入情况： 船只（顾客）到达时间和服务时间。 系统状态： 排队等候的船只（顾客）数目（队 　　　　长）L(t) 　　　　　卸载设备（服务员)是否在工作或服务效率等；; 排队系统中，顾客到达时刻数据如何收集？对每个顾客的服务时间如何？ X： x1，x2，…，xn;随机数的产生;3、(非)常见分布随机数如何产生？（离散）;由均匀随机数R∈[0,1]，决定X的抽样值。 当0R＜0.3，抽样值xi∈（0, 3.5]； xi=3.5 当0.3≤R＜0.8，抽样值xi∈（3.5, 5.5]； xi=5.5 当R≥0.8，抽样值xi∈（5.5, 8]； xi=8;r=rand(10); lisanrnd1.m y=[]; for i=1:10 if 0r(i)r(i)0.3 y(i)=3.5; elseif 0.3=r(i)r(i)0.8 y(i)=5.5; else y(i)=8; end end y % 产生随机变量数组。 ;MATLAB中的随机数生成函数;卸载货物（排队服务）系统的模拟;1、模拟目的;假设由概率分布产生了如下的数据：;3、系统的状态;4、初始条件和终止条件;5、模拟过程;ATDT;6、系统的性能指标;7、与分析模型的比较;分析模型基本假设： 1、到达间隔服从参数λ的指数分布； 2、服务时间服从参数为μ的指数分布； 3、排队规则是“先到先服务”，顾客数目和排队长度都无限制。 排队论中将这样的单服务员系统记作M/M/1。;服务强度：;