rainbowthink
2010-02-04, 17:40
各位高手，请教我想出等高线图，却做不出来，麻烦哪位帮帮忙啊！！
clc,clear
Q=83.5*1000;%d=410mm,质量流量，g/s;
u=1;%风速，m/s;
H=0;
xmin=0;xmax=350;ymin=-100;ymax=100;z=0;n=1000;
[x,y]=meshgrid(xmin:((xmax-xmin)/(n-1)):xmax,ymin:((ymax-ymin)/(n-1)):ymax);%网格划分
y1=0.22*x.*(1+0.0001*x).^(-1/2);
z1=0.2*x;
c=Q./(2*pi*u*y1.*z1).*exp(-y.^2./(2*y1.^2)).*(exp(-(z-H)^2./(2*z1.^2))+exp(-(z+H)^2./(2*z1.^2)))*1000; %高斯方程计算扩散浓度
%绘制多条等值浓度曲线程序
cm=6700:100:113800;
for i=1:length(cm)
a(i)=contour(x,y,c,[cm(i) cm(i)]);
clabel(a(i))
hold on
end
xlabel('下风向距离x（m）'),ylabel('横风向距离y（m）'),title('特定浓度的等值线');
最好出这样的图
http://www.chinavib.com/forum/attachments/month_0903/20090324_6f387074816aefe7fbc9BYVzlFNPdKlE.jpg
clc,clear
Q=83.5*1000;%d=410mm,质量流量，g/s;
u=1;%风速，m/s;
H=0;
xmin=0;xmax=350;ymin=-100;ymax=100;z=0;n=1000;
[x,y]=meshgrid(xmin:((xmax-xmin)/(n-1)):xmax,ymin:((ymax-ymin)/(n-1)):ymax);%网格划分
y1=0.22*x.*(1+0.0001*x).^(-1/2);
z1=0.2*x;
c=Q./(2*pi*u*y1.*z1).*exp(-y.^2./(2*y1.^2)).*(exp(-(z-H)^2./(2*z1.^2))+exp(-(z+H)^2./(2*z1.^2)))*1000; %高斯方程计算扩散浓度
%绘制多条等值浓度曲线程序
cm=6700:100:113800;
for i=1:length(cm)
a(i)=contour(x,y,c,[cm(i) cm(i)]);
clabel(a(i))
hold on
end
xlabel('下风向距离x（m）'),ylabel('横风向距离y（m）'),title('特定浓度的等值线');
最好出这样的图
http://www.chinavib.com/forum/attachments/month_0903/20090324_6f387074816aefe7fbc9BYVzlFNPdKlE.jpg