MATLAB实验-8计算机模拟.docVIP

PAGE  PAGE 4 新乡学院 数学与信息科学系实验报告 实验名称 计算机模拟 所属课程 数学软件与实验 实验类型 综合型实验 专 业 信息与计算科学 班 级 2011级1班 学 号 11111021012 姓 名 李欢丽 指导教师 朱耀生 一、实验概述【实验目的】 学会用模拟法和建立微分方程的方法计算每隔一段时间流动容器内某种物质的含量. 并用图像表示. 学会用rand()函数来模拟满足均匀分布的随机试验. 【实验原理】 rand(m,n)生成一个满足均匀分布的m*n随机试验，矩阵的每个元素都在(0,1)之间。 subplot(m,n,p)函数把一个图形窗口分割成m*n个子区域，可通过参数p调用各个绘图区域进行操作，子绘图区域的编号为按行从左至右编号。 ezplot(x,y,[tmin,tmax])符号作图命令函数，x,y为关于参数的符号函数，[tmin,tmax]为的取值范围 【实验环境】 MatlabR2010b二、实验内容问题1 容器内有100L盐水，含盐10kg。现以每分钟3L的速度注入含有0.01kg的盐水（假定两种盐水立即均匀混合），同时又以每分钟2L的均匀速度从容器内抽出盐水.按要求完成下面题目： 1． 用模拟法每隔10min计算一次容器内盐水的体积、含盐量和含盐率，列出一张表，从表中观测60min后盐水的含盐量； 2． 建立微分方程；计算60min后容器内盐水的含盐量； 3． 做出用上述两种方法得到的容器内盐水含盐量随时间的变化曲线. 1.分析问题 根据题意分析出隔一时间t后，容积内盐水体积，含盐量，含盐率，然后推导出通式后用matlab实现。2.问题求解 v=100; s=10; r=0.1; a=[0;100;10;0.1]; for i=10:10:60 v=v+10; s=s+0.01*10; s=s*(1-2/v); v=v-2; t=i/10+1; a(1,t)=i; a(2,t)=v; a(3,t)=s; a(4,t)=r; end a subplot(2,1,1) plot(a(1,:),a(3,:)) s=dsolve(Ds=0.01-(s+0.01*t)/(100+3*t)*2,s(0)=10,t); subplot(2,1,2) ezplot(s,[0,60])3.结果 0 10.0000 20.0000 30.0000 40.0000 50.0000 60.0000 100.0000 108.0000 116.0000 124.0000 132.0000 140.0000 148.0000 10.0000 9.9164 9.8466 9.7887 9.7411 9.7025 9.6718 0.1000 0.0918 0.0849 0.0789 0.0738 0.0693 0.06544.结论及分析 经多次实验，推知分析无误，问题2 两人约定于12点到1点到某地会面，先到者等20分钟后离去，试求两人能会面的概率？ 1.分析问题 用rand()函数来模拟满足均匀分布的随机试验.2.问题求解 n=100000; m=0; for i=1:n x=rand(1); y=rand(1); if abs(x-y)=1/3 m=m+1; end end m/n3.结果 ans = 0.44324.结论及分析 经多次试验证明分析无误。三、实验小结  subplot(m,n,p)函数把一个图形窗口分割成m*n个子区域，可通过参数p调用各个绘图区域进行操作，子绘图区域的编号为按行从左至右编号。