下面来介绍一下在连续的污染源的情况下,如何提高MATLAB编程来实现模拟四维空间分布:
本程序需要注意一下几个问题:
1.本程序运行需要1-3分钟左右时间(视个人电脑配置情况而定);
2.两个M文件都放在MATLAB同一目录下;
3.在命令窗口中输入lianxuwuran(t),回车即可运行,其中t为时间,例如lianxuwuran(20),表示经过20秒扩散的仿真图像。
4.运行出来的图形有等高线数值,但是某些等高线可能没有等高线数值,这是因为数值过于密集,为了避免重叠不便展示,故而进行了编辑。
首先在MATLAB的主界面编辑器中分别写入下列代码:
主程序代码:
%这是主程序,命名为lianxuwuran.m
function Value=lianxuwuran(t)
%现在开始编写连续污染源模型
global xx yy zz;
xmin=-500;dx=10;xmax=500;ymin=-100;dy=5;ymax=100;zmin=-10;dz=1;zmax=10;
Cxyz_t((xmax-xmin)/dx+1,(ymax-ymin)/dy+1,(zmax-zmin)/dz+1)=0;
ii=0;jj=0;kk=0;
for zz=zmin:dz:zmax;kk=kk+1;
for yy=ymin:dy:ymax;jj=jj+1;
for xx=xmin:dx:xmax;ii=ii+1;
Cxyz_t(ii,jj,kk)=quadl(@fun3D,1,t);
end