实验课(六).pptVIP

数学模型实验  （六） 仿真模型 例 1. 池水含盐 池中有水 2000 m3，含盐 2 kg， 以 6m3 / 分 的速率向池中注入浓度 0.5 kg / m3 的盐水 又以 4 m3 / 分的速率从池中流出混合后的盐水 问欲使池中盐水浓度达到 0.2 kg / m3，需要多长时间 V(t+1)=V(t)+(rI-rO) =V(t)+2 S(t+1)=S(t)+[rIpI-rOp(t)]=S(t)+[6*0.5-4*P(t)] P(t+1)=S(t+1)/V(t+1) t=0;v=[2000];s=[2];p=[0.001];x=[0]; while p(end)=0.2; t=t+1;v1=v(t)+2;s1=s(t)+3-4*p(t);p1=s1/v1; v=[v,v1];s=[s,s1];p=[p;p1];x=[x,t]; end; n=x(end)+1;m=floor(n/10);Z=zeros(m+1,n); for k=1:m;Z(k,(k-1)*10+1)=1;end; Z(m+1,n)=1;A=[x’,((10^(-3)).*v)’,s’,p’);B=Z*A 问题 1. 在池水含盐的问题中令 rO=rI=6m3/分 10. 池中盐水的浓度如何变化？ 20. 若当p(t) = 0.3kg/m3 时令pI=0, 需要多少时间达到 p*= 0.2kg / m3? t=0;v=[2000];s=[2];p=[0.001];x=[0]; while p(end)=0.3; t=t+1;v1=v(t);s1=s(t)+3-6*p(t);p1=s1/v1; v=[v,v1];s=[s,s1];p=[p;p1];x=[x,t]; end; while p(end)=0.2; t=t+1;v1=v(t);s1=s(t)-6*p(t);p1=s1/v1; v=[v,v1];s=[s,s1];p=[p;p1];x=[x,t]; end; n=x(end)+1;m=floor(n/10);Z=zeros(m+1,n); for k=1:m;Z(k,(k-1)*10+1)=1;end; Z(m+1,n)=1;A=[x’,v’,s’,p’);B=Z*A 问题2. 若池中盐水的初始浓度为 p*=0.2， 对于不同的初始体积V0和流速，当pI=0时， 计算池中盐水浓度降低一半（p(t)=0.1)所用的时间。 V0=1000, 2000, 3000, 4000, 5000 r=4,5,6,7,8,10 t=0;v=[2000];s=[400];p=[0.2];x=[0]; while p(end)=0.1; t=t+1;v1=v(t);s1=s(t)-6*p(t);p1=s1/v1; v=[v,v1];s=[s,s1];p=[p;p1];x=[x,t]; end; n=x(end)+1;m=floor(n/10);Z=zeros(m+1,n); for k=1:m;Z(k,(k-1)*10+1)=1;end; Z(m+1,n)=1;A=[x’,v’,s’,p’);B=Z*A 例2. 战争问题. 两支军队，红军(R)和蓝军(B)，进行战斗。 在这场常规战中，伤亡是由于直接交火（步兵）和火炮射击（炮兵） 。 假设直接交火的伤亡数与敌军步兵数成正比。 由炮火造成的伤亡数与敌军的炮兵数和友军的密度两者都有关系。 红军聚集了五个师袭击两个师的蓝军。 蓝军具有防御的和武器精良的优势。 蓝军为赢得战斗该尽多大的努力？ 我们将用两个状态变量： x1=R，红军部队的兵力单位数量 x2=B，蓝军部队的兵力单位数量 战斗问题的离散时间动态模型： x1(t+1)=x1(t) - a1x2 - b1x1x2 x2(t+1)=x2(t) - a2x1 - b2x1x2 初始x1(0)=5和 x2(0)=2。时间步长?t=1小时。 令 a2 =0.05, b2=0.005, 假设 a1=?a2 , b1 =? b2, ?1。 对：?=1、1.5、2、3 和5模拟这场战斗的结局。 改变a2的数值，改变初始军力的对比，影响如何？ 考虑战争问题中天气对战争的