计算机模拟版本3.pptVIP

数学建模培训;概述; 状态随时间连续变化的系统称为连续系统。通常该系统的模型一般可以用微分方程的形式表达，通过一些物理机理推导出来。模拟结果往往是近似的。 例如，追逐问题，浓度问题。;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(θ);3. 选取足够小的△t ，模拟到任意两人的距 离小于v △t 为止。 4. 连接四人在各时刻下的位置，就得到所求的 运动轨迹。 Matlab (simu2.m);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：; x=dsolve(D2x*(120-y)/sqrt((Dx)^2+1)=90/450,Dx(0)=0,x(0)=0,y);如何决定导弹位置 P2(x2, y2)?;P(x(t), y(t));function [num,y_j,L,T]=simu1(x,y,t,eps) k=0; while k1000 (simu1.m) 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;;x=[0,0]; (huatu.m) axis(equal); axis([0,140,0,140]);grid on hold on t=0.001; %步长 while (abs(x(2)-120)0.1) % 终止条件() for k=1:280 p=90*k*t-x(1); q=120-x(2); d1=p/(p^2+q^2)^0.5; d2=q/(p^2+q^2)^0.5; x(1)=x(1)+450*t*d1; x(2)=x(2)+450*t*d2; x1(1)=90*k*t; x1(2)=120; h1=line(color,[0,0.2,0.4],linewidth,2); h2=line(color,[0,0.6,0.9],linewidth,3); set(h1,xdata,x1(1),ydata,x1(2)); set(h2,xdata,x(1),ydata,x(2)); end end hold on;（2）如果当基地发射导弹的同时，敌艇立即由仪器发觉。假定敌艇为高速快艇，它即刻以135千米/小时的速度与导弹方向成垂直的方向逃逸，问导弹何时何地击中敌艇？ 不停地改变逃跑策略，运动轨迹如何？ 修正敌艇速度：90 → 135千米/小时;二、离散系统;　在我方某前沿防守地域，敌人以一个炮排（含两门火炮）为单位对我方进行干扰和破坏．为躲避我方打击，敌方对其阵地进行了伪装并经常