预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共34页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

计算机仿真引论简单的仿真举例计算机仿真(computersimulation)是系统科学、运筹学和计算机应用的一个分支,是科学研究和解决问题的一种有效方法,在各个领域有着广泛的应用。长江三峡工程进行仿真的理由:例1混凝土搅拌中心的位置解1:数学计算(解析解) 记第K个工地的位置为(xk,yk), 中心的位置为(x0,y0), 则目标函数为解2:计算机的模拟。 将平面区域划分成许多小方格,利用计算机循环计算 搅拌中心设立在每个交叉点的费用,再经过比较就可以得到最佳位置。仿真的常用术语:仿真时钟:1。系统分析3。模型的运行与改进仿真步骤图:时间步长法选取系统的一个初始的状态作为仿真时钟的零点,仿真时钟每步进一次,就对系统的所有的实体和属性以及活动进行一次全部的扫描考察,按照预定的计划和目标进行分析,计算和记录系统状态的变化,这个过程一直进行到仿真的时钟结束为止。其流程图为:例2池水含盐量问题池水含盐量随时间变化的示意图理论分析法解和图形计算机仿真计算机仿真程序 h=1;s(1)=2;w(1)=2000;r(1)=s(1)/w(1);t(1)=h; y(1)=(2000000+3000000*h+3000*h^2+h^3)/(1000+h)^2; fori=2:200 t(i)=i*h;s(i)=s(i-1)+0.5*6*h-4*h*r(i-1); w(i)=w(i-1)+2*h;r(i)=s(i)/w(i); y(i)=(2000000+3000000*t(i)+3000*t(i)^2+t(i)^3)/(1000+t(i))^2; m=floor(i/10); ifi/10-m<0.1 tm(m)=m;wm(m)=w(i);sm(m)=s(i);rm(m)=r(i); end ifr(i)>0.2t02=i*h;r02=r(i);break;end end [10*tm',sm',rm'] [t02,r02] plot(t,s,'blue',t,y,‘red')计算机仿真的的结果结果比较池水含盐量问题的推广池水含盐量推广问题的示意图池水含盐量推广问题的仿真流程计算机仿真程序 h=1;t(1)=h; s1(1)=0;w1(1)=1000;r1(1)=s1(1)/w1(1); s2(1)=0;w2(1)=1000;r2(1)=s2(1)/w2(1); s3(1)=0;w3(1)=1000;r3(1)=s3(1)/w3(1); fori=2:200 t(i)=i*h;s1(i)=s1(i-1)+4*h-3*h*r1(i-1); w1(i)=w1(i-1)+h;r1(i)=s1(i)/w1(i); s2(i)=s2(i-1)+r1(i-1)*3*h+2*(1-cos(t(i)*3.14/180))*h-2*h*r2(i-1); w2(i)=w2(i-1)+h+2*(1-cos(t(i)*3.14/180))*h; r2(i)=s2(i)/w2(i); s3(i)=s3(i-1)+r2(i-1)*2*h+2*(1-sin(t(i)*3.14/180))*h-h*r3(i-1); w3(i)=w3(i-1)+h+2*(1-sin(t(i)*3.14/180))*h; r3(i)=s3(i)/w3(i); end figure(1);plot(t,s1,'blue',t,s2,'red',t,s3,'green') figure(2);plot(t,r1,'blue',t,r2,'red',t,r3,'green') figure(3);plot(t,w1,'blue',t,w2,'red',t,w3,'green')水的变化池1---兰,池2---红,池3---绿含盐量的变化池1---兰,池2---红,池3---绿含盐率的变化池1---兰,池2---红,池3---绿讨论问题此区域形状复杂,理论分析困难,可以用计算机仿真实现。将可能的区域等分,考察每个小区域是否在此区域中,将在此区域中的小面积想加即可其仿真图如下:Matlab程序为 x=-2:0.01:6; y=-2:0.01:2; s=0; h=0.01; fori=1:800 forj=1:400 xx=-2+i*h; yy=-2+j*h; ifxx^2+yy^2<=16 ifxx^2/36+yy^2<=1 if(xx-2)^2+(yy+1)^2<=9 s=s+h^2; end end end end end s 运行后给出面积的值8.8310。