用MATLAB绘制南昌地区室外气象焓湿频图程序.docx 15页

  • 108
  • 0
  • 0
  • 约6.32千字
  • 2018-05-22 发布

用MATLAB绘制南昌地区室外气象焓湿频图程序.docx

文档工具:
    1. 1、本文档共15页,可阅读全部内容。
    2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,可选择认领,认领后既往收益都归您。
    3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形,可联系本站下载客服投诉处理。
    4. 文档侵权举报电话:19940600175。
    附录一Vaporpress子程序:function [p]=Vaporpress(t)%定义子函数名c1=-5674.5359;c2=6.3925247;c3=-0.9677843*10^-2;c4=010^-6;c5=010^-18;c6=-0.9484024*10^-12;c7=4.1635019;c8=-5800.2206;c9=1.3914993;c10=-0.048640239;c11=010^-4;c12=-010^-7;c13=6.5459673;B=101005;%南昌地区大气压力T=t+273.15;if t<0 p=exp(c1/T+c2+c3*T+c4*T^2+c5*T^3+c6*T^4+c7*log(T));elsep=exp(c8/T+c9+c10*T+c11*T^2+c12*T^3+c13*log(T));endend附录二焓湿图:c1=-5674.5359;c2=6.3925247;c3=-0.9677843*10^(-2);c4=010^(-6);c5=010^-18;c6=-0.9484024*10^-12;c7=4.1635019;c8=-5800.2206;c9=1.3914993;c10=-0.048640239;c11=010^-4;c12=-010^-7;c13=6.5459673;B=101005;xlim([-1 45]);ylim([-30 65]);%等温线t=-30;while(t<=65) [p]=Vaporpress(t);%引用饱和蒸汽压子程序 d=622*p/(B-p);if d<=42 d=d;else d=42;endx1=0; y1=t;x2=d; y2=((1.01*t+0.001*(2501+1.84*t)*d)/1.01-2.5*d);line([x1 x2],[y1 y2]); t=t+1;end%等相对湿度线绘制w=1;while(w>0.04) t=-30;while(t<=65) [p]=Vaporpress(t); d=w*622*p/(B-w*p);if d<=42 d=d;else d=42;endx1=d; y1=(1.01*t+d*0.001*(2501+1.84*t))/1.01-2.5*d; t=t+0.1; [p]=Vaporpress(t); d=w*622*p/(B-w*p);if d<=42 d=d;else d=42;end y2=(1.01*t+d*0.001*(2501+1.84*t))/1.01-2.5*d;x2=d;line([x1 x2],[y1 y2]);end w=w-0.05;end%等含湿量线绘制t=-30;d=0;while(d<43) [p]=Vaporpress(t); d1=622*p/(B-p);if d1<d t1=t-0.05;while(d1<d) t1=t1+0.05; T1=t1+273.15;if t1<0p1=exp(c1/T1+c2+c3*T1+c4*T1^2+c5*T1^3+c6*T1^4+c7*log(T1));elsep1=exp(c8/T1+c9+c10*T1+c11*T1^2+c12*T1^3+c13*log(T1));end d1=622*p1/(B-p1);end%得到大于d的d1 t2=t1-0.05; t3=(t1+t2)/2;else t1=t+0.05;while(d1>d) t1=t1-0.05; T1=t1+273.15;if t1<0p1=exp(c1/T1+c2+c3*T1+c4*T1^2+c5*T1^3+c6*T1^4+c7*log(T1));else p1=exp(c8/T1+c9+c10*T1+c11*T1^2+c12*T1^3+c13*log(T1));end d1=622*p1/(B-p1);end%小于d的d1 t2=t1+0.05; t3=(t2+t1)/2;end [x1,x

    文档评论(0)

    • 内容提供方:moon8888
    • 审核时间:2018-05-22
    • 审核编号:6014225002001152

    相似文档