容器内有100L盐水,含盐10kg.现以每分钟3L的速度注入含有0.01kg的盐水（假定两种盐水立即均匀混合）,同时又以每分钟2L的均匀速度从容器内抽出盐水.按要求完成下面题目：⑴用模拟法求出60min

Question

容器内有100L盐水,含盐10kg.现以每分钟3L的速度注入含有0.01kg的盐水（假定两种盐水立即均匀混合）,同时又以每分钟2L的均匀速度从容器内抽出盐水.按要求完成下面题目：⑴ 用模拟法求出60min（每隔5min计算一次）时容器内盐水的含盐率,；⑵ 建立微分方程；计算60min后容器内盐水的含盐量；⑶ 做出用上述两种方法得到的容器内盐水含盐量随时间的变化曲线.

Answer 1

V=100; W=10; for i=1:60 V=V+3; W=W+3*0.01; W=W*(1-2/V); V=V-2; if mod(i,10)==0 data(i/10,1)=i; data(i/10,2)=V; data(i/10,3)=W; data(i/10,4)=W/V; end end figure(1) subplot(2,1,1) plot(data(:,1),data(:,3)) %设含盐量为x,时间为t,其微分方程为dx/dt=3*0.01-2*x/(100+t) s=dsolve('Dx=3*0.01-2*x/(100+t)','x(0)=10') subplot(2,1,2) ezplot(s,[10,60])