用Matlab作三维图谢谢!

用Matlab作三维图谢谢!

clear all;[x,y]=meshgrid(0:50:25000,-2000:50:2000); %设定网格,下风向排放点距离x米,烟气中心轴水平距离y米Q=30;%源强,释放速度gam1=0.110726;gam2=0.104534;alp1=0.929418;alp2=0.826212;gam3=0.146669;gam4=0.400167;alp3=0.888723;alp4=0.632023; sigy=gam1*x.^alp1.*(x>=0&x<=1000)+gam3*x.^alp3.*(x>1000);%水平扩散系数sigz=gam2*x.^alp2.*(x>=0&x<=1000)+gam4*x.^alp4.*(x>1000);%垂直扩散系数u=3.24;%风速He=183.51;z=0;%高度 Qpi=2*Q./(2*pi*u.*sigy.*sigz+eps);%计算公式第一部分ex1=exp(-0.5*(y./(sigy+eps)).^2);%计算公式第二部分ex2=exp(-0.5*((z-He)./(sigz+eps)).^2)+exp(-0.5*((z+He)./(sigz+eps)).^2);%计算公式第三部分C=Qpi.*ex1.*ex2;%计算扩散气体浓度 mesh(x,y,C);xlabel('x轴下风向距离(m)'),ylabel('y轴向距离(m)'),zlabel('气体扩散浓度'),[Cmax,I]=max(C(:));%计算c的最大值（Cmax）,并得到下标xmax=x(I);%利用下标得到x的最大值（xmax）

我学会了