搜狗做题网用Matlab作三维图
来源：学生作业帮 编辑：搜狗做题网作业帮 分类：综合作业 时间：2024/04/28 04:03:08
用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）
转到另两个方向：
在workspace中查看cmax和xmax的值为3.4226e-005和8100,为了检验是否正确直接在图中读取,证明代码正确：
[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）
转到另两个方向：在workspace中查看cmax和xmax的值为3.4226e-005和8100,为了检验是否正确直接在图中读取,证明代码正确：