连续污染源的四维模拟程序

下面来介绍一下在连续的污染源的情况下，如何提高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