Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2009-06-18
年龄: 40
帖子: 9
声望力: 0 ![]() |
![]()
各位高手,请教我想出等高线图,却做不出来,麻烦哪位帮帮忙啊!!
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 ![]() ![]() 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/attach...VzlFNPdKlE.jpg |
![]() |
![]() |